ALBERT

Description: Phase-locked to austral winter and spring, canonical positive Indian Ocean Dipoles (pIODs) generally peak in spring. In recent decades, there has been an increase in unseasonable pIODs which, different from canonical pIODs, peak and decay by September. Distinguishing unseasonable pIODs from canonical pIODs is important, as conditions leading to more frequent unseasonable events are projected to persist in a warming climate. Here using superimposition of the first two seasonally evolving dominant modes of tropical Indian Ocean rainfall variability, we differentiate these types of pIODs. The first mode reflects characteristics of canonical pIODs, in which anomalies intensify with seasonal evolution. However, the second mode, with cool and dry anomalies extending from the eastern pole, reverses from winter to spring, signifying the demise of unseasonable pIODs. Processes embedded in the second mode reflect timing of propagation in equatorial Kelvin waves, and their relative importance to the first mode, in generating different pIODs.

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: 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: 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: Ocean wave activity excites seismic waves that propagate through the solid earth, known as microseismic noise. Here we use a network of 57 ocean-bottom seismometers (OBS) deployed around La Réunion Island in the southwest Indian Ocean to investigate the noise generated in the secondary microseismic band as a tropical cyclone moved over the network. Spectral and polarization analyses show that microseisms strongly increase in the 0.1 - 0.35 Hz frequency band as the cyclone approaches, and that this noise is composed of both compressional and surface waves, confirming theoretical predictions. We infer the location of maximum noise amplitude in space and time, and show that it coincides roughly with the location of maximum ocean-wave interactions. Although this analysis was performed retrospectively, microseisms recorded on the seafloor can be considered a novel source of information for future real-time tracking and monitoring of major storms, complementing atmospheric, oceanographic and satellite observations.

Description: Anthropogenic CO 2 is causing warming and ocean acidification. Coral reefs are being severely impacted, yet confusion lingers regarding how reefs will respond to these stressors over this century. Since the 1982-83 El Niño-Southern Oscillation warming event, the persistence of reefs around the Galápagos Islands has differed across an acidification gradient. Reefs disappeared where pH 〈 8.0 and aragonite saturation state (Ω arag ) ≤ 3 and have not recovered, whereas one reef has persisted where pH 〉 8.0 and Ω arag  〉 3. Where upwelling is greatest, calcification by massive Porites is higher than predicted by a published relationship with temperature despite high CO 2 , possibly due to elevated nutrients. However, skeletal P/Ca, a proxy for phosphate exposure, negatively correlates with density ( R  = -0.822, p  〈 0.0001). We propose that elevated nutrients have the potential to exacerbate acidification by depressing coral skeletal densities and further increasing bioerosion already accelerated by low pH.

Description: Changes to the temporal distribution of daily precipitation were investigated using a dataset of 12,513 land-based stations from the Global Historical Climatology Network (GHCN). The distribution of precipitation was measured using the Gini index (which describes how uniformly precipitation is distributed throughout a year) and the annual number of wet days. The Mann-Kendall test and a regression analysis were used to assess the direction and rate of change to both indices. Over the period of 1976-2000, East Asia, central America and Brazil exhibited a decrease in the number of both wet and light precipitation days, and eastern Europe exhibited a decrease in the number of both wet and moderate precipitation days. In contrast, the USA, southern South America, western Europe and Australia exhibited an increase in the number of both wet and light precipitation days. Trends in both directions were field significant at the global scale.

Description: In the lowermost mantle, seismic velocity variations beneath Pacific margins have been related to the perovskite (pv) to post-perovskite (pPv) phase transition. We investigate the influence of this phase transformation on the geoid using 3D spherical mantle circulation models based on a seismic tomography model and strong lateral viscosity variations in the lower mantle. We demonstrate that the geoid anomalies are strongly affected by the presence of pPv because of phase-dependent viscosity changes relative to the surrounding mantle. Whereas geoid heights above subduction zones are increased for high-viscosity pPv, the presence of weak pPv reduces them, thereby improving the fit to the observed geoid. An investigation using two different tomography models, different pPv density contrasts, and the presence or absence of a global thermal boundary layer and of lateral viscosity variations in the lower mantle demonstrate the various effects of weak pPv on the geoid.

Description: The Interdecadal Pacific Oscillation (IPO) influences multidecadal drought risk across the Pacific, but there are no millennial-length, high resolution IPO reconstructions for quantifying long-term drought risk. In Australia, drought risk increases in positive phases of the IPO, yet few suitable rainfall proxies and short (~100 y) instrumental records mean large uncertainties remain around drought frequency and duration. Likewise, it is unknown whether mega-droughts have occurred in Australia's past. In this study, an atmospheric teleconnection in the Indian Ocean mid-latitudes linking East Antarctica and Australia is exploited to produce the first accurate, annually dated millennial-length IPO reconstruction from the Law Dome (East Antarctica) ice core. Combined with an eastern Australian rainfall proxy from Law Dome, the first millennial-length Australian mega-drought (〉5 y duration) reconstruction is presented. Eight mega-droughts are identified including one 39 y drought (AD 1174–1212), which occurred during an unprecedented century of aridity (AD 1102–1212).

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: Many details of how thunderstorms generate terrestrial gamma-ray flashes (TGFs) and other forms of high energy radiation remain uncertain, including the basic question of where they are produced. We exploit the association of distinct low frequency radio emissions with generation of terrestrial gamma-ray flashes (TGFs) to directly measure for the first time the TGF source altitude. Analysis of two events reveals source altitudes of 11.8 ± 0.4 km and 11.9 ± 0.9 km. This places the source region in the interior of the thunderstorm between the two main charge layers, and implies an intrinsic TGF brightness of approximately 10 18 runaway electrons. The electric current in this non-traditional lightning process is found to be strong enough to drive nonlinear effects in the ionosphere, and in one case is comparable to the highest peak current lightning processes on the planet.

Description: A large volcanic eruption might constitute a climate emergency, significantly altering global temperature and precipitation for several years. Major future eruptions will occur, but their size or timing cannot be predicted. We show, for the first time, that it may be possible to counteract these climate effects through deliberate emissions of short-lived greenhouse gases, dampening the abrupt impact of an eruption. We estimate an emission pathway countering a hypothetical eruption three times the size of Mt Pinatubo in 1991. We use a global climate model to evaluate global and regional responses to the eruption, with and without counter emissions. We then raise practical, financial and ethical question related to such a strategy. Unlike the more commonly-discussed geoengineering to mitigate warming from long-lived greenhouse gases, designed emissions to counter temporary cooling would not have the disadvantage of needing to be sustained over long periods. Nevertheless, implementation would still face significant challenges.

