En période hivernale, le golfe du Lion est sujet à des conditions de vents régionaux forts (Mistral et/ou Tramontane) qui transportent des masses d'air continentales froides et sèches au dessus de la mer. Ces événements, les Cold Air Outbreaks (CAO) , conduisent à d'intenses échanges air-mer et donc à un pompage de chaleur qui favorise la formation d'eaux denses et le déclenchement de la convection océanique profonde. La bonne représentation de ces échanges air-mer intenses dans les modèles de climat et de prévision numérique du temps reste à l'heure actuelle une problématique majeure. Elle est au cœur du projet ASICS-MED centré sur compréhension des mécanismes de formation d'eaux denses en Méditerranée et qui s'inscrit dans le cadre de la thématique " Échanges air-mer intenses " du programme HyMeX dédié à l'étude du cycle de l'eau en Méditerranée. Les processus qui s'opèrent au sein de la couche limite atmosphérique marine (CLAM) et de la couche de mélange océanique (CMO) interagissent entre eux à différentes échelles spatiales et temporelles. La compréhension de l'évolution globale de la CLAM mais également des mécanismes locaux nécessitent la prise en compte de l'ensemble des processus. L'étude présentée ici est consacrée à la structure moyenne et turbulente de la CLAM en conditions de vents forts. L'objectif est de déterminer comment l'organisation du champ turbulent est impactée lors d'épisodes de CAO et d'estimer les flux de surface associés à ces conditions de vents forts. La méthodologie adoptée est basée sur l'utilisation conjointe d'observations aéroportées collectées lors de la campagne de mesure SOP2 d'HyMeX et de simulations numériques. La campagne de mesure SOP2 d'HyMeX qui a eu lieu au cours de l'hiver 2013 dans le golfe du Lion a permis de documenter grâce à l'avion de recherche ATR42 la structure moyenne et turbulente de la CLAM lors de 11 épisodes de CAO. Une analyse spectrale s'appuyant sur un modèle analytique a été réalisé sur 181 paliers (i.e. segments de vol rectilignes et stabilisés en altitude). Les profils verticaux des échelles turbulentes caractéristiques ainsi que la forme du spectre de la vitesse verticale ont permis de mettre en évidence un allongement des structures énergétiques dans l'axe du vent moyen associé à l'organisation du champ turbulent sous la forme de rouleaux longitudinaux. Une description unidirectionnelle du champ turbulent tridimensionnel peut conduire à une représentativité limitée des structures cohérentes au sein des échantillons. Cependant, la connaissance des profils de flux sur toute l'épaisseur de la CLAM est nécessaire pour l'estimation des échanges air-mer. Une méthode de correction des flux turbulents calculés par eddy correlation a été appliqué afin de prendre en compte les erreurs systématique et aléatoire relatives à la mesure et au traitement de données. Cette correction a permis de déterminer les meilleures estimations possibles des flux extrapolés à la surface avec une marge d'incertitude pour les 11 épisodes de CAO documentés lors de la campagne SOP2 d'HyMeX. La comparaison de ces estimations aéroportées aux autres sources d'information dérivées de paramétrisations des flux a permis de mettre en évidence une sous-estimation systématique du flux de chaleur latente en conditions de vents forts. Une approche numérique a permis de compléter l'analyse de la structure moyenne et turbulente de la CLAM lors d'épisodes de CAO. / During winter, local strong winds (Mistral or Tramontana) occurred in the Gulf of Lion which bring cold and dry continental air over a warmer sea. Those events, the cold air outbreaks, can lead to intense air-sea interactions which favour dense water formation and deep oceanic convection. The representation of air-sea exchanges is a fundamental aspect of of climate modelling and numerical weather forecasting. The ASICS-MED project aims to identify fine-scale processes leading to dense water formation and is a part of the "Intense air-sea exchanges" topic of the HyMeX program devoted to hydrological cycle in the Mediterranean. The processes occurring within the marine atmospheric boundary layer (MABL) and the oceanic mixing layer (ML) interact with one another at different spatial and temporal scales. Understanding the overall evolution of the MABL but also the local mechanisms requires taking into account all the processes. The study presented here is devoted to the mean and turbulent structure of the MABL under strong wind conditions. The objective is to determine how the organization of the turbulent field is impacted during CAO events and to estimate the surface fluxes associated with these strong wind conditions. The methodology adopted is based on the joint use of airborne observations collected during the HyMeX-SOP2 field campaign and numerical simulations. The HyMeX-SOP2 field campaign took place during the winter of 2013 in the Gulf of Lion. The research aircraft ATR42 was operated to document the mean and turbulent structure of the MABL during 11 CAO events. A spectral analysis based on an analytic model was carried out on 181 legs (i.e. stacked straight and level runs stabilized in altitude). The vertical profiles of the turbulent characteristic scales as well as the shape of the vertical velocity spectrum revealed an elongation of the energy structures in the mean wind direction associated with the organization of the turbulent field into longitudinal rolls. A unidirectional sampling of the three-dimensional turbulent field may lead to a limited representativeness of the coherent structures within the samples. However, knowledge of kinematic fluxes profiles over the entire thickness of the CLAM is necessary to estimate air-sea exchanges. A correction method was applied to turbulent fluxes calculated by eddy correlation in order to take into account systematic and random errors related to measurement and data processing. This correction made it possible to determine the best possible estimates of the extrapolated surface fluxes with a margin of uncertainty for the 11 CAO events documented during the HyMeX-SOP2 field campaign. The comparison of these airborne estimates with the other sources of information derived from bulk parameterizations show a systematic underestimation of the latent heat flux under strong wind conditions. A numerical approach allowed to complete the analysis of the mean and turbulent structure of the MABL during CAO events. The numerical study, based on the non-hydrostatic Meso-NH model, focuses on an episode of strong Tramontana with winds greater than 25m/s documented during the HyMeX-SOP2 field campaign. In a first step, a one-dimensional framework made it possible to understand the forcing terms necessary to reproduce in a realistic way the development of the observed MABL. This reference configuration allowed, in a second time, a Large-Eddy Simulation of the CAO event. This simulation has been validated using airborne data and has allowed to deepen the description of the turbulent field as well as the evolution of the coherent structures oriented in the axis of the mean wind.

