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: 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: Herein, we report on the crystal structures of Nb 2 AlC and TiNbAlC—actual composition (Ti 0.45 ,Nb 0.55 ) 2 AlC—compounds determined from Rietveld analysis of neutron diffraction patterns in the 300–1173 K temperature range. The average linear thermal expansion coefficients of a Nb 2 AlC sample in the a and c directions are, respectively, 7.9(5) × 10 −6 and 7.7(5) × 10 −6  K −1 on one neutron diffractometer and 7.3(3) × 10 −6 and 7.0(2) × 10 −6  K −1 on a second diffractometer. The respective values for the (Ti 0.45 ,Nb 0.55 ) 2 AlC composition—only tested on one diffractometer—are 8.5(3) × 10 −6 and 7.5(5) × 10 −6  K −1 . These values are relatively low compared to other MAX phases. Like other MAX phases, however, the atomic displacement parameters (APDs) show that the Al atoms vibrate with higher amplitudes than the Ti and C atoms, and more along the basal planes than normal to them. When the predictions of the APDs obtained from density functional theory are compared to the experimental results, good quantitative agreement is found for the Al atoms. In case of the Nb and C atoms, the agreement was more qualitative.

Description: In this study, lithium disilicate (LS 2 ) glass samples with different particle sizes ranging from less than 105 to 850 μm were prepared. These specimens were inserted in a Pt-Rh DSC crucible and heated to 850°C at different rates (ϕ = 0.5–30 K/min) to identify their crystallization peaks. The activation energies for the overall crystallization ( E ) and the Avrami coefficient ( n ) were evaluated using different nonisothermal models. Specifically, n was evaluated using the Augis–Benett model and the Ozawa method, and E was evaluated using the Kissinger and Ligero methods. As expected, the coarse particles mainly crystallized in the volume, while surface crystallization was predominant in the samples with particle sizes of less than 350 μm. This result was confirmed through SEM analysis of the double stage heat-treated samples. In contrast with previous studies, our results demonstrated that the activation energy decreased as the particle size increased. In addition, no clear correlation between the peak intensity (δ T p ) and the particle size was observed.

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: 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: ABSTRACT With the objective to understand the generation, propagation and nonlinear evolution of ion cyclotron waves (ICWs) in the corona and solar wind, we use electromagnetic hybrid (kinetic ions, fluid electrons) simulations with a non-uniform magnetic field. ICWs are generated by the temperature anisotropy of O 5+ ions as minority species in a proton-electron plasma with uniform density. A number of magnetic field models are used including radial and spiral with field strength decreasing linearly or with the square of the radial distance. O 5+ ions with perpendicular temperature larger than parallel are initially placed in the high magnetic field regions. These ions are found to expand outward along the magnetic field. Associated with this expansion, ion cyclotron waves propagating along the magnetic field are also seen to expand outward. These waves are generated at frequencies below the local gyro-frequency of O 5+ ions propagating parallel and anti-parallel to the magnetic field. Through analysis of the simulation results we demonstrate that wave generation and absorption takes place at all radial distances. Comparing the simulation results to observations of ICWs in the solar wind shows some of the observed wave characteristics may be explained by the mechanism discussed in this paper.

Description: The characteristics of nighttime medium-scale travelling ionospheric disturbance (MSTID) features observed over Yonaguni (24.5 o N, 123.0 o E; 19.3 o N dip latitude), Japan are studied using all-sky imaging of OI 630.0 nm airglow emission. The uniqueness of these observations is that the area observed by the imager covers the transition region between low to middle latitudes in the ionosphere. Typical low latitude limit of mid-latitude type nighttime MSTIDs possessing phase front alignments along the northwest to the southeast occurs in this region. These MSTID features are rarely sighted at dip latitudes below 15 o . We selected two year period for analysis in which one year corresponded to the solar minimum conditions and another year to the solar maximum conditions. The MSTIDs were observed to extend to farther lower latitudes during the solar minimum conditions than during the solar maximum periods. Their observed range of wavelengths, phase velocities, phase front alignment and propagation directions are similar to those observed at typical mid-latitude sites. However, on many occasions the phase fronts of the observed MSTIDs did not extend over the whole field of view of the imager indicating that some process inhibits their extension to further lower latitudes. Detailed investigation suggests that the poleward propagating enhancement of airglow intensity, probably associated with the midnight pressure bulge, causes the MSTID features to disappear when they reach lower latitudes later in the night. When the MSTIDs reach lower latitudes well before midnight, they are found to be inhibited by the equatorial ionization anomaly crest region.