Description: Vertical component data from 206 broadband seismometer stations from Korean networks KIGAM and KMA, the Japanese F-net network, and the Chinese IC and NECESSArray networks are collected for the year 2011 and the ambient seismic noise is analyzed. Rayleigh wave group velocity distribution maps are created in the period range 10 to 70 seconds. Our results are largely consistent with previous studies of the area, but provide greater detail in the Korean peninsula and the East Sea (Sea of Japan). Low group velocities are observed in the Ulleung basin, and the Chubu-Kanto and Kyushu regions in Japan. At 10 s period, sediment basins in the East Sea appear as low group velocity regions relative to higher group velocity continental regions. At periods longer than 40 s, a low group velocity region emerges in the Ulleung basin region, and is bounded by the Korean peninsula.

Description: After more than 4.5 years in orbit, the Gravity field and steady-state Ocean Circulation Explorer (GOCE) mission ended with the reentry of the satellite on 11 November 2013. This publication serves as a reference for the 5th gravity field model based on the time-wise approach (EGM_TIM_RL05), a global model only determined from GOCE observations. Due to its independence of any other gravity data, a consistent and homogeneous set of spherical harmonic coefficients up to degree and order 280 (corresponding to spatial resolution of 71.5 km on ground) is provided including a full covariance matrix characterizing the uncertainties of the model. The associated covariance matrix realistically describes the model quality. It is the first model which is purely based on GOCE including all observations collected during the entire mission. The achieved mean global accuracy is 2.4 cm in terms of geoid heights and 0.7 mGal for gravity anomalies at a spatial resolution of 100 km.

Description: It has been suggested that an increased melting of continental ice in the Amundsen Sea (AS) and Bellinghausen Sea (BS) is a likely source of the observed freshening of Ross Sea (RS) water. To test this hypothesis, we simulate the spreading of glacial melt water using the Finite-Element Sea-ice/ice-shelf/Ocean Model. Based on the spatial distribution of simulated passive tracers, most of the basal melt water from AS ice shelves flows towards the RS with more than half of the melt originating from the Getz Ice Shelf. Further, the model results show that a slight increase of the basal mass loss can substantially intensify the transport of melt water into the RS due to a strengthening of the melt-driven shelf circulation and the westward flowing coastal current. This supports the idea that the basal melting of AS and BS ice shelves is one of the main sources for the RS freshening.

Description: We report on the development of an easily deployable LF near-field Interferometric-TOA 3D Lightning Mapping Array applied to imaging of entire lightning flashes. An interferometric cross-correlation technique is applied in our system to compute windowed two-sensor time differences with sub-microsecond time resolution before TOA is used for source location. Compared to previously reported LF lightning location systems, our system captures many more LF sources. This is due mainly to the improved mapping of continuous lightning processes by using this type of hybrid interferometry/TOA processing method. We show with five station measurements that the array detects and maps different lightning processes, such as stepped and dart leaders, during both in-cloud and cloud-to-ground flashes. Lightning images mapped by our LF system are remarkably similar to those created by VHF mapping systems, which may suggest some special links between LF and VHF emission during lightning processes.

Description: We estimate the snow albedo forcing and direct radiative forcing (DRF) of black carbon (BC) in the Tibetan Plateau using a global chemical transport model in conjunction with a stochastic snow model and a radiative transfer model. The annual mean BC snow albedo forcing is 2.9 W m −2 averaged over snow-covered Plateau regions, which is a factor of three larger than the value over global land snowpack. BC-snow internal mixing increases the albedo forcing by 40-60% compared with external mixing and coated BC increases the forcing by 30-50% compared with uncoated BC aggregates, whereas Koch snowflakes reduce the forcing by 20-40% relative to spherical snow grains. The annual BC DRF at the top of the atmosphere is 2.3 W m −2 with uncertainties of −70-85% in the Plateau after scaling the modeled BC absorption optical depth to Aerosol Robotic Network (AERONET) observations. The BC forcings are attributed to emissions from different regions.

Description: ABSTRACT Basal melt is a major cause of ice shelf thinning affecting the stability of the ice shelf and reducing its buttressing effect on the inland ice. The Fimbul ice shelf (FIS) in Dronning Maud Land (DML), East Antarctica, is fed by the fast-flowing Jutulstraumen glacier, responsible for 10% of ice discharge from the DML sector of the ice sheet. Current estimates of the basal melt rates of the FIS come from regional ocean models, autosub measurements and satellite observations, which vary considerably. This discrepancy hampers evaluation of the stability of the Jutulstraumen catchment. Here, we present estimates of basal melt rates of the FIS using ground based interferometric radar. We find a low average basal melt rate on the order of 1 m/yr, with the highest rates located at the ice shelf front, which extends beyond the continental shelf break. Furthermore, our results provide evidence for a significant seasonal variability.

Description: A significant source of ozone in the troposphere is transport from the stratosphere. The stratospheric contribution has been estimated mainly using global models that attribute the transport process largely to the global scale Brewer-Dobson circulation and synoptic scale dynamics associated with upper tropospheric jet streams. We report observations from research aircraft that reveal additional transport of ozone-rich stratospheric air downward into the upper troposphere by a leading-line-trailing-stratiform (LLTS) mesoscale convective system (MCS) with convection overshooting the tropopause altitude. The fine-scale transport demonstrated by these observations poses a significant challenge to global models that currently do not resolve storm scale dynamics. Thus the upper tropospheric ozone budget simulated by global chemistry-climate models where large-scale dynamics and photochemical production from lightning-produced NO are the controlling factors may require modification.

Description: Near simultaneous periodic dispersive features of fast magnetosonic mode emissions are observed by both Van Allen Probes spacecraft while separated in magnetic local time by ~5 hours: Probe A at 15 and Probe B at 9–11 hours. Both spacecraft see similar frequency features, characterized by a periodic repetition at ~180 s. Each repetition is characterized by a rising frequency. Since no modulation is observed in the proton shell distribution, the plasma density, or in the background magnetic field at either spacecraft we conclude that these waves are not generated near the spacecraft but external to both spacecraft locations. Probe A while outside the plasmapause sees the start of each repetition ~40 s before probe B while deep inside the plasmasphere. We can qualitatively reproduce the dispersive features, but not the quantitative details. The cause for this phenomena remains to be identified.

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: Gas hydrates, pervasive in continental margin sediments, are expected to release methane in response to ocean warming, but the geographic range of dissociation and subsequent flux of methane to the ocean are not well constrained. Sediment column thermal models based on observed water column warming trends offshore Washington (USA) show that a substantial volume of gas hydrate along the entire Cascadia upper continental slope is vulnerable to modern climate change. Dissociation along the Washington sector of the Cascadia margin alone has the potential to release 45–80 Tg of methane by 2100. These results highlight the importance of lower latitude warming to global gas hydrate dynamics and suggest that contemporary warming and downslope retreat of the gas hydrate reservoir occurs along a larger fraction of continental margins worldwide than previously recognized.