Understanding the interactions between ocean and atmosphere in regions of oceanographic fronts is of vital importance for the improvement of numerical models for weather and climate forecasting. In the South Atlantic Ocean (SAO) the meeting between the warm waters of the Brazil Current (BC) and the cold waters of the Malvinas (Falkland) Current (MC) in the region known as the Brazil-Malvinas Confluence (BMC), results in intense mesoscale oceanic activity and, for this reason, this region is considered one of the most energetic of the Global Ocean. The interactions resulting from the thermal contrast in regions oceanographic fronts of the OAS are investigated in this work through estimates of heat fluxes based on data collected in situ and by satellite. The results of this study show that the response to the thermal contrasts found in the ocean is in the form of heat fluxes and these fluxes are critical in modulating the atmospheric boundary layer (ABL). Estimation based on data collected in situ show that in the warm side (north) of the oceanographic front the fluxes are more intense (latent heat: 62 W/m² and sensible heat: 0.6 W/m²) than in the cold side (south) (latent heat: 5.8 W/m² and sensible heat: -13.8 W/m²). In the South Atlantic Current (SAC) along the 30° S parallel, heat fluxes are directly related to the meandering characteristic of the current. The data collected in situ, in addition to allow heat flux estimates at a better spatial resolution, were used to develop a new method for estimating the heat energy exchanged between the atmosphere and the ocean caused by the presence of mesoscale oceanic structures. This methodology consists in the comparison of a radiosonde profile taken over waters of the structure of interest and another taken over waters which do not belong to this structure. The methodology was used to estimate the heat energy transfer between the atmosphere and the ocean over the top of three structures sampled in the OAS. The estimation of the heat energy transferred by a warm eddy detached from the BC points to an energy in the latent (sensible) form of 1.6 1017 J (-2.8 1016 J) which corresponds to approximately 0.011 % of the total heat energy of the eddy transferred to the atmosphere during the field experiment and 0.78 % transferred during the supposed lifetime of the eddy (3 months). Along the CSA two oceanic structures were studied: (i) a cold meander that receives from the atmosphere energy in the latent (sensible) form of 1.4 106 J/m2 (5.4 105 J/m2), and (ii) warmer waters associated with a detached eddy from the Agulhas Current (AC) that transfer to the atmosphe heat energy of approximately 4 106 J/m2 an 5.7 106 J/m2 in the latent and sensible forms, respectively. The estimation of heat energy transfer on top of mesoscale oceanic structures clearly demonstrate the importance of these structures for the heat exchanges between the ocean and the atmosphere and must be taken into account in future works about this subject in the SAO. / A compreensão das interações entre oceano e atmosfera em regiões de frentes oceanográficas é de vital importância para o melhoramento de modelos numéricos de previsão do tempo e clima. No Oceano Atlântico Sul (OAS) o encontro entre as águas quentes da Corrente do Brasil (CB) com as águas frias da Corrente das Malvinas (CM), na região denominada Confluência Brasil-Malvinas (CBM), resulta em intensa atividade oceânica de mesoescala e, por esse motivo, essa região é considerada uma das mais energéticas do Oceano Global. As interações resultantes do contraste termal ao longo de regiões de frentes oceanográficas no OAS são investigadas neste trabalho através de estimativas de fluxos de calor baseadas em dados de satélite e dados coletados in situ. Os resultados do trabalho demonstram que a resposta aos contrastes termais encontrados no oceano se dá na forma de fluxos de calor e que esses fluxos são fundamentais na modulação da Camada Limite Atmosférica (CLA). As estimativas com base em dados coletados in situ demonstram que no lado quente (norte) da frente oceanográfica os fluxos são mais intensos (calor latente: 62 W/m² e calor sensível: 0,6 W/m²) que nos lado frio (sul) (calor latente: 5,8 W/m² e calor sensível: -13,8 W/m²). Na Corrente Sul Atlântica (CSA), ao longo do paralelo de 30° S, os fluxos de calor estão diretamente relacionados a característica meandrante da corrente. Os dados coletados in situ, além de possibilitarem estimativas de fluxo de calor com uma melhor resolução espacial, foram usados no desenvolvimento de uma nova metodologia para estimativa da energia calorífica trocada entre oceano e atmosfera em virtude da presença de estruturas oceânicas de mesoescala. Essa metodologia consiste na comparação entre um perfil de radiossonda tomado sobre águas da estrutura de interesse e outro tomado sobre águas que não pertencem a essa estrutura. A metodologia desenvolvida foi utilizada para determinar a transferência de energia calorífica entre oceano e atmosfera em três estruturas amostradas no OAS. A estimativa da energia calorífica transferida por um vórtice quente desprendido da CB aponta para uma energia na forma latente (sensível) de 1,6 1017 J (-2,8 1016 J) que corresponde a aproximadamente 0,011 % da energia calorífica total do vórtice transferida durante o experimento de campo e de 0,78 % da energia do vórtice transferidos durante o tempo suposto de vida do vórtice (3 meses). Ao longo da CSA, duas estruturas oceânicas foram estudadas: (i) um meandro frio que recebe da atmosfera uma energia na forma latente (sensível) de 1,4 106 J/m2 (5,4 105 J/m2) e (ii) águas mais quentes associadas a um vórtice desprendido da Corrente das Agulhas (CA) que transferem para a atmosfera uma energia calorífica de aproximadamente 4 106 J/m2 e 5,7 106 J/m2 nas formas latente e sensível, respectivamente. As estimativas da transferência de energia calorífica sobre estruturas oceânicas de mesoescala demonstram claramente a importância destas nas trocas de calor entre o oceano e a atmosfera e devem ser levadas em consideração em trabalhos futuros sobre o tema no OAS.

