The atmospheric stable boundary layer: Data analysis and comparison with similarity theories

The objective of this dissertation is to study the structure and behavior of the Atmospheric Boundary Layer (ABL) in stable conditions. This type of boundary layer is not completely well understood yet, although it is very important for many practical uses, from forecast modeling to atmospheric dispersion of pollutants. We analyzed data from the SABLES98 experiment (Stable Atmospheric Boundary Layer Experiment in Spain, 1998), and compared the behaviour of this data using Monin-Obukhov's similarity functions for wind speed and potential temperature. Analyzing the vertical profiles of various variables, in particular the thermal and momentum fluxes, we identified two main contrasting structures describing two different states of the SBL, a traditional and an upside-down boundary layer. We were able to determine the main features of these two states of the boundary layer in terms of vertical profiles of potential temperature and wind speed, turbulent kinetic energy and fluxes, studying the time series and vertical structure of the atmosphere for two separate nights in the dataset, taken as case studies. We also developed an original classification of the SBL, in order to separate the influence of mesoscale phenomena from turbulent behavior, using as parameters the wind speed and the gradient Richardson number. We then compared these two formulations, using the SABLES98 dataset, verifying their validity for different variables (wind speed and potential temperature, and their difference, at different heights) and with different stability parameters (zita or Rg). Despite these two classifications having completely different physical origins, we were able to find some common behavior, in particular under weak stability conditions.

Turbulence models in the atmospheric boundary layer under convective conditions

In this work a modelization of the turbulence in the atmospheric boundary layer, under convective condition, is made. For this aim, the equations that describe the atmospheric motion are expressed through Reynolds averages and, then, they need closures. This work consists in modifying the TKE-l closure used in the BOLAM (Bologna Limited Area Model) forecast model. In particular, the single column model extracted from BOLAM is used, which is modified to obtain other three different closure schemes: a non-local term is added to the flux- gradient relations used to close the second order moments present in the evolution equation of the turbulent kinetic energy, so that the flux-gradient relations become more suitable for simulating an unstable boundary layer. Furthermore, a comparison among the results obtained from the single column model, the ones obtained from the three new schemes and the observations provided by the known case in literature ”GABLS2” is made.

The wind profile in the atmospheric stable boundary layer over complex terrain and heterogeneous surfaces: limitations of local similarity theory

L'indagine ha riguardato il profilo del vento nei primi 30 metri dello strato limite atmosferico stabile nell'ambito della teoria di similarità locale. Ad oggi, diversi esperimenti hanno confermato la validità della teoria per strati-limite su terreni livellati e superfici omogenee. Tali condizioni ideali sono però infrequenti nella realtà ed è perciò importante capire quali siano i limiti della similarità locale per strati-limite su terreni complessi e superfici disomogenee. Entrambe le condizioni sono presenti a Ny-Alesund (Svalbard, Norvegia) dove il Consiglio Nazionale delle Ricerche (CNR), nel 2009, ha installato una torre di 30 m, la Amudsen-Nobile Climate Change Tower (CCT), per lo studio dello strato-limite artico. Il lavoro di tesi ha riguardato misure di vento e turbolenza acquisite sulla CCT da maggio 2012 a maggio 2014. Il confronto tra le velocità del vento misurate dagli anemometri installati sulla CCT, ha rivelato criticità nel dato sonico manifestatesi con sovrastime sistematiche e maggiore erraticità rispetto alle misure provenienti dagli anemometri a elica. Un test condotto fra diversi metodi per il calcolo dei gradienti verticali della velocità del vento ha rivelato scarsa sensibilità dei risultati ottenuti al particolare metodo utilizzato. Lo studio ha riguardato i gradienti verticali adimensionali della velocità del vento nei primi 30-m dello strato limite stabile. Deviazioni significative tra i tra le osservazioni e i valori predetti dalla similarità locale sono state osservate in particolare per i livelli più distanti dal suolo e per valori crescenti del parametro di stabilità z/L (L, lunghezza di Obukhov locale). In particolare, si sono osservati gradienti adimensionali inferiori a quelli predetti dalle più usate relazioni di flusso-gradiente. Tali deviazioni, presenti perlopiù per z/L>0.1, sono state associate ad un effetto di accentuazione della turbolenza da parte delle irregolarità del terreno. Per condizioni meno stabili, z/L<0.1, scarti positivi tra i gradienti osservati e quelli attesi sono stati attribuiti alla formazione di strati limite interni in condizioni di vento dal mare verso la costa. Sono stati proposti diversi metodi per la stima dell'effetto della self-correlazione nella derivazione delle relazioni di flusso-gradiente, dovuta alla condivisione della variabile u*. La formula per il coefficiente lineare di self correlazione e le sue distribuzioni di probabilità empiriche sono state derivate e hanno permesso di stimare il livello di self-correlazione presente nel dataset considerato.

Stabilization of subcritical bypass transition in the leading-edge boundary layer by suction