Description: The jet stream over the eastern North Pacific (ENPJ) as a core of the atmospheric flow is known to strongly fluctuate meridionally, and its meridional displacement directly influences adjacent regional climate. Here, we investigate how this jet will be changed due to global warming. By analyzing the future scenario experiments of CMIP3 and 5, and found that both ENPJ and the eastern tropical Pacific ITCZ tend to move southward, which are closely related to the tropical eastern Pacific warming trend. Tropical eastern Pacific warming leads to not only the southward migration of ITCZ by southward-shifting the off-equatorial eastern Pacific warm pool, but also the southward shift of ENPJ by increasing baroclinic instability of the atmosphere in subtropical region through intensifying the meridional SST gradient. Not primary but yet secondly the southward shift of ITCZ contributes to the southward shift of ENPJ through a kinematic connection bridged by local Hadley circulation.

Description: Groundwater recharge affects water budgets and groundwater quality on the deltas and floodplains of South and Southeast Asia. Rain and flooding rivers recharge groundwater during the monsoon; irrigated rice fields and surface-water bodies recharge aquifers during the dry season. Groundwater throughout the region is severely contaminated by arsenic and recent research suggests that quantifying and characterizing recharge is important to understand whether recharge flushes or mobilizes arsenic from aquifers. At a field site in Bangladesh, we found that burrows of terrestrial crabs short-circuit low-permeability surface sediments, providing the primary conduit for recharge. We combine field observations along with a model that couples isotope and water balances to quantify the effect of crab burrows on aquifer recharge. Given the wide distribution of burrowing crabs and the surficial geology, we suggest that crab burrows provide widespread conduits for groundwater recharge.

Description: We demonstrate that broadband low frequency electromagnetic field fluctuations embedded within fast flows throughout the Earth's plasma sheet may drive significant ion heating. This heating is nearly entirely in the direction perpendicular to the background magnetic field and is estimated to occur at an average rate of ~1 eV/s with rates in excess of 10 eV/s within one standard deviation of the average value over all observed events. For an Earthward flow the total change in temperature along a flow path may exceed one keV and for ‘wave-rich’ flows can be comparable to that expected due to conservation of the first adiabatic invariant. The consequent increase in plasma pressure and flux tube entropy may lead to braking of inward motion and the suppression of plasma interchange.

Description: Along the continental margins, rivers and submarine groundwater supply nutrients, trace elements, and radionuclides to the coastal ocean, supporting coastal ecosystems and, increasingly, causing harmful algal blooms and eutrophication. While the global magnitude of gauged riverine water discharge is well known, the magnitude of submarine groundwater discharge (SGD) is poorly constrained. Using an inverse model combined with a global compilation of 228 Ra observations, we show that the SGD integrated over the Atlantic and Indo-Pacific Oceans between 60°S and 70°N is (12 ± 3) x 10 13  m 3  yr -1 , which is 3 to 4 times greater than the freshwater fluxes into the oceans by rivers. Unlike the rivers, where more than half of the total flux is discharged into the Atlantic, about 70% of SGD flows into the Indo-Pacific Oceans. We suggest that SGD is the dominant pathway for dissolved terrestrial materials to the global ocean, and this necessitates revisions for the budgets of chemical elements including carbon.

Description: Fast chemical reactions in geophysical flows are controlled by fluid mixing, which perturbs local chemical equilibria and thus triggers chemical reactions. Spatial fluctuations in the flow velocity lead to deformation of material fluid elements, which form the support volumes of transported chemical species. We develop an approach based on a lamellar representation of fluid mixing that provides a direct link between fluid deformation, the distribution of concentration gradiens, and the upscaled reaction rates for fast reversible reactions. The temporal evolution of effective reaction rates are determined by the flow topology and the distribution of local velocity gradients. This leads to a significant increase of the reaction efficiency, which turns out to be orders of magnitude larger than in homogeneous flow. This approach allows for the systematic evaluation of the temporal evolution of equilibrium reaction rates, and establishes a direct link between the reaction efficiency and the spatial characteristics of the underlying flow field as quantified by the deformation of material fluid elements.

Description: The 1 April 2014 Iquique M w 8.1 earthquake ruptured a segment of the megathrust fault offshore of northern Chile and generated a moderate-size tsunami across the Pacific. Tide gauges in Hawaii recorded over 1 m of wave height despite the long distance from the source and position away from the main radiated energy lobe. Inversion of global teleseismic body wave observations combined with forward modeling of the tsunami at four near-field DART stations arrives iteratively at a self-consistent finite-fault model with very compact dimensions. The slip distribution produces a NNE-SSW trending seafloor uplift patch that enhances the tsunami directionality in the WNW, resulting in good matches to observed DART and tide gauge records around the Hawaiian Islands. The relatively large waves at selected locations in Hawaii can be attributed to a combination of the spatial slip distribution and the resulting short-period waves that triggered localized resonance over the insular shelves. This event highlights the importance of characterizing detailed slip distributions in analysis or forecasting of tsunamis even for a compact source.

Description: The first decade of the twenty-first century was characterised by a hiatus in global surface warming. Using ocean model hindcasts and reanalyses we show that heat uptake between the 1990s and 2000s increased by 0.7 ± 0.3 Wm −2 . Approximately 30% of the increase is associated with colder sea surface temperatures in the eastern Pacific. Other basins contribute via reduced heat loss to the atmosphere, in particular the Southern and subtropical Indian Oceans (30%), and the subpolar North Atlantic (40%). A different mechanism is important at longer timescales (1960s-present) over which the Southern Annular Mode trended upwards. In this period, increased ocean heat uptake has largely arisen from reduced heat loss associated with reduced winds over the Agulhas Return Current and southward displacement of Southern Ocean westerlies.

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: The EPICA Dome C (EDC) ice core has allowed for the reconstruction of atmospheric CO 2 concentrations for the last 800,000 years. Here, we revisit the oldest part of the EDC CO 2 record using different air extraction methods and sections of the core. For our established cracker system, we found an analytical artifact, which increases over the deepest 200 m and reaches 10.1 ± 2.4 ppm in the oldest/deepest part. The governing mechanism is not yet fully understood, but it is related to insufficient gas extraction in combination with ice relaxation during storage and ice structure. The corrected record presented here resolves partly the issue with a different correlation between CO 2 and Antarctic temperatures found in this oldest part of the records, however, this anomaly still persists. In addition we provide here an update of 800,000 years atmospheric CO 2 history including recent studies covering the last glacial cycle.