Doppler radar wind profilers are already widely used to measure atmospheric winds throughout the free troposphere and stratosphere. Several methods have been developed to quantify atmospheric turbulence with such radars, but to date they have remained largely un-tested; this thesis presents the first comprehensive validation of one such method. Conventional in-situ measurements of turbulence have been concentrated in the surface layer, with some aircraft and balloon platforms measuring at higher altitudes on a case study basis. Radars offer the opportunity to measure turbulence near continuously, and at a range of altitudes, to provide the first long term observations of atmospheric turbulence above the surface layer. Two radars were used in this study, a Mesosphere-Stratosphere-Troposphere (MST) radar, at Capel Dewi, West Wales, and the Facility for Ground Based Atmospheric Measurements (FGAM) mobile boundary layer profiler. In-situ measurements were made using aircraft and tethered-balloon borne turbulence probes. The spectral width method was chosen for detailed testing, which uses the width of a radar's Doppler spectrum as a measure of atmospheric velocity variance. Broader Doppler spectra indicate stronger turbulence. To obtain Gaussian Doppler spectra (a requirement of the spectral width method), combination of between five and seven consecutive spectra was required. Individual MST spectra were particularly non-Gaussian, because of the sparse nature of turbulence at its observation altitudes. The width of Gaussian fits to the Doppler spectrum were compared to those from the `raw' spectrum, to ensure that non-atmospheric signals were not measured. Corrections for non-turbulent broadening, such as beam broadening, and signal processing, were investigated. Shear broadening was found to be small, and the errors in its calculation large, so no corrections for wind shear were applied. Beam broadening was found to be the dominant broadening contribution, and also contributed the largest uncertainty to spectral widths. Corrected spectral widths were found to correlate with aircraft measurements for both radars. Observing spectral widths over time periods of 40 and 60 minutes for the boundary layer profiler and MST radar respectively, gave the best measure of turbulence intensity and variability. Median spectral widths gave the best average over that period, with two-sigma limits (where sigma is the standard deviation of spectral widths) giving the best representation of the variability in turbulence. Turbulent kinetic energies were derived from spectral widths; typical boundary layer values were 0.13 m 2.s (-2) with a two-sigma range of 0.04-0.25 m 2.s (-2), and peaked at 0.21 m 2.s (-2) with a two-sigma range of 0.08-0.61 m 2.s (-2). Turbulent kinetic energy dissipation rates were also calculated from spectral widths, requiring radiosonde measurements of atmospheric stability. Dissipation rates compared well width aircraft measurements, reaching peaks of 1x10 (-3) m 2.s (-3) within 200 m of the ground, and decreasing to 1-2x10 (-5) m 2.s (-3) near the boundary layer capping inversion. Typical boundary layer values were between 1-3x10 (-4) m 2.s (-3). Those values are in close agreement with dissipation rates from previous studies.

