Investigating microscale patchiness of motile microbes driven by the interaction of turbulence and gyrotaxis in a 3D simulated convective mixed layer.

Microbes play a primary role in aquatic ecosystems and biogeochemical cycles. Spatial patchiness is a critical factor underlying these activities, influencing biological productivity, nutrient cycling and dynamics across trophic levels. Incorporating spatial dynamics into microbial models is a long-standing challenge, particularly where small-scale turbulence is involved. Here, we combine a fully 3D direct numerical simulation of convective mixed layer turbulence, with an individual-based microbial model to test the key hypothesis that the coupling of gyrotactic motility and turbulence drives intense microscale patchiness. The fluid model simulates turbulent convection caused by heat loss through the fluid surface, for example during the night, during autumnal or winter cooling or during a cold-air outbreak. We find that under such conditions, turbulence-driven patchiness is depth-structured and requires high motility: Near the fluid surface, intense convective turbulence overpowers motility, homogenising motile and non-motile microbes approximately equally. At greater depth, in conditions analogous to a thermocline, highly motile microbes can be over twice as patch-concentrated as non-motile microbes, and can substantially amplify their swimming velocity by efficiently exploiting fast-moving packets of fluid. Our results substantiate the predictions of earlier studies, and demonstrate that turbulence-driven patchiness is not a ubiquitous consequence of motility but rather a delicate balance of motility and turbulent intensity.

Evidence for an unidentified non-photochemical ground-level source of formaldehyde in the Po Valley with potential implications for ozone production

Ozone concentrations in the Po Valley of northern Italy often exceed international regulations. As both a source of radicals and an intermediate in the oxidation of most volatile organic compounds (VOCs), formaldehyde (HCHO) is a useful tracer for the oxidative processing of hydrocarbons that leads to ozone production. We investigate the sources of HCHO in the Po Valley using vertical profile measurements acquired from the airship Zeppelin NT over an agricultural region during the PEGASOS 2012 campaign. Using a 1-D model, the total VOC oxidation rate is examined and discussed in the context of formaldehyde and ozone production in the early morning. While model and measurement discrepancies in OH reactivity are small (on average 3.4 ± 13%), HCHO concentrations are underestimated by as much as 1.5 ppb (45%) in the convective mixed layer. A similar underestimate in HCHO was seen in the 2002–2003 FORMAT Po Valley measurements, though the additional source of HCHO was not identified. Oxidation of unmeasured VOC precursors cannot explain the missing HCHO source, as measured OH reactivity is explained by measured VOCs and their calculated oxidation products. We conclude that local direct emissions from agricultural land are the most likely source of missing HCHO. Model calculations demonstrate that radicals from degradation of this non-photochemical HCHO source increase model ozone production rates by as much as 0.6 ppb h−1 (12%) before noon.

Horizontal turbulent diffusion in a convective mixed layer

A large eddy simulation (LES) is used to estimate a reliable horizontal turbulent diffusion coefficient, $\def \xmlpi #1{}\def \mathsfbi #1{\boldsymbol {\mathsf {#1}}}\let \le =\leqslant \let \leq =\leqslant \let \ge =\geqslant \let \geq =\geqslant \def \Pr {\mathit {Pr}}\def \Fr {\mathit {Fr}}\def \Rey {\mathit {Re}}K_{{h}}$, in a convective mixed layer (CML). The introduction of a passive scalar field with a fixed horizontal gradient at a given time enables $K_{{h}}$ estimation as a function of height, based on the simulated turbulent horizontal scalar flux. Here $K_{{h}}$ is found to be of the order of $100\ {\mathrm{m}}^2\ {\mathrm{s}}^{-1}$ for a typical terrestrial atmospheric CML. It is shown to scale by the product of the CML convective velocity, $w_{*}$, and its depth, $h$. Here $K_{{h}}$ is characterized by a vertical profile in the CML: it is large near both the bottom and top of the CML, where horizontal flows associated with convection are large. The equation pertaining to the temporal rate of change of a horizontal scalar flux suggests that $K_{{h}}$ is determined by a balance between production and pressure correlation at a fully developed stage. Pressure correlation near the bottom of the CML is localized in convergence zones near the boundaries of convective cells and becomes large within an eddy turnover time, $h/w_{*}$, after the introduction of the passive scalar field.