Description: Equivalent effective stratospheric chlorine (EESC) construct of ozone regression models attributes ozone changes to EESC changes using a single value of the sensitivity of ozone to EESC over the whole period. Using space-based total column ozone (TCO) measurements, and a synthetic TCO time series constructed such that EESC does not fall below its late 1990s maximum, we demonstrate that the EESC-based estimates of ozone changes in the polar regions (70–90°) after 2000 may, falsely, suggest an EESC-driven increase in ozone over this period. An EESC-based regression of our synthetic “failed Montreal Protocol with constant EESC" time series suggests a positive TCO trend that is statistically significantly different from zero over 2001–2012 when in fact no recovery has taken place. Our analysis demonstrates that caution needs to be exercised when using explanatory variables, with a single fit coefficient, fitted to the entire data record, to interpret changes in only part of the record.

Description: We present a statistical study of the temporal and spatial scale characteristics of different field-aligned current (FAC) types derived with the Swarm satellite formation. We divide FACs into two classes: small-scale, up to some ten kilometer, which arecarried predominantly by kinetic Alfvén waves, and large-scale FACs with sizes of more than 150 km. For determining temporal variability we consider measurements at the same point, the orbital crossovers near the poles, but at different times. Fromcorrelation analysis we obtain a persistent period of small-scale FACs of order 10 s, while large-scale FACs can be regarded stationary for more than 60 s. For the first time we investigate the longitudinal scales. Large-scale FACs are different on dayside and nightside. On the nightside the longitudinal extension is on average 4 times the latitudinal width, while on the dayside, in particular in the cusp region, latitudinal and longitudinal scales are comparable.

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: Although the relative velocity of a single crystal or bubble in a quiescent fluid (melt) is well-characterized, the interplay of crystals/bubbles in multiparticle systems and its effect on their settling/rising velocity is poorly quantified. We propose a theoretical model for the hindered velocity of non-Brownian emulsions and suspensions of non-deformable fluid and solid particles in the creeping flow regime. The model is based on three sets of correction; two on the drag coefficient experienced by each particle to account for both return flow and Smoluchowski effects, and a correction on the rheology to account for non-local interactions introduced as a mean-field effective viscosity. Our model is tested against new and published experimental data over a wide range of particle volume fraction and viscosity ratio between the fluids. We find an excellent agreement between our model and experiments. The model is then applied to show that hindered settling can increase mineral residence time by up to an order of magnitude in convecting magma chambers.

Description: 14 C content in tree rings and 10 Be concentration records in polar ice core provide information about past cosmic ray intensities. The AD 774–775 cosmic ray event has been identified by 14 C measurement in several tree rings from all over the world. Although the quasi-decadal 10 Be Dome Fuji data in the Antarctic ice core also shows a sharp peak around AD 775, annual 10 Be variations in the Dome Fuji core or in other cores have not been revealed. We have measured quasi-annual 10 Be concentrations from approximately AD 763–794 in the Dome Fuji ice core, and detected a clear increase (~80% above the baseline) in 10 Be concentration around AD 775. However, an accurate height of this increase is not straightforwardly estimated due to the background variation in 10 Be concentration. The 10 Be increase can be due to the same cosmic ray event as shown in the 14 C content in AD 774–775.

Description: For the first time a mesoscale-resolving whole atmosphere general circulation model (GCM) has been developed, using the NCAR Whole Atmosphere Community Climate Model (WACCM) with ~0.25° horizontal resolution and 0.1 scale height vertical resolution above the middle stratosphere (higher resolution below). This is made possible by the high accuracy and high scalability of the spectral element dynamical core from the High-Order Method Modeling Environment (HOMME). For the simulated January-February period, the latitude-height structure and the magnitudes of the temperature variance compare well with those deduced from SABER observations. The simulation reveals the increasing dominance of gravity waves (GWs) at higher altitudes through both the height dependence of the kinetic energy spectra, which display a steeper slope (~-3) in the stratosphere and an increasingly shallower slope above, and the increasing spatial extent of GWs (including a planetary-scale extent of a concentric GW excited by a tropical cyclone) at higher altitudes. GW impacts on the large-scale flow is evaluated in terms of zonal mean zonal wind and tides: with no GW drag parameterized in the simulations, forcing by resolved GWs does reverse the summer mesospheric wind, albeit at an altitude higher than climatology, and only slows down the winter mesospheric wind without closing it. The hemispheric structures and magnitudes of diurnal and semidiurnal migrating tides compare favorably with observations.

Description: We demonstrate via numerical simulations that the impact of a ~ lunar-sized body with Mars is capable of creating a hemispherical magma ocean that upon cooling and solidification resulted in the formation of the southern highlands and thus the Martian dichotomy. The giant impact may have contributed a significant amount of iron to the Martian core and generated a deep thermal anomaly that led to the onset and development of the volcanism in the southern highlands. Our model predicts several mantle plumes converging to the South Pole from the equatorial regions as well as new plumes forming in the equatorial region and also an absence of significant large-scale volcanism in the northern lowlands. The core heat flux evolution obtained from our numerical models is consistent with the decline of the magnetic field. We argue that such a scenario is more consistent with a range of observations than a northern giant impact (excavating the Borealis basin) for the formation of the Martian dichotomy.

Description: We present ice velocities observed with global positioning systems and TerraSAR-X/TanDEM-X in a land-terminating region of the south-west Greenland ice sheet (GrIS) during the melt-year 2012–2013, to examine the spatial pattern of seasonal and annual ice motion. We find that whilst spatial variability in the configuration of the subglacial drainage system controls ice motion at short timescales, this configuration has negligible impact on the spatial pattern of the proportion of annual motion which occurs during summer. Whilst absolute annual velocities vary substantially, the proportional contribution of summer motion to annual motion does not. These observations suggest that in land-terminating margins of the GrIS, subglacial hydrology does not significantly influence spatial variations in net summer speedup. Furthermore, our findings imply that not every feature of the subglacial drainage system needs to be resolved in ice sheet models.