Title: Urban Ventilation Dependence on Geometric Configuration Author: RNDr. Ing. Libor Kukačka Department: Department of Atmospheric Physics Supervisor: prof. RNDr. Zbyněk Jaňour, DrSc., Academy of Sciences of the Czech Republic, Institute of Thermomechanics, v. v. i. Abstract: The main goal of the thesis is to investigate the impact of urban geometry on the urban ventilation using wind-tunnel modelling. To measure the pollutant transport, both advective and turbulent, within complex urban geometries with a high temporal resolution a special measurement method was developed. At first, the pollution of a simplified urban area was simulated by a ground-level point source and the ventilation of the intersection with respect to four wind directions was studied. Later, the pollution of other simplified and complex urban areas was simulated by a ground-level line source and the ventilation of three different street canyons with respect to wind direction perpendicular and oblique to their along-canyon axis was investigated. The clear impact of urban complexity and wind direction on street canyon ventilation is demonstrated at lateral and top openings of all investigated canyons and the intersection. Whilst the dominance of the pollutant advection is demonstrated at the eaves of pitched roofs, at the roof ridges...

In this thesis, we have evaluated a stochastic Lagrangian model for computing particle deposition rates with prospects to use for atmospheric deposition applications.  The model is one-dimensional and models the particle dynamics in the boundary layers near walls and obstacles by simulating the coherent turbulent structures and Brownian motion governing the wall-normal transport. The deposition model is used with a hybrid deterministic/stochastic particle dispersion model governing the dynamics in the turbulent bulk flow. We used a steady-state RANS k-ϵ turbulence model to simulate the turbulent fluid flow in a neutral atmospheric boundary layer (ABL) using the with inflow boundary conditions by Richards & Hoxey (1993). The turbulence model is solved with the SIMPLE algorithm using the OpenFOAM software. The mean-field characteristic of the turbulent flow in the computational domain is exported and used for the particle model. The particle model is a Lagrangian Langevin-type model, consisting of a system of stochastic differential equations. The particle model was solved using a weakly first order a-stable scheme. We evaluated the deposition model by computing the deposition rate for a range of particle sizes and compared our results with collected experimental wind tunnel data. The numerical experiment was done in a computational domain based on the ABL model by Hargreaves & Wright (2007), a rectangular domain with a logarithmic wind profile. We used a particle source near the inflow boundary with an instantaneously release at the initial time. Results showed disagreement with the experimental data and was only valid for medium sized particles. However, time restrictions led to the analysis being cut short and only a single simulation was conducted. A definite conclusion on the suitability of the method could not be made based solely on this single results. Some uncertainties were identified and discussed for further potential work on the evaluation of the method. However, one conclusion was drawn on the performance of the method. The computational cost was concluded to be too high with the first order particle scheme used and higher order schemes is required for any practical use of the method for atmospheric deposition applications.

