Vertical Mass Transport Research Articles

Multiscale Atmospheric Overturning of the Indian Summer Monsoon as Seen through Isentropic Analysis

Abstract This study investigates multiscale atmospheric overturning during the 2009 Indian summer monsoon (ISM) using a cloud-permitting numerical model. The isentropic analysis technique adopted here sorts vertical mass fluxes in terms of the equivalent potential temperature of air parcels, which is capable of delineating the atmospheric overturning between ascending air parcels with high entropy and subsiding air parcels with low entropy. The monsoonal overturning is further decomposed into contributions from three characteristic scales: the basinwide ascent over the Indian monsoon domain, the regional-scale overturning associated with synoptic and mesoscale systems, and the convective-scale overturning. Results show that the convective-scale component dominates the upward mass transport in the lower troposphere while the region-scale component plays an important role by deepening the monsoonal overturning. The spatial variability of the convective-scale overturning is analyzed, showing intense convection over the Western Ghats and the Bay of Bengal while the deepest overturning is localized over northern India and the Himalayan foothills. The equivalent potential temperature in convective updrafts is higher over land than over the ocean or coastal regions. There is also substantial variability in the atmospheric overturning associated with the intraseasonal variability. The upward mass and energy transport increase considerably during the active phases of the ISM. A clear northeastward propagation in the peak isentropic vertical mass and energy transport over different characteristic regions can be found during the ISM, which corresponds to the intraseasonal oscillations of the ISM. Altogether, this study further demonstrates the utility of the isentropic analysis technique to characterize the spatiotemporal variations of convective activities in complex atmospheric flows.

General hydro-geological impact of cleats on underground coal gasification

Potential aquifer contamination has been a huge environmental concern for the implementation of underground coal gasification (UCG) as coal seams are frequently overlaid by an aquifer. To explore the general impact of cleats on the mass transfer problems in UCG, particularly in the vertical direction, a dual-permeability (DP) cleat model was integrated into a 3-D numerical model which is dimensionally extensive to account for the hydrostatic effect of an aquifer. It was found that cleats favour UCG production, and the prime operational parameter, i.e., the injection pressure, plays a dominant role in governing vertical mass transport in UCG. Syngas production retains some constancy with variable geological settings unless a surge in the injection pressure occurs. Higher pressure injection generates higher peak rates but delayed production. At near-hydrostatic injection, gas species attain a weak presence in the overburden contrary to expectations based on key geological parameters. However, gas saturation in the cap rock increases appreciably under above-hydrostatic injection. The role of higher injection pressure becomes so overwhelming that variations or heterogeneity in geological parameters produce little difference in the transfer problems.

Atmospheric Overturning across Multiple Scales of an MJO Event during the CINDY/DYNAMO Campaign

Abstract This study investigates the atmospheric overturning of the October 2011 MJO event observed during the Cooperative Indian Ocean Experiment on Intraseasonal Variability in the Year 2011 (CINDY)/DYNAMO field experiment using a cloud-permitting numerical model. The isentropic analysis is used to sort the vertical mass transport in terms of the equivalent potential temperature of the air parcels, which naturally decomposes the atmospheric overturning between ascending air with high entropy and subsiding air with low entropy. The circulation is further decomposed into contributions of four main scales: basinwide ascent, meridional overturning, regional overturning, and convection. Results show that the convective scale dominates the upward mass transport while larger scales play an important role both by allowing a deeper overturning and by modulating convective activity. There are substantial changes in the atmospheric overturning during different phases of this MJO event. Increased convective activity at low levels precedes the onset of the MJO by several days. The initiation of the MJO itself is associated with a substantial increase in the atmospheric overturning over the Indian Ocean. The subsequent eastward propagation of the MJO event can be clearly captured by the evolutions of convective-scale vertical mass fluxes at different altitudes. The equivalent potential temperatures of the rising and subsiding air parcels in the convective-scale overturning are also increased in the troposphere during the active phase of the MJO.

Isentropic Analysis on the Intensification of Hurricane Edouard (2014)