This work investigates the subcritical spatial transition in the swept Hiemenz boundary layer by means of direct numerical simulations (DNS). A pair of steady co-rotating vortices located at the attachment line is enforced as a primary disturbance leading to streaks which are stable. A small secondary, time-dependent disturbance interacts with these streaks such that instability and breakdown to turbulence may occur. The instability only occurs for a certain band of secondary disturbance frequencies. Positive secondary instability growth rates could be observed for Reynolds numbers as low as , whereas the linear critical Reynolds number is. Uniform wall suction is shown to stabilise this transition mechanism, analogously to results from linear stability theory. The effects of suction on the formation of primary streaks and on the secondary growth rate are decoupled. For streaks of different suction whose amplitude is held constant by adjusting the Reynolds number, the suction is shown to increase the growth rate of the secondary instability. The stabilising influence of wall suction consists in decreasing the streak amplitude only. Depending on the Reynolds number and the suction strength, breakdown may either occur locally and may be convected along the far-field streamlines, or occur globally and cover broad regions in the downstream direction.

Study of mineral dust entrainment in the planetary boundary layer by lidar depolarisation technique

Measurements on 27 June 2011 were performed over the Southern Iberian Peninsula at Granada EARLINET station, using active and passive remote sensing and airborne and surface in-situ data in order to study the entrainment processes between aerosols in the free troposphere and those in the planetary boundary layer (PBL). To this aim the temporal evolution of the lidar depolarisation, backscatter-related Angström exponent and potential temperature profiles were used in combination with the PBL contribution to the aerosol optical depth (AOD). Our results show that the mineral dust entrainment in the PBL was caused by the convective processes which ‘trapped’ the lofted mineral dust layer, distributing the mineral dust particles within the PBL. The temporal evolution of ground-based in-situ data evidenced the impact of this process at surface level. Finally, the amount of mineral dust in the atmospheric column available to be dispersed into the PBL was estimated by means of POLIPHON (Polarizing Lidar Photometer Networking). The dust mass concentration derived from POLIPHON was compared with the coarse-mode mass concentration retrieved with airborne in-situ measurements. Comparison shows differences below 50 µg/m³ (30% relative difference) indicating a relative good agreement between both techniques.

Observations of the Sensitivity of Beam Attenuation to Particle Size in a Coastal Bottom Boundary Layer

The goal of this study was to test the hypothesis that the aggregated state of natural marine particles constrains the sensitivity of optical beam attenuation to particle size. An instrumented bottom tripod was deployed at the 12-m node of the Martha's Vineyard Coastal Observatory to monitor particle size distributions, particle size-versus-settling-velocity relationships, and the beam attenuation coefficient (c(p)) in the bottom boundary layer in September 2007. An automated in situ filtration system on the tripod collected 24 direct estimates of suspended particulate mass (SPM) during each of five deployments. On a sampling interval of 5 min, data from a Sequoia Scientific LISST 100x Type B were merged with data from a digital floc camera to generate suspended particle volume size distributions spanning diameters from approximately 2 mu m to 4 cm. Diameter-dependent densities were calculated from size-versus-settling-velocity data, allowing conversion of the volume size distributions to mass distributions, which were used to estimate SPM every 5 min. Estimated SPM and measured c(p) from the LISST 100x were linearly correlated throughout the experiment, despite wide variations in particle size. The slope of the line, which is the ratio of c(p) to SPM, was 0.22 g m(-2). Individual estimates of c(p):SPM were between 0.2 and 0.4 g m(-2) for volumetric median particle diameters ranging from 10 to 150 mu m. The wide range of values in c(p):SPM in the literature likely results from three factors capable of producing factor-of-two variability in the ratio: particle size, particle composition, and the finite acceptance angle of commercial beam-transmissometers.

Atmospheric nitrate concentrations and isotope composition for samples collected in the marine boundary layer in 2006 and 2007

The comprehensive isotopic composition of atmospheric nitrate (i.e., the simultaneous measurement of all its stable isotope ratios: 15N/14N, 17O/16O and 18O/16O) has been determined for aerosol samples collected in the marine boundary layer (MBL) over the Atlantic Ocean from 65°S (Weddell Sea) to 79°N (Svalbard), along a ship-borne latitudinal transect. In nonpolar areas, the d15N of nitrate mostly deriving from anthropogenically emitted NOx is found to be significantly different (from 0 to 6 per mil) from nitrate sampled in locations influenced by natural NOx sources (-4 ± 2) per mil. The effects on d15N(NO3-) of different NOx sources and nitrate removal processes associated with its atmospheric transport are discussed. Measurements of the oxygen isotope anomaly (D17O = d17O - 0.52 × d18O) of nitrate suggest that nocturnal processes involving the nitrate radical play a major role in terms of NOx sinks. Different D17O between aerosol size fractions indicate different proportions between nitrate formation pathways as a function of the size and composition of the particles. Extremely low d15N values (down to -40 per mil) are found in air masses exposed to snow-covered areas, showing that snowpack emissions of NOx from upwind regions can have a significant impact on the local surface budget of reactive nitrogen, in conjunction with interactions with active halogen chemistry. The implications of the results are discussed in light of the potential use of the stable isotopic composition of nitrate to infer atmospherically relevant information from nitrate preserved in ice cores.