Accurate forecasts of wind speed and power production are of great value for wind power producers. In Southwest Iceland, wind power installations are being planned by various entities. This study aims to create optimal wind speed and wind power production forecasts for wind power production in Southwest Iceland by applying statistical post-processing methods to a deterministic HARMONIE-AROME forecast at a single point in space. Three such methods were implemented for a 22 month-long set of forecast-observation samples in 1h resolution: Temporal Smoothing (TS), Observational Distributions on Discrete Intervals (ODDI - a relatively simple classification algorithm) and Quantile Regression Forest (QRF - a relatively complicated Machine Learning Algorithm). Wind power forecasts were derived directly from forecasts of wind speed using an idealized power curve. Four different metrics were given equal weight in the evaluation of the methods: Root Mean Square Error (RMSE), Miss Rate of the 95-percent forecast interval (MR95), Mean Median Forecast Interval Width (MMFIW - a metric to measure the forecast sharpness) and Continuous Ranked Probability Score (CRPS). Of the three methods, TS performed inadequately while ODDI and QRF performed significantly better, and similarly to each other. Both ODDI and QRF predict wind speed and power production slightly more accurately than deterministic AROME in terms of their Root Mean Square Error. In addition to an overall evaluation of all three methods, ODDI and QRF were evaluated conditionally. The results indicate that QRF performs significantly better  than ODDI at forecasting wind speed and wind power at wind speeds above 13 m/s. Else, no strong discrepancies were found between their conditional performance. The results of this study are limited by a relatively scarce data set and correspondingly short time series. The results indicate that applying statistical post-processing methods of varying complexity to deterministic wind speed forecasts is a viable approach to gaining a probabilistic insight into the wind power potential at a given location.

The effect of building geometry on the wind environment of cities is such that it can cause problems like wind danger, discomfort and poor ventilation of airborne pollutants. Computational fluid dynamics (CFD) can play a role in assessing changes in wind environment caused by building projects before realisation at little cost. However, the current state-of-the-art methods, RANS and LES, force a steep trade-off between accuracy and computational cost, and neither method is truly predictive. Time-resolved adaptive direct finite element simulation (DFS) is a method for CFD that is predictive and automatically optimises the mesh for a goal quantity, making it both efficient and accurate. In this thesis, DFS was implemented in FEniCS and used on basic validation cases to provide a proof of concept for the use of this method in the building aerodynamics, on resources freely available to anyone. The results show that the method is accurate to within 10% of the validation data with respect to the goal quantity. Visually, the expected flow features are clearly identifiable. DFS was successfully applied to a relatively complicated building geometry, with a total computation time of about 120 core-hours. We conclude that DFS has significant potential as a method for evaluating urban wind environments. Furthermore, because of its ease of use and lack of parameters, DFS can play an important role in helping architects, designers and students understand the effect of urban geometries on the wind environment. This report provides a basis for further research on DFS for building aerodynamics, as validation on more diverse urban geometries is still necessary. / Effekten av byggnaders form och geometri är så viktig att den kan ge problem för ventilation av t.ex. föroreningar, för energieffektivitet, och för vindfaror med t.ex. hög vindhastihet som kan vara farligt eller skapa obehag. Beräkningsströmningsdynamik (CFD) kan ha en roll i bedömningen av byggnadsprojekt i ett tidigt skede till liten kostnad. Dock är de etablerade och ledande metodikerna, RANS och LES, inte prediktiva och tvingar fram en kompromiss mellan beräkningskosnad och noggrannhet. Vår metodik “Time-resolved adaptive direct finite element simulation” (DFS) är en metod för CFD som är prediktiv och automatiskt optimerar beräkningsnätet (och därmed beräkningskostnaden) för en given målkvantitet, som ger både effektivitet och noggrannhet. I denna avhandling implementerades DFS i FEniCS och användes i grundläggande valideringsfall för att ge ett proof of conceptför användning av denna metod i byggnadsaerodynamik, på resurser som är fritt tillgängliga för alla. Resultaten visar att metoden är korrekt inom 10% av valideringsdata med avseende på målkvantiteten. Visuellt är de förväntade flödesfunktionerna tydligt identifierbara. DFS applicerades framgångsrikt på en relativt komplicerad byggnadsgeometri med en total beräkningstid på cirka 120 kärntimmar, vilket är en försumbar kostnad. Vi drar slutsatsen att DFS har en betydande potential som metod för utvärdering av stadsvindmiljöer. Dessutom, på grund av dess användarvänlighet och frihet från parametrar, kan DFS spela en viktig roll för att hjälpa arkitekter, designers och studenter att förstå effekterna av stadsgeometrier på vindmiljön. Denna rapport ger en grund för vidare forskning om DFS för aerodynamik, eftersom validering av mer olika stadsgeometrier fortfarande är nödvändig.