Abstract An isentropic analysis technique is adopted in this study to investigate the intensification of Hurricane Edouard (2014) predicted by an experimental real-time convection-permitting hurricane analysis and forecast system. This technique separates the vertical mass transport in terms of equivalent potential temperature θe for the rising air parcels at high entropy from the subsiding air at low entropy. It is found that as Edouard intensifies the vertical circulation becomes wider via the expansion of upward (downward) mass flux to higher (lower) θe. In the early developing stages, the asymmetric convection dominates the vertical circulation and leads to a remarkable upward mass flux maximum center in the upper troposphere. When Edouard becomes intense, the axisymmetric convection becomes important to the upper-level vertical mass transport while the asymmetric convection still dominates the low-level vertical mass transport. Development of the warm core in the eye leads to double maxima along the θe axis for both the isentropic-mean relative humidity and tangential velocity. The isentropic-mean properties such as the mid- to upper-level relative humidity, vertical velocity, and radial outflow decrease considerably while the mid- to upper-level vorticity enhances on the high-θe side before the onset of rapid intensification. The isentropic analysis also reveals that as Edouard intensifies the eye characterized by warm and dry core first forms in the low to middle troposphere and then gradually expands upward. The abovementioned results indicate that the isentropic framework may have the advantages of binning common variables with θe that could reflect the changes of the tropical cyclone structure in the inner-core region without a prior specification of the location of the storm center.

Sources and processes affecting the distribution of dissolved Nd isotopes and concentrations in the West Pacific

In the Atlantic, where deep circulation is vigorous, the dissolved neodymium (Nd) isotopic composition (expressed as ɛNd) is largely controlled by water mass mixing. In contrast, the factors influencing the ɛNd distribution in the Pacific, marked by sluggish circulation, is not clear yet. Indication for regional overprints in the Pacific is given based on its bordering volcanic islands. Our study aims to clarify the impact and relative importance of different Nd sources (rivers, volcanic islands), vertical (bio)geochemical processes and lateral water mass transport in controlling dissolved ɛNd and Nd concentration ([Nd]) distributions in the West Pacific between South Korea and Fiji. We find indication for unradiogenic continental input from South Korean and Chinese rivers to the East China Sea. In the tropical West Pacific, volcanic islands supply Nd to surface and subsurface waters and modify their ɛNd to radiogenic values of up to +0.7. These radiogenic signatures allow detailed tracing of currents flowing to the east and differentiation from westward currents with open ocean Pacific ɛNd composition in the complex tropical Pacific zonal current system. Modified radiogenic ɛNd of West Pacific intermediate to bottom waters upstream or within our section also indicates non-conservative behavior of ɛNd due to boundary exchange at volcanic island margins, submarine ridges, and with hydrothermal particles. Only subsurface to deep waters (3000 m) in the open Northwest Pacific show conservative behavior of ɛNd. In contrast, we find a striking correlation of extremely low (down to 2.77 pmol/kg Nd) and laterally constant [Nd] with the high-salinity North and South Pacific Tropical Water, indicating lateral transport of preformed [Nd] from the North and South Pacific subtropical gyres into the study area. This observation also explains the previously observed low subsurface [Nd] in the tropical West Pacific. Similarly, Western South Pacific Central Water, Antarctic Intermediate Water, and Lower Circumpolar Deep Water in the southern and equatorial West Pacific are marked by vertically and laterally almost invariant [Nd] indicating a dominance of conservative behavior of [Nd]. In contrast, Central and Intermediate Water in the North West Pacific are characterized by increasing [Nd] with depth reflecting Nd release from particles.Overall, our data demonstrate a dominant lateral transport control on [Nd] distributions and clear non-conservative modification of ɛNd in the West Pacific. The latter affords tracing of surface and subsurface zonal transport in the tropical Pacific, but prevents the use of ɛNd as strictly conservative tracer of the major meridionally circulating water masses in the West Pacific between 15°S and 28°N.

The role of fluid migration in the occurrence of shallow gas and gas hydrates in the south of the Pearl River Mouth Basin, South China Sea

