ALBERT

Description: This investigation analyzes the effect of vortex wakes on the Lagrangian displacement of particles induced by the passage of an obstacle in a two-dimensional incompressible and inviscid fluid. In addition to the trajectories of individual particles, we also study their drift and the corresponding total drift areas in the Föppl and Kirchhoff potential flow models. Our findings, which are obtained numerically and in some regimes are also supported by asymptotic analysis, are compared to the wakeless potential flow which serves as a reference. We show that in the presence of the Föppl vortex wake, some of the particles follow more complicated trajectories featuring a second loop. The appearance of an additional stagnation point in the Föppl flow is identified as a source of this effect. It is also demonstrated that, while the total drift area increases with the size of the wake for large vortex strengths, it is actually decreased for small circulation values. On the other hand, the Kirchhoff flow model is shown to have an unbounded total drift area. By providing a systematic account of the wake effects on the drift, the results of this study will allow for more accurate modeling of hydrodynamic stirring.

Description: Stratospheric ozone depletion from future nitrous oxide increases Atmospheric Chemistry and Physics, 14, 12967-12982, 2014 Author(s): W. Wang, W. Tian, S. Dhomse, F. Xie, J. Shu, and J. Austin We have investigated the impact of the assumed nitrous oxide (N 2 O) increases on stratospheric chemistry and dynamics using a series of idealized simulations with a coupled chemistry-climate model (CCM). In a future cooler stratosphere the net yield of NO y from N 2 O is shown to decrease in a reference run following the IPCC A1B scenario, but NO y can still be significantly increased by extra increases of N 2 O over 2001–2050. Over the last decade of simulations, 50% increases in N 2 O result in a maximal 6% reduction in ozone mixing ratios in the middle stratosphere at around 10 hPa and an average 2% decrease in the total ozone column (TCO) compared with the control run. This enhanced destruction could cause an ozone decline in the first half of this century in the middle stratosphere around 10 hPa, while global TCO still shows an increase at the same time. The results from a multiple linear regression analysis and sensitivity simulations with different forcings show that the chemical effect of N 2 O increases dominates the N 2 O-induced ozone depletion in the stratosphere, while the dynamical and radiative effects of N 2 O increases are overall insignificant. The analysis of the results reveals that the ozone depleting potential of N 2 O varies with the time period and is influenced by the environmental conditions. For example, carbon dioxide (CO 2 ) increases can strongly offset the ozone depletion effect of N 2 O.

Description: Corrigendum to "Air-sea exchange and gas-particle partitioning of polycyclic aromatic hydrocarbons in the Mediterranean" published in Atmos. Chem. Phys., 14, 8905–8915, 2014 Atmospheric Chemistry and Physics, 14, 12965-12965, 2014 Author(s): M. D. Mulder, A. Heil, P. Kukučka, J. Klánová, J. Kuta, R. Prokeš, F. Sprovieri, and G. Lammel No abstract available.

Description: Mechanistic study of secondary organic aerosol components formed from nucleophilic addition reactions of methacrylic acid epoxide Atmospheric Chemistry and Physics, 14, 12951-12964, 2014 Author(s): A. W. Birdsall, C. R. Miner, L. E. Mael, and M. J. Elrod Recently, methacrylic acid epoxide (MAE) has been proposed as a precursor to an important class of isoprene-derived compounds found in secondary organic aerosol (SOA): 2-methylglyceric acid (2-MG) and a set of oligomers, nitric acid esters, and sulfuric acid esters related to 2-MG. However, the specific chemical mechanisms by which MAE could form these compounds have not been previously studied with experimental methods. In order to determine the relevance of these processes to atmospheric aerosol, MAE and 2-MG have been synthesized and a series of bulk solution-phase experiments aimed at studying the reactivity of MAE using nuclear magnetic resonance (NMR) spectroscopy have been performed. The present results indicate that the acid-catalyzed MAE reaction is more than 600 times slower than a similar reaction of an important isoprene-derived epoxide, but is still expected to be kinetically feasible in the atmosphere on more acidic SOA. The specific mechanism by which MAE leads to oligomers was identified, and the reactions of MAE with a number of atmospherically relevant nucleophiles were also investigated. Because the nucleophilic strengths of water, sulfate, alcohols (including 2-MG), and acids (including MAE and 2-MG) in their reactions with MAE were found to be of similar magnitudes, it is expected that a diverse variety of MAE + nucleophile product species may be formed on ambient SOA. Thus, the results indicate that epoxide chain reaction oligomerization will be limited by the presence of high concentrations of non-epoxide nucleophiles (such as water); this finding is consistent with previous environmental chamber investigations of the relative humidity dependence of 2-MG-derived oligomerization processes and suggests that extensive oligomerization may not be likely on ambient SOA because of other competitive MAE reaction mechanisms.

Description: Megacity emission plume characteristics in summer and winter investigated by mobile aerosol and trace gas measurements: the Paris metropolitan area Atmospheric Chemistry and Physics, 14, 12931-12950, 2014 Author(s): S.-L. von der Weiden-Reinmüller, F. Drewnick, Q. J. Zhang, F. Freutel, M. Beekmann, and S. Borrmann For the investigation of megacity emission plume characteristics mobile aerosol and trace gas measurements were carried out in the greater Paris region in July 2009 and January–February 2010 within the EU FP7 MEGAPOLI project (Megacities: Emissions, urban, regional and Global Atmospheric POLlution and climate effects, and Integrated tools for assessment and mitigation). The deployed instruments measured physical and chemical properties of sub-micron aerosol particles, gas phase constituents of relevance for urban air pollution studies and meteorological parameters. The emission plume was identified based on fresh pollutant (e.g., particle-bound polycyclic aromatic hydrocarbons, black carbon, CO 2 and NO x ) concentration changes in combination with wind direction data. The classification into megacity influenced and background air masses allowed a characterization of the emission plume during summer and winter environmental conditions. On average, a clear increase of fresh pollutant concentrations in plume compared to background air masses was found for both seasons. For example, an average increase of 190% (+ 8.8 ng m −3 ) in summer and of 130% (+ 18.1 ng m −3 ) in winter was found for particle-bound polycyclic aromatic hydrocarbons in plume air masses. The aerosol particle size distribution in plume air masses was influenced by nucleation and growth due to coagulation and condensation in summer, while in winter only the latter process (i.e., particle growth) seemed to be initiated by urban pollution. The observed distribution of fresh pollutants in the emission plume – its cross sectional Gaussian-like profile and the exponential decrease of pollutant concentrations with increasing distance to the megacity – are in agreement with model results. Differences between model and measurements were found for plume center location, plume width and axial plume extent. In general, dilution was identified as the dominant process determining the axial variations within the Paris emission plume. For in-depth analysis of transformation processes occurring in the advected plume, simultaneous measurements at a suburban measurement site and a stationary site outside the metropolitan area using the mobile laboratory have proven to be most useful. Organic aerosol oxidation was observed in summer, while in winter transformation processes seemed to occur at a slower rate.

Description: A regional CO 2 observing system simulation experiment for the ASCENDS satellite mission Atmospheric Chemistry and Physics, 14, 12897-12914, 2014 Author(s): J. S. Wang, S. R. Kawa, J. Eluszkiewicz, D. F. Baker, M. Mountain, J. Henderson, T. Nehrkorn, and T. S. Zaccheo Top–down estimates of the spatiotemporal variations in emissions and uptake of CO 2 will benefit from the increasing measurement density brought by recent and future additions to the suite of in situ and remote CO 2 measurement platforms. In particular, the planned NASA Active Sensing of CO 2 Emissions over Nights, Days, and Seasons (ASCENDS) satellite mission will provide greater coverage in cloudy regions, at high latitudes, and at night than passive satellite systems, as well as high precision and accuracy. In a novel approach to quantifying the ability of satellite column measurements to constrain CO 2 fluxes, we use a portable library of footprints (surface influence functions) generated by the Stochastic Time-Inverted Lagrangian Transport (STILT) model in combination with the Weather Research and Forecasting (WRF) model in a regional Bayesian synthesis inversion. The regional Lagrangian particle dispersion model framework is well suited to make use of ASCENDS observations to constrain weekly fluxes in North America at a high resolution, in this case at 1° latitude × 1° longitude. We consider random measurement errors only, modeled as a function of the mission and instrument design specifications along with realistic atmospheric and surface conditions. We find that the ASCENDS observations could potentially reduce flux uncertainties substantially at biome and finer scales. At the grid scale and weekly resolution, the largest uncertainty reductions, on the order of 50%, occur where and when there is good coverage by observations with low measurement errors and the a priori uncertainties are large. Uncertainty reductions are smaller for a 1.57 μm candidate wavelength than for a 2.05 μm wavelength, and are smaller for the higher of the two measurement error levels that we consider (1.0 ppm vs. 0.5 ppm clear-sky error at Railroad Valley, Nevada). Uncertainty reductions at the annual biome scale range from ~40% to ~75% across our four instrument design cases and from ~65% to ~85% for the continent as a whole. Tests suggest that the quantitative results are moderately sensitive to assumptions regarding a priori uncertainties and boundary conditions. The a posteriori flux uncertainties we obtain, ranging from 0.01 to 0.06 Pg C yr −1 across the biomes, would meet requirements for improved understanding of long-term carbon sinks suggested by a previous study.

Description: Estimating regional fluxes of CO 2 and CH 4 using space-borne observations of XCH 4 : XCO 2 Atmospheric Chemistry and Physics, 14, 12883-12895, 2014 Author(s): A. Fraser, P. I. Palmer, L. Feng, H. Bösch, R. Parker, E. J. Dlugokencky, P. B. Krummel, and R. L. Langenfelds We use the GEOS-Chem global 3-D atmospheric chemistry transport model to interpret XCH 4 :XCO 2 column ratios retrieved from the Japanese Greenhouse Gases Observing Satellite (GOSAT). The advantage of these data over CO 2 and CH 4 columns retrieved independently using a full physics optimal estimation algorithm is that they are less prone to scattering-related regional biases. We show that the model is able to reproduce observed global and regional spatial (mean bias =0.7%) and temporal variations (global r 2 =0.92) of this ratio with a model bias 〈 2.5%. We also show that these variations are driven by emissions of CO 2 and CH 4 that are typically 6 months out of phase, which may reduce the sensitivity of the ratio to changes in either gas. To simultaneously estimate fluxes of CO 2 and CH 4 we use a maximum likelihood estimation approach. We use two approaches to resolve independent flux estimates of these two gases using GOSAT observations of XCH 4 :XCO 2 : (1) the a priori error covariance between CO 2 and CH 4 describing common source from biomass burning; and (2) also fitting independent surface atmospheric measurements of CH 4 and CO 2 mole fraction that provide additional constraints, improving the effectiveness of the observed GOSAT ratio to constrain flux estimates. We demonstrate the impact of these two approaches using numerical experiments. A posteriori flux estimates inferred using only the GOSAT ratios and taking advantage of the error covariance due to biomass burning are not consistent with the true fluxes in our experiments, as the inversion system cannot judge which species' fluxes to adjust. This reflects the weak dependence of XCH 4 :XCO 2 on biomass burning. We find that adding the surface data effectively provides an "anchor" to the inversion that dramatically improves the ability of the GOSAT ratios to infer both CH 4 and CO 2 fluxes. We show that the regional flux estimates inferred from GOSAT XCH 4 :XCO 2 ratios together with the surface mole fraction data during 2010 are typically consistent with or better than the corresponding values inferred from fitting XCH 4 or the full-physics XCO 2 data products, as judged by a posteriori uncertainties. We show that the fluxes inferred from the ratio measurements perform best over regions where there is a large seasonal cycle such as Tropical South America, for which we report a small but significant annual source of CO 2 compared to a small annual sink inferred from the XCO 2 data. We argue that given that the ratio measurements are less compromised by systematic error than the full physics data products, the resulting a~posteriori estimates and uncertainties provide a more faithful description of the truth. Based on our analysis we also argue that by using the ratios we may be reaching the current limits on the precision of these observed space-based data.

Description: The slow motion of a circular cylinder in a plane Poiseuille flow in a microchannel is analyzed for a wide range of cylinder radii and positions across the channel. The cylinder translates parallel to the channel walls and rotates about its axis. The Stokes approximation is used and the problem is solved analytically using the Papkovich-Fadle eigenfunction expansion and the least-squares method. The stream function and the pressure distribution of the flow field are obtained as results. The force and moment exerted on the cylinder, and the pressure change far from the cylinder, are calculated and shown as functions of the size and location of the cylinder. The results confirm some reciprocal relations exactly. In particular, the translational and rotational velocities of the drifting cylinder in the existing Poiseuille flow are determined. The induced pressure change, when the cylinder drifts in the Poiseuille flow, is also calculated. Some typical streamline patterns, depending on the size and location of the cylinder, are shown and discussed. When the cylinder translates and/or rotates in the channel blocked at infinity, a series of Moffatt eddies appears far from the cylinder in the channel, as expected.

Description: Interactions between capillary and elastic effects are relevant to a variety of applications from micro- and nano-scale manufacturing to biological systems. In this work, we investigate capillary flows in flexible, millimeter-scale cylindrical elastic tubes. We demonstrate that surface tension can cause sufficiently flexible tubes to collapse and coalesce spontaneously through non-axisymmetric buckling, and develop criteria for the initial deformation and complete collapse of a circular tube. Experimental results for capillary rise and evaporation of a liquid in a flexible tube are presented, and several regimes are seen for the equilibrium state of a flexible tube deforming under capillary pressure. Deformations of the tube walls are measured in different regimes and compared with a shell theory model. Analysis and experimental results show that despite the complex and non-axisymmetric deformed shapes of cylindrical structures, the elastocapillary length used in previous literature for flat plates and sheets can be used to predict the behavior of flexible tubes.

Description: Meteorological factors controlling low-level continental pollutant outflow across a coast Atmospheric Chemistry and Physics, 14, 13295-13312, 2014 Author(s): D. L. Peake, H. F. Dacre, J. Methven, and O. Coceal Coastal outflow describes the horizontal advection of pollutants from the continental boundary layer (BL) across a coastline. The outflow can ventilate polluted continental BLs and thus regulate air quality in highly populated coastal regions. This paper investigates the factors controlling coastal outflow and quantifies their importance as a ventilation mechanism. Tracers in the Met Office Unified Model (MetUM) are used to examine the magnitude and variability of coastal outflow over the eastern United States during summer 2004. Over the 4 week period examined, ventilation of tracer from the continental BL via coastal outflow occurs with the same magnitude as vertical ventilation via convection and advection. The relative importance of tracer decay rate, cross-coastal advection rate, and a parameter based on the relative continental and marine BL heights on coastal outflow is assessed by reducing the problem to a time-dependent box model. The ratio of the advection rate and decay rate is a dimensionless parameter which determines whether tracers are long-lived or short-lived. Long- and short-lived tracers exhibit different behaviours with respect to coastal outflow. Short-lived tracers exhibit large diurnal variability in coastal outflow but long-lived tracers do not. For short-lived tracers, increasing the advection rate increases the diurnally averaged magnitude of coastal outflow, but this has the opposite effect for very long-lived tracers. By using the box-model solutions to interpret the MetUM simulations, a land width is determined which represents the distance inland over which emissions contribute significantly to coastal outflow. A land width of between 100 and 400 km is found to be representative for a tracer with a lifetime of 24 h.

Description: We present laboratory experimental results demonstrating that librational forcing of an ellipsoidal container of water can produce intense motions through the mechanism of a libration driven elliptical instability (LDEI). These libration studies are conducted using an ellipsoidal acrylic container filled with water. A particle image velocimetry method is used to measure the 2D velocity field in the equatorial plane over hundreds libration cycles for a fixed Ekman number, E = 2 × 10 −5 . In doing so, we recover the libration induced base flow and a time averaged zonal flow. Further, we show that LDEI in non-axisymmetric container geometries is capable of driving both intermittent and saturated turbulent motions in the bulk fluid. Additionally, we measure the growth rate and amplitude of the LDEI induced excited flow in a fully ellipsoidal container at more extreme parameters than previously studied [Noir et al. , “Experimental study of libration-driven flows in nonaxisymmetric containers,” Phys. Earth Planet. Inter. 204-205 , 1 (2012); Cébron et al. , Phys. Fluids 24 , 061703, “Libration driven elliptical instability,” (2012)]. Excitation of bulk filling turbulence by librational forcing provides a mechanism for transferring rotational energy into turbulent fluid motion and thus can play an important role in the thermal evolution, interior dynamics, and magneto-hydrodynamics of librating bodies, as appear to be common in solar system settings [e.g., Comstock and Bills, “A solar system survey of forced librations in longitude,” J. Geophys. Res. Planets 108 , 1 (2003)].

Description: We numerically study the displacement flow of two iso-viscous Newtonian fluids in an inclined two-dimensional channel, formed by two parallel plates. The results are complementary to our previous studies on displacement flows in pipes and channels. The heavier displacing fluid moves the lighter displaced fluid in the downward direction. Three dimensionless groups largely describe these flows: the densimetric Froude number ( Fr ), the Reynolds number ( Re ), and the duct inclination (β). As a first order approximation, we are able to classify different flow regimes phenomenologically in a two-dimensional ( Fr ; Re cosβ/ Fr )-plane and provide leading order expressions for the transitions between different regimes. The stabilizing and/or de-stabilizing effects of the imposed mean flow on buoyant exchange flows (zero imposed velocity) are described for a broad range of dimensionless parameters.