Description: We study a magnetosphere-ionosphere coupling at low-latitudes during a moderate (CIR/HSS-driven) geomagnetic storm on 22 July 2009. Recently, it has been shown that during major (CME-driven) storms, quasi-trapped 〉30 keV electrons largely enhance below the radiation belt in the forbidden zone and produce an additional ionization in the topside ionosphere. In this work, we examine a case of the recurrent storm when the magnetosphere-ionosphere coupling through the quasi-trapped electrons also may take place. Data from NOAA/POES and Japanese GOSAT satellites were used to identify the forbidden electron enhancement (FEE). We find a positive vertical gradient of the electron fluxes that indicates to the radiation belt as a source of FEE. Using global ionospheric maps (GIM), radiotomography reconstructions from beacon data and COSMIC/FS3 radio occultation measurements, we have observed an unusually large area in the night-time ionosphere with increased total electron content (TEC) and prominent elevation of the F-layer at low-latitudes that coincides with FEEs spatially and temporarily. Ionizing particles are considered as an addition source of ionization along with generally accepted mechanisms for storm time TEC increase (a positive ionospheric storm). We discuss relative contributions of the FEE and disturbance dynamo electric field in the TEC increases during the storm recovery phase.

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: Bursty bulk flow (BBF) events, frequently observed in the magnetotail, carry significant energy and mass from the tail region at distances that are often greater than 20 R E into the near-Earth plasma sheet at ~10 R E where the flow is slowed and/or diverted. This region at ~10 R E is referred to as the BBF braking region. A number of possible channels are available for the transfer or dissipation of energy in BBF events including adiabatic heating of particles, the propagation of Alfvén waves out of the BBF braking region and into the auroral region, diverted flow out of the braking region, and energy dissipation within the braking region itself. This study investigates the generation of intense high-frequency electric field activity observed within the braking region. When present, these intense electric fields have power above the ion cyclotron frequency and almost always contain nonlinear structures such as electron phase space holes and double layers, which are often associated with field-aligned currents. A hypothesis in which the observed high-frequency electric field activity is generated by field-aligned currents resulting from turbulence in the BBF braking region is considered. Although linear Alfvén waves can generate field-aligned currents, based on theoretical calculations, the required currents are likely not the result of linear waves. Observations from the THEMIS satellites support the picture of a turbulent plasma leading to the generation of nonlinear kinetic structures. This work provides a possible mechanism for energy dissipation in turbulent plasmas.

Description: Dipolarization fronts (DFs) are often associated with the leading edge of Earthward bursty bulk flows in the magnetotail plasma sheet. Here multi-spacecraft THEMIS observations are used to show that a spatially limited region of counter-propagating ion beams, whose existence is not evident in either the plasma moments or the electric field, is observed on the low density side of DFs. The THEMIS magnetic field data are used to establish appropriate comparison cuts through a particle-in-cell (PIC) simulation of reconnection, and very good agreement is found between the observed and simulated ion distributions on both sides of the DF. Self-consistent back-tracing shows that the ion beams originate from the thermal component of the pre-existing high density plasma into which the DF is propagating; they do not originate from the inflow region in the traditional sense. Forward tracing shows that some of these ions can subsequently overtake the DF and pass back into the high density pre-existing plasma sheet with an order-of-magnitude increase in energy; this process is distinct from other ion reflection processes that occur directly at the DF. The interaction of the reconnection jet with the pre-existing plasma sheet therefore occurs over a macroscopic region, rather than simply being limited to the thin DF interface. A more general consequence of this study is the conclusion that reconnection jets are not simply fed by plasma inflow across the separatrices, but are also fed by plasma from the region into which the jet is propagating; the implications of this finding are discussed.

Description: Low energy (1-10 MeV) neutrons emanating from the Sun provide unique information about accelerated ions with steep energy spectra that may be produced in weak solar flares. However, observation of these solar neutrons can only be made in the inner heliosphere where measurement is difficult due to high background rates from neutrons produced by energetic ions interacting in the spacecraft. These ions can be from solar energetic particle events or produced in passing shocks associated with fast coronal mass ejections. Therefore, it is of the utmost importance that investigators rule out these secondary neutrons before making claims about detecting neutrons from the Sun. The MErcury Surface, Space ENvironment, GEochemistry, and Ranging ( tect MESSENGER ) Neutron Spectrometer recorded an hour-long neutron transient beginning at 15:45 UTC on 2011 June 4 for which [13] claim there is “strong evidence" that the neutrons were produced by the interaction of ions in the solar atmosphere. We studied this event in detail using data from the MESSENGER neutron spectrometer, gamma-ray spectrometer, X-ray Spectrometer, and Energetic Particle Spectrometer, and from the particle spectrometers on STEREO A . We demonstrate that the transient neutrons were secondaries produced by energetic ions, probably accelerated by a passing shock, that interacted in the spacecraft. We also identify significant faults with the authors’ arguments in favor of a solar neutron origin for the transient.