High-amplitude anomalies (HAAs), bottom-simulating reflectors (BSRs), vertical acoustic blanking features, polygonal faults, and mass transport complexes are observed in high-resolution 3D seismic data from the south of the Pearl River Mouth Basin, South China Sea. We have focused on the role of fluid migration in the occurrence of shallow gas and gas hydrates. HAAs are interpreted to be shallow gas accumulations or gas hydrates, depending on their polarity and the presence of BSRs. Acoustic blanking features are probably caused by the presence of gas chimneys and overlying free gas due to features including chaotic and weak seismic reflections, pull-downs, and long distances of vertical extension. These gas chimneys root into the underlying late Oligocene gas reservoir. There is a close spatial relationship among the gas reservoir, gas chimneys, shallow gas accumulations, and BSRs, suggesting a genetic link. We have developed a schematic model to explain the processes of fluid migration, which supplied gas to the shallow gas accumulations and the gas hydrate system. When the gas reservoir was critically pressured, the focused fluid migration pathways may be formed or be reactivated, and gas was transported vertically along these conduits to the shallow strata. Then, there were two destinies for the migrated gas at shallow depths. One is that gas was trapped below low-permeability seals, such as mudstone or gas hydrates, forming shallow gas accumulations or free gas zones. The other one is that gas entered the gas hydrate stability zone and formed gas hydrates. In addition, this research indicates that the shallow hydrocarbon resource in this area may not be very prolific probably due to limited fluid migration.

Physically driven patchy O2 changes in the North Atlantic Ocean simulated by the CMIP5 Earth system models

AbstractThe subpolar North Atlantic is a key region for the oceanic uptake of heat, oxygen, and carbon dioxide. Centennial oxygen (O2) changes are investigated in the upper 700 m of the North Atlantic Ocean using a subset of Earth system models (ESMs) included in the Coupled Model Intercomparison Project phase 5. The climatological distributions of dissolved O2 averaged for the recent past period (1975–2005) are generally well captured, although the convective activity differs among the models in space and strength, and most models show a cold bias south of Greenland. By the end of the twenty‐first century, all models predict an increase in depth‐integrated temperature of 2–3°C, resultant solubility decrease, weakened vertical mass transport, decreased nutrient supply into the euphotic layer, and weakened export production. Despite an overall tendency of the North Atlantic to lose oxygen, patchy regions of O2 increase are observed due to the weakening of the North Atlantic Current (NAC) causing a regional solubility increase (the warming hole effect) and a decrease in the advection of subtropical, low‐O2 waters into the subpolar regions (the nutrient stream effect). Additionally, a shift in the NAC position contributes to localized O2 changes near the boundaries of water masses. The net O2 change reflects the combination of multiple factors leading to highly heterogeneous and model‐dependent patterns. Our results imply that changes in the strength and position of the NAC will likely play crucial roles in setting the pattern of O2 changes in future projections.

Role of convective mixing and evaporative cooling in shallow convection

AbstractLarge‐eddy simulations of shallow convection are used to evaluate the role of convective mixing and evaporative cooling in the vertical transport of mass, heat, and moisture in nonprecipitating shallow convection. Evaporative cooling is found to increase mass flux and the magnitude of heat and moisture fluxes, comparing twin large‐eddy simulations with either suppressed or active evaporative cooling. Nonetheless, subsiding shells transport mass downward even when evaporative cooling is suppressed, emphasizing that evaporative cooling is not the primary cause of existence of subsiding shells and accompanied buoyancy reversal. Instead, vertical convective mixing is found to be the primary reason of buoyancy reversal. Evaporative cooling yet accelerates downdrafts (updrafts) in the shell (cloudy) regions as well as increases the cloud cover in the lower cloud layer. The cloudy regions are more humid, and the liquid water potential temperature is lower compared to the evaporative‐cooling‐suppressed experiment. The primary effect of evaporative cooling is thus to increase the updraft core anomalies, thus enhancing vertical turbulent fluxes.

Nonlinear Etch Rate of Au-Assisted Chemical Etching of Silicon.

We demonstrated time-dependent mass transport mechanisms of Au-assisted chemical etching of Si substrates. Variations in the etch rate and surface topology were correlated with catalyst features and etching duration. Nonlinear etching characteristics were associated with the formation of pinholes and whiskers. Variable rates of mass transport as a function of whisker density accounted for the nonlinear etch rates of Si. Nanopinholes on Au catalysts facilitated the vertical mass transport of reactants and byproducts, which dramatically changed the etch rate, surface topology, and porosity of Si. The suggested transport models describe the transient mass transport and the corresponding chemical reactions.

Three-dimensional intracellular transport in neuron bodies and neurites investigated by label-free dispersion-relation phase spectroscopy.

