Detection and Modeling of Boundary Layer Height

This thesis tackles the problem of the automated detection of the atmospheric boundary layer (BL) height, h, from aerosol lidar/ceilometer observations. A new method, the Bayesian Selective Method (BSM), is presented. It implements a Bayesian statistical inference procedure which combines in an statistically optimal way different sources of information. Firstly atmospheric stratification boundaries are located from discontinuities in the ceilometer back-scattered signal. The BSM then identifies the discontinuity edge that has the highest probability to effectively mark the BL height. Information from the contemporaneus physical boundary layer model simulations and a climatological dataset of BL height evolution are combined in the assimilation framework to assist this choice. The BSM algorithm has been tested for four months of continuous ceilometer measurements collected during the BASE:ALFA project and is shown to realistically diagnose the BL depth evolution in many different weather conditions. Then the BASE:ALFA dataset is used to investigate the boundary layer structure in stable conditions. Functions from the Obukhov similarity theory are used as regression curves to fit observed velocity and temperature profiles in the lower half of the stable boundary layer. Surface fluxes of heat and momentum are best-fitting parameters in this exercise and are compared with what measured by a sonic anemometer. The comparison shows remarkable discrepancies, more evident in cases for which the bulk Richardson number turns out to be quite large. This analysis supports earlier results, that surface turbulent fluxes are not the appropriate scaling parameters for profiles of mean quantities in very stable conditions. One of the practical consequences is that boundary layer height diagnostic formulations which mainly rely on surface fluxes are in disagreement to what obtained by inspecting co-located radiosounding profiles.

Concentrations of organic carbon, nitrogen and phosphorus on vertical profiles in waters of the tropical North Atlantic in March 2002

In the literature, an inconsistency exists between estimates of biotically-effected carbon export inferred from large-scale geochemical studies (Jenkins 1982; 47 gC m-2 a-1) and local measurements of turbulent nutrient supply (Lewis et al. 1986; 4 gC m-2 a-1) in the eastern subtropical North Atlantic. Nutrient supply to the upper ocean by turbulent mixing is reexamined using local standard oceanographic measurements and high-resolution vertical profiles of nutrients averaged over a large region directly comparable to that investigated by Jenkins (1982). Turbulent fluxes induced by internal waves and salt fingering, respectively, are separated according to Gregg (1989) and Zhang et al. (1998). Nutrient transport into the nutrient-consuming surface layer by salt fingering is more than fivefold higher than transport due to internal-wave induced turbulence. Still, this cannot resolve the above- mentioned apparent inconsistency, even if additional physical transport mechanisms such as eddy pumping, advection and horizontal diffusion are accounted for. Estimated nitrate fluxes due to vertical turbulent diffusion are 0.05-0.15 mol m-2 a-1, corresponding to 4-11 gC m-2 a-1. Observed NO3/PO4 turbulent flux ratios of up to 23 are interpreted as the imprint of N2 fixation.

Implementation of non-local boundary layer schemes in the Regional Atmospheric Modeling System and its impact on simulated mesoscale circulations

This paper proposes the implementation of different non-local Planetary Boundary Layer schemes within the Regional Atmospheric Modeling System (RAMS) model. The two selected PBL parameterizations are the Medium-Range Forecast (MRF) PBL and its updated version, known as the Yonsei University (YSU) PBL. YSU is a first-order scheme that uses non-local eddy diffusivity coefficients to compute turbulent fluxes. It is based on the MRF, and improves it with an explicit treatment of the entrainment. With the aim of evaluating the RAMS results for these PBL parameterizations, a series of numerical simulations have been performed and contrasted with the results obtained using the Mellor and Yamada (MY) scheme, also widely used, and the standard PBL scheme in the RAMS model. The numerical study carried out here is focused on mesoscale circulation events during the summer, as these meteorological situations dominate this season of the year in the Western Mediterranean coast. In addition, the sensitivity of these PBL parameterizations to the initial soil moisture content is also evaluated. The results show a warmer and moister PBL for the YSU scheme compared to both MRF and MY. The model presents as well a tendency to overestimate the observed temperature and to underestimate the observed humidity, considering all PBL schemes and a low initial soil moisture content. In addition, the bias between the model and the observations is significantly reduced moistening the initial soil moisture of the corresponding run. Thus, varying this parameter has a positive effect and improves the simulated results in relation to the observations. However, there is still a significant overestimation of the wind speed over flatter terrain, independently of the PBL scheme and the initial soil moisture used, even though a different degree of accuracy is reproduced by RAMS taking into account the different sensitivity tests.

First-year and second-year snow properties on sea ice in the Weddell Sea during POLARSTERN cruise ANT-XXII/2 (ISPOL) during the drift