Description: Compressible granular materials are involved in many applications, some of them being related to energetic porous media. Gas permeation effects are important during their compaction stage, as well as their eventual chemical decomposition. Also, many situations involve porous media separated from pure fluids through two-phase interfaces. It is thus important to develop theoretical and numerical formulations to deal with granular materials in the presence of both two-phase interfaces and gas permeation effects. Similar topic was addressed for fluid mixtures and interfaces with the Discrete Equations Method (DEM) [R. Abgrall and R. Saurel, “Discrete equations for physical and numerical compressible multiphase mixtures,” J. Comput. Phys. 186 (2), 361-396 (2003)] but it seemed impossible to extend this approach to granular media as intergranular stress [K. K. Kuo, V. Yang, and B. B. Moore, “Intragranular stress, particle-wall friction and speed of sound in granular propellant beds,” J. Ballist. 4 (1), 697-730 (1980)] and associated configuration energy [J. B. Bdzil, R. Menikoff, S. F. Son, A. K. Kapila, and D. S. Stewart, “Two-phase modeling of deflagration-to-detonation transition in granular materials: A critical examination of modeling issues,” Phys. Fluids 11 , 378 (1999)] were present with significant effects. An approach to deal with fluid-porous media interfaces was derived in Saurel et al. [“Modelling dynamic and irreversible powder compaction,” J. Fluid Mech. 664 , 348-396 (2010)] but its validity was restricted to weak velocity disequilibrium only. Thanks to a deeper analysis, the DEM is successfully extended to granular media modelling in the present paper. It results in an enhanced version of the Baer and Nunziato [“A two-phase mixture theory for the deflagration-to-detonation transition (DDT) in reactive granular materials,” Int. J. Multiphase Flow 12 (6), 861-889 (1986)] model as symmetry of the formulation is now preserved. Several computational examples are shown to validate and illustrate method’s capabilities.

Description: We perform a theoretical and numerical study of the Coulomb-driven electroconvection flow of a dielectric liquid between two coaxial cylinders. The specific case, where the inner to outer diameter ratio is 0.5, is analyzed. A strong unipolar injection of ions either from the inner or outer cylinder is considered to introduce free charge carriers into the system. A finite volume method is used to solve all governing equations including Navier-Stokes equations and a simplified set of Maxwell’s equations. The flow is characterized by a subcritical bifurcation in the finite amplitude regime. A linear stability criterion and a nonlinear one that correspond to the onset and stop of the flow motion, respectively, are linked with a hysteresis loop. In addition, we also explore the behavior of the system for higher values of the stability parameter. For inner injection, we observe a transition between the patterns made of 7 and 8 cells, before an oscillatory regime is attained. Such a transition leads to a second finite amplitude stability criterion. A simple modal analysis reveals that the competition of different modes is at the origin of this behavior. The charge density, as well as velocity field distributions is provided to help understand the bifurcation behavior.

Description: Reevaluation of stratospheric ozone trends from SAGE II data using a simultaneous temporal and spatial analysis Atmospheric Chemistry and Physics, 14, 13455-13470, 2014 Author(s): R. P. Damadeo, J. M. Zawodny, and L. W. Thomason This paper details a new method of regression for sparsely sampled data sets for use with time-series analysis, in particular the Stratospheric Aerosol and Gas Experiment (SAGE) II ozone data set. Non-uniform spatial, temporal, and diurnal sampling present in the data set result in biased values for the long-term trend if not accounted for. This new method is performed close to the native resolution of measurements and is a simultaneous temporal and spatial analysis that accounts for potential diurnal ozone variation. Results show biases, introduced by the way data are prepared for use with traditional methods, can be as high as 10%. Derived long-term changes show declines in ozone similar to other studies but very different trends in the presumed recovery period, with differences up to 2% per decade. The regression model allows for a variable turnaround time and reveals a hemispheric asymmetry in derived trends in the middle to upper stratosphere. Similar methodology is also applied to SAGE II aerosol optical depth data to create a new volcanic proxy that covers the SAGE II mission period. Ultimately this technique may be extensible towards the inclusion of multiple data sets without the need for homogenization.

Description: Phytoplankton patchiness, namely the heterogeneous distribution of microalgae over multiple spatial scales, dramatically impacts marine ecology. A spectacular example of such heterogeneity occurs in thin phytoplankton layers (TPLs), where large numbers of photosynthetic microorganisms are found within a small depth interval. Some species of motile phytoplankton can form TPLs by gyrotactic trapping due to the interplay of their particular swimming style (directed motion biased against gravity) and the transport by a flow with shear along the direction of gravity. Here we consider gyrotactic swimmers in numerical simulations of the Kolmogorov shear flow, both in laminar and turbulent regimes. In the laminar case, we show that the swimmer motion is integrable and the formation of TPLs can be fully characterized by means of dynamical systems tools. We then study the effects of rotational Brownian motion or turbulent fluctuations (appearing when the Reynolds number is large enough) on TPLs. In both cases, we show that TPLs become transient, and we characterize their persistence.

Description: In this study, rock friction ‘stick-slip’ experiments are used to develop constraints on models of earthquake recurrence. Constant-rate loading of bare rock surfaces in high quality experiments produces stick-slip recurrence that is periodic at least to second order. When the loading rate is varied, recurrence is approximately inversely proportional to loading rate. These laboratory events initiate due to a slip rate-dependent process that also determines the size of the stress drop and as a consequence, stress drop varies weakly but systematically with loading rate. This is especially evident in experiments where the loading rate is changed by orders of magnitude, as is thought to be the loading condition of naturally occurring, small repeating earthquakes driven by afterslip, or low-frequency earthquakes loaded by episodic slip. As follows from the previous studies referred to above, experimentally observed stress drops are well described by a logarithmic dependence on recurrence interval that can be cast as a non-linear slip-predictable model. The fault's rate dependence of strength is the key physical parameter. Additionally, even at constant loading rate the most reproducible laboratory recurrence is not exactly periodic, unlike existing friction recurrence models. We present example laboratory catalogs that document the variance and show that in large catalogs, even at constant loading rate, stress drop and recurrence co-vary systematically. The origin of this covariance is largely consistent with variability of the dependence of fault strength on slip rate. Laboratory catalogs show aspects of both slip and time predictability and successive stress drops are strongly correlated indicating a ‘memory’ of prior slip history that extends over at least one recurrence cycle.

Description: The model of gas bubble growth in high-viscous gas-saturated magmatic melt, subjected to rapid decompression, is presented in the current study. It is shown that consideration of unsteady character of the process is extremely important in a wide range of supersaturation. The analytical solution is found for the profile of dissolved gas concentration and the rate of bubble growth. The model of kinetics of overall degassing is developed. This model is based on distinguishing the so-called “forbidden” zone in the melt volume with suppressed formation of the new nucleation sites. The simple analytical dependences of the number of nucleating bubbles and typical nucleation time on the value of initial decompression were derived together with time dependence of volumetric concentration of the gas phase. Our results match the available experimental data.

Description: We investigate the effect of viscosity contrast on the stability of gravitationally unstable, diffusive layers in porous media. Our analysis helps evaluate experimental observations of various diffusive (boundary) layer models that are commonly used to study the sequestration of CO 2 in brine aquifers. We evaluate the effect of viscosity contrast for two basic models that are characterized with respect to whether or not the interface between CO 2 and brine is allowed to move. We find that diffusive layers are in general more unstable when viscosity decreases with depth within the layer compared to when viscosity increases with depth. This behavior is in contrast to the one associated with the classical displacement problem of gravitationally unstable diffusive layers that are subject to mean flow. For the classical problem, a greater instability is associated with the displacement of a more viscous, lighter fluid along the direction of gravity by a less viscous, heavier fluid. We show that the contrasting behavior highlighted in this study is a special case of the classical displacement problem that depends on the relative strength of the displacement and buoyancy velocities. We demonstrate the existence of a critical viscosity ratio that determines whether the flow is buoyancy dominated or displacement dominated. We explain the new behaviors in terms of the interaction of vorticity components related to gravitational and viscous effects.

Description: The development of a round liquid jet under the influence of a confined coaxial flow of an immiscible liquid of comparable density (central to annular flow density ratio of 8:10) was investigated in the vicinity of the nozzle exit. Two flow regimes were considered; one where the annular flow is faster than the central jet, so the central liquid jet is accelerated and one where the annular flow is slower, so the central liquid jet is decelerated. The central jet was visualised by high speed photography. Three modes of jet development were identified and classified in terms of the Reynolds number, Re, of the central jet which was in the range of 525 〈 Re 〈 2725, a modified definition of the Weber number, We, which allows the distinction between accelerating and deceleration flows and was in the range of −22 〈 We 〈 67 and the annular to central Momentum Ratio, MR, of the two streams which was in the range of 3.6 〈 MR 〈 91. By processing the time resolved jet images using Proper Orthogonal Decomposition (POD), it was possible to reduce the description of jet morphology to a small number of spatial modes, which isolated the most significant morphologies of the jet development. In this way, the temporal and spatial characteristics of the instabilities on the interface were clearly identified which highlights the advantages of POD over direct observation of the images. Relationships between the flow parameters and the interfacial waves were established. The wavelength of the interfacial instability was found to depend on the velocity of the fastest moving stream, which is contrary to findings for fluids with large density differences.

Description: Interaction of a vortex ring impinging on multiple permeable screens orthogonal to the ring axis was studied to experimentally investigate the persistence and decay of vortical structures inside the screen array using digital particle image velocimetry in a refractive index matched environment. The permeable screens had porosities (open area ratios) of 83.8%, 69.0%, and 55.7% and were held by a transparent frame that allowed the screen spacing to be changed. Vortex rings were generated using a piston-cylinder mechanism at nominal jet Reynolds numbers of 1000, 2000, and 3000 with piston stroke length-to-diameter ratios of 2 and 3. The interaction of vortex rings with the porous medium showed a strong dependence of the overall flow evolution on the screen porosity, with a central flow being preserved and vortex ring-like structures (with smaller diameter than the primary vortex ring) being generated near the centerline. Due to the large rod size used in the screens, immediate reformation of the transmitted vortex ring with size comparable to the primary ring (as has been observed with thin screens) was not observed in most cases. Since the screens have lower complexity and high open area ratios, centerline vortex ring-like flow structures formed with comparable size to the screen pore size and penetrated through the screens. In the case of low porosity screens (55.7%) with large screen spacing, re-emergence of large scale (large separation), weak vortical structures/pairs (analogous to a transmitted vortex ring) was observed downstream of the first screen. Additional smaller scale vortical structures were generated by the interaction of the vortex ring with subsequent screens. The size distribution of the generated vortical structures were shown to be strongly affected by porosity, with smaller vortical structures playing a stronger role as porosity decreased. Finally, porosity significantly affected the decay of total energy, but the effect of screen spacing decreased as porosity decreased.

Description: Vortex cavitation forming in the leading-edge vortices of a delta wing was examined to determine how the individual cavitation bubbles incepted, grew, interacted with the underlying vortical flow and produced acoustic tones. The non-cavitating vortical flow over the delta wing was chosen to be similar to those previously reported in the literature. It was found that vortex breakdown was unaffected by the presence of incipient and developed vortex cavitation bubbles in the vortex core. While some cavitation bubbles incepted, grew, and collapsed relatively quickly, others reached an equilibrium position wherein the bubble tip was stationary in the laboratory frame at a particular location along the vortex axis. For a given attack angle, the equilibrium location moved upstream with a reduction in free stream cavitation number. It is shown that the existence of these stationary vortex bubbles is possible when there is a balance between the axial growth of the bubble along the vortex axis and the opposite motion of the axial jetting flow in the vortex core, and only a single equilibrium position is possible along the axially evolving vortex for a given free stream cavitation number. These transient and stationary vortex bubbles emit significant cavitation noise upon inception, growth, and collapse. The spectral content of the noise produced was expected to be related to the interaction of the bubble with the surrounding vortical flow in a manner similar to that reported in previous studies, where sustained tones were similar to the underlying vortex frequency. However, in the present study, the dominant frequency and higher harmonics of the tones occur at a higher frequency than that of the underlying vortex. Hence, it is likely that the highly elongated stationary bubbles have higher-order volume oscillations compared to the two-dimensional radial mode of the vortex cores of vortex cavitation bubbles with much smaller diameter-to-length ratios.

Description: Emission of iodine-containing volatiles by selected microalgae species Atmospheric Chemistry and Physics, 14, 13327-13335, 2014 Author(s): U. R. Thorenz, L. J. Carpenter, R.-J. Huang, M. Kundel, J. Bosle, and T. Hoffmann In this study we present the results of an emission study of different phytoplankton samples in aqueous media treated with elevated ozone levels. Halocarbon measurements show that the samples tested released bromoform and different iodocarbons, including iodomethane, iodochloromethane and diiodomethane. Iodide and iodate levels in the liquid phase were representative of concentrations of surface water in a natural environment. Measurement of volatile iodine (I 2 ) emissions from two diatom samples ( Mediopyxis helysia and Porosira glacialis ) and the background sample (F/2 medium from filtered natural seawater) showed that the quantity of evolved I 2 depends on the ozone concentration in the air. This behaviour was assumed to be caused by the oxidation reaction mechanism of iodide with ozone. The I 2 emission flux agrees with model calculations at different iodide concentrations. The I 2 emission of a natural plankton concentrate sample was, however, very low compared to other samples and showed no dependence on ozone. The reason for this was shown to be the low iodide concentration in the algal suspension, which seems to be the limiting factor in the oxidative formation of I 2 .

Description: We derive equations for HTI and orthorhombic symmetries to analyze fluid substitution effects in porous fractured media. The derivations are based on the anisotropic Gassmann equation and linear slip theory. We assess the influence of fluid substitution (gas, brine, and oil), on elastic moduli, velocities, anisotropy, and azimuthal amplitude variations. We find that in the direction normal to fractures, P-wave moduli increase as much as 56% and P-wave velocity increases up to 19% for gas-to-brine substitution. For the direction parallel to fractures, P-wave velocity remains almost constant when porosity is low (5%), but can increase up to 4% if porosity is high (25%). Since P-waves in two different directions have different sensitivities to fluids and fractures, the Thomsen's parameters (defined for HTI and orthorhombic symmetries), ε and δ , are sensitive to fluid types and fractures. We also found that δ is sensitive to porosity for liquid saturation, but insensitive to porosity for the case of gas saturation. Gassmann assumes (and as has been observed) that shear modulus does not depend on fluids. And we observe no changes in shear-wave splitting ( γ ) for different fluids. The azimuthal amplitude variation is dependent on fluid types, fractures and porosity. We observe up to 12% increase in azimuthal amplitude variation for low porosity gas sands after brine saturation, and 6% decrease for high porosity gas sands. We find that the percentage changes in gas-to-oil substitution are about half that of the gas-to-brine case. The equations we have derived provide a useful tool to quantitatively evaluate the effects of fluid substitution on seismic anisotropy.

Description: In the central part of Fennoscandia the crust is currently rising, because of the delayed response of the viscous mantle to melting of the Late Pleistocene ice sheet. This process, called Glacial Isostatic Adjustment (GIA), causes a negative anomaly in the present-day static gravity field as isostatic equilibrium has not been reached yet. Several studies have tried to use this anomaly as a constraint on models of GIA, but the uncertainty in crustal and upper mantle structures has not been fully taken intoaccount. Therefore, our aim is to revisit this using improved crustal models and compensation techniques. We find that, in contrast with other studies, the effect of crustal anomalies on the gravity field cannot be effectively removed, because of uncertainties in the crustal and upper mantle density models. Our second aim is to estimate the effects on geophysical models, which assume isostatic equilibrium, after correcting the observed gravity field with numerical models for GIA. We show that correcting for GIA in geophysical modelling can give changes of several km in the thickness of structural layers of modeled lithosphere, which is a small but significant correction. Correcting the gravity field for GIA prior to assuming isostatic equilibrium and inferring density anomalies might be relevant in other areas with ongoing post-glacial rebound such as North America and the polar regions.

Description: Large variability of earthquake stress drops and scaled energy has been commonly reported in the literature, but it is difficult to assess how much of this variability is caused by underlying physical source processes rather than simply observational uncertainties. Here, we examine a variety of dynamically realistic rupture scenarios for circular and elliptical faults and investigate to what extent the variability in seismically estimated stress drops and scaled energy comes from differences in source geometry, rupture directivity, and rupture speeds. We numerically simulate earthquake source scenarios using a cohesive-zone model with the small-scale yielding limit, where the solution approaches a singular crack model with spontaneous healing of slip. Compared to symmetrical circular source models, asymmetrical models result in larger variability of estimated corner frequencies and scaled energy over the focal sphere. The general behavior of the spherical averages of corner frequencies and scaled energy in the subshear regime extends to the supershear regime, although shear Mach waves generated by the propagation of supershear rupture lead to much higher corner-frequency and scaled-energy estimates locally. Our results suggest that at least a factor of two difference in the spherical average of corner frequencies is expected in observational studies simply from variability in source characteristics almost independent of the actual stress drops, translating into a factor ofeight difference in estimated stress drops. Furthermore, radiation efficiency estimates derived from observed seismic spectra should not be directly interpreted as describing rupture properties unless there are independent constraints on rupture speed and geometry.

Description: Laminar flow over a periodic array of cylindrical surface roughness elements is simulated with an immersed boundary spectral method both to validate the method for subsequent studies and to examine how persistent streamwise vortices are introduced by a low Reynolds number roughness element. Direct comparisons are made with prior studies at a roughness-based Reynolds number Re k (= U ( k ) k / ν ) of 205 and a diameter to spanwise spacing ratio d / λ of 1/3. Downstream velocity contours match present and past experiments very well. The shear layer developed over the top of the roughness element produces the downstream velocity deficit. Upstream of the roughness element, the vortex topology is found to be consistent with juncture flow experiments, creating three cores along the recirculation line. Streamtraces stemming from these upstream cores, however, have unexpectedly little effect on the downstream flowfield as lateral divergence of the boundary layer quickly dissipates their vorticity. Long physical relaxation time of the recirculating wake behind the roughness remains a prominent issue for simulating this type of flowfield.