Description: The response of the D-region low latitude ionosphere has been examined for extreme space weather event of 14-16 December 2006 associated with a X1.5 solar flare and an intense geomagnetic storm ( Dst  = -146 nT) using VLF signals from NWC (19.8 kHz) and NPM (21.4 kHz) transmitters monitored at Suva (Geog. 18.10 ο S, 178.40 ο E), Fiji. Modeling of flare associated amplitude and phase enhancements of NWC (3.6 dB, 223 o ) and NPM (5 dB, 153 o ) using Long Wave Propagation Capability code shows reduction in the D-region reflection height ( H ') by 11.1 km and 9.4 km, and enhancement in ionization gradients described by increases in the exponential sharpness factor ( β ) by 0.122 and 0.126 km -1 , for the NWC and NPM paths, respectively. During the storm the daytime signal strengths of the NWC and NPM signals were reduced by 3.2 dB on 15 and 16 December (for about 46 hrs) and recovered by 17 December. Modelling for the NWC path shows that storm-time values of H ' and β were reduced by 1.2 km and 0.06 km -1 , respectively. Morlet wavelet analysis of signals amplitudes shows no clearly strong signatures of gravity wave propagation to low latitudes during the main and recovery phases. The reduction in VLF signal strength is due to increased signal attenuation and absorption by the Earth-ionosphere waveguide due to storm-induced D-region ionization changes and hence changes in D-region parameters. The long duration of the storm effect results from the slow diffusion of changed composition/ionization at D-region altitudes compared with higher altitudes in the ionosphere.

Description: We study the ion density and temperature in the pre- and post-dipolarization plasma sheets in the Earth's magnetotail, using 9 years (2001–2009) of Cluster data. For our study we selected cases when Cluster observed dipolarization fronts (DFs) with an earthward plasma flow greater than 150 km/s. We perform a statistical study of the temperature and density variations during the DF crossings. Earlier studies concluded that on average the temperature increases while the densitydecreases across the DF. Our statistical results show a more diverse picture: While ~53% of the DFs follow this pattern (category A), for ~28% the temperature decreases while the density increases across the DF (category B). We found an overall decrease in thermal pressure for category A DFs with a more pronounced decrease at the DF duskside, while DFs of the category B showed no clear pattern in the pressure change. Both categories are associated with earthward plasma flows, but with some difference: (1) Category A flows are faster than category B flows. (2) The observations indicate that category B flows are directed perpendicular to the current in the near-Earth current sheet while category A flows are tilted slightly duskward from this direction. (3) The background B z of category B is higher than that of category A. Based on these results we hypothesize that after reconnection takes place, a BBF emerges with category A characteristics, and as it travels earthward it further evolves into category B characteristics, which is in a more dipolarized region with slower plasma flow (closer to the flow braking region).

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: In this letter, 2-μm Pb 0.97 La 0.02 (Zr 0.75 Sn 0.18 Ti 0.07 )O 3 antiferroelectric thick film with tetragonal structure was prepared. The effects of operating electric field, temperature, and frequency on the thermal–electrical energy harvesting capacity of the film were studied by using the Olsen cycle. The results demonstrated that giant energy harvesting effect could be realized in the antiferroelectric thick film. The maximum harvestable energy density per cycle of the film was about 7.8 J/cm 3 at 1 kHz, which was the largest reported value to date. The corresponding energy harvesting efficiency was 0.53%. Moreover, the film had a low leakage current density (about 7.3 × 10 −7 and 3.9 × 10 −5  A/cm 2 at 25 and 200°C, respectively), which was favorable for its application in the devices of the thermal–electrical energy harvesting.

Description: Discharged energy properties of PbO–SrO–Na 2 O–Nb 2 O 5 –SiO 2 glass-ceramics with crystallization time from 1 to 1000 min were investigated by measuring their hysteresis loops (described as quasi-static measuring method) and pulse-discharge current-time curves (described as dynamic measuring method). The results show the same trend for both measuring methods: With the increment of crystallization time, the discharged energy density increases gradually, while the energy efficiency decreases. The highest energy efficiencies were obtained in the sample with crystallization time of 1 min, which are 96.3% and 82.4%, corresponding quasi-static and dynamic measurement, respectively. The reduction of energy efficiency with crystallization time is attributed to combined effect of ferroelectric polarization and interfacial polarization, and part of the corresponding energy could not release in the pulse-discharge process.

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: Ice cores are archives of climate change and possibly large solar proton events (SPEs). Wolff et al . [2012] used a single event, a nitrate peak in the GISP2-H core, which McCracken et al . [2001a] time associated with the poorly quantified 1859 Carrington event, to discredit SPE-produced, impulsive nitrate deposition in polar ice. This is not the ideal test case. We critique the Wolff et al. analysis and demonstrate that the data they used cannot detect impulsive nitrate events because of resolution limitations. We suggest re-examination of the top of the Greenland ice sheet at key intervals over the last two millennia with attention to fine resolution and replicate sampling of multiple species. This will allow further insight into polar depositional processes on a sub-seasonal scale, including atmospheric sources, transport mechanisms to the ice sheet, post-depositional interactions, and a potential SPE association.

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: New lead-free perovskite solid solution ceramics of (1  − x )( Bi 1/2 Na 1/2 ) TiO 3 – x Ba ( Ni 1/2 Nb 1/2 ) O 3 [(1− x )BNT– x BNN, x  =   0.02–0.06) were prepared and their dielectric, ferroelectric, piezoelectric, and electromechanical properties were investigated as a function of the BNN content. The X-ray diffraction results indicated that the addition of BNN has induced a morphotropic phase transformation from rhombohedral to pseudocubic symmetry approximately at x  =   0.045, accompanying an evolution of dielectric relaxor behavior as characterized by enhanced dielectric diffuseness and frequency dispersion. In the proximity of the ferroelectric rhombohedral and pseudocubic phase coexistence zone, the x  =   0.045 ceramics exhibited optimal piezoelectric and electromechanical coupling properties of d 33 ~121 pC/N and k p ~0.27 owing to decreased energy barriers for polarization switching. However, further addition of BNN could cause a decrease in freezing temperatures of polar nanoregions till the coexistence of nonergodic and ergodic relaxor phases occurred near room temperature, especially for the x  =   0.05 sample which has negligible negative strains and thus show the maximum electrostrain of 0.3% under an external electric field of 7 kV/mm, but almost vanished piezoelectric properties. This was attributed to the fact that the induced long-range ferroelectric order could reversibly switch back to its original ergodic state upon removal of external electric fields.

Description: Light transmission in polycrystalline magnesium fluoride was studied as a function of the mean grain size at different wavelengths. The mean grain size was varied by annealing hot-pressed billets in argon atmosphere at temperatures ranging from 600°C to 800°C for 1 h. The grain-size and grain-orientation distributions were characterized by electron back scatter diffraction. The scattering coefficients were calculated from the in-line transmittance measured at various wavelengths. The scattering coefficient of polycrystalline magnesium fluoride increased linearly with the mean grain size and inversely with the square of the wavelength of light. It is shown that these trends are consistent with theoretical models based on both a limiting form of the Raleigh–Gans–Debye (RGD) theory of particle scattering and light retardation theories that take refractive index variations along the light path. Quantitative predictions of the theories are, however, subject to uncertainly due to the restrictive assumptions made in the theories and difficulties in representing the microstructure in the theoretical models. In particular, grain-size distribution has a significant influence on the scattering coefficient calculated using particle scattering models.

Description: Determining the absolute chronology of ceramic artifacts has significant implications for archeological and historical research. Wilson, Hall et al . recently suggested a new technique for direct absolute dating of archeological ceramics based on a moisture-induced chemical reaction, called rehydroxylation (RHX) dating. RHX dating proceeds by measuring the mass of chemically combined water in the ceramics in the form of OH hydroxyls, and the mass gain rate at the Effective Lifetime Temperature (ELT) that the ceramics experienced over its lifetime. To date, ELT determinations have been based on estimates of the ceramic's lifetime temperature history; taking into account weather and climate data and the depth at which the artifact was found. The uncertainty in determining the ELT can be a major component of the overall dating uncertainty. Here, we propose an alternative method which relies minimally on weather and climate data, and provides more precise determinations of the ELT and the ceramic age . The proposed method (SAS: Same Age Samples) involves a minimum of four measurements of the RHX mass gain rate constant for two ceramic samples of the same age at two temperatures. We show via simulations that the proposed SAS method can determine the ELT with a precision of 0.2 K which is comparable to the best ELT determination based on lifetime temperature history, and also comparable to available microbalance temperature resolutions of around 0.1 K. The corresponding percent age error is then 1.4%, or 43 yr for a 3000-yr-old ceramic. The proposed SAS method should be tested with ceramic samples of different ages, whose ELT are well-known.

Description: Glasses in the Na 2 O–CaO–SrO–ZnO–SiO 2 system have previously been investigated for suitability as a reagent in Al-free glass polyalkenoate cements (GPCs). These materials have many properties that offer potential in orthopedics. However, their applicability has been limited, to date, because of their poor strength. This study was undertaken with the aim of increasing the mechanical properties of a series of these Zn-based GPC glasses by doping with nitrogen to give overall compositions of: 10Na 2 O–10CaO–20SrO–20ZnO–(40−3 x )SiO 2 – x Si 3 N 4 ( x is the no. of moles of Si 3 N 4 ). The density, glass-transition temperature, hardness, and elastic modulus of each glass were found to increase fairly linearly with nitrogen content. Indentation fracture resistance also increases with nitrogen content according to a power law relationship. These increases are consistent with the incorporation of N into the glass structure in threefold coordination with silicon resulting in extra cross-linking of the glass network. This was confirmed using 29 Si MAS-NMR which showed that an increasing number of Q 2 units and some Q 3 units with extra bridging anions are formed as nitrogen content increases at the expense of Q 1 units. A small proportion of Zn ions are found to be in tetrahedral coordination in the base oxide glass and the proportion of these increases with the presence of nitrogen.

Description: Severe geomagnetic storms have a strong impact on space communication and satellite navigation systems. Forecasting the appearance of geomagnetically induced disturbances in the ionosphere is one of the urgent goals of the space weather community. The challenge is that the processes governing the distribution of the crucial ionospheric parameters has a rather poor quantitative description and the models, built using the empirical parameterisations, have limited capabilities for operational purposes. On the other hand, data assimilation techniques are becoming more and more popular for nowcasting the state of the large-scale geophysical systems. We present an example of an ionospheric data assimilation system performance assessment during a strong geomagnetic event, which took place on 26 September 2011. The first-principle model has assimilated slant total electron content measurements from a dense network of ground stations, provided by the Norwegian Mapping Authority. The results have shown satisfactory agreement with independent data and demonstrate that the assimilation model is accurate to about 2–4 TEC units, and can be used for operational purposes in high-latitude regions. The operational system performance assessment is the subject of future work.

Description: [1]  The last solar minimum period was anomalously extended and low in EUV irradiance compared with previous solar minima. It can readily be expected that the thermosphere and ionosphere must be correspondingly affected by this low solar activity. While there have been unanimous reports on the thermospheric changes, being cooler and lower in its density as expected, the ionospheric responses to low solar activity in previous studies were not consistent with each other, probably due to the limited ionospheric observations used for them. In this study, we utilized the measurements of total electron content (TEC) from TOPEX and JASON-1 satellites during the periods of 1992 to 2010, which includes both the last two solar minimum periods, in order to investigate how the ionosphere responded to the extremely low solar activity during the last solar minimum compared with previous solar minimum. Although the global daily mean TECs show negligible differences between the two solar minimum periods, the global TEC maps reveal that there are significant systematic differences ranging from about -30% to +50% depending on local time, latitude and season. The systematic variations of the ionospheric responses seem to mainly result from the relative effects of reduced solar EUV production and reduced recombination rate due to thermospheric changes during the last solar minimum period.

Description: [1]  Some of the potentially most destructive effects of severe space weather storms are caused by the Geomagnetically Induced Currents. GICs can cause failures of electric transformers and result in wide-spread blackouts. GICs are induced by the time variability of the magnetic field, and are closely related to the time derivative of the local magnetic field perturbation. Predicting dB/dt is rather challenging, since the local magnetic perturbations and their time derivatives are both highly fluctuating quantities, especially during geomagnetic storms. The currently available first-principles based and empirical models cannot predict the detailed minute-scale or even faster time variation of the local magnetic field. On the other hand, Pulkkinen et al. [2013] demonstrated recently that several models can predict with positive skill scores whether the horizontal component of dB/dt at a given magnetometer station will exceed some threshold value in a 20-minute time interval. In this paper we investigate if one can improve the efficiency of the prediction further. We find that the Space Weather Modeling Framework, the best performing among the five models compared by Pulkkinen et al. [2013], shows significantly better skill scores in predicting the magnetic perturbation than predicting its time derivative, especially for large deviations. We also find that there is a strong correlation between the magnitude of dB/dt and the magnitude of the horizontal magnetic perturbation itself. Combining these two results one can devise an algorithm that gives better skill scores for predicting dB/dt exceeding various thresholds in 20-minute time intervals than the direct approach.

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: [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: [1]  As a weakly magnetized planet, Mars ionosphere/atmosphere interacts directly with the shocked solar wind plasma flow. Even though many numerical studies have been successful in reproducing numerous features of the interaction process, these earlier studies focused mainly on interaction under steady solar wind conditions. Recent observations suggest that plasma escape fluxes are significantly enhanced in response to solar wind dynamic pressure pulses. In this study, we focus on the response of the ionosphere to pressure enhancements in the solar wind. Through modeling of two idealized events using a magnetohydrodynamics model, we find that the upper ionosphere of Mars responds almost instantaneously to solar wind pressure enhancements, while the collision dominated lower ionosphere (below ~150 km) does not have noticeable changes in density. We also find that ionospheric perturbations in density, magnetic field and velocity can last more than an hour after the solar wind returns to the quiet conditions. The topside ionosphere forms complicated transient shapes in response, which may explain unexpected ionospheric behaviors in recent observations. We also find that ionospheric escape fluxes do not correlate directly with simultaneous solar wind dynamic pressure. Rather, their intensities also depend on the earlier solar wind conditions. It takes a few hours for the ionospheric/atmospheric system to reach a new quasi-equilibrium state.

Description: [1]  Combining THEMIS wave and particle observations and a quantitative calculation of linear wave growth rate, we demonstrate that magnetosonic (MS) waves can be locally excited by ion ring distributions in the Earth's magnetosphere when the ion ring energy is comparable to the local Alfven energy. MS waves in association with ion ring distributions were observed by THEMIS A on 24 November 2010 in the afternoon sector, both outside the plasmapause where the wave spectrum varied with f LHR and inside the plasmapause where the wave frequency band remained nearly constant. Our plasma instability analysis in three different regions shows that higher and narrow frequency band MS waves are excited locally outside the plasmapause, and lower and broad frequency band MS waves are excited in the region where the density slightly increases. However, there is no evidence for wave excitation inside the plasmapause, and wave propagation from a distant source is needed to explain their existence. The simulation of the MS wave growth rate spectra during this event agrees reasonably well with the observed wave magnetic field power spectra. We also simulated a MS wave event on 19 October 2011 in the dusk sector, and found that the ion ring distribution with an ion ring energy slightly higher than the local Alfven energy can excite the typical broad band MS waves outside the plasmapause.

Description: [1]  The second Radio Aurora Explorer (RAX-2) satellite has completed more than 30 conjunction experiments with the AMISR chain of incoherent scatter radars in Alaska, and Resolute Bay, Canada. Coherent radar echoing occurred during four of the passes: three when E region electron drifts exceeded the ion acoustic speed threshold and one during HF heating of the ionosphere by the HAARP heater. In this paper, we present the results for the first three passes associated with backscatter from natural irregularities. We analyze, in detail, the largest drift case because the plasma turbulence was the most intense and because the corresponding ground-to-space bi-static scattering geometry was the most favorable for magnetic aspect sensitivity analysis. A set of data analysis procedures including interference removal, autocorrelation analysis, and the application of a radar beam deconvolution algorithm mapped the distribution of E region backscatter with 3 km resolution in altitude and ∼ 0.1 ∘ in magnetic aspect angle. To our knowledge, these are the highest resolution altitude-resolved magnetic aspect sensitivity measurements made at UHF frequencies in the auroral region. In this paper, we show that, despite the large electron drift speed of ∼ 1500 m/s, the magnetic aspect sensitivity of sub-meter scale irregularities is much higher than previously reported. The root-mean-square of the aspect angle distribution varied monotonically between 0.5-0.1 ∘ for the altitude range 100-110 km. Findings from this single but compelling event suggest that sub-meter scale waves propagating at larger angles from the main E  ×  B flow direction (secondary waves) have parallel electric fields that are too small to contribute to E region electron heating. It is possible anomalous electron heating in the auroral electrojet can be explained by (a) the dynamics of those sub-meter scale waves propagating in the E  ×  B direction (primary waves) or (b) the dynamics of longer wavelengths.

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: [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: [1]  Working toward a physical understanding of how solar-wind/magnetosphere coupling works, four arguments are presented indicating that the solar-wind electric field v sw  × B sw does not control the rate of reconnection between the solar wind and the magnetosphere. Those four arguments are (1) that the derived rate of dayside reconnection is not equal to solar-wind electric field, (2) that electric-field driver functions can be improved by a simple modification that disallows their interpretation as the solar-wind electric field, (3) that the electric field in the magnetosheath is not equal to the electric field in the solar wind, and (4) that the magnetosphere can mass load and reduce the dayside reconnection rate without regard for the solar-wind electric field. The data is more consistent with a coupling function based on local control of the reconnection rate than the Axford conjecture that reconnection is controlled by boundary conditions irrespective of local parameters. Physical arguments that the solar-wind electric field controls dayside reconnection are absent; it is speculated that it is a coincidence that the electric field does so well at correlations with geomagnetic indices.

Description: [1]  This paper presents a case study from a single, six-hour observing period to illustrate the application of techniques developed for interferometric radio telescopes to the spectral analysis of observations of ionospheric fluctuations with sparse arrays. We have adapted the deconvolution methods used for making high dynamic range images of cosmic sources with radio arrays to making comparably high dynamic range maps of spectral power of wavelike ionospheric phenomena. In the example presented here, we have used observations of the total electron content (TEC) gradient derived from Very Large Array (VLA) observations of synchrotron emission from two galaxy clusters at 330 MHz as well as GPS-based TEC measurements from a sparse array of 33 receivers located within New Mexico near the VLA. We show that these techniques provide a significant improvement in signal to noise (S/N) of detected wavelike structures by correcting for both measurement inaccuracies and wavefront distortions. This is especially true for the GPS data when combining all available satellite/receiver pairs, which probe a larger physical area and likely have a wider variety of measurement errors than in the single-satellite case. In this instance, we found the peak S/N of the detected waves was improved by more than an order of magnitude. The data products generated by the deconvolution procedure also allow for a reconstruction of the fluctuations as a two-dimensional waveform/phase screen that can be used to correct for their effects.

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]  In this paper, we test whether time periods with hot proton temperature anisotropy are associated with EMIC waves, and whether the plasma conditions during the observed waves satisfy the linear theory threshold condition. We identify 865 events observed by the Composition DIstribution Function (CODIF) instrument onboard Cluster spacecraft 4 (SC4) during 1 January 2001 – 1 January 2011 that exhibit a positive temperature anisotropy ( A hp  =  T ⊥  h / T ∥  h  − 1) in the 10-40 keV protons. The events occur over an L range from 4 to 10 in all magnetic local times and at magnetic latitudes (MLAT) within ±50°. Of these Hot Proton Temperature Anisotropy (HPTA) events, only 68 events have electromagnetic ion cyclotron (EMIC) waves. In these 68 HPTA events, for those at 3.8〈 L  ≤ 5 and |MLAT| ≤ 10 ° , the EMIC waves with powers 〉1.0 nT 2 /Hz mainly appear in the region with f EMIC / f H , eq  〈 0.8. Two stop bands are present, one near the region with f EMIC / f H , eq  ≈ 0.33, the other in the region with 0.8 〈  f EMIC / f H , eq  〈 0.9. Most of the EMIC waves in the He, H, and 〉 H bands satisfy A hp /( A hp  + 1) 〉  f EMIC / f H , lo , A hp /( A hp  + 1) 〉 0.45 *  f EMIC / f H , lo , and A hp /( A hp  + 1) 〈 0.45 *  f EMIC / f H , lo .   f EMIC , f H , eq   and f H , lo are the EMIC wave frequency, the magnetic equatorial and the local proton gyrofrequencies. We also find that the EMIC waves predominantly occur with A hp  〉 0.25. By testing a threshold equation for the EMIC instability based on linear theory, we find that for EMIC waves with |MLAT| ≤ 10 ° in the He, H and 〉 H bands the percentages that satisfy the predicted conditions for wave growth by the threshold equation are 15.2%, 24.6% and 25.6%. For the EMIC waves with |MLAT| 〉 10 ° the percentages that satisfy the wave growth predicted conditions are only 2.8%, 2.6% and 0.0%. Finally, possible reasons for the low forecast accuracies of EMIC waves are suggested.

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: [1]  The simultaneous onset of the preliminary impulse (PI) of the geomagnetic sudden commencement at high latitude and dayside dip equator is explained by means of the TM 0 mode waves propagating at the speed of light in the Earth-ionosphere waveguide (EIW) [Kikuchi et al., 1978]. A couple of issues remain to be addressed in the EIW model: (1) How is the TM 0 mode wave is excited by the field-aligned currents in the polar region? (2) How are the quasi-steady ionospheric currents are achieved by the TM 0 mode waves? (3) How simultaneous or delayed are the onset and peak of the equatorial PI with respect to the high latitude PI? To address these issues, we examine the TEM (TM 0 ) mode wave propagation in the finite-length transmission lines replacing the pair of FACs (magnetosphere-ionosphere (MI) transmission line) and the Earth-ionosphere waveguide (ionosphere-ground (IG) transmission line). The issue (1) is addressed by showing that a fraction of the TEM mode wave is transmitted from the MI to IG transmission lines through the polar ionosphere. To address the issues (2) and (3), we examine the properties of the finite-length IG transmission line with finite ionospheric conductivity. It is shown that the ionospheric currents start to grow instantaneously and continue to grow gradually with time constants of 1-10 sec depending on the ionospheric conductivity. The MIG transmission line enables us to explain the instantaneous onset and delayed peak time of the equatorial PI and quick electric field response of the low latitude ionosphere and inner magnetosphere.

Description: This study presents a thermodynamic analysis to predict the type of initial, amorphous oxide overgrowth (i.e., am - Al 2 O 3 or am - SiO 2 ) on bare Al – Si alloy substrates. This analysis have taken into account the energies associated with both its interfaces (interface between the Al – Si alloy substrate and the thin oxide film and interface between the thin oxide film and vacuum) along with the bulk Gibbs free energy of oxide formation. This developed analysis is then applied for various parameters, such as, Si alloying element content at the substrate/oxide interface, the growth temperature, the oxide film thickness (up to 1 nm), and various low-index crystallographic surfaces of the substrate. It is found that am - SiO 2 overgrowth is thermodynamically preferred for a combination of lower oxide film thickness, lower growth temperature, and lower Si alloying content at the alloy/oxide interface. This is because of the overcompensation of the lower energies of both the interfaces over the bulk Gibbs free energy. Furthermore, it is found that for all cases, am - Al 2 O 3 forms a more stable interface with Al – Si alloy than am - SiO 2.

Description: [1]  Plasmoids and other reconnection-related signatures have been observed in Jupiter's magnetotail through analysis of magnetic field and energetic particle data. Previous studies have established the spatial distribution and recurrence period of tail reconnection events, and identified the location of a statistical x-line separating inward and outward flow. Here we present new analysis focusing specifically on 43 plasmoid signatures observed in magnetometer data in order to establish the average properties and internal structure of Jovian plasmoids. We present statistics on the observed plasmoid length scale, duration, radial position, and local time distribution. On average, the observed plasmoids have a ~3 R J radial extent and ~7 minute duration, and result in the closure of ~4-8 GWb of open flux from reconnection of open field lines in the post-plasmoid plasma sheet. We also determine the amount of mass released and the magnetic flux closed in order to understand the role of tail reconnection in the transport of mass and flux in Jupiter's magnetosphere. The observed plasmoid properties are consistent with a mass loss rate of ~0.7-120 kg/s and a flux closure rate of ~7-70 GWb/day. We conclude that tail reconnection and plasmoid release is an important method of flux transport at Jupiter but likely cannot account for the mass input from Io, suggesting that additional mass loss mechanisms may be significant. Finally, we examine the plasmoid interior structure through minimum variance analysis and find that most plasmoids lack a core field and are better described by magnetic loops rather than flux ropes.

Description: [1]  Hot flow anomalies (HFAs) represent a subset of solar wind discontinuities interacting with collisionless bow shocks. They are typically formed when the normal component of the motional (convective) electric field points toward the embedded current sheet on at least one of its sides. The core region of an HFA contains hot and highly deflected ion flows and rather low and turbulent magnetic field. In this paper, we report observations of possible HFA-like events at Mercury identified over a course of two planetary years. Using data from the orbital phase of the MErcury Surface, Space ENvironment, GEochemistry, and Ranging (MESSENGER) mission, we identify a representative ensemble of active current sheets magnetically connected to Mercury's bow shock. We show that some of these events exhibit magnetic and particle signatures of HFAs similar to those observed at other planets, and present their key physical characteristics. Our analysis suggests that Mercury's bow shock does not only mediate the flow of supersonic solar wind plasma but also provides conditions for local particle acceleration and heating as predicted by previous numerical simulations. Together with earlier observations of HFA activity at Earth, Venus, Mars, and Saturn, our results suggest that hot flow anomalies could be a common property of planetary bow shocks, and show that the characteristic size of these events is controlled by the bow shock standoff distance and/or local solar wind conditions.

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]  This paper investigates the influence of two solar eclipses on the ionosphere complexity measures: Tsallis entropy, Renyi entropy, Hurst exponent, beta exponent, fractal dimension. The study used GPS TEC measured at 3 locations in Japan during the solar eclipses of 22 July 2009 and 21 May 2012.This is the first effort to compare the complexity measures by comparing TEC time series of the eclipse day with those from the day before and day after the eclipse. It was found from analysis of the TEC observations that there were no abnormal variations of the complexity parameters from their expected values for either eclipse. Model calculations also show that TEC deviations during the eclipses are small.

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: 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]  We have analyzed the data of the world neutron monitor network for the first ground level enhancement of solar cycle 24, the GLE on May 17, 2012. A newly computed neutron monitor yield function and an inverse method are applied to estimate the energy spectrum, anisotropy axis direction and pitch-angle distribution of the high-energy solar particles in interplanetary space. The method includes the determination of the asymptotic viewing cones of neutron monitor stations through computations of trajectories of cosmic rays in a model magnetosphere. The cosmic ray particle trajectories are determined with the MAGNETOCOSMICS code using Tsyganenko 1989 and IGRF models. Subsequent calculation of the neutron monitor responses with the model function is carried out, that represents an initial guess of the inverse problem. Derivation of the solar energetic particle characteristics is fulfilled by fitting the data of the global neutron monitor network using the Levenberg-Marquardt method over the nine-dimensional parameter space. The pitch-angle distribution and rigidity spectrum of high-energy protons are obtained as function of time in the course of the GLE. The angular distribution appears quite complicated. It comprises a focused beam along the interplanetary magnetic field line from the Sun and a loss-cone feature around the opposite direction, possibly indicative of the particle transport in interplanetary magnetic field structures associated with previous coronal mass ejections.

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]  Geomagnetic activity is strongly controlled by solar wind and Interplanetary Magnetic Field (IMF) conditions, especially the southward component of IMF (IMF Bs). We analyze the statistical properties of IMF Bs at 1 AU using in situ observations for more than a solar cycle (1995 - 2010). IMF Bs-events are defined as continuous IMF Bs intervals with varying thresholds of Bs magnitude and duration, and categorized by different solar wind structures, such as magnetic cloud (MC), interplanetary small-scale magnetic flux rope (ISMFR), interplanetary coronal mass ejection (ICME) without MC signature (ejecta), stream interacting region (SIR), and shock, as well as events unrelated with well-defined solar wind structures. The statistical properties of IMF Bs-events and their geoeffectiveness are investigated in detail based on satellite and ground measurements. We find that the integrated duration and number of Bs-events follow the sunspot number when Bz 〈 -5 nT. We also find that in extreme Bs-events (t 〉 6 hours, Bz 〈 -10 nT), a majority (53 %) are related to MC and 10 % are related with ejecta, but nearly a quarter are not associated with any well-defined solar wind structure. We find different geomagnetic responsesfor Bs-events with comparable duration and magnitude depending on what type of solar wind structures they are associated with. We also find that great Bs-events (t 〉 3 hours, Bz 〈 -10 nT) do not always trigger magnetic storms.