Description: Measurements of particulate δ 15 N in coastal marine laminated sediments provide a high-resolution proxy for fluctuations in the intensity of denitrification in the water column. In the Eastern Tropical North Pacific oxygen minimum zone, this denitrification signal is transported northward by the California Undercurrent, thus serving as a tracer of ocean circulation. This is verified through comparisons between salinity in the thermocline off Southern California (Santa Monica Basin) and the difference between δ 15 N sed within age equivalent sediments from a southern (Pescadero Slope) and northern (Santa Monica Basin) site. Trends in this parameter, Δ δ 15 N sed , relate to Pacific Decadal Oscillation (PDO) phase changes between 1900 and 1990. We hypothesize that the decline in Δ δ 15 N sed during warm PDO phases is due to a strengthening of the California Undercurrent transporting 15  N-enriched nitrate from the ETNP northward. The deviation from this trend after 1990 suggests recent changes in circulation and/or California Current water nutrient biogeochemistry.

Description: Natural climate variations in the United States wind resource are assessed by using cyclostationary empirical orthogonal functions (CSEOFs) to decompose wind reanalysis data. Compared to approaches that average climate signals or assume stationarity of the wind resource on interannual timescales, the CSEOF analysis isolates variability associated with specific climate oscillations, as well as their modulation from year to year. Contributions to wind speed variability from the modulated annual cycle (MAC) and the El Niño-Southern Oscillation (ENSO) are quantified, and information provided by the CSEOF analysis further allows the spatial variability of these effects to be determined. The impacts of the MAC and ENSO on the wind resource are calculated at existing wind turbine locations in the United States, revealing variations in the wind speed of up to 30% at individual sites. The results presented here have important implications for predictions of wind plant power output and siting.

Description: Understanding the spatial variation of anisotropy in the upper mantle is important for characterizing the lithospheric deformation and mantle flow dynamics. In this study, we apply a full-wave approach to image the upper-mantle anisotropy in Southern California using 5954 SKS splitting data. Three-dimensional sensitivity kernels combined with a wavelet-based model parameterization are adopted in a multiscale inversion. Spatial resolution lengths are estimated based on a statistical resolution matrix approach, showing a finest resolution length of ~25 km in regions with densely-distributed stations. The anisotropic model displays structural fabric in relation to surface geologic features such as the Salton Trough, the Transverse Ranges and the San Andreas Fault. The depth variation of anisotropy does not suggest a lithosphere-asthenosphere decoupling. At long-wavelengths, the fast directions of anisotropy are aligned with the absolute plate motion inside the Pacific and North American plates.

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: Climate extremes such as droughts and heatwaves affect terrestrial ecosystems and may alter local carbon budgets. However, it still remains uncertain to what degree extreme impacts in the carbon cycle influence the carbon cycle-climate feedback both today and the near future. Here, we analyze spatiotemporally contiguous negative extreme anomalies in gross primary production (GPP) and net ecosystem production (NEP) in model output of the Coupled Model Intercomparison Project Phase 5 (CMIP5) ensemble and investigate their future development and attribution to climatic drivers. We find that relative to the overall increase in global carbon uptake, negative extremes in GPP and NEP lose importance towards the end of the 21st century. This effect can be related to elevated CO 2 concentrations and higher amounts of available water at the global scale, partially mitigating the impacts of droughts and heatwaves, respectively. Overall, based on CMIP5 models we hypothesize that terrestrial ecosystems might be more resilient against future climate extremes than previously thought. Future work will have to further scrutinize these results considering that various biological and biogeochemical feedbacks are not yet integrated within earth system models.

Description: The mass eruption rate feeding a volcanic plume is commonly estimated from its maximum height. Winds are known to affect the column dynamics causing bending and hence reducing the maximum plume height for a given mass eruption rate. However, the quantitative predictions including wind effects on mass eruption rate estimates are not well constrained. To fill this gap, we present a series of new laboratory experiments on forced plumes rising in a density-stratified crossflow. We identify three dynamical regimes corresponding to increasing effect of wind on the plume rise. The transition from one regime to another is governed by two dimensionless velocity scales defined as a function of source and environmental parameters. The results are found consistent with the conditions of historical eruptions and provide new empirical relationships to estimate mass eruption rate from plume height in windy conditions, leading to valuable tools for eruption risk assessment.

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: [1]  Whistler mode chorus is an important magnetospheric emission, playing a dual role in the acceleration and loss of relativistic electrons in the Earth's outer radiation belt. Chorus is typically generated in the equatorial region in the frequencyrange 0.1 - 0.8 f ce , where f ce is the local electron gyrofrequency. However, as the waves propagate to higher latitudes, significant wave power can occur at frequencies below 0.1 f ce . Since this wave power is largely omitted in current radiation belt models we construct a global model of low frequency chorus, f LHR  〈  f  〈 0.1 f ce , using data from six satellites. We find that low frequency chorus is strongest, with an average intensity of 200 pT 2 , in the pre-noon sector during active conditions at mid latitudes (20 o  〈 | λ m | 〈 50 o ) from 4 〈  L ∗  〈 8. Such mid-latitude, low frequency chorus wave power will contribute to the acceleration and loss of relativistic electrons and should be taken into account in radiation belt models.

Description: [1]  Megaripples are distinguished from regular ripples by their larger size and bimodal sediment distribution. The interplay between wind, grain size and morphology controls their development, but the exact mechanisms that limit the size of megaripples have been unclear. Using wind tunnel experiments we found two main mechanisms that limit the height of megaripples. The first mechanism is megaripple flattening due to strong enough winds that drive the coarse grains into saltation; the second mechanism is megaripple deflation by impacts of faster saltation grains. In this latter mechanism, the coarse grains are propelled by the impacts of fine saltating grains. The occurrence of both these mechanisms depends on the grain size distribution, and increases with both megaripple height and wind speed. Thus, for a given wind environment and grain size distribution, there exists a limit on the size of megaripples, which is determined by these two mechanisms.