Description: Long-term variability of dust events in Iceland (1949–2011) Atmospheric Chemistry and Physics, 14, 13411-13422, 2014 Author(s): P. Dagsson-Waldhauserova, O. Arnalds, and H. Olafsson The long-term frequency of atmospheric dust observations was investigated for the southern part of Iceland and interpreted together with earlier results obtained from northeastern (NE) Iceland (Dagsson-Waldhauserova et al., 2013). In total, over 34 dust days per year on average occurred in Iceland based on conventionally used synoptic codes for dust observations. However, frequent volcanic eruptions, with the re-suspension of volcanic materials and dust haze, increased the number of dust events fourfold (135 dust days annually). The position of the Icelandic Low determined whether dust events occurred in the NE (16.4 dust days annually) or in the southern (S) part of Iceland (about 18 dust days annually). The decade with the most frequent dust days in S Iceland was the 1960s, but the 2000s in NE Iceland. A total of 32 severe dust storms (visibility 〈 500 m) were observed in Iceland with the highest frequency of events during the 2000s in S Iceland. The Arctic dust events (NE Iceland) were typically warm, occurring during summer/autumn (May–September) and during mild southwesterly winds, while the subarctic dust events (S Iceland) were mainly cold, occurring during winter/spring (March–May) and during strong northeasterly winds. About half of the dust events in S Iceland occurred in winter or at sub-zero temperatures. A good correlation was found between particulate matter (PM 10 ) concentrations and visibility during dust observations at the stations Vík and Stórhöfði. This study shows that Iceland is among the dustiest areas of the world and that dust is emitted year-round.

Description: On the origin of the occasional spring nitrate peak in Greenland snow Atmospheric Chemistry and Physics, 14, 13361-13376, 2014 Author(s): L. Geng, J. Cole-Dai, B. Alexander, J. Erbland, J. Savarino, A. J. Schauer, E. J. Steig, P. Lin, Q. Fu, and M. C. Zatko Ice core nitrate concentrations peak in the summer in both Greenland and Antarctica. Two nitrate concentration peaks in one annual layer have been observed some years in ice cores in Greenland from samples dating post-1900, with the additional nitrate peak occurring in the spring. The origin of the spring nitrate peak was hypothesized to be pollution transport from the mid-latitudes in the industrial era. We performed a case study on the origin of a spring nitrate peak in 2005 measured from a snowpit at Summit, Greenland, covering 3 years of snow accumulation. The effect of long-range transport of nitrate on this spring peak was excluded by using sulfate as a pollution tracer. The isotopic composition of nitrate (δ 15 N, δ 18 O and Δ 17 O) combined with photochemical calculations suggest that the occurrence of this spring peak is linked to a significantly weakened stratospheric ozone (O 3 ) layer. The weakened O 3 layer resulted in elevated UVB (ultraviolet-B) radiation on the snow surface, where the production of OH and NO x from the photolysis of their precursors was enhanced. Elevated NO x and OH concentrations resulted in enhanced nitrate production mainly through the NO 2 + OH formation pathway, as indicated by decreases in δ 18 O and Δ 17 O of nitrate associated with the spring peak. We further examined the nitrate concentration record from a shallow ice core covering the period from 1772 to 2006 and found 19 years with double nitrate peaks after the 1950s. Out of these 19 years, 14 of the secondary nitrate peaks were accompanied by sulfate peaks, suggesting long-range transport of nitrate as their source. In the other 5 years, low springtime O 3 column density was observed, suggesting enhanced local production of nitrate as their source. The results suggest that, in addition to direct transport of nitrate from polluted regions, enhanced local photochemistry can also lead to a spring nitrate peak. The enhanced local photochemistry is probably associated with the interannual variability of O 3 column density in the Arctic, which leads to elevated surface UV radiation in some years. In this scenario, enhanced photochemistry caused increased local nitrate production under the condition of elevated local NO x abundance in the industrial era.

Description: Aerosol hygroscopicity and cloud condensation nuclei activity during the AC 3 Exp campaign: implications for cloud condensation nuclei parameterization Atmospheric Chemistry and Physics, 14, 13423-13437, 2014 Author(s): F. Zhang, Y. Li, Z. Li, L. Sun, R. Li, C. Zhao, P. Wang, Y. Sun, X. Liu, J. Li, P. Li, G. Ren, and T. Fan Aerosol hygroscopicity and cloud condensation nuclei (CCN) activity under background conditions and during pollution events are investigated during the Aerosol-CCN-Cloud Closure Experiment (AC 3 Exp) campaign conducted at Xianghe, China in summer 2013. A gradual increase in size-resolved activation ratio (AR) with particle diameter ( D p ) suggests that aerosol particles have different hygroscopicities. During pollution events, the activation diameter ( D a ) measured at low supersaturation (SS) was significantly increased compared to background conditions. An increase was not observed when SS was 〉 0.4%. The hygroscopicity parameter (κ) was ~ 0.31–0.38 for particles in accumulation mode under background conditions. This range in magnitude of κ was ~ 20%, higher than κ derived under polluted conditions. For particles in nucleation or Aitken mode, κ ranged from 0.20–0.34 for background and polluted cases. Larger particles were on average more hygroscopic than smaller particles. The situation was more complex for heavy pollution particles because of the diversity in particle composition and mixing state. A non-parallel observation CCN closure test showed that uncertainties in CCN number concentration estimates ranged from 30–40%, which are associated with changes in particle composition as well as measurement uncertainties associated with bulk and size-resolved CCN methods. A case study showed that bulk CCN activation ratios increased as total condensation nuclei (CN) number concentrations ( N CN ) increased on background days. The background case also showed that bulk AR correlated well with the hygroscopicity parameter calculated from chemical volume fractions. On the contrary, bulk AR decreased with increasing total N CN during pollution events, but was closely related to the fraction of the total organic mass signal at m/z 44 ( f 44 ), which is usually associated with the particle's organic oxidation level. Our study highlights the importance of chemical composition in determining particle activation properties and underlines the significance of long-term observations of CCN under different atmospheric environments, especially regions with heavy pollution.

Description: Petrophysical properties of rocks and their applicability at larger scale are a challenging topic in Earth sciences. Petrophysical properties of rocks are severely affected by: boundary conditions, rock fabric/microstructure, and tectonics that require a multi-scale approach to be properly defined. Here we: (1) report laboratory measurements of density, porosity, permeability and P-wave velocities at increasing confining pressure conducted on Miocene foredeep sandstones (Frosinone Fm.); (2) compare the laboratory results with larger-scale geophysical investigations; (3) discuss the effect of thrusting on the properties of sandstones. At ambient pressure, laboratory porosity varied from 2.2% to 13.8% and P-wave velocities (Vp) from 1.5 km/s to 2.7 km/s. The P-wave velocity increased with confining pressure, reaching between 3.3 km/s to 4.7 km/s at 100 MPa. In situ Vp profiles, measured using sonic logs, matched the ultrasonic laboratory measurement well. The permeability varied between 1.4 × 10 -15  m 2 to 3.9 × 10 -15  m 2 and was positively correlated with porosity. The porosity and permeability of samples taken at various distances to the Olevano-Antrodoco fault plane progressively decreased with distance while P-wave velocity increased. At about 1 km from the fault plane, the relative variations reached 43%, 65% and 20% for porosity, permeability and P-wave velocity, respectively. This suggests that tectonic loading changed the petrophysical properties inherited from sedimentation and diagenesis. Using field constraints and assuming overburden-related inelastic compaction in the proximity of the fault plane, we conclude that the fault reached the mechanical condition for rupture in compression at differential stress of 64.8 MPa at a depth of 1500 m.

Description: Atmospheric submicron aerosol composition and particulate organic nitrate formation in a boreal forestland–urban mixed region Atmospheric Chemistry and Physics, 14, 13483-13495, 2014 Author(s): L. Q. Hao, A. Kortelainen, S. Romakkaniemi, H. Portin, A. Jaatinen, A. Leskinen, M. Komppula, P. Miettinen, D. Sueper, A. Pajunoja, J. N. Smith, K. E. J. Lehtinen, D. R. Worsnop, A. Laaksonen, and A. Virtanen The Puijo aerosol–cloud observation station is a unique measurement site for its location in the mixed region between the boreal forestland and the municipality of Kuopio, Finland. A measurement campaign was carried out at the station during fall 2010. An Aerodyne high-resolution time-of-flight aerosol mass spectrometer (HR-Tof-AMS) was deployed to characterize the atmospheric submicron aerosols. Positive matrix factorization (PMF) was applied to the unified high-resolution mass spectra organic species with NO + and NO 2 + ions to discover the intrinsic relationships between the organic and inorganic species and their daily cycles. On average, the submicron aerosols in this study were dominated by organic and sulfate species, composing 48.2 and 28.7% of total observed aerosol mass, respectively, with smaller contributions from ammonium (9.3%), nitrate (4.9%), chloride (0.8%) and BC (8.1%). The sources of these species included the primary emissions originating from the city area, secondary formation from both natural and anthropogenic emissions and regional transport. The PMF analysis succeeded in separating the mixed organic and inorganic spectra into three distinct organic and one inorganic factors. For organic factors, the semi-volatile oxygenated organic aerosol (SVOOA) and low-volatility oxygenated OA (LVOOA) accounted for 54.8 and 36.3% of total organic masses, respectively, while the hydrocarbon-like organic aerosol (HOA) accounted for 8.9% of total organics, with its main source from urban emissions. The inorganic factor is identified as NH 4 NO 3 , comprising 6.9% of the fitted aerosol mass by PMF. Based on the PMF results, the nitrate species were separated into organic and inorganic components, with the organic nitrates contributing one-third of the total nitrate mass. The results highlight both anthropogenic and biogenic emissions as important atmospheric aerosol sources in a forest–urban mixed region.

Description: Dynamical analysis of sea-breeze hodograph rotation in Sardinia Atmospheric Chemistry and Physics, 14, 13471-13481, 2014 Author(s): N. Moisseeva and D. G. Steyn This study investigates the diurnal evolution of sea-breeze (SB) rotation over an island at the middle latitudes. Earlier research on sea breezes in Sardinia shows that the onshore winds around various coasts of the island exhibit both the theoretically predicted clockwise rotation as well as seemingly anomalous anticlockwise rotation. A non-hydrostatic fully compressible numerical model (WRF) is used to simulate wind fields on and around the island on previously studied sea-breeze days, and is shown to capture the circulation on all coasts accurately. Diurnal rotation of wind is examined, and patterns of clockwise and anticlockwise rotation are identified. A dynamical analysis is performed by extracting individual forcing terms from the horizontal momentum equations. Analysis of several regions around the island shows that the direction of rotation is a result of a complex interaction between near-surface and synoptic pressure gradient, Coriolis and advection forcings. An idealized simulation is performed over an artificial island with dramatically simplified topography yet similar dimensions and latitude to Sardinia. Dynamical analysis of the idealized case reveals a rather different pattern of hodograph rotation to the real Sardinia, yet similar underlying dynamics. The research provides new insights into the dynamics underlying sea-breeze hodograph rotation, especially in coastal zones with a complex topography and/or coastline.

Description: Atmospheric measurements of ratios between CO 2 and co-emitted species from traffic: a tunnel study in the Paris megacity Atmospheric Chemistry and Physics, 14, 12871-12882, 2014 Author(s): L. Ammoura, I. Xueref-Remy, V. Gros, A. Baudic, B. Bonsang, J.-E. Petit, O. Perrussel, N. Bonnaire, J. Sciare, and F. Chevallier Measurements of CO 2 , CO, NO x and selected Volatile Organic Compounds (VOCs) mole fractions were performed continuously during a 10-day period in the Guy Môquet tunnel in Thiais, a peri-urban area about 15 km south of the centre of Paris, between 28 September and 8 October 2012. This data set is used here to identify the characteristics of traffic-emitted CO 2 by evaluating its ratios to co-emitted species for the first time in the Paris region. High coefficients of determination ( r 2 〉 0.7) are observed between CO 2 and certain compounds that are characteristic of the traffic source (CO, NO x , benzene, xylenes and acetylene). Weak correlations ( r 2 〈 0.2) are found with species such as propane, n-butane and i -butane that are associated with fuel evaporation, an insignificant source for CO 2 . To better characterise the traffic signal we focus only on species that are well-correlated with CO 2 and on rush-hour periods characterised by the highest traffic-related mole fractions. From those mole fractions we remove the nighttime-average weekday mole fraction obtained for each species that we infer to be the most appropriate background signal for our study. Then we calculate observed Δspecies / ΔCO 2 ratios, which we compare with the ones provided by the 2010 bottom–up high-resolved regional emission inventory from Airparif (the association in charge of monitoring the air quality in Île-de-France), focusing on local emission data for the specific road of the tunnel. We find an excellent agreement (2%) between the local inventory emission CO / CO 2 ratio and our observed ΔCO / ΔCO 2 ratio. Former tunnel experiments carried out elsewhere in the world provided observed ΔCO / ΔCO 2 ratios that differ from 49 to 592% to ours. This variability can be related to technological improvement of vehicles, differences in driving conditions, and fleet composition. We also find a satisfactory agreement with the Airparif inventory for n-propylbenzene, n-pentane and xylenes to CO 2 ratios. For most of the other species, the ratios obtained from the local emission inventory overestimate the observed ratios to CO 2 by 34 to more than 300%. However, the emission ratios of NO x , o -xylene and i -pentane are underestimated by 30 to 79%. One main cause of such high differences between the inventory and our observations is likely the obsolete feature of the VOCs speciation matrix of the inventory that has not been updated since 1998, although law regulations on some VOCs have occurred since that time. Our study bears important consequences, discussed in the conclusion, for the characterisation of the urban CO 2 plume and for atmospheric inverse modelling of urban CO 2 emissions.

Description: The use of SMILES data to study ozone loss in the Arctic winter 2009/2010 and comparison with Odin/SMR data using assimilation techniques Atmospheric Chemistry and Physics, 14, 12855-12869, 2014 Author(s): K. Sagi, D. Murtagh, J. Urban, H. Sagawa, and Y. Kasai The Superconducting Submillimeter-Wave Limb-Emission Sounder (SMILES) on board the International Space Station observed ozone in the stratosphere with high precision from October 2009 to April 2010. Although SMILES measurements only cover latitudes from 38° S to 65° N, the combination of data assimilation methods and an isentropic advection model allows us to quantify the ozone depletion in the 2009/2010 Arctic polar winter by making use of the instability of the polar vortex in the northern hemisphere. Ozone data from both SMILES and Odin/SMR (Sub-Millimetre Radiometer) for the winter were assimilated into the Dynamical Isentropic Assimilation Model for OdiN Data (DIAMOND). DIAMOND is an off-line wind-driven transport model on isentropic surfaces. Wind data from the operational analyses of the European Centre for Medium- Range Weather Forecasts (ECMWF) were used to drive the model. In this study, particular attention is paid to the cross isentropic transport of the tracer in order to accurately assess the ozone loss. The assimilated SMILES ozone fields agree well with the limitation of noise induced variability within the SMR fields despite the limited latitude coverage of the SMILES observations. Ozone depletion has been derived by comparing the ozone field acquired by sequential assimilation with a passively transported ozone field initialized on 1 December 2009. Significant ozone loss was found in different periods and altitudes from using both SMILES and SMR data: The initial depletion occurred at the end of January below 550 K with an accumulated loss of 0.6–1.0 ppmv (approximately 20%) by 1 April. The ensuing loss started from the end of February between 575 K and 650 K. Our estimation shows that 0.8–1.3 ppmv (20–25 %) of O 3 has been removed at the 600 K isentropic level by 1 April in volume mixing ratio (VMR).

Description: Henry's law constants of polyols Atmospheric Chemistry and Physics, 14, 12815-12837, 2014 Author(s): S. Compernolle and J.-F. Müller Henry's law constants (HLC) are derived for several polyols bearing between 2 and 6 hydroxyl groups, based on literature data for water activity, vapour pressure and/or solubility. While deriving HLC and depending on the case, also infinite dilution activity coefficients (IDACs), solid state vapour pressures or activity coefficient ratios are obtained as intermediate results. An error analysis on the intermediate quantities and the obtained HLC is included. For most compounds, these are the first values reported, while others compare favourably with literature data in most cases. Using these values and those from a previous work (Compernolle and Müller, 2014), an assessment is made on the partitioning of polyols, diacids and hydroxy acids to droplet and aqueous aerosol.

Description: Nitrous oxide emissions from a commercial cornfield ( Zea mays ) measured using the eddy covariance technique Atmospheric Chemistry and Physics, 14, 12839-12854, 2014 Author(s): H. Huang, J. Wang, D. Hui, D. R. Miller, S. Bhattarai, S. Dennis, D. Smart, T. Sammis, and K. C. Reddy Increases in observed atmospheric concentrations of the long-lived greenhouse gas nitrous oxide (N 2 O) have been well documented. However, information on event-related instantaneous emissions during fertilizer applications is lacking. With the development of fast-response N 2 O analyzers, the eddy covariance (EC) technique can be used to gather instantaneous measurements of N 2 O concentrations to quantify the exchange of nitrogen between the soil and atmosphere. The objectives of this study were to evaluate the performance of a new EC system, to measure the N 2 O flux with the system, and finally to examine relationships of the N 2 O flux with soil temperature, soil moisture, precipitation, and fertilization events. An EC system was assembled with a sonic anemometer and a fast-response N 2 O analyzer (quantum cascade laser spectrometer) and applied in a cornfield in Nolensville, Tennessee during the 2012 corn growing season (4 April–8 August). Fertilizer amounts totaling 217 kg N ha −1 were applied to the experimental site. Results showed that this N 2 O EC system provided reliable N 2 O flux measurements. The cumulative emitted N 2 O amount for the entire growing season was 6.87 kg N 2 O-N ha −1 . Seasonal fluxes were highly dependent on soil moisture rather than soil temperature. This study was one of the few experiments that continuously measured instantaneous, high-frequency N 2 O emissions in crop fields over a growing season of more than 100 days.