Observations of snow properties, superimposed ice, and atmospheric heat fluxes have been performed on first-year and second-year sea ice in the western Weddell Sea, Antarctica. Snow in this region is particular as it does usually survive summer ablation. Measurements were performed during Ice Station Polarstern (ISPOL), a 5-week drift station of the German icebreaker RV Polarstern. Net heat flux to the snowpack was 8 W/m**2, causing only 0.1 to 0.2 m of thinning of both snow cover types, thinner first-year and thicker second-year snow. Snow thinning was dominated by compaction and evaporation, whereas melt was of minor importance and occurred only internally at or close to the surface. Characteristic differences between snow on first-year and second-year ice were found in snow thickness, temperature, and stratigraphy. Snow on second-year ice was thicker, colder, denser, and more layered than on first-year ice. Metamorphism and ablation, and thus mass balance, were similar between both regimes, because they depend more on surface heat fluxes and less on underground properties. Ice freeboard was mostly negative, but flooding occurred mainly on first-year ice. Snow and ice interface temperature did not reach the melting point during the observation period. Nevertheless, formation of discontinuous superimposed ice was observed. Color tracer experiments suggest considerable meltwater percolation within the snow, despite below-melting temperatures of lower layers. Strong meridional gradients of snow and sea-ice properties were found in this region. They suggest similar gradients in atmospheric and oceanographic conditions and implicate their importance for melt processes and the location of the summer ice edge.

Global surface net-radiation at 5 km from MODIS Terra

Reliable and fine resolution estimates of surface net-radiation are required for estimating latent and sensible heat fluxes between the land surface and the atmosphere. However, currently, fine resolution estimates of net-radiation are not available and consequently it is challenging to develop multi-year estimates of evapotranspiration at scales that can capture land surface heterogeneity and are relevant for policy and decision-making. We developed and evaluated a global net-radiation product at 5 km and 8-day resolution by combining mutually consistent atmosphere and land data from the Moderate Resolution Imaging Spectroradiometer (MODIS) on board Terra. Comparison with net-radiation measurements from 154 globally distributed sites (414 site-years) from the FLUXNET and Surface Radiation budget network (SURFRAD) showed that the net-radiation product agreed well with measurements across seasons and climate types in the extratropics (Wilmott’s index ranged from 0.74 for boreal to 0.63 for Mediterranean sites). Mean absolute deviation between the MODIS and measured net-radiation ranged from 38.0 ± 1.8 W∙m−2 in boreal to 72.0 ± 4.1 W∙m−2 in the tropical climates. The mean bias was small and constituted only 11%, 0.7%, 8.4%, 4.2%, 13.3%, and 5.4% of the mean absolute error in daytime net-radiation in boreal, Mediterranean, temperate-continental, temperate, semi-arid, and tropical climate, respectively. To assess the accuracy of the broader spatiotemporal patterns, we upscaled error-quantified MODIS net-radiation and compared it with the net-radiation estimates from the coarse spatial (1° × 1°) but high temporal resolution gridded net-radiation product from the Clouds and Earth’s Radiant Energy System (CERES). Our estimates agreed closely with the net-radiation estimates from the CERES. Difference between the two was less than 10 W•m−2 in 94% of the total land area. MODIS net-radiation product will be a valuable resource for the science community studying turbulent fluxes and energy budget at the Earth’s surface.

Influence of the sea state on Mediterranean heavy precipitation: a case-study from HyMeX SOP1

Sea state can influence the turbulent air–sea exchanges, especially the momentum flux, by modifying the sea-surface roughness. The high-resolution non-hydrostatic convection-permitting model MESO-NH is used here to investigate the impact of a more realistic representation of the waves on heavy precipitation during the Intense Observation Period (IOP) 16a of the first HyMeX Special Observation Period (SOP1). Several quasi-stationary mesoscale convective systems developed over the western Mediterranean region, two of them over the sea, and resulted in heavy precipitation on the French and Italian coasts on 26 October 2012. Three different bulk parametrizations are tested in this study: a reference case (NOWAV) without any wave effect, a parametrization taking into account theoretical wave effects (WAV) and a last one with realistic wave characteristics from the MFWAM analyses (WAM). Using a realistic wave representation in WAM significantly increases the roughness length and the friction velocity with respect to NOWAV and WAV. The three MESO-NH sensitivity experiments of the IOP16a show that this surface-roughness increase in WAM generates higher momentum fluxes and directly impacts the low-level dynamics of the atmosphere, with a slowdown of the 10 m wind, when and where the wind speed exceeds 10 m s−1 and the sea state differs from the idealized one. The turbulent heat fluxes are not significantly influenced by the waves, these fluxes being controlled by the moisture content rather than by the wind speed in the simulations. Although the convective activity is globally well reproduced by all the simulations, the difference in the low-level dynamics of the atmosphere influences the localization of the simulated heavy precipitation. Objective evaluation of the daily rainfall amount and of the 10 m wind speed against the observations confirms the positive impact of the realistic wave representation on this simulation of heavy precipitation.

Scaling arguments for the fluxes in turbulent miscible fountains