Description: [1]  This technique paper describes a novel method for quantitatively and routinely identifying auroral breakup following substorm onset using the Time History of Events and Macroscale Interactions During Substorms (THEMIS) all-sky imagers (ASIs). Substorm onset is characterised by a brightening of the aurora that is followed by auroral poleward expansion and auroral breakup. This breakup can be identified by a sharp increase in the auroral intensity i(t) and the time derivative of auroral intensity i'(t) . Utilising both i(t) and i'(t) we have developed an algorithm for identifying the time interval and spatial location of auroral breakup during the substorm expansion phase within the field of view of ASI data based solely on quantifiable characteristics of the optical auroral emissions. We compare the time interval determined by the algorithm to independently identified auroral onset times from three previously published studies. In each case the time interval determined by the algorithm is within error of the onset independently identified by the prior studies. We further show the utility of the algorithm by comparing the breakup intervals determined using the automated algorithm to an independent list of substorm onset times. We demonstrate that up to 50% of the breakup intervals characterised by the algorithm are within the uncertainty of the times identified in the independent list. The quantitative description and routine identification of an interval of auroral brightening during the substorm expansion phase provides a foundation for unbiased statistical analysis of the aurora to probe the physics of the auroral substorm as a new scientific tool for aiding the identification of the processes leading to auroral substorm onset.