Description: On the hiatus in the acceleration of tropical upwelling since the beginning of the 21st century Atmospheric Chemistry and Physics, 14, 12803-12814, 2014 Author(s): J. Aschmann, J. P. Burrows, C. Gebhardt, A. Rozanov, R. Hommel, M. Weber, and A. M. Thompson Chemistry–climate models predict an acceleration of the upwelling branch of the Brewer–Dobson circulation as a consequence of increasing global surface temperatures, resulting from elevated levels of atmospheric greenhouse gases. The observed decrease of ozone in the tropical lower stratosphere during the last decades of the 20th century is consistent with the anticipated acceleration of upwelling. However, more recent satellite observations of ozone reveal that this decrease has unexpectedly stopped in the first decade of the 21st century, challenging the implicit assumption of a continuous acceleration of tropical upwelling. In this study we use three decades of chemistry-transport-model simulations (1980–2013) to investigate this phenomenon and resolve this apparent contradiction. Aside from a high-bias between 1985–1990, our model is able to reproduce the observed tropical lower stratosphere ozone record. A regression analysis identifies a significant decrease in the early period followed by a statistically robust trend-change after 2002, in qualitative agreement with the observations. We demonstrate that this trend-change is correlated with structural changes in the vertical transport, represented in the model by diabatic heating rates taken from the reanalysis product Era-Interim. These changes lead to a hiatus in the acceleration of tropical upwelling between 70–30 hPa and a southward shift of the tropical pipe at 30 and 100 hPa during the past decade, which appear to be the primary causes for the observed trend-change in ozone.

Description: Emissions of biogenic volatile organic compounds and subsequent photochemical production of secondary organic aerosol in mesocosm studies of temperate and tropical plant species Atmospheric Chemistry and Physics, 14, 12781-12801, 2014 Author(s): K. P. Wyche, A. C. Ryan, C. N. Hewitt, M. R. Alfarra, G. McFiggans, T. Carr, P. S. Monks, K. L. Smallbone, G. Capes, J. F. Hamilton, T. A. M. Pugh, and A. R. MacKenzie Silver birch ( Betula pendula ) and three Southeast Asian tropical plant species ( Ficus cyathistipula , Ficus benjamina and Caryota millis ) from the pantropical fig and palm genera were grown in a purpose-built and environment-controlled whole-tree chamber. The volatile organic compounds emitted from these trees were characterised and fed into a linked photochemical reaction chamber where they underwent photo-oxidation under a range of controlled conditions (relative humidity or RH ~65–89%, volatile organic compound-to-NO x or VOC / NO x ~3–9 and NO x ~2 ppbV). Both the gas phase and the aerosol phase of the reaction chamber were monitored in detail using a comprehensive suite of on-line and off-line chemical and physical measurement techniques. Silver birch was found to be a high monoterpene and sesquiterpene but low isoprene emitter, and its emissions were observed to produce measurable amounts of secondary organic aerosol (SOA) via both nucleation and condensation onto pre-existing seed aerosol ( Y SOA 26–39%). In contrast, all three tropical species were found to be high isoprene emitters with trace emissions of monoterpenes and sesquiterpenes. In tropical plant experiments without seed aerosol there was no measurable SOA nucleation, but aerosol mass was shown to increase when seed aerosol was present. Although principally isoprene emitting, the aerosol mass produced from tropical fig was mostly consistent (i.e. in 78 out of 120 aerosol mass calculations using plausible parameter sets of various precursor specific yields) with condensation of photo-oxidation products of the minor volatile organic compounds (VOCs) co-emitted; no significant aerosol yield from condensation of isoprene oxidation products was required in the interpretations of the experimental results. This finding is in line with previous reports of organic aerosol loadings consistent with production from minor biogenic VOCs co-emitted with isoprene in principally isoprene-emitting landscapes in Southeast Asia. Moreover, in general the amount of aerosol mass produced from the emissions of the principally isoprene-emitting plants was less than would be expected from published single-VOC experiments, if co-emitted species were solely responsible for the final SOA mass. Interpretation of the results obtained from the fig data sets leaves room for a potential role for isoprene in inhibiting SOA formation under certain ambient atmospheric conditions, although instrumental and experimental constraints impose a level of caution in the interpretation of the results. Concomitant gas- and aerosol-phase composition measurements also provide a detailed overview of numerous key oxidation mechanisms at work within the systems studied, and their combined analysis provides insight into the nature of the SOA formed.

Description: Surface-to-mountaintop transport characterised by radon observations at the Jungfraujoch Atmospheric Chemistry and Physics, 14, 12763-12779, 2014 Author(s): A. D. Griffiths, F. Conen, E. Weingartner, L. Zimmermann, S. D. Chambers, A. G. Williams, and M. Steinbacher Atmospheric composition measurements at Jungfraujoch are affected intermittently by boundary-layer air which is brought to the station by processes including thermally driven (anabatic) mountain winds. Using observations of radon-222, and a new objective analysis method, we quantify the land-surface influence at Jungfraujoch hour by hour and detect the presence of anabatic winds on a daily basis. During 2010–2011, anabatic winds occurred on 40% of days, but only from April to September. Anabatic wind days were associated with warmer air temperatures over a large fraction of Europe and with a shift in air-mass properties, even when comparing days with a similar mean radon concentration. Excluding days with anabatic winds, however, did not lead to a better definition of the unperturbed aerosol background than a definition based on radon alone. This implies that a radon threshold reliably excludes local influences from both anabatic and non-anabatic vertical-transport processes.

Description: Rare temperature histories and cirrus ice number density in a parcel and a one-dimensional model Atmospheric Chemistry and Physics, 14, 13013-13022, 2014 Author(s): D. M. Murphy A parcel and a one-dimensional model are used to investigate the temperature dependence of ice crystal number density. The number of ice crystals initially formed in a cold cirrus cloud is very sensitive to the nucleation mechanism and the detailed history of cooling rates during nucleation. A possible small spread in the homogeneous freezing threshold due to varying particle composition is identified as a sensitive nucleation parameter. In a parcel model, the slow growth rate of ice crystals at low temperatures inherently leads to a strong increase in ice number density at low temperatures. This temperature dependence is not observed. The model temperature dependence occurs for a wide range of assumptions and for either homogeneous or, less strongly, heterogeneous freezing. However, the parcel model also shows that random temperature fluctuations result in an extremely wide range of ice number densities. A one-dimensional model is used to show that the rare temperature trajectories resulting in the lowest number densities are disproportionately important. Low number density ice crystals sediment and influence a large volume of air. When such fall streaks are included, the ice number becomes less sensitive to the details of nucleation than it is in a parcel model. The one-dimensional simulations have a more realistic temperature dependence than the parcel mode. The one-dimensional model also produces layers with vertical dimensions of meters even if the temperature forcing has a much broader vertical wavelength. Unlike warm clouds, cirrus clouds are frequently surrounded by supersaturated air. Sedimentation through supersaturated air increases the importance of any process that produces small numbers of ice crystals. This paper emphasizes the relatively rare temperature trajectories that produce the fewest crystals. Other processes are heterogeneous nucleation, sedimentation from the very bottom of clouds, annealing of disordered to hexagonal ice, and entrainment.

Description: An improved dust emission model – Part 1: Model description and comparison against measurements Atmospheric Chemistry and Physics, 14, 13023-13041, 2014 Author(s): J. F. Kok, N. M. Mahowald, G. Fratini, J. A. Gillies, M. Ishizuka, J. F. Leys, M. Mikami, M.-S. Park, S.-U. Park, R. S. Van Pelt, and T. M. Zobeck Simulations of the dust cycle and its interactions with the changing Earth system are hindered by the empirical nature of dust emission parameterizations in weather and climate models. Here we take a step towards improving dust cycle simulations by using a combination of theory and numerical simulations to derive a physically based dust emission parameterization. Our parameterization is straightforward to implement into large-scale models, as it depends only on the wind friction velocity and the soil's threshold friction velocity. Moreover, it accounts for two processes missing from most existing parameterizations: a soil's increased ability to produce dust under saltation bombardment as it becomes more erodible, and the increased scaling of the dust flux with wind speed as a soil becomes less erodible. Our treatment of both these processes is supported by a compilation of quality-controlled vertical dust flux measurements. Furthermore, our scheme reproduces this measurement compilation with substantially less error than the existing dust flux parameterizations we were able to compare against. A critical insight from both our theory and the measurement compilation is that dust fluxes are substantially more sensitive to the soil's threshold friction velocity than most current schemes account for.

Description: The impact of volcanic aerosol on the Northern Hemisphere stratospheric polar vortex: mechanisms and sensitivity to forcing structure Atmospheric Chemistry and Physics, 14, 13063-13079, 2014 Author(s): M. Toohey, K. Krüger, M. Bittner, C. Timmreck, and H. Schmidt Observations and simple theoretical arguments suggest that the Northern Hemisphere (NH) stratospheric polar vortex is stronger in winters following major volcanic eruptions. However, recent studies show that climate models forced by prescribed volcanic aerosol fields fail to reproduce this effect. We investigate the impact of volcanic aerosol forcing on stratospheric dynamics, including the strength of the NH polar vortex, in ensemble simulations with the Max Planck Institute Earth System Model. The model is forced by four different prescribed forcing sets representing the radiative properties of stratospheric aerosol following the 1991 eruption of Mt. Pinatubo: two forcing sets are based on observations, and are commonly used in climate model simulations, and two forcing sets are constructed based on coupled aerosol–climate model simulations. For all forcings, we find that simulated temperature and zonal wind anomalies in the NH high latitudes are not directly impacted by anomalous volcanic aerosol heating. Instead, high-latitude effects result from enhancements in stratospheric residual circulation, which in turn result, at least in part, from enhanced stratospheric wave activity. High-latitude effects are therefore much less robust than would be expected if they were the direct result of aerosol heating. Both observation-based forcing sets result in insignificant changes in vortex strength. For the model-based forcing sets, the vortex response is found to be sensitive to the structure of the forcing, with one forcing set leading to significant strengthening of the polar vortex in rough agreement with observation-based expectations. Differences in the dynamical response to the forcing sets imply that reproducing the polar vortex responses to past eruptions, or predicting the response to future eruptions, depends on accurate representation of the space–time structure of the volcanic aerosol forcing.

Description: A model study on changes of European and Swiss particulate matter, ozone and nitrogen deposition between 1990 and 2020 due to the revised Gothenburg protocol Atmospheric Chemistry and Physics, 14, 13081-13095, 2014 Author(s): S. Aksoyoglu, J. Keller, G. Ciarelli, A. S. H. Prévôt, and U. Baltensperger We report a study of changes in air quality due to emission reductions using the chemical transport model CAMx. The model domain includes all of Europe with a nested domain over Switzerland. The model simulations were performed with emissions for 1990 (the reference year for the Gothenburg Protocol), 2005 (the reference year for the revised Gothenburg Protocol), 2006 (for model validation) and 2020 (the target year for the revised Gothenburg Protocol) using three emission scenarios prepared by IIASA/GAINS. Changes in ozone, particulate matter and nitrogen deposition are the central theme of the study. The modelled relative changes in the annual average PM 2.5 concentrations between 1990 and 2005 look reasonable based on various PM 10 and PM 2.5 observations in the past. The results obtained in this study suggest that annual mean concentrations of PM 2.5 decreased by about 20–50% in Europe. Simulations using the baseline scenario (BL 2020) suggest that PM 2.5 concentrations in 2020 will be about 30% lower than those in 2005. The largest predicted decrease in PM 2.5 , based on the MTFR (maximum technically feasible reduction) scenario, was about 60% and was located mainly in the eastern part of Europe. In the case of ozone, both model results and measurements show an increase in the mean ozone mixing ratios between 1990 and 2005. The observations, however, suggest a larger increase, indicating the importance of background ozone levels. Although emission reductions caused a decrease in peak ozone values, average ozone levels in polluted regions increased due to reduced titration with nitric oxide (NO). This caused a change in the frequency distribution of ozone. Model simulations using emission scenarios for 2020 suggest that annual average ozone mixing ratios will continue to increase. Changes in the levels of the damage indicators AOT40 for forests and SOMO35 are reported as well. The model results suggest that nitrogen deposition has decreased by 10–30% in the eastern part of Europe since 1990, while it has increased by about 20% in the Iberian Peninsula. The decrease is mainly due to the deposition of oxidized nitrogen species, whereas deposition of reduced nitrogen compounds increased. In Switzerland, nitrogen deposition is larger in the northern part of the Alps, where ammonia emissions are the highest. Applying the baseline scenario, we found that the deposition of oxidized nitrogen compounds will have decreased by a further 40% by 2020, whereas deposition of reduced species will continue to increase. This will lead to a 10–20% decrease in the total nitrogen deposition in most of the model domain, with a 10% increase in the eastern part of Europe.

Description: An improved dust emission model – Part 2: Evaluation in the Community Earth System Model, with implications for the use of dust source functions Atmospheric Chemistry and Physics, 14, 13043-13061, 2014 Author(s): J. F. Kok, S. Albani, N. M. Mahowald, and D. S. Ward The complex nature of mineral dust aerosol emission makes it a difficult process to represent accurately in weather and climate models. Indeed, results in the companion paper indicate that many large-scale models underestimate the dust flux's sensitivity to the soil's threshold friction velocity for erosion. We hypothesize that this finding explains why many dust cycle simulations are improved by using an empirical dust source function that shifts emissions towards the world's most erodible regions. Here, we both test this hypothesis and evaluate the performance of the new dust emission parameterization presented in the companion paper. We do so by implementing the new emission scheme into the Community Earth System Model (CESM) and comparing the resulting dust cycle simulations against an array of measurements. We find that the new scheme shifts emissions towards the world's most erodible regions in a manner that is strikingly similar to the effect of implementing a widely used source function based on satellite observations of dust source regions. Furthermore, model comparisons against aerosol optical depth measurements show that the new physically based scheme produces a statistically significant improvement in CESM's representation of dust emission, which exceeds the improvement produced by implementing a source function. These results indicate that the need to use an empirical source function is eliminated, at least in CESM, by the additional physics in the new scheme, and in particular by its increased sensitivity to the soil's threshold friction velocity. Since the threshold friction velocity is affected by climate changes, our results further suggest that many large-scale models underestimate the global dust cycle's climate sensitivity.

Description: Solutal Marangoni instability (SMI) is investigated both in 2D and 3D using a combined Cahn-Hilliard and Navier-Stokes model in a finite system. Fe-Sn is chosen as a representative alloy system since the phase diagram reveals a region with a miscibility gap, where two liquid phases, namely, the Fe-rich phase L 1 and the Sn-rich phase L 2 , are in chemical equilibrium. In 3D, considering a perturbed liquid cylinder ( L 2 phase) with a length of λ and a radius of R 0 embedded in the middle of a simulation box of λ × H × H (length × width × height) surrounded by the phase L 1 , we find that the perturbation induced Marangoni flow is either clockwise or anti-clockwise depending on the mean curvature difference between the convex and concave regions which is affected by the ratio of λ/ R 0 . The critical ratio of λ/ R 0 for SMI is shown to be invariant for different Marangoni numbers as well as independent of the geometrical properties of the L 1 phase. In 2D, a perturbed liquid pipe with a length of λ and a radius of R 0 embedded in the middle of a simulation box of λ × H (length × height) is taken into account. Due to different curvature constitution, the critical ratio of λ/ R 0 for SMI depends on the height of the L 1 phase.

Description: For several years, a promising Plasma Synthetic Jet actuator for high-speed flow control has been under development at ONERA. So far, its confined geometry and small space-time scales at play have prevented its full experimental characterization. Complementary accurate numerical simulations are then considered in this study in order to provide a complete aerothermodynamic description of the actuator. Two major obstacles have to be overcome with this approach: the modeling of the energy deposited by the electric arc and the accurate computation of the transient response of the cavity generating the pulsed jet. To solve the first problem, an Euler solver coupled with an electric circuit model was used to evaluate the energy deposition in the cavity. Such a coupling is performed by considering the electric field between the two electrodes. The second issue was then addressed by injecting these source terms in large Eddy simulations of the entire actuator. Aerodynamic results were finally compared with Schlieren visualizations. Using the proposed methodology, the temporal evolution of the jet front is remarkably well predicted.

Description: Using cloud ice flux to parametrise large-scale lightning Atmospheric Chemistry and Physics, 14, 12665-12682, 2014 Author(s): D. L. Finney, R. M. Doherty, O. Wild, H. Huntrieser, H. C. Pumphrey, and A. M. Blyth Lightning is an important natural source of nitrogen oxide especially in the middle and upper troposphere. Hence, it is essential to represent lightning in chemistry transport and coupled chemistry–climate models. Using ERA-Interim meteorological reanalysis data we compare the lightning flash density distributions produced using several existing lightning parametrisations, as well as a new parametrisation developed on the basis of upward cloud ice flux at 440 hPa. The use of ice flux forms a link to the non-inductive charging mechanism of thunderstorms. Spatial and temporal distributions of lightning flash density are compared to tropical and subtropical observations for 2007–2011 from the Lightning Imaging Sensor (LIS) on the Tropical Rainfall Measuring Mission (TRMM) satellite. The well-used lightning flash parametrisation based on cloud-top height has large biases but the derived annual total flash density has a better spatial correlation with the LIS observations than other existing parametrisations. A comparison of flash density simulated by the different schemes shows that the cloud-top height parametrisation has many more instances of moderate flash densities and fewer low and high extremes compared to the other parametrisations. Other studies in the literature have shown that this feature of the cloud-top height parametrisation is in contrast to lightning observations over certain regions. Our new ice flux parametrisation shows a clear improvement over all the existing parametrisations with lower root mean square errors (RMSEs) and better spatial correlations with the observations for distributions of annual total, and seasonal and interannual variations. The greatest improvement with the new parametrisation is a more realistic representation of the zonal distribution with a better balance between tropical and subtropical lightning flash estimates. The new parametrisation is appropriate for testing in chemistry transport and chemistry–climate models that use a lightning parametrisation.