Description: [1]  Very low frequency (VLF) remote sensing observations of multifaceted local and conjugate ionospheric perturbations from geographically identified and well characterized oceanic lightning discharges are presented for the first time. Lightning-induced electron precipitation (LEP) events are shown to produce disturbances first in the conjugate hemisphere and subsequently in the hemisphere of the causative lightning discharge in agreement with theoretical predictions. A rough threshold peak current of ~100 kA is identified for lightning discharges to generate LEP events for the geomagnetic conditions present during observations. The occurrence of early VLF events and the subsequent duration of their recovery does not seem to fit any simple metric of lightning discharge peak current or proximity to great circle path. Knowledge of the full spectral density of the lightning EMP not just its peak current and the subionospheric mode structure are likely necessary to determine if a specific lightning discharge will generate an early VLF perturbation.

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 have investigated the influence that megathrust earthquake slip has on the activation of splay faults using a 2D Finite Element Method (FEM), taking into account the effects of gravity and variations in the frictional strength properties of splay faults. We simulated both landward-dipping and seaward-dipping splay fault geometries, and imposed depth-variable slip distributions of subduction events. Our results indicate that the two types of splay fault exhibit a similar behavior, with variations in frictional properties along the faults affecting only the seismic magnitude. The triggering process is controlled by a critical depth. Megathrust slip concentrated at depths shallower than the critical depth will favor normal displacement, while megathrust slip concentrated at depths deeper than the critical depth is likely to result in reverse motion. Our results thus provide a useful tool for predicting the activation of secondary faults and may have direct implications for tsunami hazard research.

Description: [1]  Atmospheric low frequency sound, i.e., infrasound, from underwater events has not been considered thus far, due to the high impedance contrast of the water-air interface making it almost fully reflective. Here, we report for the first time on atmospheric infrasound from a large underwater earthquake (Mw 8.1) near the Macquarie Ridge, which was recorded at 1,325 km from the epicenter. Seismic waves coupled to hydro-acoustic waves at the ocean floor, after which the energy entered the SOund Fixing And Ranging (SOFAR) channel and was detected on a hydrophone array. The energy was diffracted by a sea mount and an oceanic ridge, which acted as a secondary source, into the water column followed by coupling into the atmosphere. The latter results from evanescent wave coupling and the attendant anomalous transparency of the sea surface for very low frequency acoustic waves.

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: [1]  We present surface velocity measurements from a high-elevation site located 140 km from the western margin of the Greenland Ice Sheet, and ~ 50 km into its accumulation area. Annual velocity increased each year from 51.78 ± 0.01 m yr -1 in 2009 to 52.92 ± 0.01 m yr -1 in 2012 – a net increase of 2.2 %. These data also reveal a strong seasonal velocity cycle of up to 8.1 % above the winter mean, driven by seasonal melt and supraglacial lake drainage. Sole et al. (2013) recently argued that ice motion in the ablation area is mediated by reduced winter flow following the development of efficient subglacial drainage during warmer, faster, summers. Our data complement this analysis and reveal a year-on-year increase in annual velocity above the ELA, where despite surface melt increasing, it is still sufficiently low to hinder the development of efficient drainage under thick ice.

Description: [1]  The dual-spacecraft Van Allen Probes mission has provided a new window into megaelectron Volt (MeV) particle dynamics in the Earth's radiation belts. Observations (up to E ~10 MeV) show clearly the behavior of the outer electron radiation belt at different time scales: months-long periods of gradual inward radial diffusive transport and weak loss being punctuated by dramatic flux changes driven by strong solar wind transient events. We present analysis of multi-MeV electron flux and phase space density (PSD) changes during March 2013 in the context of the first year of Van Allen Probes operation. This March period demonstrates the classic signatures both of inward radial diffusive energization as well as abrupt localized acceleration deep within the outer Van Allen zone (L ~4.0 ± 0.5). This reveals graphically that both “competing” mechanisms of multi-MeV electron energization are at play in the radiation belts, often acting almost concurrently or at least in rapid succession.

Description: [1]  High-resolution physical, geochemical and geochronological analyses on the sedimentary sequence of Yeniçağa Lake, located in a fault-bounded basin along the North Anatolian Fault, reveal fingerprints of paleoearthquakes. A robust sediment chronology, spanning the last 3400 years, is constructed by radiocarbon dating and time-stratigraphical correlation with the precisely dated Sofular Cave speleothem record. Yeniçağa sedimentary sequence contains eleven seismically induced event deposits characterized by siliciclastic-enriched intervals. Some of the event deposits are also associated with implications of sudden lake deepening, which may be related to co-seismic subsidence. The paleoearthquake series having an average recurrence interval of ca. 260 years are interrupted by two possible seismic gaps of ca. 420 and 540 years.

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]  The annual runoff from the melting of large glaciers and snow fields along the northern perimeter of the Gulf of Alaska is a critical component of marine physical and biological systems; yet, most of this freshwater is not measured. Here we show estimates of melt for the watershed that contains the largest and longest glacier in North America, the Bering Glacier. The procedure combines in-situ observations of snow and ice melt acquired by a long-term monitoring program, multi-spectral satellite observations and nearby temperature measurements. The estimated melt is 40 km 3 per melt season, ± 3.0 km 3 , observed over the decadal period, 2002-2012. As a result of climate change these estimates could increase to 60 km 3 per year by 2050. This technique and the derived melt coefficients can be applied to estimate melt from Alaska to Washington glaciers.

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: [1]  This study emphasizes the importance of rainstorm events in mobilizing carbon at the soil-stream interface from tropical rainforests. Half-hourly geochemical/isotopic records over a 13.5 hour period from a 20 km 2 tropical rainforest headwater in Guyana show an order of magnitude increase in dissolved organic carbon (DOC) concentration in less than 30 minutes (10.6-114 mg/L). The composition of DOC varies significantly and includes optically invisible dissolved organic matter (iDOM) that accounts for a large proportion (4-89%) of the total DOC, quantified using size exclusion chromatography (SEC). SEC suggests that iDOM is comprised of low molecular weight organic moieties, which are likely sourced from fresh leaf litter and/or topsoil, as shown in soils from the surrounding environment. Although poorly constrained at present, the presence of iDOM further downstream during the wet season suggests that this organic matter fraction may represent an un-quantified source of riverine CO 2 outgassing in tropical headwaters, requiring further consideration.

Description: [1]  The 2004 Mw 8.1 event on 2004 December 23 near the Macquarie Ridge is a very large intraplate event that has been overshadowed by the Mw 9.3 Sumatra-Andaman event only three days later. We are able to track the progress of source evolution by estimating the progression of the points of energy emission, exploiting the good azimuthal distribution of available stations. The results indicate that this event ruptured on two nearby fault systems reactivating fossil fracture zones, with the second sub-event to the west triggered by the first. The total duration of high-frequency radiation is quite short, about 60 s, for such a large event. Much of the high-frequency radiation occurs on a fault sub-parallel to that inferred from long-period studies. This composite fault behaviour with displaced triggered failure appears to be a characteristic of large intraplate events in the oceans.