For established axisymmetric turbulent miscible Boussinesq fountains in quiescent uniform environments, expressions are developed for the fluxes of volume, momentum and buoyancy at the outflow from the fountain: the outflow referring to the counterflow at the horizontal plane of the source. The fluxes are expressed in terms of the fountain source conditions and two dimensionless functions of the source Froude number, Fr0: a radial function (relating a horizontal scale of the outflow to the source radius) and a volume flux function (relating the outflow and source volume fluxes). The forms taken by these two functions at low Fr0 and high Fr0 are deduced, thereby providing the outflow fluxes and outflow Froude number solely in terms of the source conditions. For high Fr0, the outflow Froude number, Frout, is shown to be invariant, indicating (by analogy with plumes for which the 'far-field' Froude number is invariant with source Froude number) that the outflow may be regarded as 'far-field' since the fluxes within the fountain have adjusted to attain a balance which is independent of the source conditions. Based on Frout, the fluxes in the plume that forms beyond the fountain outflow are deduced. Finally, from the results of previously published studies, we show that the scalings deduced for fountains are valid for 0.0025 ≲ Fr0 ≲ 1.0 for low Fr0 and Fr0≳ 3.0 for high Fr0. © 2014 Cambridge University Press.

Determination of turbulent heat fluxes using a large aperture scintillometer over undulating mixed agricultural terrain

Scintillometry is an established technique for determining large areal average sensible heat fluxes. The scintillometer measurement is related to sensible heat flux via Monin–Obukhov similarity theory, which was developed for ideal homogeneous land surfaces. In this study it is shown that judicious application of scintillometry over heterogeneous mixed agriculture on undulating topography yields valid results when compared to eddy covariance (EC). A large aperture scintillometer (LAS) over a 2.4 km path was compared with four EC stations measuring sensible (H) and latent (LvE) heat fluxes over different vegetation (cereals and grass) which when aggregated were representative of the LAS source area. The partitioning of available energy into H and LvE varied strongly for different vegetation types, with H varying by a factor of three between senesced winter wheat and grass pasture. The LAS derived H agrees (one-to-one within the experimental uncertainty) with H aggregated from EC with a high coefficient of determination of 0.94. Chronological analysis shows individual fields may have a varying contribution to the areal average sensible heat flux on short (weekly) time scales due to phenological development and changing soil moisture conditions. Using spatially aggregated measurements of net radiation and soil heat flux with H from the LAS, the areal averaged latent heat flux (LvELAS) was calculated as the residual of the surface energy balance. The regression of LvELAS against aggregated LvE from the EC stations has a slope of 0.94, close to ideal, and demonstrates that this is an accurate method for the landscape-scale estimation of evaporation over heterogeneous complex topography.

Case study modeling of turbulent and mesoscale fluxes over the BOREAS region

Results from aircraft and surface observations provided evidence for the existence of mesoscale circulations over the Boreal Ecosystem-Atmosphere Study (BOREAS) domain. Using an integrated approach that included the use of analytical modeling, numerical modeling, and data analysis, we have found that there are substantial contributions to the total budgets of heat over the BOREAS domain generated by mesoscale circulations. This effect is largest when the synoptic flow is relatively weak, yet it is present under less favorable conditions, as shown by the case study presented here. While further analysis is warranted to document this effect, the existence of mesoscale flow is not surprising, since it is related to the presence of landscape patches, including lakes, which are of a size on the order of the local Rossby radius and which have spatial differences in maximum sensible heat flux of about 300 W m−2. We have also analyzed the vertical temperature profile simulated in our case study as well as high-resolution soundings and we have found vertical profiles of temperature change above the boundary layer height, which we attribute in part to mesoscale contributions. Our conclusion is that in regions with organized landscapes, such as BOREAS, even with relatively strong synoptic winds, dynamical scaling criteria should be used to assess whether mesoscale effects should be parameterized or explicitly resolved in numerical models of the atmosphere.

On the parameterization of turbulent surface fluxes over heterogeneous sea ice surfaces

Turbulent surface fluxes of momentum and sensible and latent heat as well as surface temperature, air temperature, air humidity, and wind speed were measured by the German Falcon research aircraft over the marginal ice zone (MIZ) of the northern Baltic Sea and the Fram Strait. Applying the bulk formulas and the stability functions to the measurements, the roughness lengths for momentum z0, sensible heat zT, and latent heat zq were calculated. As mean values over a wide range of sea ice conditions, we obtain z0 = 5 � 10�4 m, zT = 1 � 10�8 m, and zq = 1 � 10�7 m. These correspond to the following mean values (± standard deviations) of neutral transfer coefficients reduced to 10 m height, CDN10 = (1.9 ± 0.8) � 10�3, CHN10 = (0.9 ± 0.3) � 10�3, and CEN10 = (1.0 ± 0.2) � 10�3. An average ratio of z0/zT � 104 was observed over the range of 10�6 m < z0 < 10�2 m and differs from previously published results over compact sea ice (10�1 < z0/zT < 103). Other observational results over heterogeneous sea ice do not exist. However, our z0/zT ratio approximately agrees with observations over heterogeneous land surfaces. Flux parameterizations based on commonly used roughness lengths ratios (z0 = zT = zq) overestimate the surface heat fluxes compared to our measurements by more than 100%.