Description: Tropospheric carbon monoxide over the Pacific during HIPPO: two-way coupled simulation of GEOS-Chem and its multiple nested models Atmospheric Chemistry and Physics, 14, 12649-12663, 2014 Author(s): Y.-Y. Yan, J.-T. Lin, Y. Kuang, D. Yang, and L. Zhang Global chemical transport models (CTMs) are used extensively to study air pollution and transport at a global scale. These models are limited by coarse horizontal resolutions that do not allow for a detailed representation of small-scale nonlinear processes over the pollutant source regions. Here we couple the global GEOS-Chem CTM and its three high-resolution nested models to simulate the tropospheric carbon monoxide (CO) over the Pacific Ocean during five High-performance Instrumented Airborne Platform for Environmental Research (HIAPER) Pole-to-Pole Observations (HIPPO) campaigns between 2009 and 2011. We develop a two-way coupler, the PeKing University CouPLer (PKUCPL), allowing for the exchange and interaction of chemical constituents between the global model (at 2.5° long. × 2° lat.) and the three nested models (at 0.667° long. × 0.5° lat.) covering Asia, North America, and Europe. The coupler obtains nested model results to modify the global model simulation within the respective nested domains, and simultaneously acquires global model results to provide lateral boundary conditions (LBCs) for the nested models. Compared to the global model alone, the two-way coupled simulation results in enhanced CO concentrations in the nested domains. Sensitivity tests suggest the enhancement to be a result of improved representation of the spatial distributions of CO, nitrogen oxides, and non-methane volatile organic compounds, the meteorological dependence of natural emissions, and other resolution-dependent processes. The relatively long lifetime of CO allows for the enhancement to be accumulated and carried across the globe. We found that the two-way coupled simulation increased the global tropospheric mean CO concentrations in 2009 by 10.4%, with a greater enhancement at 13.3% in the Northern Hemisphere. Coincidently, the global tropospheric mean hydroxyl radical (OH) was reduced by 4.2%, resulting in a 4.2% enhancement in the methyl chloroform lifetime (MCF; via reaction with tropospheric OH). The resulting CO and OH contents and MCF lifetime are closer to observation-based estimates. Both the global and the two-way coupled models capture the general spatiotemporal patterns of HIPPO CO over the Pacific. The two-way coupled simulation is much closer to HIPPO CO, with a mean bias of 1.1 ppb (1.4%) below 9 km compared to the bias at −7.2 ppb (−9.2%) for the global model alone. The improvement is most apparent over the North Pacific. Our test simulations show that the global model alone could resemble the two-way coupled simulation (especially below 4 km) by increasing its global CO emissions by 15% for HIPPO-1 and HIPPO-3, by 25% for HIPPO-2 and HIPPO-4, and by 35% for HIPPO-5. This has important implications for using the global model alone to constrain CO emissions. Thus, the two-way coupled simulation is a significantly improved model tool for studying the global impacts of air pollutants from major anthropogenic source regions.

Description: In a subduction system the force and the energy required to deform the overriding plate are generally thought to come from the negative buoyancy of the subducted slab and its potential energy, respectively. Such deformation might involve extension and backarc basin formation, or shortening and mountain building. How much of the slab's potential energy is consumed during overriding plate deformation remains unknown. In this work, we present dynamic three-dimensional laboratory experiments of progressive subduction with an overriding plate to quantify the force ( F OPD ) that drives overriding plate deformation and the associated energy dissipation rate (Φ OPD ), and we compare them with the negative buoyancy ( F BU ) of the subducted slab and its total potential energy release rate (Φ BU ), respectively. We varied the viscosity ratio between the plates and the sub-lithospheric upper mantle with η SP /η UM  = 157-560, and the thickness of the overriding plate with T OP  = 0.5-2.5 cm (scaling to 25-125 km in nature). The results show that F OPD / F BU has average values of 0.5-2.0%, with a maximum of 5.3% and Φ OPD /Φ BU has average values of 0.05-0.30%, with a maximum of 0.41%. The results indicate that only a small portion of the negative buoyancy of the slab and its potential energy are used to deform the overriding plate. Our models also suggest that the force required to deform the overriding plate is of comparable magnitude to the ridge push force. Furthermore, we show that in subduction models with an overriding plate bending dissipation at the subduction zone hinge remains low (3-15% during steady-state subduction).

Description: We investigated the influence of earthquake source complexity on the extent of inundation caused by the resulting tsunami. We simulated 100 scenarios with collocated sources of variable slip on the Hikurangi subduction-interface in the vicinity of Hawke's Bay and Poverty Bay in New Zealand and investigated the tsunami effects on the cities of Napier and Gisborne. Rupture complexity was found to have a first order effect on flow depth and inundation extent for local tsunami sources. The position of individual asperities in the slip distribution on the rupture interface control to some extent how severe inundation will be. However, predicting inundation extent in detail from investigating the distribution of slip on the rupture interface proves difficult. Assuming uniform slip on the rupture interface in tsunami models can underestimate the potential impact and extent of inundation. For example, simulation of an M W 8.7 to M W 8.8 earthquake with uniform slip reproduced the area that could potentially be inundated by equivalent non-uniform slip events of M W 8.4. De-aggregation, to establish the contribution of different sources with different slip distributions to the probabilistic hazard, cannot be performed based on magnitude considerations alone. We propose two predictors for inundation severity based on the offshore tsunami wave field using the linear wave equations in an attempt to keep costly simulations of full inundation to a minimum.

Description: Teleseismic transfer function analysis is used to investigate crust and upper mantle velocity structure in the vicinity of the active eastern Tennessee seismic zone (ETSZ). The ETSZ is associated with the New York – Alabama (NY-AL) magnetic lineament, a prominent aeromagnetic anomaly indicative of Grenville-age, basement structure. Radial component, P-wave transfer functions for ten short-period stations operated by the Center for Earthquake Research and Information (CERI) are inverted for velocity structure. Velocity profiles are also determined for three broadband stations by converting the instrument response to that of an S-13 short-period seismometer. Distinct differences in the velocity profiles are found for stations located on either side of the NY-AL magnetic lineament; velocities west of the lineament are lower than velocities to the east of the lineament in the upper 10 km and in the depth range 30 to 50 km. A gradational Moho boundary is found beneath several stations located in the Valley and Ridge province. A Moho boundary is absent at four Valley and Ridge stations located east of the magnetic lineament and south of 35.5°N.

Description: This integrated study provides significant insight into parameters controlling the acoustic and reservoir properties of microporous limestones, improving the knowledge of the relationships among petrophysic and microstructural content. Petrophysical properties measured from laboratory and logging tools (porosity, permeability, electrical conductivity and acoustic properties) have been coupled with thin section and SEM observations on the EST205 borehole from the Oxfordian limestone aquifer of the Eastern part of the Paris Basin. A major achievement is the establishment of the link between micrite microtexture types (particle morphology and nature of inter-crystal contacts) and the physical response, introducing a new effective and interesting rock-typing approach for microporous reservoirs. Fluid-flow properties are enhanced by the progressive augmentation of intercrystalline microporosity and associated pore throat diameter, as the coalescence of micrite particles decreases. Concerning acoustic properties, the slow increase of P-wave velocity can be seen as a reflection of crystal size and growing contact cementation leading to a more cohesive and stiffer micrite microtexture. By applying poroelasticity theory on our samples, we show that velocity dispersion can be a very useful tool for data discrimination in carbonates. This dispersion analysis highlights the presence of microcracks in the rocks, and their overall effect on acoustic and transport properties. The presence of microcracks is also confirmed with observations and permeability measurements under high confining pressure. Finally, a possible origin of high porous levels in neritic limestones is a mineralogical transformation of carbonates through freshwater-related diagenesis during subaerial exposure time. Finally, by applying poroelasticity theory on our samples, we show that velocity dispersion can be a very useful tool for data discrimination in carbonates.

Description: The “local porosity theory” proposed by Hilfer was revisited to develop a “local clay theory” (LCT) that establishes a quantitative relationship between the effective electrical conductivity and clay distribution in clay-rocks. This theory is primarily based on a “local simplicity” assumption; under this assumption, the complexity of spatial clay distribution can be captured by two local functions, namely, the local clay distribution and the local percolation probability, which are calculated from a partitioning of a mineral map. The local clay distribution provides information about spatial clay fluctuations and the local percolation probability describes the spatial fluctuations in the clay connectivity. This LCT was applied to (a) a mineral map made from a Callovo-Oxfordian mudstone sample and (b) (macroscopic) electrical conductivity measurements performed on the same sample. The direct and inverse modeling show two results. First, the textural and classical model assuming that the electrical anisotropy of clay-rock is mainly controlled by the anisotropy of the sole clay matrix provides inconsistent inverted values. An other textural effect, the anisotropy induced by elongated and oriented non-clayey grains, should be considered. Second, the effective conductivity values depend primarily on the choice of the inclusion-based models used in the LCT. The impact of local fluctuations of clay content and connectivity on the calculated effective conductivity is lower.

Description: Volcanic ash is often deposited in a hot state. Volcanic ash containing glass, deposited above the glass transition interval, has the potential to sinter viscously both to itself (particle-particle) and to exposed surfaces. Here, we constrain the kinetics of this process experimentally under non-isothermal conditions using standard glasses. In the absence of external load, this process is dominantly driven by surface relaxation. In such cases the sintering process is rate-limited by the melt viscosity, the size of the particles and the melt-vapour interfacial tension. We propose a polydisperse continuum model that describes the transition from a packing of particles to a dense pore-free melt and evaluate its efficacy in describing the kinetics of volcanic viscous sintering. We apply our model to viscous sintering scenarios for cooling crystal-poor rhyolitic ash using the 2008 eruption of Chaitén volcano as a case example. We predict that moderate linear cooling rates of 〉10 -1  °C.min -1 can result in the common observation of incomplete sintering and the preservation of pore networks.

Description: Snapshot and classical proper orthogonal decomposition (POD) are used to examine the large-scale, energetic motions in fully developed turbulent pipe flow at Re D = 47,000 and 93,000. The snapshot POD modes come in pairs, representing the same azimuthal mode number but with a simple phase shift. The first 10 snapshot POD modes, associated with the very large scale motions (VLSMs), contribute 43% of the average Reynolds shear stress, and for first 80 modes u ′ and v ′ are anti-correlated so that they all contribute to positive shear stress events. The attached motions are contained in the lower order modes, and detached motions do not appear until snapshot POD mode numbers ≥15. We find that snapshot POD can introduce mode mixing, which is avoided in classical POD. Classical POD also gives frequency information, confirming that the low order modes capture well the behavior of the very large scale motions.

Description: We explore the instabilities developed in a fluid in which viscosity depends on temperature. In particular, we consider a dependency that models a very viscous (and thus rather rigid) lithosphere over a convecting mantle. To this end, we study a 2D convection problem in which viscosity depends on temperature by abruptly changing its value by a factor of 400 within a narrow temperature gap. We conduct a study which combines bifurcation analysis and time-dependent simulations. Solutions such as limit cycles are found that are fundamentally related to the presence of symmetry. Spontaneous plate-like behaviors that rapidly evolve towards a stagnant lid regime emerge sporadically through abrupt bursts during these cycles. The plate-like evolution alternates motions towards either the right or the left, thereby introducing temporary asymmetries on the convecting styles. Further time-dependent regimes with stagnant and plate-like lids are found and described.

Description: We study the shock wave structure in a rarefied polyatomic gas based on a simplified model of extended thermodynamics in which the dissipation is due only to the dynamic pressure. In this case the differential system is very simple because it is a variant of Euler system with a new scalar equation for the dynamic pressure [T. Arima, S. Taniguchi, T. Ruggeri, and M. Sugiyama, Phys. Lett. A376, 2799–2803 (2012)]. It is shown that this theory is able to describe the three types of the shock wave structure observed in experiments: the nearly symmetric shock wave structure (Type A, small Mach number), the asymmetric structure (Type B, moderate Mach number), and the structure composed of thin and thick layers (Type C, large Mach number).

Description: In this paper, the scaling property of the inverse energy cascade and forward enstrophy cascade of the vorticity filed ω( x , y ) in two-dimensional (2D) turbulence is analyzed. This is accomplished by applying a Hilbert-based technique, namely Hilbert-Huang transform, to a vorticity field obtained from a 8192 2 grid-points direct numerical simulation of the 2D turbulence with a forcing scale k f = 100 and an Ekman friction. The measured joint probability density function p ( C , k ) of mode C i ( x ) of the vorticity ω and instantaneous wavenumber k ( x ) is separated by the forcing scale k f into two parts, which correspond to the inverse energy cascade and the forward enstrophy cascade. It is found that all conditional probability density function p ( C | k ) at given wavenumber k has an exponential tail. In the inverse energy cascade, the shape of p ( C | k ) does collapse with each other, indicating a nonintermittent cascade. The measured scaling exponent ζ ω I ( q ) is linear with the statistical order q , i.e., ζ ω I ( q ) = − q / 3 , confirming the nonintermittent cascade process. In the forward enstrophy cascade, the core part of p ( C | k ) is changing with wavenumber k , indicating an intermittent forward cascade. The measured scaling exponent ζ ω F ( q ) is nonlinear with q and can be described very well by a log-Poisson fitting: ζ ω F ( q ) = 1 3 q + 0.45 1 − 0 . 43 q . However, the extracted vorticity scaling exponents ζ ω ( q ) for both inverse energy cascade and forward enstrophy cascade are not consistent with Kraichnan's theory prediction. New theory for the vorticity field in 2D turbulence is required to interpret the observed scaling behavior.

Description: [1]  We analyse 1980-2010 ground displacements, to discern similarities or differences between Campi Flegrei (CF) inflations and deflations and highlight possible anomalies in particular areas. We show that the deformation pattern can be decomposed into two stationary (constant over time, except for a mere scaling factor) parts; both of them are satisfied by simple deformation sources. A quasi-horizontal elongated crack (oriented NW to SE, and embedded in an elastic layered half-space at a depth of about 3600 m) satisfies large-scale deformation. All source parameters but potency (volume change) are constant over time. Residual deformation is confined to the area of the Solfatara fumarolic field and satisfied by a small spheroid located at about 1900 m in depth. Again, all source parameters but potency are constant over time. The histories of the two sources are somewhat similar but not equal, supporting the existence of a genuine local deformation source at Solfatara against the emergence of a mere distortion of large-scale deformation. Although reality is probably much more complex, our simple model explains 1980-2010 CF deformation within ground-displacement data errors, and is consistent with Solfatara geochemical conceptual models, fumarolic geochemical data, and seismic attenuation imaging of CF. The observation that the CF deformation pattern can be decomposed into two stationary parts is hardly compatible with several recent works which proposed multiple sources with different features acting in different periods, fluid injections implying ample changes of large-scale deformation pattern over time, complex spatial and temporal patterns of distributions of volumetric sources.

Description: We report experimental observations of the controlled deformation of a dielectric liquid jet subjected to a local high-voltage electrostatic field in the direction normal to the jet. The jet deforms to the shape of an elliptic cylinder upon application of a normal electrostatic field. As the applied electric field strength is increased, the elliptic cylindrical jet deforms permanently into a flat sheet, and eventually breaks-up into droplets. We interpret this observation—the stretch of the jet is in the normal direction to the applied electric field—qualitatively using the Taylor-Melcher leaky dielectric theory, and develop a simple scaling model that predicts the critical electric field strength for the jet-to-sheet transition. Our model shows a good agreement with experimental results, and has a form that is consistent with the classical drop deformation criterion in the Taylor-Melcher theory. Finally, we statistically analyze the resultant droplets from sheet breakup, and find that increasing the applied electric field strength improves droplet uniformity and reduces droplet size.

Description: Normal impingement of a single droplet on a thin liquid film is investigated numerically solving the axisymmetric Navier-Stokes equations. Gravity and viscosity are taken into account whereas compressibility effects are neglected. Two phases are tracked by means of volume of fluid method and adaptive mesh refinement is used to increase accuracy of the interface. Numerical results are validated both qualitatively and quantitatively using experimental measurements. Effects of gas density, gas viscosity, and film thickness on the crown behavior are studied. Influence of droplet deviation from spherical shape on the crown behavior is investigated. It is shown that increasing the gas density leads to reduction of crown radius evolution rate, while gas viscosity does not affect the rate of crown radius evolution. Development rate of crown height decreases by increasing the gas density. Reynolds number and splashing regime can change the effect of gas viscosity on the crown height evolution. Deviation of droplet from sphere can change behavior of crown completely as result of change in droplet mass center position. Difference between numerical results and experimental ones is justified using different droplet shapes.

Description: We present an analytical study of peak mode isotachophoresis (ITP), and provide closed form solutions for sample distribution and electric field, as well as for leading-, trailing-, and counter-ion concentration profiles. Importantly, the solution we present is valid not only for the case of fully ionized species, but also for systems of weak electrolytes which better represent real buffer systems and for multivalent analytes such as proteins and DNA. The model reveals two major scales which govern the electric field and buffer distributions, and an additional length scale governing analyte distribution. Using well-controlled experiments, and numerical simulations, we verify and validate the model and highlight its key merits as well as its limitations. We demonstrate the use of the model for determining the peak concentration of focused sample based on known buffer and analyte properties, and show it differs significantly from commonly used approximations based on the interface width alone. We further apply our model for studying reactions between multiple species having different effective mobilities yet co-focused at a single ITP interface. We find a closed form expression for an effective-on rate which depends on reactants distributions, and derive the conditions for optimizing such reactions. Interestingly, the model reveals that maximum reaction rate is not necessarily obtained when the concentration profiles of the reacting species perfectly overlap. In addition to the exact solutions, we derive throughout several closed form engineering approximations which are based on elementary functions and are simple to implement, yet maintain the interplay between the important scales. Both the exact and approximate solutions provide insight into sample focusing and can be used to design and optimize ITP-based assays.