Description: [1]  The Brazilian Northeast (NE) is strongly affected by ENSO. During La Niña events, the precipitation over the NE is generally above average. However, during the last La Niña event in 2011-12, the NE went through its worst drought in the last 30 years. In this study, observations and numerical simulations are used to determine what made the 2011-12 event different from other events. We find that eastern Pacific (canonical) La Niña events cause a cooling of the tropical North Atlantic and warming of the tropical South Atlantic that lead to a southward migration of the ITCZ, which in turn brings rain to the NE. On the other hand, La Niña events with the cooling concentrated in the central Pacific cause the opposite meridional SST gradient in the tropical Atlantic, leading to droughts over the NE. The 2011-12 event was of the latter type. This study also shows that it is possible to predict the sign of the NE rainfall anomaly during ENSO events using a simple SST index.

Description: [1]  The impact of ice clouds on weather and climate is a function of ice particle shape through light scattering properties and cloud lifetime through ice particle sedimentation rates. Many weather forecast and climate models use two categories to represent ice cloud particles: cloud ice, and snow, though the distinction between particle categories is generally without observational justification. Improved characterization of cloud ice and snow as well as the transition between them will make models more realistic. An analysis of particle imagery data from high resolution aircraft particle imaging probes indicates that atmospheric ice particles can easily be separated by particle complexity. In this work, a technique is described which enables the clear separation of vapor grown particles from aggregates of particles. When applied to two example datasets, the technique shows that the separation between these categories occurs at 150 and 250 microns, for two example datasets.

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: 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: 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: 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: [1]  We present the first satellite-detected perturbations of the outgoing longwave radiation (OLR) associated with blowing snow events over the Antarctic ice sheet using data from CALIOP and CERES. Significant cloud-free OLR differences are observed between the clear and blowing snow sky, with the sign and magnitude depending on season and time of the day. During nighttime, OLRs are usually larger when blowing snow is present; the average difference in OLRs between without and with blowing snow over the East Antarctica ice sheet is about -5.2 W/m 2 for the winter months of 2009. During daytime, in contrast, the OLR perturbation is usually smaller or even has the opposite sign. The observed seasonal variations and day-night differences in the OLR perturbation are consistent with theoretical calculations of the influence of blowing snow on OLR. Detailed atmospheric profiles are needed to quantify the radiative effect of blowing snow from the satellite observations.

Description: [1]  Data from the Coupled Model Intercomparison Project Phase 5 (CMIP5) for historical and future climate scenarios are examined for changes in the energy cycle component of land surface feedbacks on the atmosphere, namely through the linkages from soil moisture to sensible heat flux to the height of the lifting condensation level marking the cloud base. Climate models project heightened sensitivity in both of these segments of the feedback pathway over most of the globe. This is in agreement with studies showing similar increases in land-atmosphere feedbacks through the water cycle, despite different physical processes, and may contribute to prevalent droughts and floods found in most climate change forecasts.

Description: [1]  The El Niño-Southern Oscillation (ENSO) is the leading mode of interannual variability, with global impacts on weather and climate that have seasonal predictability. Research on the link between interannual ENSO variability and the leading mode of intraseasonal variability, the Madden-Julian Oscillation (MJO), has focused mainly on the role of MJO initiating or terminating ENSO. We use observational analysis and modeling to show that the MJO has an important simultaneous link to ENSO: strong MJO activity significantly weakens the atmospheric branch of ENSO. For weak MJO conditions relative to strong MJO conditions, the average magnitude of ENSO-associated tropical precipitation anomalies increases by 63% and the strength of hemispheric teleconnections increases by 58%. Since the MJO has predictability beyond three weeks, the relationships shown here suggest that there may be subseasonal predictability of the ENSO teleconnections to continental circulation and precipitation.

Description: [1]  Measurements taken by the Gravity Recovery And Climate Experiment (GRACE) satellites indicated a continued water storage increase over the Missouri River Basin (MRB) prior to the 2011 flood event. An analysis of the major hydrologic variables in the MRB, i.e. those of soil moisture, streamflow, groundwater storage and precipitation, show a marked variability at the 10-15 year timescale coincident with the water storage increase. A climate diagnostic analysis was conducted to determine what climate forcing conditions preceded the long-term changes in these variables. It was found that precipitation over the MRB undergoes a profound modulation during the transition points of the Pacific Quasi-Decadal Oscillation (QDO) and associated teleconnections. The results infer a prominent teleconnection forcing in driving the wet/dry spells in the MRB, and this connection implies the potential for climate prediction of future wet/dry extreme events.

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: Freshwater in the Arctic Ocean plays an important role in the regional ocean circulation, sea ice and global climate. From salinity observed by a variety of platforms we are able, for the first time, to estimate a statistically reliable liquid freshwater trend from monthly gridded fields over all upper Arctic Ocean basins. From 1992 to 2012 this trend was 600 ± 300  km 3 yr − 1 . A numerical model agrees very well with the observed freshwater changes. A decrease in salinity made up about 2/3 of the freshwater trend and a thickening of the upper layer 1/3. The Arctic Ocean Oscillation index, a measure for the regional wind stress curl, correlated well with our freshwater timeseries. No clear relation to Arctic Oscillation or Arctic Dipole indices could be found. Following other observational studies, an increased Bering Strait freshwater import to the Arctic Ocean, a decreased Davis Strait export and enhanced net sea ice melt could have played an important role in the freshwater trend we observed.

Description: [1]  The timing and evolution of Jabat Island, Marshall Islands, was investigated using morphostratigraphic analysis and radiometric dating. Results show the first evidence of island building in the Pacific during latter stages of Holocene sea-level rise. A three-phase model of development of Jabat is presented. Initially, rapid accumulation of coarse sediments on Jabat occurred 4,800-4,000 yrs BP across a reef flat higher than present level, as sea level continued to rise. During the highstand, island margins and particularly the western margin, accreted vertically to 2.5-3.0 m above contemporary ridge elevations. This accumulation phase was dominated by sand-size sediments. Phase three involved deposition of gravel ridges on the northern reef, as sea level fell to present position. Jabat has remained geomorphically stable for the past 2,000 years. Findings suggest reef platforms may accommodate the oldest reef islands in atoll systems, which may have profound implications for questions of prehistoric migration through Pacific archipelagos.