La physique de la couche limite atmosphérique en domaine océanique est principalement régie par les processus couplés liés au vent, à l’état de mer local, et à des effets de flottabilité. Leur compréhension reste néanmoins parcellaire et leurs descriptions théoriques et stochastiques sont pour le moins lacunaires, lorsqu’elles ne sont tout simplement pas mises à mal par les rares observations. Dans un contexte d’exploitation croissante de la ressource éolienne offshore, la mise en place de méthodes numériques visant à une description plus fine des propriétés turbulentes de cette couche limite sera une étape déterminante dans la réduction des coûts et l’optimisation des structures pour des rendements de récupération d’énergie améliorés. Ainsi, un outil numérique a été mis en place afin d’étudier le couplage entre un écoulement atmosphérique et l’état de mer. Un code Large-Eddy Simulation massivement parallèle pour la simulation des écoulements atmosphériques incompressibles développé par P. Sullivan au National Center for Atmospheric Research est couplé à un code spectral d’états de mer non-linéaires développé au Laboratoire de recherche en Hydrodynamique, Energétique et Environnement Atmosphérique. De nombreuses configurations de vents et d’états de mer sont modélisées. On montre que les lois semi empiriques souvent utilisées pour représenter la distribution verticale de la vitesse moyenne du vent sont une bonne approximation dans les situations où un petit état de mer est soumis à un fort vent. Néanmoins, dans le cas de houles très rapides se propageant dans des zones de faible vent, la création d’un jet de vent par la houle invalide ces lois semi-empiriques. / Modelling the dynamic coupling of ocean-atmosphere systems requires a fundamental and quantitative understanding of the mechanisms governing the windwave interactions: despite numerous studies, our current understanding remains quite incomplete and, in certain conditions, sparse field observations contradict the usual theoretical and stochastic models. Within the context of a growing exploitation of the offshore wind energy and the development of met ocean models, a fine description of this resource is a key issue. Field experiments and numerical modelling have revealed that atmospheric stability and wave effects, including the dynamic sea surface roughness, are two major factors affecting the wind field over oceans. A numerical tool has been implemented in order to study the coupling between an atmospheric flow and the seastate. A massively parallel large-eddy simulation developed by P. Sullivan at the National Center for Atmospheric Research is then coupled to a High-Order Spectral wave model developed at the Hydrodynamics,Energetics & Atmospheric Environment Laboratory in Ecole Centrale de Nantes. Numerous configurations of wind and sea states are investigated. It appears that, under strongly forced wind conditions above a small sea state, the semi-empirical laws referred to as standards in the international guidelines are a good approximation for the vertical profile of the mean wind speed. However, for light winds overlying fast-moving swell, the presence of a wave induced wind jet is observed, invalidating the use of such logarithmic laws.