Description: A novel waveform modified from the standard-sinusoidal function is adopted to enhance the virtual aeroshaping effect of the synthetic jets positioned at the front stagnation point of a circular cylinder. The waveform is characterized by a control parameter, namely, the suction duty cycle factor k , which is the ratio of the time duration of the suction cycle to that of the blowing cycle. The strength of the synthetic jet vortex pair is enhanced by increasing the suction duty cycle factor. The periodic closed envelope forms upstream of the circular cylinder for k ≤ 1.00, while the quasi-steady open envelope forms for k ≥ 2.00, acting the virtual aeroshaping effect. As a result, both the statistical characteristics and the vortex dynamics of the near-wake flow field change with the suction duty cycle factor. The recirculation region downstream of the circular cylinder becomes smaller or even disappears, and thus, the drag coefficient over the circular cylinder is reduced by increasing the suction duty cycle factor to k ≥ 1.00. The statistical mean and fluctuating velocities show corresponding changes in the near wake with the different wake patterns. For k ≤ 0.50, the wake vortex shows the antisymmetric shedding mode which is similar with the natural case. For 1.00 ≤ k ≤ 2.00, the wake vortex shows the bistable state mode, where vortex sheds with symmetric or antisymmetric mode; the antisymmetric shedding mode dominates the global flow field for k = 1.00, while it is the symmetric shedding mode that dominates the flow field for k = 2.00. For k = 4.00, it shows the antisymmetric shedding mode with a shorter vortex formation length than the natural case. The above findings indicate that the virtual aeroshaping effect of the synthetic jets can be enhanced by increasing the suction duty cycle factor so as to increase the momentum coefficient while keeping other control parameters unchanged, providing us another way for effective flow control.

Description: The dynamics of vapor-liquid interface are important because interfacial instability determines bubble growth, detachment frequency, waiting time, shape of bubbles, and the interrelationship between bubble formation sites. In this study, a detailed numerical simulation has been performed to understand the transition in bubble release pattern and multimode bubble formation in saturated pool boiling. The interfaces drop down alternatively at the nodes and antinodes of the wavelengths dictated by Rayleigh-Taylor instability and Taylor-Helmholtz instability. Due to higher degrees of superheat, vapor jets emanate from nodes and antinodes. An attempt has been made to predict the maximum and minimum heat fluxes during saturated pool boiling.

Description: We present a combination of experiment, theory, and modelling on laminar mixing at large Péclet number. The flow is produced by oscillating electromagnetic forces in a thin electrolytic fluid layer, leading to oscillating dipoles, quadrupoles, octopoles, and disordered flows. The numerical simulations are based on the Diffusive Strip Method (DSM) which was recently introduced (P. Meunier and E. Villermaux, “The diffusive strip method for scalar mixing in two-dimensions,” J. Fluid Mech.662, 134–172 (2010)) to solve the advection-diffusion problem by combining Lagrangian techniques and theoretical modelling of the diffusion. Numerical simulations obtained with the DSM are in reasonable agreement with quantitative dye visualization experiments of the scalar fields. A theoretical model based on log-normal Probability Density Functions (PDFs) of stretching factors, characteristic of homogeneous turbulence in the Batchelor regime, allows to predict the PDFs of scalar in agreement with numerical and experimental results. This model also indicates that the PDFs of scalar are asymptotically close to log-normal at late stages, except for the large concentration levels which correspond to low stretching factors.

Description: Janus droplets are compound droplets that consist of two adhering drops of different fluids that are suspended in a third fluid. We use the Shan-Chen lattice Boltzmann method for multicomponent mixtures to simulate Janus droplets at rest and in shear. In this simulation model, interfacial tensions are not known a priori from the model parameters and must be determined using numerical experiments. We show that interfacial tensions obtained with the Young-Laplace law are consistent with those measured from the equilibrium geometry. The regimes of adhering, separated, and engulfing droplets were explored. Two different adhesion geometries were considered for two-dimensional simulations of Janus droplets in shear. The first geometry resembles two adhering circles with small overlap. In the second geometry, the two halves are semicircular. For both geometries, the rotation rate of the droplet depends on its orientation. The width of the periodic simulation domain also affects the rotation rate of both droplet types up to an aspect ratio of 6:1 (width:height). While the droplets with the first geometry oscillated about the middle of the domain, the droplets of the second geometry did not translate while rotating. A four-pole vortex structure inside droplets of the second geometry was found. These simulations of single Janus droplets reveal complex behaviour that implies a rich range of possibilities for the rheology of Janus emulsions.

Description: [1]  Detailed three-dimensional (3-D) P - and S -wave attenuation ( Qp and Qs ) models of the crust and upper mantle under the entire Northeast Japan (Tohoku) arc from the Japan Trench to the Japan Sea coast are determined, for the first time, using a large number of high-quality t * data measured precisely from P - and S -wave spectra of local earthquakes. The suboceanic earthquakes used in this work are relocated precisely using sP depth phases. Our results reveal a prominent landward dipping high- Q zone representing the subducting Pacific slab, a landward dipping intermediate to high Q zone in the mantle wedge between the Pacific coast and the volcanic front, and significant low- Q anomalies in the crust and mantle wedge between the volcanic front and the Japan Sea coast. Prominent high- Q patches surrounded by low- Q anomalies are revealed in the interplate megathrust zone under the Tohoku forearc where the great 2011 Tohoku-oki earthquake ( Mw 9.0) occurred. The high- Q patches in the megathrust zone generally exhibit large coseismic slips of megathrust earthquakes and large slip deficit on the plate interface. We think that these high- Q patches represent asperities in the megathrust zone, whereas the low- Q anomalies reflect weakly coupled areas. We also find that the hypocenters of the 2011 Tohoku-oki earthquake and its large foreshock ( Mw 7.3) and two large megathrust aftershocks ( Mw 7.4, 7.7) are located in areas where Qp , Qs and Qp / Qs change abruptly. These results suggest that structural heterogeneities in the megathrust zone control the interplate seismic coupling and the nucleation of megathrust earthquakes.

Description: Understanding the physics of water evaporation from saline porous media is important in many natural and engineering applications such as durability of building materials and preservation of monuments, water quality, and mineral-fluid interactions. We applied synchrotron x-ray micro-tomography to investigate the pore-scale dynamics of dissolved salt distribution in a three dimensional drying saline porous media using a cylindrical plastic column (15 mm in height and 8 mm in diameter) packed with sand particles saturated with CaI 2 solution (5% concentration by mass) with a spatial and temporal resolution of 12 μ m and 30 min, respectively. Every time the drying sand column was set to be imaged, two different images were recorded using distinct synchrotron x-rays energies immediately above and below the K-edge value of Iodine. Taking the difference between pixel gray values enabled us to delineate the spatial and temporal distribution of CaI 2 concentration at pore scale. Results indicate that during early stages of evaporation, air preferentially invades large pores at the surface while finer pores remain saturated and connected to the wet zone at bottom via capillary-induced liquid flow acting as evaporating spots. Consequently, the salt concentration increases preferentially in finer pores where evaporation occurs. Higher salt concentration was observed close to the evaporating surface indicating a convection-driven process. The obtained salt profiles were used to evaluate the numerical solution of the convection-diffusion equation (CDE). Results show that the macro-scale CDE could capture the overall trend of the measured salt profiles but fail to produce the exact slope of the profiles. Our results shed new insight on the physics of salt transport and its complex dynamics in drying porous media and establish synchrotron x-ray tomography as an effective tool to investigate the dynamics of salt transport in porous media at high spatial and temporal resolution.

Description: Pressure dependent isotopic fractionation in the photolysis of formaldehyde-d 2 Atmospheric Chemistry and Physics, 14, 551-558, 2014 Author(s): E. J. K. Nilsson, J. A. Schmidt, and M. S. Johnson The isotope effects in formaldehyde photolysis are the key link between the δD of methane emissions and the δD of atmospheric in situ hydrogen production. A few recent studies have suggested that a pressure dependence in the isotopic fractionation can partly explain enrichment of deuterium with altitude in the atmosphere. The mechanism and the extent of this pressure dependency is, however, not adequately described. In the present work D 2 CO and H 2 CO were photolyzed in a static reaction chamber at bath gas pressures of 50, 200, 400, 600 and 1000 mbar; these experiments compliment and extend our earlier work with HDCO vs. H 2 CO. The UV lamps used for photolysis emit light at wavelengths that primarily dissociate formaldehyde into molecular products, CO and H 2 or D 2 . The isotope effect k (H 2 CO)/ k (D 2 CO) = 3.16 ± 0.03 at 1000 mbar is in good agreement with results from previous studies. Similarly to what was previously shown for k (H 2 CO)/ k (HDCO), the isotope effect decreased as pressure decreased. In addition, a model was constructed using RRKM theory to calculate the lifetime of excited formaldehyde on the S 0 surface, to investigate its role in the observed pressure dependent photolytic fractionation of deuterium. The model shows that part of the fractionation is a result of competition between the isotopologue dependent rates of unimolecular dissociation and collisional relaxation. We suggest that the remaining fractionation is due to isotope effects in the rate of the non-radiative transition from S 1 to S 0 , which are not considered in the present model.

Description: On the consistency between global and regional methane emissions inferred from SCIAMACHY, TANSO-FTS, IASI and surface measurements Atmospheric Chemistry and Physics, 14, 577-592, 2014 Author(s): C. Cressot, F. Chevallier, P. Bousquet, C. Crevoisier, E. J. Dlugokencky, A. Fortems-Cheiney, C. Frankenberg, R. Parker, I. Pison, R. A. Scheepmaker, S. A. Montzka, P. B. Krummel, L. P. Steele, and R. L. Langenfelds Satellite retrievals of methane weighted atmospheric columns are assimilated within a Bayesian inversion system to infer the global and regional methane emissions and sinks for the period August 2009 to July 2010. Inversions are independently computed from three different space-borne observing systems and one surface observing system under several hypotheses for prior-flux and observation errors. Posterior methane emissions are compared and evaluated against surface mole fraction observations via a chemistry-transport model. Apart from SCIAMACHY (SCanning Imaging Absorption spectroMeter for Atmospheric CartograpHY), the simulations agree fairly well with the surface mole fractions. The most consistent configurations of this study using TANSO-FTS (Thermal And Near infrared Sensor for carbon Observation – Fourier Transform Spectrometer), IASI (Infrared Atmospheric Sounding Interferometer) or surface measurements induce posterior methane global emissions of, respectively, 565 ± 21 Tg yr −1 , 549 ± 36 Tg yr −1 and 538 ± 15 Tg yr −1 over the one-year period August 2009–July 2010. This consistency between the satellite retrievals (apart from SCIAMACHY) and independent surface measurements is promising for future improvement of CH 4 emission estimates by atmospheric inversions.

Description: Two hundred fifty years of aerosols and climate: the end of the age of aerosols Atmospheric Chemistry and Physics, 14, 537-549, 2014 Author(s): S. J. Smith and T. C. Bond Carbonaceous and sulfur aerosols have a substantial global and regional influence on climate, resulting in a net cooling to date, in addition to their impact on health and ecosystems. The magnitude of this influence has changed substantially over the past and is expected to continue to change into the future. An integrated picture of the changing climatic influence of black carbon, organic carbon and sulfate over the period 1850 through 2100, focusing on uncertainty, is presented using updated historical inventories and a coordinated set of emission projections. We describe, in detail, the aerosol emissions from the RCP4.5 scenario and its associated reference scenario. While aerosols have had a substantial impact on climate over the past century, we show that, by the end of the 21st century, aerosols will likely be only a minor contributor to radiative forcing due to increases in greenhouse gas forcing and a net global decrease in pollutant emissions. This outcome is even more certain under a successful implementation of a policy to limit greenhouse gas emissions as low-carbon energy technologies that do not emit appreciable aerosol or SO 2 are deployed.

Description: The role of horizontal model resolution in assessing the transport of CO in a middle latitude cyclone using WRF-Chem Atmospheric Chemistry and Physics, 14, 609-627, 2014 Author(s): C. A. Klich and H. E. Fuelberg We use the Weather Research and Forecasting with Chemistry (WRF-Chem) online chemical transport model to simulate a middle latitude cyclone in East Asia at three different horizontal resolutions (45, 15, and 5 km grid spacing). The cyclone contains a typical warm conveyor belt (WCB) with an embedded squall line that passes through an area having large surface concentrations (〉 400 ppbv) of carbon monoxide (CO). Model output from WRF-Chem is used to compare differences between the large-scale CO vertical transport by the WCB (the 45 km simulation) with the smaller-scale transport due to its convection (the 5 km simulation). Forward trajectories are calculated from WRF-Chem output using the Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) model. At 45 km grid spacing, the WCB exhibits gradual ascent, lofting surface CO to 6–7 km. Upon reaching the warm front, the WCB and associated CO ascend more rapidly and later turn eastward over the Pacific Ocean. Convective transport at 5 km resolution with explicitly resolved convection occurs much more rapidly, with surface CO lofted to altitudes greater than 10 km in 1 h or less. We also compute CO vertical mass fluxes over specified areas and times to compare differences in transport due to the different grid spacings. Upward CO flux exceeds 110 000 t h −1 in the domain with explicit convection when the squall line is at peak intensity, while fluxes from the two coarser resolutions are an order of magnitude smaller. Specific areas of interest within the 5 km domain are defined to compare the magnitude of convective transport to that within the entire 5 km region. Although convection encompasses only a small portion of the 5 km domain, it is responsible for ~40% of the upward CO transport. We also examine the vertical transport due to a short wave trough and its associated area of convection, not related to the cyclone, that lofts CO to the upper troposphere. Results indicate that fine-scale resolution with explicitly resolved convection is important when assessing the vertical transport of surface emissions in areas of deep convection.

Description: Comparing ECMWF AOD with AERONET observations at visible and UV wavelengths Atmospheric Chemistry and Physics, 14, 593-608, 2014 Author(s): V. Cesnulyte, A. V. Lindfors, M. R. A. Pitkänen, K. E. J. Lehtinen, J.-J. Morcrette, and A. Arola This paper presents validation results of the European Centre for Medium-Range Weather Forecasts (ECMWF) Integrated Forecasting System MACC (Monitoring Atmospheric Composition and Climate) re-analysis aerosol optical depth (AOD) for the period 2003–2006. We evaluate the MACC AOD at a UV wavelength (340 nm) and at mid-visible (500 and 550 nm) by comparing against ground-based AERONET measurements at 12 sites. The AERONET sites cover various parts of the globe and are categorized in three groups: urban/anthropogenic, biomass burning and dust, depending on the typically dominating aerosol type. This is the first time a global model such as the ECMWF has been evaluated for the performance of AOD at a UV wavelength. The results show that the MACC system generally provides a good representation of the AOD on a monthly basis, showing a realistic seasonal cycle. The model is mostly able to capture major dust load events and also the peak months of biomass burning correctly. For Kanpur and Solar Village, however, the model overestimates the AOD during the monsoon period when the aerosol load is generally low. When comparing hourly AOD values, the model–measurement agreement is better for biomass burning and dust sites than for urban sites, with an average correlation coefficient around 0.90 for biomass burning sites, around 0.77 for dust sites, and below 0.70 for urban sites. The AOD at 500 nm averaged over all sites shows only a small systematic difference between modeled and measured values, with a relative mean bias of 0.02. However, for the AOD at 340 nm the relative mean bias is −0.2. All sites included in the study show a relative mean bias at 340 nm smaller (or more negative) than that at 500 nm, indicating a strong wavelength dependence in the performance of the AOD in the MACC system. A comparison against fine and coarse mode AOD of the AERONET indicates that this has to do with the size distribution of the model: generally, the ECMWF model overestimates the contribution by coarse mode particles.

Description: The behaviour of low Reynolds number, non-Boussinesq fountains from four different nozzle geometries (one circular and three rectangular nozzles) are studied. High speed laser schlieren imaging is used to study the fountain behaviour (frequency and penetration height). Bi-orthogonal decomposition and dynamic mode decomposition (DMD) are used to understand the unsteady characteristics of fountains. The flow regimes of fountains are classified as steady, flapping, and flapping-bobbing type. The DMD technique successfully separates the bobbing oscillation from the combined flapping-bobbing oscillation of the fountain. The frequency of the flapping oscillation, and the frequency of the bobbing oscillation in the flapping-bobbing regime scales as St h Fr h = C 1 and S t h F r h 2 = C 2 , respectively, where the characteristic length scale is the smallest dimension ( h ) of the nozzle. The mean steady state penetration heights ( Z s / h ) of “forced” low Reynolds number non-Boussinesq fountains are independent of nozzle shape (circular and rectangular), and scales linearly with the Froude number.

Description: The large-scale properties of self-similar unstably stratified homogeneous (USH) turbulence are investigated using an eddy-damped quasi-normal markovianized approximation of the nonlinear term. This analysis shows that a special role is played by the wave vectors contained in the equatorial plane, i.e., the plane perpendicular to gravity. It is indeed in this plane that turbulent spectra reach their maxima and evolve linearly from their initial condition when their initial infrared exponent is smaller than 4. At other angles, this property is not satisfied and turbulent spectra eventually undergo an evolution dominated by nonlinear backscattering processes. The self-similar evolution of USH turbulence is also shown to be related to the properties of large scales. In particular, the asymptotic growth rate of the mixing length depends on the initial infrared exponent in the equatorial plane. Besides, the self-similar asymptotic values of the concentration and velocity correlations also depend on the properties of large scales. This allows to derive relations between the correlations and the growth rate parameter.