Evidence for an unidentified ground-level source of formaldehyde in the Po Valley with potential implications for ozone production

Ozone concentrations in the Po Valley of Northern Italy often exceed international regulations. As both a source of radicals and an intermediate in the oxidation of most volatile organic compounds (VOCs), formaldehyde (HCHO) is a useful tracer for the oxidative processing of hydrocarbons that leads to ozone production. We investigate the sources of HCHO in the Po Valley using vertical profile measurements acquired from the airship Zeppelin NT over an agricultural region during the PEGASOS 2012 campaign. Using a 1-D model, the total VOC oxidation rate is examined and discussed in the context of formaldehyde and ozone production in the early morning. While model and measurement discrepancies in OH reactivity are small (on average 3.4±11%), HCHO concentrations are underestimated by as much as 1.5 ppb (45%) in the convective mixed layer. A similar underestimate in HCHO was seen in the 2002–2003 FORMAT Po-Valley measurements, though the additional source of HCHO was not identified. Oxidation of unmeasured VOC precursors cannot explain the missing HCHO source, as measured OH reactivity is explained by measured VOCs and their calculated oxidation products. We conclude that local direct emissions from agricultural land are the most likely source of missing HCHO. Model calculations demonstrate that radicals from degradation of this non-photochemical HCHO source increase model ozone production rates by as much as 0.7 ppb h−1 (10%) before noon.

Factors controlling temporal variability of near-ground atmospheric ²²²Rn concentration over central Europe

Concentration of radon (222Rn) in the near-ground atmosphere has been measured quasi-continuously from January 2005 to December 2009 at two continental sites in Europe: Heidelberg (south-west Germany) and Krakow (southern Poland). The atmosphere was sampled at ca. 30 and 20 m above the local ground. Both stations were equipped with identical instruments. Regular observations of 222Rn were supplemented by measurements of surface fluxes of this gas in the Krakow urban area, using two different approaches. The measured concentrations of 222Rn varied at both sites in a wide range, from less than 2.0 Bq m−3 to approximately 40 Bq m−3 in Krakow and 35 Bq m−3 in Heidelberg. The mean 222Rn content in Krakow, when averaged over the entire observation period, was 30% higher than in Heidelberg (5.86 ± 0.09 and 4.50 ± 0.07 Bq m−3, respectively). Distinct seasonality of 222Rn signal is visible in the obtained time series of 222Rn concentration, with higher values recorded generally during late summer and autumn. The surface 222Rn fluxes measured in Krakow also revealed a distinct seasonality, with broad maximum observed during summer and early autumn and minimum during the winter. When averaged over a 5-year observation period, the night-time surface 222Rn flux was equal to 46.8 ± 2.4 Bq m−2 h−1. Although the atmospheric 222Rn levels at Heidelberg and Krakow appeared to be controlled primarily by local factors, it was possible to evaluate the "continental effect" in atmospheric 222Rn content between both sites, related to gradual build-up of 222Rn concentration in the air masses travelling between Heidelberg and Krakow. The mean value of this build-up was equal to 0.78 ± 0.12 Bq m−3. The measured minimum 222Rn concentrations at both sites and the difference between them was interpreted in the framework of a simple box model coupled with HYSPLIT (Hybrid Single Particle Lagrangian Integrated Trajectory) analysis of air mass trajectories. The best fit of experimental data was obtained for the mean 222Rn flux over the European continent equal to 52 Bq m−2 h−1, the mean transport velocity of the air masses within the convective mixed layer of the planetary boundary layer (PBL) on their route from the Atlantic coast to Heidelberg and Krakow equal to 3.5 m s−1, the mean rate constant of 222Rn removal across the top of the PBL equal to the 222Rn decay constant and the mean height of the convective mixed layer equal to 1600 m.