Description: [1]  Variable supply of iron to the ocean is often invoked to explain part of past changes in atmospheric CO 2 (CO 2atm ). Using model simulations we find that CO 2atm is sensitive on the order 15, 2 and 1 ppm to sedimentary, dust and hydrothermal iron input. CO 2atm is insensitive to dust because it is not the major iron input to the Southern Ocean. Modifications to the relative export of Si(OH) 4 to low latitudes are opposite to those predicted previously. Although hydrothermalism is the major control on the iron inventory in ~25% of the ocean, it remains restricted to the deep ocean, with minor effects on CO 2atm . Nevertheless, uncertainties regarding the iron binding ligand pool can have significant impacts on CO 2atm . Ongoing expansion of iron observations as part of GEOTRACES will be invaluable in refining these results.

Description: [1]  Hydrothermal venting often occurs at submarine volcanic calderas on island arc chains, typically at shallower depths than mid-ocean ridges. The effect of these systems on ocean biogeochemistry has been under-investigated to date. Here we show that hydrothermal effluent from an island arc caldera was rich in Fe(III) colloids (0.02 – 0.2 µm; 46% of total Fe), contributing to a fraction of hydrothermal Fe that was stable in ocean water. Iron(III) colloids from island arc calderas may be transferred into surrounding waters (generally 0-1500 m depth) by ocean currents, thereby potentially stimulating surface ocean primary productivity. Hydrothermal Fe oxy-hydroxide particles (〉0.2 µm) were also pervasive in the studied caldera and contained high concentrations of oxyanions of phosphorus (P), vanadium (V), arsenic (As), and manganese (Mn). Hydrothermal island arcs may be responsible for 〉 50% of global hydrothermal P scavenging and 〉 40% V scavenging, despite representing 〈10% of global hydrothermal fluid flow.

Description: [1]  Sea surface topography observations are deduced from an airborne reflectometry experiment. The GORS (GNSS Occultation Reflectometry Scatterometry) receiver was set up aboard the German HALO (High Altitude LOng range) research aircraft. Flights were conducted over the Mediterranean Sea about 3500 m above sea level. A signal path model divided into large and small scale contributions is used for phase altimetry. The results depict geoid undulations and resolve anomalies of the sea surface topography.For the whole experiment 65 tracks over the Mediterranean sea are retrieved and compared with a topography model. Tracks distinguish between RHCP/LHCP. The differnce, however, is not significant for this study. Precision and spatial resolution decrease disproposionately at low elevations. Eight tracks with centimeter precision are obtained between 11 ∘ and 33 ∘ of elevation. At higher elevation angles the number of tracks is significantly reduced due to surface roughness. In future such retrievals could contribute to ocean eddy detection.

Description: [1]  The electrical conductivity of single-crystal brucite was measured as a function of pressure and temperature (P-T) using impedance spectroscopy (IS). IS measurements demonstrated that electrical conductivity perpendicular to the c -axis is nearly half order of magnitude higher than that parallel to the c -axis under the same P-T conditions. Electrical conductivity increased by two orders of magnitude during compression from 3.7 to 11 GPa, irrespective of crystallographic direction. However, the conductivity increase with pressure became less significant upon further compression from 11 to 13 GPa. The pressure effect is closely related to the interactions between neighboring hydroxyls. The ratio of free protons to total hydrogen increases from 2% to 33% when pressure increases from 3.7 to 13 GPa at 950 K. This indicates that the most of protons are bound within the crystal structure at low pressures, whereas more protons become free and mobile at higher pressures.

Description: [1]  Global stacks of receiver functions clearly exhibit the upper mantle stratification. Besides the most prominent seismic discontinuities, such as the Moho and the 410 and 660 km discontinuities, a negative discontinuity is detected at a depth of ~600 km, indicating a low velocity layer at the base of the mantle transition zone. The slant slack technique helps to identify the primary conversions from the multiple reverberations. Presence of the negative 600 km discontinuity underneath both continent and ocean island stations, where the crustal thickness significantly differs, also precludes the possible cause of crustal reverberations. We conclude that the negative 600 km discontinuity could be a global feature, possibly resulted from accumulation of ancient subducted oceanic crust. The X-discontinuity at ~300 km depth is also observed in our global stacks, which can be explained by the coesite-stishovite phase transformation.

Description: [1]  Continental passive margins are characterized by a wide variety of geometries and widths. Whether these variations have an influence on subsequent dynamics of orogenesis is unresolved. To investigate, a series of upper mantle numerical experiments were performed with systematically varied continental margin widths and geometries. Results show that the vertical geometry of subducting continental margin crust controls both crustal and mantle lithosphere deformation. On both scales, deformational end-members can be identified. Namely break-off vs. delamination of continental mantle lithosphere, and double vergence vs. single vergence of crustal thrust fronts form as a direct result of passive margin geometry. We find that the subduction of upper crust to depths 〉 100 km promotes lithospheric delamination and is facilitated by an extended passive margin. Modelled orogens show decreasing upper plate deformation with increasing margin width. These results suggest that along-strike deformational variation within orogens may develop due to pre-collision passive margin geometry.

Description: [1]  Mackenzie River discharge and bathymetry effects on sea ice in the Beaufort Sea are examined in 2012 when Arctic sea ice extent hit a record low. Satellite-derived sea surface temperature revealed warmer waters closer to river mouths. By 5 July 2012, Mackenzie warm waters occupied most of an open-water area about 316,000 km 2 . Surface temperature in a common open-water area increased by 6.5 °C between 14 June and 5 July 2012, before and after the river waters broke through a recurrent landfast ice barrier formed over the shallow seafloor offshore the Mackenzie Delta. In 2012, melting by warm river waters was especially effective when the strong Beaufort Gyre fragmented sea ice into unconsolidated floes. The Mackenzie and other large rivers can transport an enormous amount of heat across immense continental watersheds into the Arctic Ocean, constituting a stark contrast to the Antarctic that has no such rivers to affect sea ice.

Description: [1]  Observed blocking trends are diagnosed to test the hypothesis that recent Arctic warming and sea ice loss has increased the likelihood of blocking over the Northern Hemisphere. To ensure robust results, we diagnose blocking using three unique blocking identification methods from the literature, each applied to four different reanalyses. No clear hemispheric increase in blocking is found for any blocking index, and while seasonal increases and decreases are found for specific isolated regions and time periods, there is no instance where all three methods agree on a significant trend. Blocking is shown to exhibit large interannual and decadal variability, highlighting the difficulty in separating any potentially forced response from natural variability.