La capa límit atmosfèrica és l'àrea directament influenciada per la presència de la superfície de la terra i la seva alçada és d'uns centenars de metres a uns pocs quilòmetres. Durant el vespre, el refredament radiatiu estratifica establement l'aire prop del sòl i es forma el que es coneix com a Capa Límit Estable (CLE). D'avui en dia, la CLE és un règim que encara no està prou ben caracteritzat. La turbulència, que no és homogènia ni isòtropa, i la gran importància dels efectes locals com l'orografia, entre d'altres factors, dificulten l'estudi d'aquest règim. Per aquest motiu, la CLE és objecte d'especial atenció, sobretot a l'hora de millorar la seva representació en models tant de temps com de clima.Aquest treball es centra en l'estudi de la CLE mitjançant 3 eines diferents: 1) simulacions explícites de grans remolins (més conegudes com a simulacions LES), per determinar el comportament dels moviments turbulents, on les resolucions són de l'ordre de metres; 2) simulacions mesoscalars, per caracteritzar els efectes locals, on les resolucions són de l'ordre de kilòmetres; 3) anàlisi de les observacions sota aquestes condicions per tal de caracteritzar i entendre millor els fenòmens observats.En primer lloc s'estudia el rang d'estabilitats a on el model LES, que considera la teoria de Kolmogorov per la dissipació de l'energia, funciona correctament. Els resultats del model són realistes tal com mostra la seva comparació amb les mesures de dues campanyes experimentals (SABLES-98 i CASES-99). Per explorar més a fons els resultats LES i per comparar-los amb les mesures s'han utilitzat les Funcions de Distribució de Probabilitat (PDF). Aquests resultats LES són també comparables als obtinguts amb altres models LES, tal com mostra la intercomparació de models LES, més coneguda com a GABLS.Un cop desenvolupades totes les eines necessàries es fa un LES d'un cas més realista, basat en les observacions d'un màxim de vent de capes baixes (més conegut com a Low-Level Jet, LLJ). L'anàlisi combinat dels resultats LES i les mesures permet entendre millor els processos de barreja que tenen lloc a través de la inversió. Finalment, la contribució dels efectes locals s'estudia mitjançant les simulacions mesoscalars, en aquest cas centrades a l'illa de Mallorca. Durant el vespre es veu com les circulacions locals es desenvolupen a les conques (de longitud al voltant de 25km), formant-se, per exemple, vents catabàtics o LLJ com l'estudiat anteriorment. En aquest cas les simulacions es verifiquen amb imatges de satèl·lit NOAA i observacions de les estacions automàtiques de mesures, donant resultats semblants. / The atmospheric boundary layer is the area directly influenced by the presence of the Earth's surface and its height is from hundreds of meters to few kilometres. During the night, the radiative cooling stratifies the layer close to the surface and it forms the Stably-stratified Atmospheric Boundary Layer (SBL). Nowadays, the SBL is a regime not well enough characterized, yet. Turbulence, which is not homogeneous either isotropic, and the great importance of the local effects, like the orography, among other factors, make the SBL be a difficult regime to study. Even so, the SBL is an object of special attention, especially when improving its representation in numerical prediction models or climate models.This work focuses on the study of the SBL through 3 different tools: 1) Large-Eddy Simulations (LES), to determine the turbulent motions, where the resolutions are about 1m; 2) Mesoscale simulations, to characterize the local effects, where resolutions are about 1km; 3) Analysis of the observations under these conditions in order to better characterize and understand the observed phenomena.In first place, it is studied the range of stabilities where the LES model, that considers the Kolmogorov theory for the dissipation of the energy, works correctly. The results are realistic as the comparison with measures from two experimental campaigns (SABLES-98 and CASES-99) shows. To explore the results more thoroughly, and to compare the LES results to the measurements, the Probability Density Functions (PDF) have been used. The LES results are also comparable to the ones obtained with other LES models, as the intercomparison of different LES models show, better known as GABLS.Then, a more realistic case is performed using the LES model, based on observations of a Low-Level Jet (LLJ). The combined inspection of the LES results and the observations allow to better understand the mixing processes that take place through the inversion layer. Finally, the contribution of the local effects is studied through a mesoscale simulation. Here the attention is focused on the Mallorca Island. During the night, the model is able to reproduce the local circulations is a basin of a characteristic size of 25km. The main features obtained previously from the LES of the LLJ are also reproduced by the mesoscale model. These runs are verified with NOAA satellite images and observations from the automatic surface weather stations, giving that the model is able to reproduce realistic results.