Description: [1]  We applied the Grad-Shafranov (GS) reconstruction technique to Martian magnetic flux ropes observed by Mars Global Surveyor in order to estimate their spatial structures. This technique can provide a magnetic field map of their cross section from single spacecraft data, under the assumption that the structure is two-dimensional, magneto-hydrostatic, and time-independent. We succeeded in recovering the spatial structure for 70 events observed between April 1999 and November 2006. The reconstruction results indicate that the flux rope axes were mostly oriented horizontal to the Martian surface, and were randomly distributed with respect to the typical plasma streamline. A subset of events with duration longer than 240 sec was observed at solar zenith angles larger than 75 deg. These events all occur downstream from strong crustal magnetic field in the southern hemisphere, indicating an association between the crustal fields and the detected flux ropes. Using the shape and size of the flux ropes obtained from the GS reconstruction, we estimate lower limits on their volume that span 2–3 orders of magnitude, with larger flux ropes observed downstream from strong crustal magnetic fields. Estimated ion escape rates associated with flux ropes are of the order of 10 22 –10 23 ion/sec, being approximately 10% of previously estimated escape rates during solar minimum.

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]  From a survey of Polar plasma waves conducted over the interval 1 April, 1996 to 4 April, 1997 (during solar minimum) at and inside the plasmasphere, magnetosonic waves were detected at all local times with a slight preference of occurrence in the midnight-postmidnight sector at L = 3 to 4. The waves occurred primarily during heightened geomagnetic (AE) activity. Wave occurrence (and intensities) peaked at ~ ±5° of the magnetic equator, with half-maxima at ~ ±10°. For other wave events, magnetosonic waves were also detected as far from the equator as +20° and -60° MLAT, but at lower intensities. An extreme magnetosonic wave intensity event of amplitude B w  = ~ ± 1 nT and E w  = ~ ± 25 mV/m was detected during the survey period. The event occurred near local midnight (0022 MLT), at the magnetic equator (MLAT = -0.5°), at the plasmapause (L = 3.5), and during an intense substorm/convection event (AE = 624 nT; SYM-H = -33 nT). If more stringent requirements (| MLAT| ≤ 5° and AE 〉 300 nT) are imposed, the wave occurrence rate approaches ~50% for the 23 to 00 MLT bin at L = 3 to 4. This strong local time anisotropy in the location of magnetosonic wave occurrence rate supports the idea of generation by protons injected from the plasmasheet into the midnight sector magnetosphere by substorm electric fields. Magnetosonic waves were also detected near late morning (1031 MLT) during relative geomagnetic quiet (low AE). We mention that one possible generation mechanism is a recovering/expanding plasmasphere engulfing preexisting energetic ions, which in turn leads to an ion instability. The wave magnetic component oscillations are aligned along B 0 , the ambient magnetic field direction, and the electric component oscillations are orthogonal to B 0 , indicating linear polarization. The magnetosonic wave amplitudes decreased at locations further from the magnetic equator, while transverse whistler mode wave amplitudes increased. We argue that intense magnetosonic waves are always present somewhere in the magnetosphere during strong substorm/convection events. We thus suggest that modelers use dynamic particle tracing codes and the maximum (rather than average) wave amplitudes to simulate wave-particle interactions.

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.