Due to the limitations of fluorescence imaging techniques, the study of intracellular cargo is typically restricted to two-dimensional analyses. To overcome low light levels and the risk of phototoxicity, we employ quantitative phase imaging, a family of full-field imaging techniques that measure the optical path length shift introduced by the specimen. Specifically, we use spatial light interference microscopy (SLIM) to study the transport of mass in whole tomographic volumes and show that a time-correlation technique, dispersion-relation phase spectroscopy (DPS), can be used to simultaneously assay the horizontal and vertical traffic of mass through a cell. To validate our method, we compare the traffic inside cell bodies and neuronal extensions, showing that the vertical transport of mass may prove a more sensitive and interesting metric than similar measurements limited to a 2D, horizontal plane. © 2017 International Society for Advancement of Cytometry.

Turbulent mixing due to the Holmboe wave instability at high Reynolds number

We consider numerically the transition to turbulence and associated mixing in stratified shear flows with initial velocity distribution $\overline{U}(z,0)\,\boldsymbol{e}_{x}=U_{0}\,\boldsymbol{e}_{x}\tanh (z/d)$ and initial density distribution $\overline{\unicode[STIX]{x1D70C}}(z,0)=\unicode[STIX]{x1D70C}_{0}[1-\tanh (z/\unicode[STIX]{x1D6FF})]$ away from a hydrostatic reference state $\unicode[STIX]{x1D70C}_{r}\gg \unicode[STIX]{x1D70C}_{0}$. When the ratio $R=d/\unicode[STIX]{x1D6FF}$ of the characteristic length scales over which the velocity and density vary is equal to one, this flow is primarily susceptible to the classic well-known Kelvin–Helmholtz instability (KHI). This instability, which typically manifests at finite amplitude as an array of elliptical vortices, strongly ‘overturns’ the density interface of strong initial gradient, which nevertheless is the location of minimum initial gradient Richardson number $Ri_{g}(0)=Ri_{b}=g\unicode[STIX]{x1D70C}_{0}d/\unicode[STIX]{x1D70C}_{r}U_{0}^{2}$, where $Ri_{g}(z)=-([g/\unicode[STIX]{x1D70C}_{r}]\,\text{d}\overline{\unicode[STIX]{x1D70C}}/\text{d}z)/(\text{d}\overline{U}/\text{d}z)^{2}$ and $Ri_{b}$ is a bulk Richardson number. As is well known, at sufficiently high Reynolds numbers ($Re$), the primary KHI induces a vigorous but inherently transient burst of turbulence and associated irreversible mixing localised in the vicinity of the density interface, leading to a relatively well-mixed region bounded by stronger density gradients above and below. We explore the qualitatively different behaviour that arises when $R\gg 1$, and so the density interface is sharp, with $Ri_{g}(z)$ now being maximum at the density interface $Ri_{g}(0)=RRi_{b}$. This flow is primarily susceptible to Holmboe wave instability (HWI) (Holmboe, Geophys. Publ., vol. 24, 1962, pp. 67–113), which manifests at finite amplitude in this symmetric flow as counter-propagating trains of elliptical vortices above and below the density interface, thus perturbing the interface so as to exhibit characteristically cusped interfacial waves which thereby ‘scour’ the density interface. Unlike previous lower-$Re$ experimental and numerical studies, when $Re$ is sufficiently high the primary HWI becomes increasingly more three-dimensional due to the emergence of shear-aligned secondary convective instabilities. As $Re$ increases, (i) the growth rate of secondary instabilities appears to saturate and (ii) the perturbation kinetic energy exhibits a $k^{-5/3}$ spectrum for sufficiently large length scales that are influenced by anisotropic buoyancy effects. Therefore, at sufficiently high $Re$, vigorous turbulence is triggered that also significantly ‘scours’ the primary density interface, leading to substantial irreversible mixing and vertical transport of mass above and below the (robust) primary density interface. Furthermore, HWI produces a markedly more long-lived turbulence event compared to KHI at a similarly high $Re$. Despite their vastly different mechanics (i.e. ‘overturning’ versus ‘scouring’) and localisation, the mixing induced by KHI and HWI is comparable in both absolute terms and relative efficiency. Our results establish that, provided the flow Reynolds number is sufficiently high, shear layers with sharp density interfaces and associated locally high values of the gradient Richardson number may still be sites of substantial and efficient irreversible mixing.