The main aim of this research is to improve our knowledge of the sea breeze circulation in the complex coastal environments, where more than one mesoscale circulations occur. Interaction of these circulations with each other and with external factors such as topographical features and large scale winds leads to pronounced changes in the thermodynamic structure of the boundary layer. The variations in sea breeze circulation also have distinct effect on the pollutant transport and dispersion mechanisms in the coastal urban areas. In this research, dynamic and thermodynamic characteristics of the sea breeze circulation and their associated air pollution potential have been investigated by utilizing observational data for two summer periods and numerical modelling techniques. Effect of some external factors such as gradient flow and terrain elevation has also been examined. Observed meteorological and air quality data was obtained from a number of monitoring sites within and around Auckland while Advanced Weather Research & Forecasting (WRF) and ‘The Air Pollution Model’ (TAPM) were employed to simulate meteorology and pollutant dispersion in Auckland. WRF is used to investigate the thermally induced mesoscale circulation while TAPM has been employed to examine the pollutant dispersion in the region. Both models were validated against observed data from six different sites within Auckland. Validation results of WRF and TAPM are also compared with surface meteorology. Validation and inter-comparison of the two models show that WRF performed better than TAPM for all the surface meteorology variables. WRF showed a positive bias in predicted winds speed and relative humidity and a cold bias in the near surface Temperature. TAPM on the other hand under-predicted surface winds, while near surface temperature and relative humidity are similar to WRF. Results show that the sea breeze occurred around 20% of the two summer periods of 2006 and 2007. Both observed data analysis and the numerical modelling results confirmed the existence of two thermally induced systems in the Auckland region. Bay breezes are initiated in the morning hours (0800 – 1000 hours) from small bodies of water (Manukau, Waitemata, and Kaipara Harbour, and along the Hauraki Gulf coastline), followed by mature sea breezes from the main bodies of water (Tasman Sea and larger Hauraki Gulf area) in the late morning. The cessation of sea breezes started after 1600 hours. Frequency of sea breeze days was the highest under coast-parallel gradient winds (southeast and northwest), with speeds < 6 m s-1. The predicted depth of the sea breeze inflow ranged between 200 and 600 m, while the depth of the return flow was in the range of 200 – 500 m. Sensible heat flux is an important control in the development of sea breeze over the region. Coastal mountain ranges helped early onset of the sea breeze, but also inhibited inland propagation. Strong jet-like westerly winds along the coastline near the Manukau Harbour are due partly to the narrow opening at the Manukau Head, reduced friction over the harbour water, and divergence of wind due to coastline shape. Gradient winds significantly affect the evolution of the sea breeze and modify many of its dynamics, such as the sea breeze inflow layer, return flow, inland penetration, sea breeze head, etc. Under northerly gradient flow northeast sea breeze lasts longer while under southerly gradient flow cessation of the westerly sea breeze was delayed. Over both east and west coasts, WRF predicted anticlockwise rotation, especially under easterly gradient wind conditions. However, inland stations near Manukau Harbour show partial and complete clockwise rotation, which is primarily due to orographic features of the region. The diurnal rotation of the sea breeze system may contribute to recirculation of pollutants in the morning hours under coast-parallel gradient wind conditions. Pollutants that are emitted during morning peak traffic hours and advected towards Manukau Harbour by the remnants of the land breeze may be returned by bay breezes in the mid morning hours. Mixed layer height over land before arrival of the sea breeze also varied a lot and ranged between 600 to 1400 m. A convective internal boundary layer (CIBL) forms in the surface layer after arrival of the sea breeze. The CIBL under coast parallel gradient winds was relatively shallow (200 – 400 m), while under coasts-normal gradient winds (southwest and northeast), the predicted depth was in the range of 400 to 500 m. However, the inland extent of the CIBL was greater under coast-normal winds, especially under south-westerly gradient winds. The ground level concentration of air pollutants thus can be increased during sea breeze inflow over the region. Both bay breeze and mature sea breeze contribute towards development, extent and strength of the sea breeze convergence zones (SBCZs). Gradient winds and terrain play an important role in the position and strength of SBCZs. Under strong south-westerly gradient flow, a SBCZ is formed along the eastern coastline, while under north-easterly gradient winds a SBCZ is formed along the west coastline. During coast-parallel gradient winds the SBCZ is formed in the middle of landmass, and is then gradually displaced eastward or westward depending on the balance between large scale PGF and surface friction effect. In addition to SBCZs, terrain and coastline-induced convergences were also evident. Higher ground level concentrations of pollutants are expected under coast-normal gradient winds, when SBCZs are formed in the middle of the land mass and the wind speed of the sea breeze inflow and the sea breeze front is relatively low. This may increase pollution concentration, especially in the evening hours, to unacceptable levels. Results of this research suggest that given the size, synoptic meteorology and specific geography of the region, significant recirculation of pollutants is not likely to happen to contribute to next day’s pollution. The pollutant concentration may increase in the SBCZs, but their ability to recirculate the pollutants requires more extensive research. A closed sea breeze circulation cell is unlikely to form in this region due to topographical influences and a strong gradient wind effect. The pollutant plume is expected to be advected in the return flow over the peaks of higher terrain and via the top of the convergence zones, but its remixing in the onshore flow is subject to many factors such as gradient wind speed and direction, direction of the return flow and nature (size and state) of the pollutant. In appropriate conditions, pollution levels may reach to unhealthy levels under coast-parallel gradient wind condition.