Description: Projections of atmospheric mercury levels and their effect on air quality in the United States Atmospheric Chemistry and Physics, 14, 783-795, 2014 Author(s): H. Lei, D. J. Wuebbles, X.-Z. Liang, Z. Tao, S. Olsen, R. Artz, X. Ren, and M. Cohen The individual and combined effects of global climate change and emissions changes from 2000 to 2050 on atmospheric mercury levels in the United States are investigated by using the global climate-chemistry model, CAM-Chem, coupled with a mercury chemistry-physics mechanism (CAM-Chem/Hg). Three future pathways from the Intergovernmental Panel on Climate Change (IPCC) Special Report on Emissions Scenarios (SRES) are considered, with the A1FI, A1B and B1 scenarios representing the upper, middle and lower bounds of potential climate warming, respectively. The anthropogenic and biomass burning emissions of mercury are projected from the energy use assumptions in the IPCC SRES report. Natural emissions from both land and ocean sources are projected by using dynamic schemes. TGM concentration increases are greater in the low latitudes than they are in the high latitudes, indicative of a larger meridional gradient than in the present day. In the A1FI scenario, TGM concentrations in 2050 are projected to increase by 2.1–4.0 ng m −3 for the eastern US and 1.4–3.0 ng m −3 for the western US. This spatial difference corresponds to potential increases in wet deposition of 10–14 μg m −2 for the eastern US and 2–4 μg m −2 for the western US. The increase in Hg(II) emissions tends to enhance wet deposition and hence increase the risk of higher mercury entering the hydrological cycle and ecosystem. In the B1 scenario, mercury concentrations in 2050 are similar to present level concentrations; this finding indicates that the domestic reduction in mercury emissions is essentially counteracted by the effects of climate warming and emissions increases in other regions. The sensitivity analyses show that changes in anthropogenic emissions contribute 32–53% of projected changes in mercury air concentration, while the independent contribution by climate change and its induced natural emissions change accounts for 47–68%.

Description: [1]  We detected 32,078 very small, local microearthquakes (average M L  = -1) during a 9-month deployment of five ocean bottom seismometers (OBS) on the periphery of the Trans-Atlantic Geotraverse (TAG) active mound. Seismicity rates were constant without any mainshock-aftershock behavior at ~243 events per day at the beginning of the experiment, 128 events per day after an instrument failed, and 97 events per day at the end of the experiment when whale calls increased background noise levels. The microearthquake seismograms are characterized by durations of 〈1 second and most have single-phase P -wave arrivals (i.e., no S -arrivals). We accurately located 6,207 of the earthquakes, with hypocenters clustered within a narrow depth interval from ~50-125 mbsf on the south and west flanks of the deposit. We model the microearthquakes as reaction-driven fracturing events caused by anhydrite deposition in the secondary circulation system of the hydrothermal mound, and show that under reasonable modeling assumptions an average event represents a volume increase of 31-58 cm 3 , yielding an annual (seismogenic) anhydrite deposition rate of 27-51 m 3 .

Description: Secondary flow cells are commonly observed in straight laboratory channels, where they are often associated with duct corners. Here, we present velocity measurements acquired with an acoustic Doppler current profiler in a straight reach of the Seine river (France). We show that a remarkably regular series of stationary flow cells spans across the entire channel. They are arranged in pairs of counter-rotating vortices aligned with the primary flow. Their existence away from the river banks contradicts the usual interpretation of these secondary flow structures, which invokes the influence of boundaries. Based on these measurements, we use a depth-averaged model to evaluate the momentum transfer by these structures, and find that it is comparable with the classical turbulent transfer.

Description: [1]  Acoustic velocity measurements were conducted during triaxial deformation tests of silty clay and clayey silt core samples from the Nankai subduction zone (IODP Expeditions 315, 316 and 333). We provide a new set of sonic wave velocity data, continuously measured during pressure increase and subsequent axial deformation. A new processing method for the experimental data was developed using seismic time series analysis. The results show that the identification of first arrivals by manual trace-by-trace picking alone can be erroneous. During axial deformation compressional wave velocities (Vp) range between about 1300 - 2200 m/s, and shear wave velocities (Vs) range between about 150 - 800 m/s. Vp slightly increases with rising effective confining pressure and effective axial stress. Samples from the accretionary prism toe show the highest Vp, while those from the forearc slope sediments show somewhat lower Vp. The samples from the incoming plate are slightly richer in clay minerals, and have the lowest values for Vp. Vs increases with higher effective confining pressures and effective axial stress, irrespective of composition and tectonic setting. Shear and bulk moduli, calculated from sonic velocities are in the range between 0.2 and 1.3 GPa and 3.85 and 8.41 GPa, respectively. The elastic moduli of the samples from the accretionary prism toe and the footwall of the megasplay fault ranging between 1.50 and 3.98 GPa are higher compared to those from the hanging wall samples and the incoming plate samples ranging between 0.59 and 0.88 Gpa. Thus, they allow differentiation between normal and over-consolidated sediments. The data show that in a tectono-sedimentary environment of only subtle compositional differences, the acoustic properties can be used to differentiate between stronger (accretionary prism toe) and weaker (forearc slope, incoming plate) sediments. Especially the Vp/Vs ratio may be instrumental in detecting zones of low effective stress and thus high pore fluid pressure.

Description: [1]  Studies of the petrology, mineral chemistry and rock-magnetic properties of nine pyroxenite xenoliths from Hannuoba basalts, northern North China Craton, have been made to determine the magnetization signature of the continental lower crust. These pyroxenites are weakly magnetic with low average susceptibility ( χ ) and saturation isothermal remanent magnetization ( M rs ) of 39.59 × 10 -8  m 3  kg -1 and 12.05 × 10 -3 Am 2 kg -1 , respectively. The magnetic minerals are mainly magnetite, pyrrhotite and Fe-rich spinel, which significantly contribute to χ and natural remanent magnetization. Magnetite occurs as interstitial micro-crystals together with zeolite aggregates, indicating a secondary origin in a supergene environment. In contrast, pyrrhotite and Fe-rich spinel were formed prior to the xenoliths’ ascent to the surface, as evidenced by their dominant occurrence as tiny inclusions and thin exsolution lamellae in pyroxene. The Fe-rich spinel has ~ 50% mole-fraction of Fe 3 O 4 and corresponds to the strongest magnetization, and its coexistence with Mg-rich spinel implies a reheating event due to the underplating of basaltic magma. Besides, armalcolite and ilmenite were found in the reaction rims between xenoliths and the basalt, but they contribute little to the whole rock magnetization. However, these pyroxenite xenoliths would be non-magnetic at in situ depths, as well as peridotite and mafic granulite xenoliths derived from the crust-mantle transition zone (~ 32-42 km). Therefore, we suggest the limiting depth of magnetization at the boundary between weakly magnetic deep-seated (lower crust and upper mantle) xenoliths and strongly magnetic Archean granulite facies rocks (~ 32 km) in Hannuoba, northern North China Craton.

Description: Analysis of fluxes across the turbulent/non-turbulent interface (TNTI) of turbulent boundary layers is performed using data from two-dimensional particle image velocimetry (PIV) obtained at high Reynolds numbers. The interface is identified with an iso-surface of kinetic energy, and the rate of change of total kinetic energy ( K ) inside a control volume with the TNTI as a bounding surface is investigated. Features of the growth of the turbulent region into the non-turbulent region by molecular diffusion of K , viscous nibbling, are examined in detail, focussing on correlations between interface orientation, viscous stress tensor elements, and local fluid velocity. At the level of the ensemble (Reynolds) averaged Navier-Stokes equations (RANS), the total kinetic energy K is shown to evolve predominantly due to the turbulent advective fluxes occurring through an average surface which differs considerably from the local, corrugated, sharp interface. The analysis is generalized to a hierarchy of length-scales by spatial filtering of the data as used commonly in Large-Eddy-Simulation (LES) analysis. For the same overall entrainment rate of total kinetic energy, the theoretical analysis shows that the sum of resolved viscous and subgrid-scale advective flux must be independent of scale. Within the experimental limitations of the PIV data, the results agree with these trends, namely that as the filter scale increases, the viscous resolved fluxes decrease while the subgrid-scale advective fluxes increase and tend towards the RANS values at large filter sizes. However, a definitive conclusion can only be made with fully resolved three-dimensional data, over and beyond the large dynamic spatial range presented here. The qualitative trends from the measurement results provide evidence that large-scale transport due to the energy-containing eddies determines the overall rate of entrainment, while viscous effects at the smallest scales provide the physical mechanism ultimately responsible for entrainment. Data spanning over a decade in Reynolds number suggest that the fluxes (or the entrainment velocity) scale with the friction velocity (or equivalently the local turbulent fluctuating velocity), whereas Taylor microscale and boundary-layer thickness are the appropriate length scales at small and large filter sizes, respectively.

Description: The direct and indirect radiative effects of biogenic secondary organic aerosol Atmospheric Chemistry and Physics, 14, 447-470, 2014 Author(s): C. E. Scott, A. Rap, D. V. Spracklen, P. M. Forster, K. S. Carslaw, G. W. Mann, K. J. Pringle, N. Kivekäs, M. Kulmala, H. Lihavainen, and P. Tunved We use a global aerosol microphysics model in combination with an offline radiative transfer model to quantify the radiative effect of biogenic secondary organic aerosol (SOA) in the present-day atmosphere. Through its role in particle growth and ageing, the presence of biogenic SOA increases the global annual mean concentration of cloud condensation nuclei (CCN; at 0.2% supersaturation) by 3.6–21.1%, depending upon the yield of SOA production from biogenic volatile organic compounds (BVOCs), and the nature and treatment of concurrent primary carbonaceous emissions. This increase in CCN causes a rise in global annual mean cloud droplet number concentration (CDNC) of 1.9–5.2%, and a global mean first aerosol indirect effect (AIE) of between +0.01 W m −2 and −0.12 W m −2 . The radiative impact of biogenic SOA is far greater when biogenic oxidation products also contribute to the very early stages of new particle formation; using two organically mediated mechanisms for new particle formation, we simulate global annual mean first AIEs of −0.22 W m −2 and −0.77 W m −2 . The inclusion of biogenic SOA substantially improves the simulated seasonal cycle in the concentration of CCN-sized particles observed at three forested sites. The best correlation is found when the organically mediated nucleation mechanisms are applied, suggesting that the first AIE of biogenic SOA could be as large as −0.77 W m −2 . The radiative impact of SOA is sensitive to the presence of anthropogenic emissions. Lower background aerosol concentrations simulated with anthropogenic emissions from 1750 give rise to a greater fractional CCN increase and a more substantial first AIE from biogenic SOA. Consequently, the anthropogenic indirect radiative forcing between 1750 and the present day is sensitive to assumptions about the amount and role of biogenic SOA. We also calculate an annual global mean direct radiative effect of between −0.08 W m −2 and −0.78 W m −2 in the present day, with uncertainty in the amount of SOA produced from the oxidation of BVOCs accounting for most of this range.

Description: [1]  The shallow velocity structure of the Lucky Strike segment of the Mid-Atlantic Ridge is investigated using seismic refraction and reflection techniques applied to downward continued multi-channel streamer data. We present a three-dimensional velocity model beneath the Lucky Strike Volcano with unprecedented spatial resolutions of a few hundred meters. These new constraints reveal large lateral variations in P-wave velocity structure beneath this feature. Throughout the study area, uppermost crustal velocities are significantly lower than those inferred from lower-resolution OBS studies, with the lowest values (1.8-2.2 km/s) found beneath the three central volcanic cones. Within the central volcano, distinct shallow units are mapped that likely represent a systematic process such as burial of older weathered surfaces. We infer that the entire upper part of the central volcano is young relative to the underlying median valley floor and that there has been little increase in the layer 2A velocities since emplacement. Layer 2A thins significantly across the axial valley bounding faults likely as the result of footwall uplift. The upper crustal velocities increase with age, on average, at a rate of ~0.875 km/s/Myr, similar to previous measurements at fast spreading ridges, suggesting hydrothermal sealing of small scale porosity is progressing at normal to enhanced rates.

Description: In a recent direct numerical simulation (DNS) study [P. K. Yeung and K. R. Sreenivasan, “ Spectrum of passive scalars of high molecular diffusivity in turbulent mixing,” J. Fluid Mech.716, R14 (2013)] with Schmidt number as low as 1/2048, we verified the essential physical content of the theory of Batchelor, Howells, and Townsend [“Small-scale variation of convected quantities like temperature in turbulent fluid. 2. The case of large conductivity,” J. Fluid Mech.5, 134 (1959)] for turbulent passive scalar fields with very strong diffusivity, decaying in the absence of any production mechanism. In particular, we confirmed the existence of the −17/3 power of the scalar spectral density in the so-called inertial-diffusive range. In the present paper, we consider the DNS of the same problem, but in the presence of a uniform mean gradient, which leads to the production of scalar fluctuations at (primarily) the large scales. For the parameters of the simulations, the presence of the mean gradient alters the physics of mixing fundamentally at low Peclet numbers. While the spectrum still follows a −17/3 power law in the inertial-diffusive range, the pre-factor is non-universal and depends on the magnitude of the mean scalar gradient. Spectral transfer is greatly reduced in comparison with those for moderately and weakly diffusive scalars, leading to several distinctive features such as the absence of dissipative anomaly and a new balance of terms in the spectral transfer equation for the scalar variance, differing from the case of zero gradient. We use the DNS results to present an alternative explanation for the observed scaling behavior, and discuss a few spectral characteristics in detail.

Description: Long-range transport of giant particles in Asian dust identified by physical, mineralogical, and meteorological analysis Atmospheric Chemistry and Physics, 14, 505-521, 2014 Author(s): G. Y. Jeong, J. Y. Kim, J. Seo, G. M. Kim, H. C. Jin, and Y. Chun Giant particles transported over long distances are generally of limited concern in atmospheric studies due to their low number concentrations in mineral dust and possible local origin. However, they can play an important role in regional circulation of earth materials due to their enormous volume concentration. Asian dust laden with giant particles was observed in Korea on 31 March 2012, after a migration of about 2000 km across the Yellow Sea from the Gobi Desert. Scanning electron microscopy (SEM) revealed that 20% of the particles exceeded 10 μm in equivalent sphere diameter, with a maximum of 60 μm. The median diameter from the number distribution was 5.7 μm, which was larger than the diameters recorded of 2.5 and 2.9 μm in Asian dust storms in 2010 and 2011, respectively, and was consistent with independent optical particle counter data. Giant particles (〉10 μm) contributed about 89% of the volume of the dust in the 2012 storm. Illite–smectite series clay minerals were the major mineral group followed by quartz, plagioclase, K-feldspar, and calcite. The total phyllosilicate content was ~52%. The direct long-range transport of giant particles was confirmed by calcite nanofibers closely associated with clays in a submicron scale identified by high-resolution SEM and transmission electron microscopy. Since giant particles consisted of clay agglomerates and clay-coated quartz, feldspars, and micas, the mineral composition varied little throughout the fine ( 20 μm) size bins. Analysis of the synoptic conditions of the 2012 dust event and its migration indicated that the mid-tropospheric strong wind belt directly stretching to Korea induced rapid transport of the dust, delivering giant particles. Giant dust particles with high settling velocity would be the major input into the terrestrial and marine sedimentary and ecological systems of East Asia and the western Pacific. Analysis of ancient aeolian deposits in Korea suggested the common deposition of giant particles from Asian dust through the late Quaternary Period. The roles of giant particles should be reviewed with regard to regional circulation of mineral particles and nutrients.

Description: A non-equilibrium wall-model based on unsteady 3D Reynolds-averaged Navier-Stokes (RANS) equations has been implemented in an unstructured mesh environment. The method is similar to that of the wall-model for structured mesh described by Wang and Moin [Phys. Fluids14, 2043–2051 (2002)], but is supplemented by a new dynamic eddy viscosity/conductivity model that corrects the effect of the resolved Reynolds stress (resolved turbulent heat flux) on the skin friction (wall heat flux). This correction is crucial in predicting the correct level of the skin friction. Unlike earlier models, this eddy viscosity/conductivity model does not have a stress-matching procedure or a tunable free parameter, and it shows consistent performance over a wide range of Reynolds numbers. The wall-model is validated against canonical (attached) transitional and fully turbulent flows at moderate to very high Reynolds numbers: a turbulent channel flow at Re τ = 2000, an H-type transitional boundary layer up to Re θ = 3300, and a high Reynolds number boundary layer at Re θ = 31 000. Application to a separated flow over a NACA4412 airfoil operating close to maximum lift is also considered to test the performance of the wall-model in complex non-equilibrium flows.

Description: Conventional shallow water theory successfully reproduces many key features of the Jovian atmosphere: a mixture of coherent vortices and stable, large-scale, zonal jets whose amplitude decreases with distance from the equator. However, both freely decaying and forced-dissipative simulations of the shallow water equations in Jovian parameter regimes invariably yield retrograde equatorial jets, while Jupiter itself has a strong prograde equatorial jet. Simulations by Scott and Polvani [“Equatorial superrotation in shallow atmospheres,” Geophys. Res. Lett.35, L24202 (2008)] have produced prograde equatorial jets through the addition of a model for radiative relaxation in the shallow water height equation. However, their model does not conserve mass or momentum in the active layer, and produces mid-latitude jets much weaker than the equatorial jet. We present the thermal shallow water equations as an alternative model for Jovian atmospheres. These equations permit horizontal variations in the thermodynamic properties of the fluid within the active layer. We incorporate a radiative relaxation term in the separate temperature equation, leaving the mass and momentum conservation equations untouched. Simulations of this model in the Jovian regime yield a strong prograde equatorial jet, and larger amplitude mid-latitude jets than the Scott and Polvani model. For both models, the slope of the non-zonal energy spectra is consistent with the classic Kolmogorov scaling, and the slope of the zonal energy spectra is consistent with the much steeper spectrum observed for Jupiter. We also perform simulations of the thermal shallow water equations for Neptunian parameter values, with a radiative relaxation time scale calculated for the same 25 mbar pressure level we used for Jupiter. These Neptunian simulations reproduce the broad, retrograde equatorial jet and prograde mid-latitude jets seen in observations. The much longer radiative time scale for the colder planet Neptune explains the transition from a prograde to a retrograde equatorial jet, while the broader jets are due to the deformation radius being a larger fraction of the planetary radius.