Spontaneous formation of three-dimensionally ordered Bi-rich nanostructures within GaAs1−xBix/GaAs quantum wells

In this work, we report on the spontaneous formation of ordered arrays of nanometer-sized Bi-rich structures due to lateral composition modulations in Ga(As,Bi)/GaAs quantum wells grown by molecular beam epitaxy. The overall microstructure and chemical distribution is investigated using transmission electron microscopy. The information is complemented by synchrotron x-ray grazing incidence diffraction, which provides insight into the in-plane arrangement. Due to the vertical inheritance of the lateral modulation, the Bi-rich nanostructures eventually shape into a three-dimensional assembly. Whereas the Bi-rich nanostructures are created via two-dimensional phase separation at the growing surface, our results suggest that the process is assisted by Bi segregation which is demonstrated to be strong and more complex than expected, implying both lateral and vertical (surface segregation) mass transport. As demonstrated here, the inherent thermodynamic miscibility gap of Ga(As,Bi) alloys can be exploited to create highly uniform Bi-rich units embedded in a quantum confinement structure.

Isentropic Analysis of a Simulated Hurricane

Abstract Hurricanes, like many other atmospheric flows, are associated with turbulent motions over a wide range of scales. Here the authors adapt a new technique based on the isentropic analysis of convective motions to study the thermodynamic structure of the overturning circulation in hurricane simulations. This approach separates the vertical mass transport in terms of the equivalent potential temperature of air parcels. In doing so, one separates the rising air parcels at high entropy from the subsiding air at low entropy. This technique filters out oscillatory motions associated with gravity waves and separates convective overturning from the secondary circulation. This approach is applied here to study the flow of an idealized hurricane simulation with the Weather Research and Forecasting (WRF) Model. The isentropic circulation for a hurricane exhibits similar characteristics to that of moist convection, with a maximum mass transport near the surface associated with a shallow convection and entrainment. There are also important differences. For instance, ascent in the eyewall can be readily identified in the isentropic analysis as an upward mass flux of air with unusually high equivalent potential temperature. The isentropic circulation is further compared here to the Eulerian secondary circulation of the simulated hurricane to show that the mass transport in the isentropic circulation is much larger than the one in secondary circulation. This difference can be directly attributed to the mass transport by convection in the outer rainband and confirms that, even for a strongly organized flow like a hurricane, most of the atmospheric overturning is tied to the smaller scales.

Coherent Structures in the Boundary and Cloud Layers: Role of Updrafts, Subsiding Shells, and Environmental Subsidence

Abstract Coherent structures, such as updrafts, downdrafts/shells, and environmental subsidence in the boundary and cloud layers of shallow convection, are investigated using a new classification method. Using large-eddy simulation data, the new method first filters out background turbulence and small-scale gravity waves from the coherent part of the flow, composed of turbulent coherent structures and large-scale transporting gravity waves. Then the algorithm divides this coherent flow into “updrafts,” “downdrafts/shells,” “subsidence,” “ascendance,” and four other flow structures using an octant analysis. The novel method can systematically track structures from the cloud-free boundary layer to the cloud layer, thus allowing systematic analysis of the fate of updrafts and downdrafts. The frequency and contribution of the coherent structures to the vertical mass flux and transport of heat and moisture can then be investigated for the first time. Updrafts, subsidence, and downdrafts/subsiding shells—to a lesser extent—are shown to be the most frequent and dominant contributors to the vertical transport of heat and moisture in the boundary layer. Contrary to previous perspective, environmental subsidence transport is shown to be weak in the cloud layer. Instead, downdrafts/shells are the main downward transport contributors, especially in the trade inversion layer. The newly developed method in this study can be used to better evaluate the entrainment and detrainment of individual—or an ensemble of—coherent structures from the unsaturated boundary layer to the cloud layer.

Comparison of Microbial and Geochemical Conditions of Lignite Coal Spoil and Overburden Areas and their Environmental Impact