Aircraft Evaluation of Ground-Based Raman Lidar Water Vapor Turbulence Profiles in Convective Mixed Layers

High temporal and vertical resolution water vapor measurements by Raman and differential absorption lidar systems have been used to characterize the turbulent fluctuations in the water vapor mixing ratio field in convective mixed layers. Since daytime Raman lidar measurements are inherently noisy (due to solar background and weak signal strengths), the analysis approach needs to quantify and remove the contribution of the instrument noise in order to derive the desired atmospheric water vapor mixing ratio variance and skewness profiles. This is done using the approach outlined by Lenschow et al.; however, an intercomparison with in situ observations was not performed.Water vapor measurements were made by a diode laser hygrometer flown on a Twin Otter aircraft during the Routine Atmospheric Radiation Measurement (ARM) Program Aerial Facility Clouds with Low Optical Water Depths Optical Radiative Observations (RACORO) field campaign over the ARM Southern Great Plains (SGP) site in 2009. Two days with Twin Otter flights were identified where the convective mixed layer was quasi stationary, and hence the 10-s, 75-m data from the SGP Raman lidar could be analyzed to provide profiles of water vapor mixing ratio variance and skewness. Airborne water vapor observations measured during level flight legs were compared to the Raman lidar data, demonstrating good agreement in both variance and skewness. The results also illustrate the challenges of comparing a point sensor making measurements over time to a moving platform making similar measurements horizontally.

Factors controlling temporal variability of near-ground atmospheric ²²²Rn concentration over Central Europe

Specific activity of 222Rn in near-ground atmosphere has been measured quasi-continuously from January 2005 to December 2009 at two continental sites in Europe: Heidelberg (south-west Germany) and Krakow (southern Poland). Atmosphere was sampled at ca. 30 m and 20 m, respectively, above the local ground. Both stations were equipped with identical instrumentation. Regular observations of 222Rn were supplemented by measurements of surface fluxes of this gas in Krakow urban area, using two entirely different approaches. Atmospheric 222Rn concentrations varied at both sites in a wide range, from less than 2 Bq m−3 to approximately 40 Bq m−3 in Krakow and ca. 35 Bq m−3 in Heidelberg. Averaged over entire observation period, the 222Rn content in Krakow was approximately 30 % higher when compared to Heidelberg (5.86 ± 0.09 Bq −3 and 4.50 ± 0.07 Bq m−3, respectively). Distinct seasonality of 222Rn signal was visible in both presented time series, with higher values recorded generally during late summer and autumn. The surface 222Rn fluxes in Krakow also revealed a distinct seasonality, with broad maximum observed during summer and early autumn and minimum during the winter. Averaged over 5 yr observation period, the night-time surface 222Rn flux was equal 46.8 ± 2.4 Bq m−2 h−1. Although the atmospheric 222Rn levels at Heidelberg and Krakow appeared to be controlled primarily by local factors, it was possible to evaluate the "continental effect" in atmospheric 222Rn content between both sites, related to the gradual build-up of 222Rn concentration in the air masses travelling between Heidelberg and Krakow. The mean value of this load was equal 0.78 ± 0.12 Bq m−3. The measured minimum 222Rn concentrations at both sites and the difference between them was interpreted in the framework of a simple box model coupled with HYSPLIT analysis of air mass trajectories. Best fit of experimental and model data was obtained for the average 222Rn flux over the European continent equal 52 Bq m−2 h−1, the mean transport velocity of the air masses within convective mixed layer of PBL on their route from the Atlantic coast to Heidelberg and Krakow equal 3.5 m s−1, the mean rate constant of 222Rn removal across the top of PBL equal to the 222Rn decay constant and the mean height of the convective mixed layer height equal 1600 m.