Description: The hydrodynamic interaction of two deformable vesicles in shear flow induces a net displacement, in most cases an increase of their distance in the transverse direction. The statistical average of these interactions leads to shear-induced diffusion in the suspension, both at the level of individual particles which experience a random walk made of successive interactions, and at the level of suspension where a nonlinear down-gradient diffusion takes place, an important ingredient in the structuring of suspension flows. We make an experimental and computational study of the interaction of a pair of lipid vesicles in shear flow by varying physical parameters, and investigate the decay of the net lateral displacement with the distance between the streamlines on which the vesicles are initially located. This decay and its dependency upon vesicle properties can be accounted for by a simple model based on the well established law for the lateral drift of a vesicle in the vicinity of a wall. In the semi-dilute regime, a determination of self-diffusion coefficients is presented.

Description: [1]  High quality streaming potential coupling coefficient measurements have been carried out using a newly designed cell with both a steady-state methodology and a new pressure transient approach. The pressure transient approach has shown itself to be particularly good at providing high quality streaming potential coefficient measurements as each transient increase or decrease allows thousands of measurements to be made at different pressures to which a good linear regression can be fitted. Nevertheless, the transient method can be up to five times as fast as the conventional measurement approaches because data from all flow rates are taken in the same transient measurement rather than separately. Test measurements have been made on samples of Berea and Boise sandstone as a function of salinity (approximately 18 salinities between 10-5 mol/dm3 and 2 mol/dm3). The data have also been inverted to obtain the zeta potential. The streaming potential coefficient becomes greater (more negative) for fluids with lower salinities, which is consistent with existing measurements. Our measurements are also consistent with the high salinity streaming potential coefficient measurements made by Vinogradov et al. (2010). Both the streaming potential coefficient and the zeta potential have also been modeled using the theoretical approach of Glover et al. (2012). This modeling allows the microstructural, electrochemical and fluid properties of the saturated rock to be taken into account in order to provide a relationship that is unique to each particular rock sample. In all cases, we found that the experimental data was a good match to the theoretical model.

Description: [1]  The present study, which is a follw-up of the JGR paper by Ji et al. (2013a), provides a new calibration for both seismic and fabric properties of antigorite serpentinites. Comparisons of the laboratory velocities of antigorite serpentinites measured at high pressures with crystallographic preferred orientation (CPO) data measured using electron backscatter diffraction (EBCD) techniques demonstrate that seismic anisotropy in high T serpentinite, which is essentially controlled by the antigorite c-axis fabric, is independent on the operating slip system, but strongly dependent on the regime and magnitude of finite strain experienced by the rock. Extrapolation of the experimental data with both pressure and temperature suggests that V p anisotropy decreases but shear wave splitting (Δ V s ) and V p / V s increase with increasing pressure in either cold or hot subduction zones. For a cold, steeply subducting slab, antigorite is most likely deformed by nearly coaxial flattening or trench-parallel movements, forming trench-parallel seismic anisotropy. For a hot, shallowly subducting slab, however, antigorite is most likely deformed by simple shear or transpression. Trench-normal seismic anisotropy can be observed when the subducting dip angle is smaller than 30°. The geophysical characteristics of the Tibetan Plateau such as strong heterogeneity in V p , V s and attenuation, shear wave splitting and electric conductivity may be explained by the presence of strongly deformed serpentinites in lithospheric shear zones reactivated along former suture zones between amalgamated blocks, hydrated zones of subducting lithospheric mantle, and the crust-mantle boundary if the temperature is below 700 °C in the region of interest.

Description: The transport of atmospheric NO x and HNO 3 over Cape Town Atmospheric Chemistry and Physics, 14, 559-575, 2014 Author(s): B. J. Abiodun, A. M. Ojumu, S. Jenner, and T. V. Ojumu Cape Town, the most popular tourist city in Africa, usually experiences air pollution with unpleasant odour in winter. Previous studies have associated the pollution with local emission of pollutants within the city. The present study examines the transport of atmospheric pollutants (NO x and HNO 3 ) over South Africa and shows how the transport of pollutants from the Mpumalanga Highveld, a major South African industrial area, may contribute to the pollution in Cape Town. The study analysed observation data (2001–2008) from the Cape Town air-quality network and simulation data (2001–2004) from a regional climate model (RegCM) over southern Africa. The simulation accounts for the influence of complex topography, atmospheric conditions, and atmospheric chemistry on emission and transport of pollutants over southern Africa. Flux budget analysis was used to examine whether Cape Town is a source or sink for NO x and HNO 3 during the extreme pollution events. The results show that extreme pollution events in Cape Town are associated with the lower level (surface – 850 hPa) transport of NO x from the Mpumalanga Highveld to Cape Town, and with a tongue of high concentration of HNO 3 that extends from the Mpumalanga Highveld to Cape Town along the south coast of South Africa. The prevailing atmospheric conditions during the extreme pollution events feature an upper-level (700 hPa) anticyclone over South Africa and a lower-level col over Cape Town. The anticyclone induces a strong subsidence motion, which prevents vertical mixing of the pollutants and caps high concentration of pollutants close to the surface as they are transported from the Mpumalanga Highveld toward Cape Town. The col accumulates the pollutants over the city. This study shows that Cape Town can be a sink for the NO x and HNO 3 during extreme pollution events and suggests that the accumulation of pollutants transported from other areas (e.g. the Mpumalanga Highveld) may contribute to the air pollution in Cape Town.

Description: [1]  We explore the concept of maximum possible earthquake magnitude, M , in a region represented by an earthquake catalog from the viewpoint of statistical testing. For this aim, we assume that earthquake magnitudes are independent events that follow a doubly-truncated Gutenberg-Richter distribution and focus on the upper truncation M . In earlier work, it has been shown that the value of M cannot be well constrained from earthquake catalogs alone. However, for two hypothesized values M and M ′, alternative statistical tests may address the question: Which value is more consistent with the data? In other words: Is it possible to reject a magnitude within reasonable errors, i.e. the error of the first and the error of the second kind? The results for realistic settings indicate that either the error of the first kind or the error of the second kind is intolerably large. We conclude that it is essentially impossible to infer M in terms of alternative testing with sufficient confidence from an earthquake catalog alone, even in regions like Japan with excellent data availability. These findings are also valid for frequency-magnitude distributions with different tail behavior, e.g. exponential tapering. Finally, we emphasize that different data may only be useful to provide additional constraints for M , if they do not correlate with the earthquake catalog, i.e. if they have not been recorded in the same observational period. In particular, longterm geological assessments might be suitable to reduce the errors, while GPS measurements provide overall the same information as the catalogs.

Description: [1]  We simulated the ascent of bubbly magma in a volcanic conduit by slow decompression experiments using syrup foam as a magma analogue. During decompression, some large voids appear in the foam. The expansion of one void deep in the foam leads to another void expansion, and the void expansion then propagates upward. The void expansion finally reaches the surface of the foam to originate outgassing. The velocity of the upward propagation of void expansions is essentially the same as the rupturing velocity of the bubble film, suggesting that the rupture of films separating each void propagates upward to create the pathway for outgassing. The calculated apparent permeability of decompressed foam can become higher than that measured for natural pumices/scoriae. The upward propagation of film ruptures thus allows for efficient outgassing. This may also appear as the mechanism for energetic gas emissions originating at a depth, such as Strombolian eruptions.

Description: Enhancement of atmospheric H 2 SO 4 / H 2 O nucleation: organic oxidation products versus amines Atmospheric Chemistry and Physics, 14, 751-764, 2014 Author(s): T. Berndt, M. Sipilä, F. Stratmann, T. Petäjä, J. Vanhanen, J. Mikkilä, J. Patokoski, R. Taipale, R. L. Mauldin III, and M. Kulmala Atmospheric H 2 SO 4 / H 2 O nucleation influencing effects have been studied in the flow tube IfT-LFT (Institute for Tropospheric Research – Laminar Flow Tube) at 293 ± 0.5 K and a pressure of 1 bar using synthetic air as the carrier gas. The presence of a possible background amine concentration in the order of 10 7 –10 8 molecule cm −3 throughout the experiments has to be taken into account. In a first set of investigations, ozonolysis of olefins (tetramethylethylene, 1-methyl-cyclohexene, α-pinene and limonene) for close to atmospheric concentrations, served as the source of OH radicals and possibly other oxidants initiating H 2 SO 4 formation starting from SO 2 . The oxidant generation is inevitably associated with the formation of organic oxidation products arising from the parent olefins. These products (first generation mainly) showed no clear effect on the number of nucleated particles within a wide range of experimental conditions for H 2 SO 4 concentrations higher than ~10 7 molecule cm −3 . Also the early growth process of the nucleated particles was not significantly influenced by the organic oxidation products in line with the expected growth by organic products using literature data. An additional, H 2 SO 4 -independent process of particle (nano-CN) formation was observed in the case of α-pinene and limonene ozonolysis for H 2 SO 4 concentrations smaller than ~10 7 molecule cm −3 . Furthermore, the findings confirm the appearance of an additional oxidant for SO 2 beside OH radicals, very likely stabilized Criegee Intermediates (sCI). A second set of experiments has been performed in the presence of added amines in the concentrations range of a few 10 7 –10 10 molecule cm −3 applying photolytic OH radical generation for H 2 SO 4 production without addition of other organics. All amines showed significant nucleation enhancement with increasing efficiency in the order pyridine 〈 aniline 〈 dimethylamine 〈 trimethylamine. This result supports the idea of H 2 SO 4 cluster stabilization by amines due to strong H 2 SO 4 ↔amine interactions. On the other hand, this study indicates that for organic oxidation products (in presence of the possible amine background as stated) a distinct H 2 SO 4 / H 2 O nucleation enhancement can be due to increased H 2 SO 4 formation caused by additional organic oxidant production (sCI) rather than by stabilization of H 2 SO 4 clusters due to H 2 SO 4 ↔organics interactions.

Description: Campholenic aldehyde ozonolysis: a mechanism leading to specific biogenic secondary organic aerosol constituents Atmospheric Chemistry and Physics, 14, 719-736, 2014 Author(s): A. Kahnt, Y. Iinuma, A. Mutzel, O. Böge, M. Claeys, and H. Herrmann In the present study, campholenic aldehyde ozonolysis was performed to investigate pathways leading to specific biogenic secondary organic aerosol (SOA) marker compounds. Campholenic aldehyde, a known α-pinene oxidation product, is suggested to be a key intermediate in the formation of terpenylic acid upon α-pinene ozonolysis. It was reacted with ozone in the presence and absence of an OH radical scavenger, leading to SOA formation with a yield of 0.75 and 0.8, respectively. The resulting oxidation products in the gas and particle phases were investigated employing a denuder/filter sampling combination. Gas-phase oxidation products bearing a carbonyl group, which were collected by the denuder, were derivatised by 2,4-dinitrophenylhydrazine (DNPH) followed by liquid chromatography/negative ion electrospray ionisation time-of-flight mass spectrometry analysis and were compared to the gas-phase compounds detected by online proton-transfer-reaction mass spectrometry. Particle-phase products were also analysed, directly or after DNPH derivatisation, to derive information about specific compounds leading to SOA formation. Among the detected compounds, the aldehydic precursor of terpenylic acid was identified and its presence was confirmed in ambient aerosol samples from the DNPH derivatisation, accurate mass data, and additional mass spectrometry (MS 2 and MS 3 fragmentation studies). Furthermore, the present investigation sheds light on a reaction pathway leading to the formation of terpenylic acid, involving α-pinene, α-pinene oxide, campholenic aldehyde, and terpenylic aldehyde. Additionally, the formation of diaterpenylic acid acetate could be connected to campholenic aldehyde oxidation. The present study also provides insights into the source of other highly functionalised oxidation products (e.g. m / z 201, C 9 H 14 O 5 and m / z 215, C 10 H 16 O 5 ), which have been observed in ambient aerosol samples and smog chamber-generated monoterpene SOA. The m / z 201 and 215 compounds were tentatively identified as a C 9 - and C 10 -carbonyl-dicarboxylic acid, respectively, based on reaction mechanisms of campholenic aldehyde and ozone, as well as detailed interpretation of mass spectral data, in conjunction with the formation of corresponding DNPH derivatives.

Description: Flow climatology for physicochemical properties of dichotomous aerosol over the western North Atlantic Ocean at Bermuda Atmospheric Chemistry and Physics, 14, 691-717, 2014 Author(s): J. L. Moody, W. C. Keene, O. R. Cooper, K. J. Voss, R. Aryal, S. Eckhardt, B. Holben, J. R. Maben, M. A. Izaguirre, and J. N. Galloway Dichotomous aerosols (nominal super- and sub-μm-diameter size fractions) in sectored on-shore flow were sampled daily from July 2006 through June 2009, at the Tudor Hill Atmospheric Observatory (THAO) on the western coast of Bermuda (32.27° N, 64.87° W) and analyzed for major chemical and physical properties. FLEXPART retroplumes were calculated for each sampling period and aerosol properties were stratified accordingly based on transport from different regions. Transport from the northeastern United States (NEUS) was associated with significantly higher (factors of 2 to 3 based on median values) concentrations of bulk particulate non-sea-salt (nss) SO 4 2- , NO 3 - , and NH 4 + and associated scattering and absorption at 530 nm, relative to transport from Africa (AFR) and the oceanic background. These differences were driven primarily by higher values associated with the sub-μm size fraction under NEUS flow. We estimate that 75(±3)% of the NEUS nss SO 4 2- was anthropogenic in origin, while only 25(±9)% of the AFR nss SO 4 2- was anthropogenic. Integrating over all transport patterns, the contribution of anthropogenic sulfate has dropped 14.6% from the early 1990s. Bulk scattering was highly correlated with bulk nss SO 4 2- in all flow regimes but the corresponding regression slopes varied significantly reflecting differential contributions to total scattering by associated aerosol components. Absorption by super-μm aerosol in transport from the NEUS versus AFR was similar although the super-μm aerosol size fraction accounted for a relatively greater contribution to total absorption in AFR flow. Significantly greater absorption Ångström exponents (AAEs) for AFR flow reflects the wavelength dependence of absorption by mineral aerosols; lower AAEs for NEUS flow is consistent with the dominance of absorption by combustion-derived aerosols. Higher AOD associated with transport from both the NEUS and AFR relative to oceanic background flow results in a top of atmosphere direct radiative forcing on the order of −1.6 to −2.5 W m −2 , respectively, showing these aerosols drive cooling. The dominance of transport from the NEUS on an annual basis coupled with the corresponding decreases in anthropogenic nss SO 4 2- aerosols since the early 1990s implies that emission reductions in the US account for a decline in atmospheric cooling over the western North Atlantic Ocean during this period.

Description: The magnetohydrodynamic Richtmyer-Meshkov instability is investigated for the case where the initial magnetic field is unperturbed and aligned with the mean interface location. For this initial condition, the magnetic field lines penetrate the perturbed density interface, forbidding a tangential velocity jump and therefore the presence of a vortex sheet. Through simulation, we find that the vorticity distribution present on the interface immediately after the shock acceleration breaks up into waves traveling parallel and anti-parallel to the magnetic field, which transport the vorticity. The interference of these waves as they propagate causes the perturbation amplitude of the interface to oscillate in time. This interface behavior is accurately predicted over a broad range of parameters by an incompressible linearized model derived presently by solving the corresponding impulse driven, linearized initial value problem. Our use of an equilibrium initial condition results in interface motion produced solely by the impulsive acceleration. Nonlinear compressible simulations are used to investigate the behavior of the transverse field magnetohydrodynamic Richtmyer-Meshkov instability, and the performance of the incompressible model, over a range of shock strengths, magnetic field strengths, perturbation amplitudes and Atwood numbers.

Description: We investigate experimentally the statistical properties of bedload transport induced by a steady, uniform, and laminar flow. We focus chiefly on lateral transport. The analysis is restricted to experiments where the flow-induced shear stress is just above the threshold for sediment transport. We find that, in this regime, the concentration of moving particles is low enough to neglect interactions between themselves. We can therefore represent bedload as a thin layer of independent walkers travelling over the bed surface. In addition to their downstream motion, the particles show significant fluctuations of their cross-stream velocity, likely due to the roughness of the underlying sediment bed. This causes particles to disperse laterally. Based on thousands of individual trajectories, we show that this lateral spreading is the manifestation of a random walk. The experiments are entirely consistent with Fickian diffusion.

Description: Near-surface meteorology during the Arctic Summer Cloud Ocean Study (ASCOS): evaluation of reanalyses and global climate models Atmospheric Chemistry and Physics, 14, 427-445, 2014 Author(s): G. de Boer, M. D. Shupe, P. M. Caldwell, S. E. Bauer, O. Persson, J. S. Boyle, M. Kelley, S. A. Klein, and M. Tjernström Atmospheric measurements from the Arctic Summer Cloud Ocean Study (ASCOS) are used to evaluate the performance of three atmospheric reanalyses (European Centre for Medium Range Weather Forecasting (ECMWF)-Interim reanalysis, National Center for Environmental Prediction (NCEP)-National Center for Atmospheric Research (NCAR) reanalysis, and NCEP-DOE (Department of Energy) reanalysis) and two global climate models (CAM5 (Community Atmosphere Model 5) and NASA GISS (Goddard Institute for Space Studies) ModelE2) in simulation of the high Arctic environment. Quantities analyzed include near surface meteorological variables such as temperature, pressure, humidity and winds, surface-based estimates of cloud and precipitation properties, the surface energy budget, and lower atmospheric temperature structure. In general, the models perform well in simulating large-scale dynamical quantities such as pressure and winds. Near-surface temperature and lower atmospheric stability, along with surface energy budget terms, are not as well represented due largely to errors in simulation of cloud occurrence, phase and altitude. Additionally, a development version of CAM5, which features improved handling of cloud macro physics, has demonstrated to improve simulation of cloud properties and liquid water amount. The ASCOS period additionally provides an excellent example of the benefits gained by evaluating individual budget terms, rather than simply evaluating the net end product, with large compensating errors between individual surface energy budget terms that result in the best net energy budget.