During the 20th century, Germany belonged to the world’s largest producers of lignite coal. Approaches for remediation of lignite spoil and overburden material were carried out since at least the last 50 years. Today, results of former remediation can be evaluated by their long-term effects. Coal spoil seepage waters can impact billions of m3 of ground water in involved areas for a very long time horizon (eternity impacts). In the last years our working group had studied especially lignite coal spoil sites in the Lusatia (East Germany) with sandy substrates and faster vertical mass transport processes. In new investigations, another coal spoil site in Middle Germany with more cohesive soil substrates and layers of drift clays was studied, and geochemical and microbial parameters were characterized. The investigated site partially had a remediated topsoil layer since 40 years, and since this time agricultural cultivation of different plants (grains, legumes, Pocaceae) was performed. In this paper, data of different microbial populations in the coal spoil substrate and geochemical background data are given, and they are compared with former data of microbial communities in sandy substrates and their impact to the environment. Finally, the success of long-term remediation of such sites, the influence onto microbial life and impact to the environments will be compared and discussed.

Mesoscale air transport at a midlatitude squall line in Europe – a numerical analysis

This study examines a fast‐propagating squall line over Western and Central Europe. A high‐resolution regional weather prediction model is used to analyze the three‐dimensional wind field and the corresponding air mass transport induced by the mesoscale convective system. The squall line can be divided into a bow‐shaped strong part with a continuous line of heavy convective precipitation, and a weaker part exhibiting several isolated multi‐cell storms. In the strong part, a deep and intense rear inflow jet occurs with relative wind speeds of 12–18 m s−1. This part is also characterized by moderate to strong wind shear, in contrast to the weaker part of the squall line where a shallow cold pool surges under the potentially unstable presquall air.To investigate the air mass transport along the squall line, in the model simulations passive fluid tracers are initialized in several vertical layers with different thermodynamic conditions. At the strong part of the squall line a very efficient vertical tracer transport occurs. Potentially warm low‐level air is lifted from the lowermost 2 km to the upper troposphere where an organized horizontal outflow is observed. The potentially coldest presquall air descends just in front of the updraught region and partly penetrates the convective precipitation area. In the weaker part of the squall line, the potential instability is not significantly reduced because the vertical air mass transport is less effective there. Additionally, passive tracers are used to analyze the composition of the cold pool which turns out to be a very heterogeneous mixture of descending post‐frontal cold air and presquall air originating at different heights.

Properties of a Mesoscale Convective System in the Context of an Isentropic Analysis

Abstract Application of an isentropic analysis of convective motions to a simulated mesoscale convective system is presented. The approach discriminates the vertical mass transport in terms of equivalent potential temperature. The scheme separates rising air at high entropy from subsiding air at low entropy. This also filters out oscillatory motions associated with gravity waves and isolates the overturning motions associated with convection and mesoscale circulation. The mesoscale convective system is additionally partitioned into stratiform and convective regions based on the radar reflectivity field. For each of the subregions, the mass transport derived in terms of height and an isentropic invariant of the flow is analyzed. The difference between the Eulerian mass flux and the isentropic counterpart is a significant and symmetric contribution of the buoyant oscillations to the upward and downward mass fluxes. Filtering out these oscillations results in substantial reduction of the diagnosed downward-to-upward convective mass flux ratio. The analysis is also applied to graupel and snow mixing ratios and number concentrations, illustrating the relationship of the particle formation process to the updrafts.

Upscale Error Growth in a High-Resolution Simulation of a Summertime Weather Event over Europe*

Abstract The growth of small-amplitude, spatially uncorrelated perturbations has been studied in a weather forecast of a 4-day period in the summer of 2007, using a large domain covering Europe and the eastern Atlantic and with explicitly resolved deep convection. The error growth follows the three-stage conceptual model of Zhang et al., with rapid initial growth (e-folding time about 0.5 h) on all scales, relaxing over about 20 h to a slow growth of the large-scale perturbations (e-folding time 12 h). The initial growth was confined to precipitating regions, with a faster growth rate where conditional instability was large. Growth in these regions saturated within 3–10 h, continuing for the longest where the precipitation rate was large. While the initial growth was mainly in the divergent part of the flow, the eventual slow growth on large scales was more in the rotational component. Spectral decomposition of the disturbance energy showed that the rapid growth in precipitating regions projected onto all Fourier components; however, the amplitude at saturation was too small to initiate the subsequent large-scale growth. Visualization of the disturbance energy showed it to expand outward from the precipitating regions at a speed corresponding to a deep tropospheric gravity wave. These results suggest a physical picture of error growth with a rapidly growing disturbance to the vertical mass transport in precipitating regions that spreads to the radius of deformation while undergoing geostrophic adjustment, eventually creating a balanced perturbation that continues to grow through baroclinic instability.

Entrainment and mixing in a laboratory model of oceanic overflow

AbstractWe present experimental measurements of a wall-bounded gravity current, motivated by characterizing natural gravity currents such as oceanic overflows. We use particle image velocimetry and planar laser-induced fluorescence to simultaneously measure the velocity and density fields as they evolve downstream of the initial injection from a turbulent channel flow onto a plane inclined at$10^\circ $with respect to horizontal. The turbulence level of the input flow is controlled by injecting velocity fluctuations upstream of the output nozzle. The initial Reynolds number based on the Taylor microscale of the flow,$R_{\lambda }$, is varied between 40 and 120, and the effects of the initial turbulence level are assessed. The bulk Richardson number$\mathit{Ri}$for the flow is${\sim }$0.3 whereas the gradient Richardson number$\mathit{Ri}_g$varies between 0.04 and 0.25, indicating that shear dominates the stabilizing effect of stratification. Kelvin–Helmholtz instability results in vigorous vertical transport of mass and momentum. We present baseline characterization of standard turbulence quantities and calculate, in several different ways, the fluid entrainment coefficient$E$, a quantity of considerable interest in mixing parameterization for ocean circulation models. We also determine the properties of mixing as represented by the flux Richardson number$\mathit{Ri}_f$as a function of$\mathit{Ri}_g$and diapycnal mixing parameter$K_{\rho }$versus the buoyancy Reynolds number$\mathit{Re}_b$. We find reasonable agreement with results from natural flows.

Aerosol pollution potential from major population centers

Abstract. Major population centers (MPCs), or megacities, represent the largest of growing urban agglomerations with major societal and environmental implications. In terms of air quality, they are seen as localized but strong emission sources of aerosols and trace gases which in turn affect air pollution levels in the city or in downwind regions. In the state-of-the-art atmospheric chemistry general circulation model EMAC, generic aerosol and gas-phase tracers with equal emission source strengths at 46 MPC locations are used to study the balance between local pollution build-up and pollution export, either vertically into the upper troposphere or horizontally in the lower troposphere. The insoluble gas-phase tracers with fixed lifetimes are transported with the atmospheric circulation, while the aerosol tracers also undergo gravitational sedimentation as well as dry and wet deposition processes. The strength of low-level tracer export depends on the location of the emission source and prevailing meteorology, in particular on atmospheric stability and the height of the boundary layer and the mixing out of this layer. In contrast, vertical transport of tracer mass depends on the tracer's solubility: the more soluble a tracer is, the less mass reaches altitudes above five kilometers. Hence, the mass of insoluble gas-phase tracer above five kilometers can be up to ten times higher than the hydrophilic aerosol mass from the same source. In the case of aerosol tracers, pollution build-up around the source is determined by meteorological factors which have only indirect effects on tracer lifetime, like surface wind, boundary layer height, and turbulent mixing, as well as those which affect the lifetime of the tracers such as precipitation. The longer a tracer stays in the atmosphere, the lower is the relative importance of the location of the source to the atmospheric mass, and thus the lower is the relative local pollution build-up. We further use aerosol deposition fields to estimate regions with high deposition, that is more than 1% or more than 5% of the corresponding tracer emission deposited in this region. In doing so, we find that the high deposition areas are larger for aerosols with diameters of 10.0 μm, and these differ less between the MPCs than for aerosols with diameters smaller than 2.5 μm due to faster deposition. Furthermore, cities in regions with high precipitation rates or unfavorable geographic locations, e.g., in a basin, suffer most of this high deposition. Most of the high deposition occurs over land, although about 50% of the MPCs are located along coastlines. By folding the aerosol deposition fields with geographical distributions of cropland, pasture, and forest, the impact on different land ecosystems is assessed. In general, forest is exposed most to deposition from MPCs while pastureland is least affected. Moreover, the impact on humans, measured with a threshold exceedance of pollutant surface mixing ratios, is more dependent on population densities than on the size of the area with a certain mixing ratio.