ALBERT

Beschreibung: 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.

Beschreibung: 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.

Beschreibung: 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.

Beschreibung: To study how an impacting plume modifies the mantle lithosphere, we analyzed the microstructures and crystal preferred orientations (CPO) of 29 peridotites and 37 pyroxenites that sample the mantle root of the Ontong Java Plateau (OJP) from 60 to 120 km depth. The peridotites show a strong compositional variability, but homogeneous coarse-granular to tabular microstructures, except for those equilibrated at the shallowest and deepest depths, which are porphyroclastic. All peridotites have clear olivine CPO, with dominant fiber-[010] patterns. Low intragranular misorientations and straight grain boundaries in olivine suggest that, above 100 km depth, annealing often followed deformation. Calculated density and P-wave velocities of the peridotites decrease weakly with depth. S-wave velocities decrease faster, resulting in increasing Vp/Vs ratio with depth. Calculated densities and seismic velocity profiles are consistent with those estimated for normal mantle compositions under a cold oceanic geotherm. Enrichment in pyroxenites may further increase seismic velocities. The calculated seismic properties cannot therefore explain the low S-waves velocities predicted by Rayleigh wave tomography and ScS data in the mantle beneath the OJP. Calculated P- and S-wave anisotropy is variable (2-12%). It is higher on average in the deeper section of the lithosphere. Because olivine has dominantly [010]-fiber CPO patterns, if foliations are horizontal, vertically propagating S-waves and Rayleigh waves will sample very weak anisotropy in the OJP mantle lithosphere. Moreover, if the orientation of the lineation changes with depth, the anisotropy-induced contrast in seismic properties might produce an intralithospheric reflector marking the stratification of the OJP mantle root.

Beschreibung: The quantification of heat and mass flow between deep reservoirs and the surface is important for understanding magmatic and hydrothermal systems. Here, we use high-resolution measurement of carbon dioxide flux (ϕCO 2 ) and heat flow at the surface to characterize the mass (CO 2 and steam) and heat released to the atmosphere from two magma-hydrothermal systems. Our soil gas and heat flow surveys at Rotokawa and White Island in the Taupō Volcanic Zone, New Zealand, include over 3,000 direct measurements of ϕCO 2 and soil temperature and 60 carbon isotopic values on soil gases. Carbon dioxide flux was separated into background and magmatic/hydrothermal populations based on the measured values and isotopic characterization. Total CO 2 emission rates (ΣCO 2 ) of 441 ± 84 t d -1 and 124 ± 18 t d -1 were calculated for Rotokawa (2.9 km 2 ) and for the crater floor at White Island (0.3 km 2 ), respectively. The total CO 2 emissions differ from previously published values by +386 t d -1 at Rotokawa and +25 t d -1 at White Island, demonstrating that earlier research underestimated emissions by 700% (Rotokawa) and 25% (White Island). These differences suggest that soil CO 2 emissions facilitate more robust estimates of the thermal energy and mass flux in geothermal systems than traditional approaches. Combining the magmatic/hydrothermal-sourced CO 2 emission (constrained using stable isotopes) with reservoir H 2 O:CO 2 mass ratios and the enthalpy of evaporation, the surface expression of thermal energy release for the Rotokawa hydrothermal system (226 MW t ) is 10 times greater than the White Island crater floor (22.5 MW t ).

Beschreibung: 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.

Beschreibung: 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.

Beschreibung: 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.

Beschreibung: 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.

Beschreibung: 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).

Beschreibung: To better constrain the mechanical behavior of sediments accreted to accretionary prism, we conducted triaxial mechanical tests on natural samples from the Miura-Boso paleo-accretionary prism (Japan) in drained conditions with confining pressures up to 200 MPa as well as post-experiments P-wave velocity (V p ) measurements. During experiments, deformation is principally non-coaxial and accommodated by two successive modes of deformation, both associated with strain-hardening and velocity-strengthening behavior: (1) compaction-assisted shearing, distributed in a several mm-wide shear zone and (2) faulting, localized within a few tens of µm-wide, dilatant fault zone. Deformation is also associated with (1) a decrease in Young's modulus all over the tests, (2) anomalously low V p in the deformed samples compared to their porosity and (3) an increase in sensitivity of V p to effective pressure. We interpret this evolution of the poroelastic properties of the material as reflecting the progressive breakage of intergrain cement and the formation of microcracks along with macroscopic deformation. When applied to natural conditions, these results suggest that the deformation style (localized vs distributed) of shallow (z 〈 a few km) sediments is mainly controlled by the variations in stress/strain rate during the seismic cycle and is therefore independent of the porosity of sediments. Finally, we show that the effect of strain, through cement breakage and microcracks formation, may lower V p for effective pressure up to 40 MPa. As a consequence, the low V p anomalies observed in Nankai accretionary prisms by seismic imaging between 2 and 4km depth [ Kitajima and Saffer , 2012] could reflect sediment deformation rather than porosity anomalies. This article is protected by copyright. All rights reserved.

Beschreibung: This study presents a new thermodynamic model for the calculation of phase relations during the melting of anhydrous spinel lherzolite at pressures of 1–2.5 GPa. The model is based on the total energy minimization algorithm for calculating phase equilibria within multicomponent systems and the thermodynamic configuration of Ueki and Iwamori [2013]. The model is based on a SiO 2 – Al 2 O 3 – FeO–Fe 3 O 4 –MgO–CaO system that includes silicate melt, olivine, clinopyroxene, orthopyroxene, and spinel as possible phases. The molar Gibbs free energy of the melt phase is modeled quasi-empirically, and the thermodynamic parameters for silicate melt end-member components are calibrated with a polybaric calibration database. The temperatures and pressures used in this newly compiled calibration dataset are 1230–1600 ∘ C and 0.9–3 GPa, corresponding to the stability range of spinel lherzolite. The modeling undertaken during this study reproduces the general features of experimentally determined melting phase relations of spinel lherzolite at 1–2.5 GPa, including the solidus temperature, the melt composition, the chemical reaction during melting and the degree of melting. This new thermodynamic modeling also reproduces phase relations of various bulk compositions from fertile to deplete spinel lherzolite and can be used in the modeling of polybaric mantle melting within various natural settings. Comparing the results derived from this new modeling with those produced using previous models indicates that the new approach outlined here, involving a combination of total energy minimization and the direct calibration of melt thermodynamic parameters at pressure and temperature conditions corresponding to mantle melting with a relatively simple melt thermodynamic equation, can accurately model polybaric melting phase relations. This article is protected by copyright. All rights reserved.

Beschreibung: 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.

Beschreibung: 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.

Beschreibung: [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.

Beschreibung: [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.

Beschreibung: Knowledge on the behavior of Mg isotopes during metamorphic dehydration is the prerequisite for applying Mg isotopes as tracers for crustal recycling. Here, we report Mg isotopic compositions of metapelites from the Onawa contact aureole, Maine. Except one sample, all metapelites across the aureole, from the wall-rock regional metamorphic rocks to the partially melted rocks adjacent to the pluton, have similar Mg isotopic compositions (δ 26 Mg = -0.09 to +0.12‰). This observation indicates limited Mg isotope fractionation during metamorphic dehydration and fluid-rock interaction, due to the low Mg concentration in fluids relative to rocks. Our results suggest that Mg isotopic compositions of metapelites can record those of their protoliths and, hence, recycled clastic sedimentary materials may preserve their low-temperature Mg isotopic signatures through subduction zones. Therefore, Mg isotopes may serve as new tracers for crustal recycling, for example, tracing components experienced weathering cycles within granite sources.

Beschreibung: [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.

Beschreibung: [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.

Beschreibung: [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.

Beschreibung: [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.

Beschreibung: [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.

Beschreibung: [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.

Beschreibung: Magnetic anomaly identifications underpin plate tectonic reconstructions and form the primary dataset from which age of the oceanic lithosphere and seafloor spreading regimes in the ocean basins can be determined. Although these identifications are an invaluable resource, their usefulness to the wider scientific community has been limited due to the lack of a central community infrastructure to organize, host and update these interpretations. We have developed an open-source, community-driven online infrastructure as a repository for quality-checked magnetic anomaly identifications from all ocean basins. We provide a global sample dataset that comprises 96,733 individually picked magnetic anomaly identifications organized by ocean basin and publication reference, and provide accompanying Hellinger-format files, where available. Our infrastructure is designed to facilitate research in plate tectonic reconstructions or research that relies on an assessment of plate reconstructions, for both experts and non-experts alike. To further enhance the existing repository and strengthen its value, we encourage others in the community to contribute to this effort.

Beschreibung: [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.

Beschreibung: [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.

Beschreibung: [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.

Beschreibung: We present an efficient method for high-volume heavy mineral separation from clay-rich rocks using an ultrasonic probe. The ultrasonic clay separator (UCS) is an easily constructed device that allows for the recovery of high-density minerals, as small as 10 microns, with a minimum of sample preparation. Heavy mineral recovery from clay-rich material with the UCS yields a greater number of small (〈100 micron) grains and approximately double the amount of material from that of gravity settling and decanting. Despite development with heavy mineral recovery in mind, the UCS should suitable for recovering small grain size geologic materials from flocculating clay-rich material.

Beschreibung: [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.

Beschreibung: [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.

Beschreibung: [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.

Beschreibung: [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.

Beschreibung: [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.

Beschreibung: In the detachment mode of slow seafloor spreading, convex-upward detachment faults take up a high proportion of the plate separation velocity exposing gabbro and serpentinized peridotite on the seafloor. Large, long-lived hydrothermal systems such as TAG are situated off-axis and may be controlled by fluid flow up a detachment fault, with the source of magmatic heat being as deep as 7 kmbsf. The consequences of such deep circulation for the evolution of fluid temperature and salinity have not previously been investigated. Microthermometry on fluid inclusions trapped in diabase, gabbro and trondjhemite, recovered at the Atlantis Massif Oceanic Core Complex (30° N, mid-Atlantic Ridge), reveals evidence for magmatic exsolution, phase separation, and mixing between hydrothermal fluids and previously phase separated fluids. Four types of fluid inclusions were identified, ranging in salinity from 1.4 wt.% NaCl to 35 wt.% NaCl, although the most common inclusions have salinities close to seawater (3.4 wt.% NaCl). Homogenization temperatures range from 160 to 〉400 °C, with the highest temperatures in hypersaline inclusions trapped in trondjhemite and the lowest temperatures in low salinity inclusions trapped in quartz veins. The fluid history of the Atlantis Massif is interpreted in the context of published thermochronometric data from the Massif, and a comparison with the inferred circulation pattern beneath the TAG hydrothermal field, to better constrain the pressure temperature conditions of trapping and when in the history of exhumation of the rocks sampled by IODP Hole U1309D fluids have been trapped.

Beschreibung: [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.

Beschreibung: [1]  On January 9, 2002 and November 14, 2001, the São Luís 30 MHz coherent backscatter radar observed unusual day-time echoes scattered from the equatorial electrojet. The electrojet echoing layers on these days, as seen in the range time intensity (RTI) maps, exhibited quasi-periodic oscillations. Time-frequency decomposition of the magnetic field perturbations ΔH , measured simultaneously by the ground-based magnetometers, also showed evidence of short period waves. The ground-based observations were aided by measurements of the brightness temperature in the water vapor and infrared bands made by the GOES 8 satellite. The GOES 8 satellite measurements indicated evidence of deep tropospheric convection activities, which are favorable for the launch of atmospheric gravity waves (AGW) near São Luís. Our multi-technique investigation, combined with an analysis of the equatorial electric field and current density, indicates that AGW forcing could have been responsible, via coupling with E-region electric fields, for the short-period electrojet oscillations observed over São Luís.

Beschreibung: [1]  The diurnal variation of the global electric circuit is investigated using the World Wide Lightning Location Network (WWLLN), which has been shown to identify nearly all thunderstorms ( [16], using WWLLN data from 2005). To create an estimate of global electric circuit activity, a clustering algorithm is applied to the WWLLN dataset to identify global thunderstorms from 2010 – 2013. Annual, seasonal, and regional thunderstorm activity is investigated in this new WWLLN thunderstorm dataset in order to estimate the source behavior of the global electric circuit. Through the clustering algorithm, the total number of active thunderstorms are counted every 30 minutes creating a measure of the global electric circuit source function. The thunderstorm clusters are compared to precipitation radar data from the Tropical Rainfall Measurement Mission satellite and with case studies of thunderstorm evolution. [2]  The clustering algorithm reveals an average of 660 ± 70 thunderstorms active at any given time with a peak-to-peak variation of 36%. The highest number of thunderstorms occurs in November (720 ± 90) and the lowest number occurs in January (610 ± 80). Thunderstorm cluster and electrified storm cloud activity are combined with thunderstorm overflight current measurements to estimate the global electric circuit thunderstorm contribution current to be 1090 ± 70 A with a variation of 24%. By utilizing the global coverage and high time resolution of WWLLN, the total active thunderstorm count and current is shown to be less than previous estimates based on compiled climatologies.

Beschreibung: [1]  Dynamics of the dayside magnetosphere and proton radiation belt was analyzed during unusual magnetic storm on 21 January 2005. We have found that during the storm from 1712 to 2400 UT, the subsolar magnetopause was continuously located inside geosynchronous orbit due to strong compression. The compression was found to be extremely strong from 1846 to 2035 UT when the dense plasma of fast erupting filament produced the solar wind dynamic pressure Pd peaked up to 〉100 nPa and, in the first time, the upstream solar wind was observed at geosynchronous orbit during almost 2 hours. Under the extreme compression, the outer magnetosphere at L  〉 5 was pushed inward and the outer radiation belt particles with energies of several tens of keV moved earthward, became adiabatically accelerated and accumulated in the inner magnetosphere at L  〈 4 that produced the intensified ring current with an exceptionally long lifetime. The observations were compared with predictions of various empirical and first principles models. All the models failed to predict the magnetospheric dynamics under the extreme compression when the minimal magnetopause distance was estimated to be ~3 Re. The inconsistencies between the model predictions and observations might result from distortions of plasma measurements by extreme heliospheric conditions consisting in very fast solar wind streams (~1000 km/s) and intense fluxes of solar energetic particles. We speculated that anomalous dynamics of the magnetosphere could be well described by the models if the He abundance in the solar wind was assumed to be 〉20%, which is well appropriate for erupting filaments and which is in agreement with the upper 27% threshold for the He/H ratio obtained from Cluster measurements.

Beschreibung: We study the elastic wave speed structure of the crust and the uppermost mantle in western Tibet using P and S-wave arrival times from regional earthquakes recorded by a temporary seismic network. We relocate the earthquakes, and subsequently invert travel time residuals for 3D distributions of wave speed. Resolution tests with a variety of input structures are used to verify the reliability of our results. The crust beneath western Tibet has low P-wave speed (5.9 - 6.3 km/s) throughout its nearly 80 km thickness, with lower values in this range concentrated within the Lhasa block. Beneath the Himalaya wave speeds are higher. Southern and western limits of the slow material beneath the Tibetan Plateau correlate with the Karakoram fault, and dip beneath the plateau at ~40° angle. We find no evidence of a sub-horizontal low velocity zone in the crust. In the uppermost mantle we find a long and narrow region of fast (up to 8.4 km/s) P-wave speed extending from the Karakoram fault in NE direction, and crossing the Bangong-Nujiang suture. In a north-south cross-section, the distribution of relatively fast P-wave speed suggests a ramp-flat geometry consistent with India underthrusting the Tibetan Plateau at least as far as 32.5°N. A plausible interpretation of the upper mantle fast feature is the formation of eclogite from the mafic lower-crustal material of India after it is underthrust beneath Tibet. Notably, in western Tibet this process only takes place in a narrow region.

Beschreibung: We have employed shear wave splitting techniques to image anisotropy beneath the D'Entrecasteaux Islands, in southeastern Papua New Guinea. Our results provide a detailed picture of the extending continent that lies immediately ahead of a propagating mid-ocean ridge tip; we image the transition from continental to oceanic extension. A dense shear wave splitting dataset from a 2010-11 passive-source seismic deployment is analyzed using single- and multi-channel methods. Splitting delay times of 1-1.5 s are observed and fast axes of anisotropy trending N-S, parallel to rifting direction, predominate the results. This trend is linked to lattice-preferred orientation of olivine, primarily in the shallow convecting mantle, driven by up to 200 km of N-S continental extension ahead of the westward-propagating Woodlark Rift. This pattern differs from several other continental rifts that evince rift-strike-parallel fast axes and is evident despite the complex recent tectonic history. We contend that across most of this rift, the unusually high rate and magnitude of extension has been sufficient to produce a regime change to a mid-ocean-ridge-like mantle fabric. Stations in the south of our array show more complex splitting that might be related to melt or to complex inherited structure at the edge of the extended region.

Beschreibung: From April 2010 through February 2011, CO 2 flux surveys were performed on Lake Rotomahana, New Zealand. The area has been hydrothermally active with fumaroles and sublacustrine hydrothermal activity before and since the eruption of Mt Tarawera in 1886. The total CO 2 emission from the lake calculated by sequential Gaussian simulation is 549 ± 72 t day -1 . Two different mechanisms of degassing, diffusion through the water-air interface and bubbling, are distinguished using a graphical statistical approach. The carbon dioxide budget calculated for the lake confirms that the main source of CO 2 to the atmosphere is by diffusion covering 94.5 % of the lake area (mean CO 2 flux 25 g m -2 day -1 ) and to a lesser extent, bubbling (mean CO 2 flux 1297 g m -2 day -1 ). Mapping of the CO 2 flux over the entire lake, including over lakefloor vents detected during the survey, correlates with eruption craters formed during the 1886 eruption. These surveys also follow regional tectonic patterns present in the southeastern sector of Lake Rotomahana suggesting a deep magmatic source (~ 10 km) for CO 2 and different pathways for the gas to escape to the surface. The values of δ 13 C CO2 (-2.88 and -2.39 ‰) confirm the magmatic origin of CO 2 .

Beschreibung: [1]  Type III radio bursts are produced near the local electron plasma frequency and/or near its harmonic by fast electrons ejected from the solar active regions and moving through the corona and solar wind. These bursts have dynamic spectra with frequency rapidly falling with time. This paper presents two new methods developed to detect type III bursts automatically in the data from High Frequency Receiver (HFR) of the STEREO/WAVES (S/WAVES) radio instrument onboard the STEREO spacecraft. The first technique is applicable to the low frequency band (HFR-1: 125 kHz to 1.975 MHz) only. This technique can possibly be implemented in on-board satellite software aimed at preliminary detection of bursts and identification of time intervals with relatively high solar activity. In the second technique the bursts are detected in both the low frequency band and the high frequency band (HFR-2: 2.025 MHz to 16.025 MHz), with the computational burden being higher by one order of magnitude as compared with that for the first technique. Preliminary tests of the method show that the performance of the first technique is quite high, P dL  = 72 % ± 3 %. The performance of the second technique is considerably higher, P dL  +  H  = 81 % ± 1%, while the number of false alarms does not exceed 10% for one daily spectrum.

Beschreibung: [1]  A dipolarizing flux bundle (DFB) is a small magnetotail flux tube (typically 〈 ~3 R E in X GSM and Y GSM ) with a significantly more dipolar magnetic field than its background. Dipolarizing flux bundles typically propagate earthward at a high speed from the near-Earth reconnection region. Knowledge of a DFB's flux transport properties leads to better understanding of near-Earth (X = -6 to -30 R E ) magnetotail flux transport and thus conversion of magnetic energy to kinetic and thermal plasma energy following magnetic reconnection. We explore DFB properties with a statistical study using data from the Time History of Events and Macroscale Interactions during Substorms (THEMIS) mission. To establish the importance of DFB flux transport, we compare it with transport by bursty bulk flows (BBFs) that typically envelop DFBs. Because DFBs coexist with flow bursts inside BBFs, they contribute 〉65% of BBF flux transport, even though they last only ~30% as long as BBFs. The rate of DFB flux transport increases with proximity to Earth and to the pre-midnight sector, as well as with geomagnetic activity and distance from the neutral sheet. Under the latter two conditions the total flux transport by a typical DFB also increases. Dipolarizing flux bundles appear more often during increased geomagnetic activity. Since BBFs have been previously shown to be the major flux transporters in the tail, we conclude that DFBs are the dominant drivers of this transport. The occurrence rate of DFBs as a function of location and geomagnetic activity informs us about processes that shape global convection and energy conversion

Beschreibung: [1]  The substorm current wedge (SCW) is a fundamental component of geomagnetic substorms. Models tend to describe the SCW as a simple line current flowing into the ionosphere towards dawn and out of the ionosphere towards dusk, linked by a westward electrojet. We use multi-spacecraft observations from perigee passes of the Cluster 1 and 4 spacecraft during a substorm on 15 Jan 2010, in conjunction with ground-based observations, to examine the spatial structuring and temporal variability of the SCW. At this time, the spacecraft travelled east-west azimuthally above the auroral region. We show that the SCW has significant azimuthal sub-structure on scales of 100 km at altitudes of 4,000-7,000 km. We identify 26 individual current sheets in the Cluster 4 data and 34 individual current sheets in the Cluster 1 data, with Cluster 1 passing through the SCW 120-240 s after Cluster 4 at 1,300-2,000 km higher altitude. Both spacecraft observed large-scale regions of net upward and downward field-aligned current, consistent with the large-scale characteristics of the SCW, although sheets of oppositely directed currents were observed within both regions. We show that the majority of these current sheets were closely aligned to a north-south direction,in contrast to the expected east-west orientation of the pre-onset aurora. Comparing our results with observations of the field-aligned current associated with bursty bulk flows (BBFs) we conclude that significant questions remain for the explanation of SCW structuring by BBF driven “wedgelets”. Our results therefore represent constraints on future modelling and theoretical frameworks on the generation of the SCW.

Beschreibung: [1]  Large horizontal winds and wind shears have been measured in the lower thermosphere by rockets, lidars, and non-specular meteor radars. This paper describes a detailed analysis of 3 multi-hour non-specular meteor radar data sets collected at the Jicamarca Radio Observatory. This provides some of the highest resolution sustained measurements in this part of the atmosphere. These show: (1) intense wind speeds, maintaining 180 m/s for half an hour and 160 m/s for another half an hour; (2) winds structured in layers that move up or, more commonly, down in the pre-dawn hours at rates of a few km/hr; (3) intense wind shears that typically persist at around 50 m/s/km but, in one instance, sustains values approaching 100 m/s/km for a few hours.

Beschreibung: [1]  We present results from an analysis of high-latitude ionosphere-thermosphere (IT) coupling to the solar wind during a moderate magnetic storm which occurred on 5-6 August 2011. During the storm, a multi-point set of observations of the ionosphere and thermosphere was available. We make use of ionospheric measurements of electromagnetic and particle energy made by the Defense Meteorological Satellite Program (DMSP), and neutral densities measured by the Gravity Recovery and Climate Experiment (GRACE) satellite to infer: (1) the energy budget and (2) timing of the energy transfer process during the storm. We conclude that the primary location for energy input to the IT system may be the extremely high latitude region. We suggest that the total energy available to the IT system is not completely captured either by observation or empirical models.

Beschreibung: [1]  The level of solar activity varies from cycle to cycle. This variability is probably caused by a combination of nonlinear and random effects. Based on surface flux transport simulations, we show that the observed inflows into active regions and towards the activity belts provide an important nonlinearity in the framework of Babcock-Leighton model for the solar dynamo. Inclusion of these inflows also leads to a reproduction of the observed proportionality between the open heliospheric flux during activity minima and the maximum sunspot number of the following cycle. A substantial component of the random variability of the cycle strength is associated with the cross-equatorial flux plumes that occur when large, highly-tilted sunspot groups emerge close to the equator. We show that the flux transported by these events is important for the amplitude of the polar fields and open flux during activity minima. The combined action of inflows and cross-equatorial flux plumes provides an explanation for the weakness of the polar fields at the end of solar cycle 23 (and hence for the relative weakness of solar cycle 24).

Beschreibung: [1]  We present the first direct measurement of neutral oxygen in the lunar exosphere, detected by the Chandrayaan-1 Energetic Neutral Analyzer (CENA). With the lunar surface consisting of about 60% of oxygen in number, the neutral oxygen detected in CENA's energy range (11 eV – 3.3 keV) is attributed to have originated from the lunar surface, where it is released through solar wind ion sputtering. We verify this proposition by comparing the measured oxygen content in two different mass spectra groups with ion sputtering theory. One group contains mass spectra that were recorded when the solar wind consisted of almost pure hydrogen and the other group contains mass spectra that were recorded when the helium content in the solar wind was very high (〉3.5%). Since helium is a much more effective sputtering agent than hydrogen (5% of alpha particles present in the solar wind typically contribute 30% of the total sputter yield), these two groups should show clear differences in the oxygen sputter yield. Fitting of CENA's mass spectra with calibration spectra from ground and in-flight data resulted in the detection of a robust oxygen signal, with a flux of 0.2 to 0.4 times the flux of backscattered hydrogen, depending, as expected, on the solar wind helium content and particle velocity. These measurements present the first in-situ detection of oxygen in the lunar exosphere. For the two solar wind types observed, we derive sub-solar surface oxygen atom densities of N 0  = (1.1 ± 0.3) ⋅ 10 7 m − 3 and (1.4 ± 0.4) ⋅ 10 7 m − 3 , respectively, which agree well with earlier model predictions and measured upper limits. From these surface densities we derive, by modeling, column densities of N C  = (1.5 ± 0.5) ⋅ 10 13 m − 2 and (1.6 ± 0.5) ⋅ 10 13 m − 2 . [2]  In addition, in the CENA mass spectra, we identified for the first time a helium component. This helium is attributed to backscattering of solar wind helium (alpha particles) from the lunar surface as neutral energetic helium atoms, which has been observed for the first time. This identification is supported by the characteristic energy of the measured helium atoms, which is roughly four times the energy of reflected solar wind hydrogen, and the correlation with solar wind helium content.

Beschreibung: [1]  The diffuse aurora is an almost permanent feature in the Earth's upper polar atmosphere, providing the major source of ionizing energy input into the high-latitude region. Previous theoretical and observational studies have demonstrated that whistler-mode chorus scattering primarily accounts for intense nightside diffuse auroral precipitation within ~ 8 R E , but what causes the dayside diffuse aurora remains poorly understood. Using conjugate satellite wave and particle observations on 13 August 2009 from the THEMIS spacecraft and ground-based all-sky imager measurements at the South Pole on the dayside, we perform a quantitative analysis of wave driven diffusion and electron precipitation. Our results demonstrate that the dayside chorus scattering was the dominant contributor to the observed dayside diffuse auroral precipitation and that the chorus wave intensity primarily controlled its brightness, indicating that dayside chorus can be the major driver of the Earth's dayside diffuse aurora. While further investigations are required to bring closure to the origin of the dayside diffuse aurora under differing solar wind conditions and geomagnetic situations, our finding is an important complement to recent work on the formation mechanism of the diffuse aurora and provides improved understanding of the roles of resonant wave-particle interactions in diffuse auroral precipitation pattern on a global scale.

Beschreibung: [1]  Ultra Low Frequency (ULF: 0.001-5 Hz) magnetic records have recently been used in the search for short term earthquake prediction methods. The separation of local and global effects in the magnetic records is the greatest challenge in this research area. Geomagnetic indices are often used to predict global ULF magnetic behavior where it is assumed that increases in a geomagnetic index correspond with an increase in ULF power. This paper examines the relationships between geomagnetic indices and ULF power, spectral polarization ratio and the relationship between the spectral polarization ratio and solar wind parameters. The power in the ULF, Pc3-5 bands (10-600 s) shows a linear correlation coefficient of  0.2 with the Kp magnetic activity index. The correlation varies with magnetic local time (MLT) and latitude. The correlation coefficient is inversely related to the integrated power in the ULF Pc3 band (10-45 s) over MLT and magnetic latitude. The ratio of spectral powers Z ( ω )/ G ( ω ) is discussed and shown to be a promising parameter in the search for earthquake precursor signals in ULF records.

Beschreibung: [1]  Recent testing of a quantitative model describing the classical (region-1-sense, referred to as the R1 current loop) substorm current wedge (SCW) revealed systematic discrepancies between the observed and predicted amplitudes, which suggested us to include additional region-2-sense currents (R2-loop) earthward of the dipolarized region (SCW2L model). In this paper we discuss alternative circuit geometries of the 3d substorm current system and interpret simultaneous observations of the magnetic field dipolarizations by NOAA Geostationary Operational Environmental Satellite (GOES) and by NASA Time History of Events and Macroscale Interactions during Substorms (THEMIS) spacecraft, to quantitatively investigate the SCW2L model parameters. During two cases of a dipole-like magnetotail configuration, the dipolarization/injection front fortuitously stopped at r  ∼ 9 Re for the entire duration of ∼ 30 min-long SCW-related dipolarization within a unique, radially-distributed multi-spacecraft constellation, which allowed us to determine the locations and total currents of both SCW2L loops. In addition, we conducted a survey of the dipolarization amplitudes in events, simultaneously observed at 6.6 Re (GOES) and 11 Re (THEMIS) under a wide range of magnetotail conditions. We infer that the ratio I 2 / I 1 varies in the range 0.2 to 0.6 (median value 0.4) and that the equatorial part of the R2 current loop stays at the distance r  〉 6.6 Re in the case of a dipole-like field geometry ( BZ 0  〉 75 nT at 6.6 Re prior to the onset), but it is located at r  〈 6.6 Re in the case of a stretched magnetic field configuration (with BZ 0  〈 60 nT). Since the ground midlatitude perturbations are sensitive to the combined effect of the R1- and R2-sense current loops with the total current roughly equal to I 1  −  I 2 , the ratio I 2 / I 1 becomes an important issue when attempting to monitor the current disruption intensity from ground observations.

Beschreibung: Regional-scale geologic structures characteristic of mantle lithosphere within cratons found in continent interiors are interpreted using geo-registered diverse data sets from the Slave craton of northwest Canada. We developed and applied a new method for mapping seismic discontinuities in three dimensions using multi-year observations at sparse, individual broadband receivers. New, fully 3-D conductivity models used all available magnetotelluric data. Discontinuity surfaces and conductivity models were geo-registered with previously published P-wave and surface wave velocity models to confirm first-order structures such as a mid-lithosphere discontinuity. Our 3-D model to 400 km depth was calibrated by ‘drill hole’ observations derived from xenolith suites extracted from kimberlites. A number of new structural discontinuities emerge from direct comparison of co-registered data sets and models. Importantly we distinguish primary mantle layers from secondary features related to younger metasomatism. Sub-horizontal Slave craton layers with tapered, wedge-shaped margins indicate construction of the craton core at 2.7 Ga by underthrusting and flat stacking of lithosphere. Mapping of conductivity and metasomatism in 3-D, the latter inferred via mineral recrystallization and resetting of isotopic ages in xenoliths, indicates overprinting of the primary layered structures. The observed distribution of relatively conductive mantle at 100–200 km depths is consistent with pervasive metasomatism; vertical ‘chimneys’ reaching to crustal depths in locations where kimberlites erupted or where Au mineralization is known.

Beschreibung: Understanding the nature of the crust has long been a goal for seismologists when imaging the Earth. This is particularly true in volcanic regions where imaging melt storage and migration can have important implications for the size and nature of an eruption. Receiver functions and the H- κ stacking (H κ ) technique are often used to constrain crustal thickness (H) and the ratio of P- to S-wave velocities ( κ ). In this paper I show that it is essential to consider anisotropy when performing H κ . I show that in a medium with horizontally transverse isotropy a strong variation in κ with back azimuth is present which characterises the anisotropic medium. In a vertically transverse isotropic medium no variation in κ with back azimuth is observed, but κ is increased across all back azimuths. Thus, estimates of κ are more difficult to relate to composition than previously thought. I extend these models to melt induced anisotropy and show that similar patterns are observed, but with more significant variations and increases in κ . Based on these observations I develop a new anisotropic H- κ stacking technique which inverts H κ data for melt fraction, aspect ratio and orientation of melt inclusions. I apply this to data for the Afar Depression and show that melt is stored in interconnected stacked sills in the lower crust, which likely supply the recent volcanic eruptions and dike intrusions. This new technique can be applied to any anisotropic medium where it can provide constraints on the average crustal anisotropy.

Beschreibung: A Fortran 90 program to visualize data on the Yin-Yang grid system is developed. The purpose of this study is to provide simulation researchers with a source code as a starting point of their own custom-made visualization tools. A basic but sufficiently diverse set of visualization methods are implemented using a Fortran 90 binding for OpenGL for scalar and vector fields defined or simulated on the Yin-Yang grid.

Beschreibung: We present Cluster observations of wave-particle interactions during an earthward-propagating dipolarization front (DF) and associated fast plasma bulk flows detected at the central current sheet in Earth's magnetotail. During this period, flux tubes behind the DF frequently contain more energetic or hotter ions than did the pre-existing flux tubes ahead of the DF. On the other hand, electrons within the DF flux tubes heat less, or are even colder, than were the pre-existing populations, and are often accompanied by superposed isolated beams. At the same time, electrostatic emissions are strongly enhanced over a wide range of frequencies (up to several times the electron cyclotron frequency) behind the DFs. This low-frequency electrostatic wave power is well correlated with ion energization. From linear theory, we find two wave modes: a high-frequency beam mode and a low-frequency whistler mode that are associated with the electron beam component. We attribute the generation of whistlers to electron beams that persist for a while before undergoing rapid thermalization. The existence of isolated beam components behind DFs detected during the 4-sec Cluster spin period indicates that DFs either provide a continuous source of electron beams or facilitate a physical process that maintains the beams against rapid thermalization. Our analysis suggests that the earthward motion of the DF flux tube, via Fermi acceleration as the magnetic field lines behind the DF shorten, can lead to the persistent electron beams that generate whistler mode waves, which in turn can heat ions. This scenario, by which free energy in electron beams generates waves that then heat ions, accounts for the Cluster observations of different energization behaviors between electrons and ions behind DFs.

Beschreibung: We examine a unique data set from seven Hubble Space Telescope (HST) ‘visits’ that imaged Saturn's northern dayside ultraviolet emissions exhibiting usual circumpolar ‘auroral oval’ morphologies, during which Cassini measured the interplanetary magnetic field (IMF) upstream of Saturn's bow shock over intervals of several hours. The auroras generally consist of a dawn arc extending towards noon centered near ~15º co-latitude, together with intermittent patchy forms at ~10º co-latitude and poleward thereof, located between noon and dusk. The dawn arc is a persistent feature, but exhibits variations in position, width, and intensity, which have no clear relationship with the concurrent IMF. However, the patchy post-noon auroras are found to relate to the (suitably lagged and averaged) IMF B z , being present during all four visits with positive B z and absent during all three visits with negative B z . The most continuous such forms occur in the case of strongest positive B z . These results suggest that the post-noon forms are associated with reconnection and open flux production at Saturn's magnetopause, related to the similarly-interpreted bifurcated auroral arc structures previously observed in this LT sector in Cassini UVIS data, whose details remain unresolved in these HST images. One of the intervals with negative IMF B z , however, exhibits a pre-noon patch of very high latitude emission extending poleward of the dawn arc to the magnetic/spin pole, suggestive of the occurrence of lobe reconnection. Overall, these data provide evidence of significant IMF-dependence in the morphology of Saturn's dayside auroras.

Beschreibung: A three-dimensional (3D) particle-in-cell (PIC) simulation of the whistler anisotropy instability is carried out for a collisionless, homogeneous, magnetized plasma with β e  = 0.10. This is the first 3D PIC simulation of the evolution of enchanced fluctuations from this growing mode driven by an anisotropic electron velocity distribution with T ⊥  e / T ‖ e  〉 1 where ⊥ and ∥ represent directions perpendicular and parallel to the background magnetic field B o , respectively. The early-time magnetic fluctuation spectrum grows with properties reflecting the predictions of linear theory with narrowband maxima at kc / ω e  ≃ 1 and k  ×  B o  = 0, and a wavevector anisotropy in the sense of k ⊥  〈 〈  k ∥ . Here ω e represents the electron plasma frequency. At later times the fluctuations undergo both a forward transfer to shorter wavelengths, also with k ⊥  〈 〈  k ∥ , and an inverse transfer to longer wavelengths with wavevector anisotropy k ⊥  〉 〉  k ∥ . The inverse transfer is consistent with a prediction of nonlinear three-wave coupling theory.

Beschreibung: Auroral ionospheric F-region density depletions observed by PFISR (Poker Flat Incoherent Scatter Radar) during the MICA (Magnetosphere-Ionosphere Coupling in the Alfvén Resonator) sounding rocket campaign are critically examined alongside complementary numerical simulations. Particular processes of interest include cavity formation due to intense frictional heating and Pedersen drifts, evolution in the presence of structured precipitation, and refilling due to impact ionization and downflows. Our analysis uses an ionospheric fluid model which solves conservation of mass, momentum, and energy equations for all major ionospheric species. These fluid equations are coupled to an electrostatic current continuity equation to self-consistently describe auroral electric fields. Energetic electron precipitation inputs for the model are specified by inverting optical data, and electric field boundary conditions are obtained from direct PFISR measurements. Thus, the model is driven in as realistic a manner as possible. Both ISR data and simulations indicate that the conversion of the F-region plasma to molecular ions and subsequent recombination is the dominant process contributing to the formation of the observed cavities, all of which occur in conjuction with electric fields exceeding ∼ 90 mV/m. Furthermore, the cavities often persist several minutes past the point when the frictional heating stops. Impact ionization and field-aligned plasma flows modulate the cavity depth in a significant way, but are of secondary importance to the molecular generation process. Informal comparisons of the ISR density and temperature fits to the model verify that the simulations reproduce most of the observed cavity features to a reasonable level of detail.

Beschreibung: Energetic particle injections in the near-Earth plasma sheet are critical for supplying particles and energy to the inner magnetosphere. Recent case studies have demonstrated a good correlation between injections and transient, narrow, fast flow channels as well as earthward reconnection (dipolarization) fronts in the magnetotail, but statistical observations beyond geosynchronous orbit (GEO) to verify the findings were lacking. By surveying trans-geosynchronous injections using THEMIS, we show that their likely origin is the earthward-traveling, dipolarizing flux bundles following near-Earth reconnection. The good correlation between injections and fast flows, reconnection fronts and impulsive, dawn-dusk electric field increases is not limited to within 12 R E , but extends out to 30 R E . Like near-Earth reconnection, both ion and electron injections are most probable in the pre-midnight sector. Similar to bursty bulk flows (BBFs), injection-time flow speeds are faster with increasing distance from Earth. With faster flows, injection intensity generally increases and extends to higher energy channels. With increased geomagnetic activity, injection occurrence rate increases (akin to that of BBFs) and spectral hardening occurs (κ decreases). The occurrence rate increase within the inner magnetosphere suggests that injections populate the radiation belts more effectively under enhanced geomagnetic activity. Our results are inconsistent with the classical concept of an azimuthally wide injection boundary moving earthward from ~9-12 R E to GEO under an enhanced cross-tail electric field. Rather, particle injection and transport occur along a large range of radial distances due to effects from earthward-penetrating, azimuthally localized, transient, strong electric fields of recently reconnected, dipolarizing flux bundles.

Beschreibung: A method for estimating the vector neutral wind profiles in the mesosphere and lower thermosphere (MLT) region of the upper atmosphere from Arecibo dual-beam incoherent scatter radar data is presented. The method yields continuous estimates of both the altitude-averaged F -region plasma drifts and all three components of the altitude-resolved neutral wind profiles in the MLT using data taken while the Arecibo feed system swings in azimuth. The problem is mixed determined, and its solution is not inherently unique. Second order Tikhonov regularization is used to find solutions consistent with the available data while being minimally structured, additional structure being unsupported by the data. The solution is found using the method of conjugate gradient least squares and sparse matrix mathematics. Example data acquired during an interval of midlatitude spread F are used to illustrate the method. The estimated wind profiles exhibit characteristics broadly consistent with gravity waves but are impulsive, with features that generally persist for less than one and a half wave periods.

Beschreibung: We report the peculiar interaction of two type III bursts observed in the solar wind. As electrons beams propagating on the same magnetic field lines cross, a spectacular depletion of the type III radio emission is observed. We combine observations from the WAVES experiment on board the STEREO mission together with kinetic plasma simulations to study the extinction of type III radio emission resulting from the interaction between two electron beams. The remote observations enable to follow the electron beams in the interplanetary medium and show that the level of radiated radio waves is recovered after the beam crossing. The in situ observations of beam-driven Langmuir waves give evidence for Langmuir decay. The density fluctuations are extracted from in situ observations. The velocity of the beams is independently evaluated from in situ observations of decaying Langmuir waves and remote radio observations. The kinetic simulations show that the level of beam-driven Langmuir waves is reduced as the two beams cross. We show that the slow beam induced a strong reduction of the quasilinear relaxation of the fast beam, limiting the amplitude of the generated Langmuir waves. Moreover, in the case of two electron beams, the lack of Langmuir waves coherence reduces the efficiency of the Langmuir parametric decay. We thus conclude that the observed depletion of the type III radio 5 is independent of the radio emission mechanism, as long as it depends on the Langmuir amplitude and coherence.

Beschreibung: The solar minimum period between solar cycles 23 and 24 was the longest since the beginning of space-based measurements, and many manifestations of solar activity were unusually low. Thermospheric neutral density was about 30% lower than during the previous solar minimum, but changes in the ionosphere between the two solar minima are more controversial. Solar radiation, geomagnetic activity, and anthropogenic increases in greenhouse gases, can all play a role in these changes. In this paper, we address the latter of these potential contributions, the degree to which secular change driven by greenhouse gases, primarily CO 2 , could be responsible for the observed changes. New 3D model simulations find a global mean density decrease at 400 km of 5.8% between the two recent solar minima, which is larger than earlier 1D model results, and in better agreement with observations. From these model simulations and from other observational work, we estimate that the contribution of secular change to global mean neutral density decrease between the two recent solar minima is less than ~6%. The contribution of secular change to the global average decrease of F-region ionosphere peak density ( N m F 2 ) and altitude ( h m F 2 ), near mid-day, is estimated to be 1.5% and 1.5 km, respectively. However, secular changes in the ionosphere exhibit large variations with local time, geographic location, and season. The mid-day change of N m F 2 seen in the model simulations ranged between +6% and -9%, and the change of h m F 2 ranged between +11 km and -11 km, depending on geographic location.

Beschreibung: In this paper we study the planetary magnetic disturbance during the magnetic storm occurring on 05 April 2010 associated with high speed solar wind stream due to a coronal hole following a CME. We separate the magnetic disturbance associated to the ionospheric disturbance dynamo (Ddyn) from the magnetic disturbance associated to the prompt penetration of magnetospheric electric field (DP2). This event exhibits different responses of ionospheric disturbance dynamo in the different longitude sectors (European-African, Asian and American). The strongest effect is observed in the European-African sector. The Ddyn disturbance reduces the amplitude of the daytime H-component at low latitudes during four consecutive days in agreement with the Blanc and Richmond's model of ionospheric disturbance dynamo. The amplitude of Ddyn decreased with time during the four days. We discuss its diverse worldwide effects. The observed signature of magnetic disturbance process in specific longitude sector is strongly dependent on which Earth's side faces the magnetic storms (i.e., there is a different response depending on which longitude sector is at noon when the SSC hits).Finally we determined an average period of 22 hours for Ddyn using wavelet analysis.

Beschreibung: Complex magnetosphere-ionosphere coupling mechanisms result in high latitude irregularities that are difficult to characterize using only Global Navigation Satellite System (GNSS) scintillation measurements. However, GNSS observations combined with physical parameters derived from modeling can be used to study the physics of these irregularities. We have developed a full three dimensional (3D) electromagnetic (EM) wave propagation model called “Satellite-beacon Ionospheric-scintillation Global Model ofthe upper Atmosphere" (SIGMA), to simulate GNSS scintillations. This model eliminates the most significant approximation made by the previous simulation approaches about the correlation length of the irregularity. Thus, for the first time, using SIGMA wecan accomplish scintillation simulations of significantly high fidelity. While the model is global, it is particularly applicable at high latitudes as it accounts for the complicated geometry of the magnetic field lines in these regions. Using SIGMA wesimulate the spatial and temporal variations in the GNSS signal phase and amplitude on the ground. In this paper, we present the model and results from a study to determine the sensitivity of the SIGMA outputs to different input parameters. We have deduced from our sensitivity study that the peak to peak (P2P) power gets most affected by the spectral index and line of sight (LOS) direction, while the P2P phase and standard deviation of the phase ( σ φ ) are more sensitive to the anisotropy of the irregularity. The sensitivity study of SIGMA narrows the parametric space to investigate when comparing the modeled results to the observations.

Beschreibung: The attenuation of VLF signals from lightning and ground-based VLF transmitters during transionospheric propagation has been the subject of recent interest, as discrepancies have been found between satellite data and model calculations. Previous modeling efforts, however, have not considered the self-absorption effect due to nonlinear heating and ionization in the lower ionosphere. A self-consistent model of ionospheric heating is presented here using a time-domain model of VLF wave propagation through the ionosphere. The model is able to estimate the attenuation of signals due to heating below ~100 km altitude. In this model, the ionospheric state is updated as the fields propagate, leading to changes in collision frequency and electron density, which in turn affect the wave propagation. We use this model for ground-based VLF transmitters at different frequencies, amplitudes, and latitudes (i.e., magnetic dip angle), and for lightning-generated sferics with different amplitudes, at different latitudes, and using a variety of ionospheric density profiles. We find that the inclusion of self-consistent heating causes a change in the transionospherically propagating wave amplitude that varies considerably with the source amplitude and other parameters. Typical values for the heating contribution to wave attenuation are 1-2 dB for VLF transmitters, but greater than 10 dB for large amplitude lightning discharges. An interesting effect is observed for VLF transmitters and low-amplitude lightning, where the signal is actually enhanced due to heating, rather than attenuated, in the direction propagating across the Earth's magnetic field.

Beschreibung: The present work has investigated the configuration of field-aligned currents (FACs) during a long period of radial interplanetary magnetic field (IMF) on 19 May 2002 by using high-resolution and precise vector magnetic field measurements of the CHAMP satellite. During the interest period IMF B y and B z are weakly positive and B x keeps pointing to the Earth for almost 10 hours. The geomagnetic indices Dst is about -40 nT and AE about 100 nT on average. The cross polar cap potential calculated from Assimilative Mapping of Ionospheric Electrodynamics and derived from DMSP observations have average values of 10-20 kV. Obvious hemispheric differences are shown in the configurations of FACs on the day and nightside. At the south pole FACs diminish in intensity to magnitudes of about 0.1  μA/m 2 , the plasma convection maintains two cell flow pattern, and the thermospheric density is quiet low. However, there are obvious activities in the northern cusp region. One pair of FACs with a downward leg toward the pole and upward leg on the equatorward side emerge in the northern cusp region, exhibiting opposite polarity to FACs typical for duskward IMF orientation. An obvious sunward plasma flow channel persists during the whole period. These ionospheric features might be manifestations of an efficient magnetic reconnection process occurring in the northern magnetospheric flanks at high latitude. The enhanced ionospheric current systems might deposit large amount of Joule heating into the thermosphere. The air densities in the cusp region get enhanced and subsequently propagate equatorward on the dayside. Although geomagnetic indices during the radial IMF indicate low level activity, the present study demonstrates that there are prevailing energy inputs from the magnetosphere to both the ionosphere and thermosphere in the northern polar cusp region.

Beschreibung: The interaction of planetary bodies with their surrounding magnetized plasma can often be described with the magneto-hydrodynamic (MHD) equations, which are commonly solved by numerical models. For these models it is necessary to define physically correctboundary conditions for the plasma mass and energy density, the plasma velocity and the magnetic field. Many planetary bodies have electrically non-conductive surfaces, which do not allow electric current to penetrate their surfaces. Magnetic boundary conditions, which consider that the associated radial electric current at the planetary surface is zero are difficult to implement because they include the curl of the magnetic field. Here we derive new boundary conditions by a decomposition of the magnetic field in poloidal and toroidal parts. We find that the toroidal part of the magnetic field needs to vanish at the surface of the insulator. For the spherical harmonics coefficients of the poloidal part we derive a Cauchy boundary condition, whichalso matches a possible intrinsic field by including its Gauss coefficients. Thus we can additionally include planetary dynamo fields as well as time-variable induction fields within electrically conductive subsurface layers. We implement the non-conducting boundary condition in the MHD simulation code ZEUS-MP using spherical geometry and provide a numerical implementation in Fortran 90 as auxiliary-material on the JGR website. We apply it to a model for Ganymede's plasma environment. Our model also includes a consistent set of boundary conditions for the other MHD variables density, velocity and energy. With this model we can describe Galileo spacecraft observations in and around Ganymede's mini-magnetosphere very well.

Beschreibung: In recent years we have benefitted greatly from the first in-orbit multi-wavelength images of Saturn's polar atmosphere from the Cassini spacecraft. Specifically, images obtained from the Cassini UltraViolet Imaging Spectrograph (UVIS) provide an excellent view of the planet's auroral emissions, which in turn give an account of the large-scale magnetosphere-ionosphere coupling and dynamics within the system. However, obtaining near-simultaneous views of the auroral regions with in situ measurements of magnetic field and plasma populations at high-latitudes is more difficult to routinely achieve. Here we present an unusual case, during Revolution 99 in January 2009, where UVIS observes the entire northern UV auroral oval during a two-hour interval while Cassini traverses the magnetic flux tubes connecting to the auroral regions near 21 LT, sampling the related magnetic field, particle, and radio and plasma wave signatures. The motion of the auroral oval evident from the UVIS images requires a careful interpretation of the associated latitudinally “oscillating” magnetic field and auroral field-aligned current signatures whereas previous interpretations have assumed a static current system. Concurrent observations of the auroral hiss (typically generated in regions of downward directed field-aligned current) support this revised interpretation of an oscillating current system. The nature of the motion of the auroral oval evident in the UVIS image sequence, and the simultaneous measured motion of the field-aligned currents (and related plasma boundary) in this interval, is shown to be related to the northern hemisphere magnetosphere oscillation phase. This is in agreement with previous observations of the auroral oval oscillatory motion.

Beschreibung: We have developed a technique by which to estimate the spatial distribution of plasmaspheric helium ions based on extreme ultraviolet (EUV) data obtained from the IMAGE satellite. The estimation is performed using a linear inversion method based on the Bayesian approach. The global imaging data from the IMAGE satellite enables us to estimate a global two-dimensional distribution of the helium ions in the plasmasphere. We applied this technique to a synthetic EUV image generated from a numerical model. This technique was confirmed to successfully reproduce the helium ion density that generated the synthetic EUV data. We also demonstrate how the proposed technique works for real data using two real EUV images.

Beschreibung: Twelve submarine mud volcanoes (MV) in the Kumano forearc basin within the Nankai Trough subduction zone were investigated for hydrocarbon origins and fluid dynamics. Gas hydrates diagnostic for methane concentrations exceeding solubilities were recovered from MVs 2, 4, 5, and 10. Molecular ratios (C 1 /C 2 〈250) and stable carbon isotopic compositions (δ 13 C-CH 4 〉−40‰ V-PDB) indicate that hydrate-bound hydrocarbons (HCs) at MVs 2, 4, and 10 are derived from thermal cracking of organic matter. Considering thermal gradients at the nearby IODP Sites C0009 and C0002, the likely formation depth of such HCs ranges between 2,300 and 4,300 m below seafloor (mbsf). With respect to basin sediment thickness and the minimum distance to the top of the plate boundary thrust we propose that the majority of HCs fueling the MVs is derived from sediments of the Cretaceous to Tertiary Shimanto belt below Plio-/Pleistocene to recent basin sediments. With respect to sizes and appearances hydrates are suggested to be relicts of higher MV activity in the past, although the sporadic presence of vesicomyid clams at MV 2 showed that fluid migration is sufficient to nourish chemosynthesis-based organisms in places. Distributions of dissolved methane at MVs 3, 4, 5, and 8 pointed at fluid supply through one or few MV conduits and effective methane oxidation in the immediate subsurface. The aged nature of the hydrates suggests that the major portion of methane immediately below the top of the methane-containing sediment interval is fueled by current hydrate dissolution rather than active migration from greater depth.

Beschreibung: 40 Ar/ 39 Ar geochronology depends critically on well calibrated standards, often traceable to first-principles K-Ar age calibrations using bulb-tracer systems. Tracer systems also provide precise standards for noble-gas studies and interlaboratory calibration. The exponential expression long-used for calculating isotope tracer concentrations in K-Ar age dating and calibration of 40 Ar/ 39 Ar age standards may provide a close approximation of those values, but is not correct. Appropriate equations are derived that accurately describe the depletion of tracer reservoirs and concentrations of sequential tracers. The true form of the expression is a power law, not exponential, and a similar expression was presented by Miiller (2006, J. Res. Natl. Inst. Stand. Technol., 111 (5), 335–360). Evaluation of the expressions demonstrates that systematic error introduced through use of the exponential approximation may be substantial where reservoir volumes are small and resulting depletion constants are large. Traditional use of large reservoir to tracer volumes and the resulting small depletion constants have kept errors well less than experimental uncertainties in most previous K-Ar and calibration studies. Use of the proper expression, however, permits use of volumes appropriate to the problems addressed.

Beschreibung: Natural hydrate-bearing sediment (HBS) predominantly exists in non-cementing habit, and its limited availability for use in laboratory studies demands a time-effective and repeatable laboratory process for forming representative samples with natural accumulation habit. This study reports on a three-step laboratory process for forming non-cementing methane hydrate in sandy sediments: (1) initial HBS formation under excess-gas conditions; (2) slow saline water (5wt % CaCl 2 ) injection under strictly controlled PT conditions; and (3) a temperature warming/cooling cycle. Changes in compressional wave velocity ( V p ) of sediment, as well as pressure-temperature (P-T) condition, were monitored throughout the tests. The evolution of V p , in good agreement with rock physics model calculations, suggested that the transition from cementing hydrate into non-cementing hydrate occurs during saline injection as well as temperature warming/cooling cycle. The proposed process appeared to be an efficient and consistent substitute for the existing methods, to form non-cementing hydrate habit in sandy sediments.

Beschreibung: Great earthquakes anticipated on the Cascadia subduction fault can potentially rupture beyond the geodetically and thermally inferred locked zone to the depths of episodic tremor and slip (ETS) or to the even deeper forearc mantle corner (FMC). To evaluate these extreme rupture limits, we map the FMC from southern Vancouver Island to central Oregon by combining published seismic velocity structures with a model of the Juan de Fuca plate. These data indicate that the FMC is somewhat shallower beneath Vancouver Island (36–38 km) and Oregon (35–40 km) and deeper beneath Washington (41–43 km). The updip edge of tremor follows the same general pattern, overlying a slightly shallower Juan de Fuca plate beneath Vancouver Island and Oregon (˜30 km) and a deeper plate beneath Washington (˜35 km). Similar to the Nankai subduction zone, the best constrained FMC depths correlate with the center of the tremor band suggesting that ETS is controlled by conditions near the FMC rather than directly by temperature or pressure. Unlike Nankai, a gap as wide as 70 km exists between the downdip limit of the inferred locked zone and the FMC. This gap also encompasses a ˜50 km wide gap between the inferred locked zones and the updip limit of tremor. The separation of these features offers a natural laboratory for determining the key controls on downdip rupture limits.

Beschreibung: We use cosmic radiation records (neutron monitor and the cosmogenic radionuclides, 10 Be and 14 C) as a proxy to compare the solar activity during the extended solar minimum 2006-9, with that during the Grand Solar Minima and Maxima that occurred between 1391-2010. The inferred cosmic ray intensities during the Spoerer, Maunder, and Dalton Grand Minima were significantly greater than those during 2006-9. The onset phases of the three Grand Minima extended over between two and five Schwabe (sunspot) cycles, the cosmic ray intensity at the Schwabe minima increasing from a value approximating that of 2006-9, to substantially higher values later in the Grand Minimum. The minimum estimated strengths of the heliospheric magnetic field near Earth during the Grand Minima were 2.4nT (Spoerer); 〈2.0nT(Maunder) and 2.6 nT (Dalton), compared to 3.9nT in 2009. We conclude that the periods of highest solar activity during the Maunder Minimum approximated those near the sunspot minima between 1954 and 1996. The average ratio of the maximum to minimum estimated HMF in the six Schwabe cycles in the Maunder Minimum is 1.54 (range 1.30-1.85) compared to 1.52 (1.31-1.63) for the modern epoch suggesting similar operation of the solar dynamo in both intervals. The onset phase of the Maunder Minimum extending over five Schwabe cycles, and the large increase in cosmic ray flux (and decrease in estimated heliospheric magnetic field), leads us to speculate that the magnetohydrodynamic amplification in the solar dynamo exhibits a relaxation time well in excess of the 11 year period of the Schwabe cycle.

Beschreibung: The temperature profile obtained from the space borne instrument ‘Sounding of Atmosphere by Broadband Emission Radiometry’ (SABER) instrument onboard ‘Thermosphere Ionosphere Mesosphere Energetics and Dynamics’ (TIMED) shows a triple layered mesospheric inversion event on the night of 23 September 2011, when there is an overpass near the low-latitude sites Gadanki (13.5°N, 79.2°E) and Tirunelveli (8.7°N, 77.8°E). The three mesospheric inversion layers (MIL) are formed in the height region around ~70 km (lower), ~80 km (middle) and ~90 km (upper) with amplitudes ~11 K, ~44 K, and ~109 K and thickness of 3.4 km, 4.9 km and 6.6 km respectively. The formation of the lower and middle MILs can only be observed in the Rayleigh lidar temperature profiles over Gadanki due to upper height limitation of the system. Nearly all the dominant causative mechanisms are examined for the occurrence of the MIL event. The lower MIL at ~70 km is inferred to be due to planetary wave dissipation, as there is a sudden decrease of planetary wave activity above 70 km. Further, it is demonstrated that the middle MIL at ~80 km is due to the turbulence generated by gravity wave breaking which is in turn due to gravity wave–semi-diurnal tidal interaction, though the height of the middle MIL descends at the rate of ~1 km/h, which is nearly the vertical phase speed of diurnal tide, whereas the upper MIL at above 90 km is due to the large chemical heating rate (~45 K/day) generated by the dominant exothermic reaction O + O + M → O 2  + M.

Beschreibung: In this paper, we analyze radar observations of E x B drift and plasma irregularities, ionosonde observations of E- and F-layer parameters including spread F, and magnetic field observations made from Indian low latitudes linked with the 2009 sudden stratospheric warming (SSW) event. E x B drift variations presented here are the first of their kind from the Indian sector as far as the effect of SSW is concerned. Difference of magnetic fields observed from the equator and low latitude (∆H) and E x B drift show linear relation and both show remarkably large positive values in the morning and negative values in the afternoon exhibiting semidiurnal behavior. Remarkable changing patterns in the critical frequency of F 2 layer (f o F 2 ) and F 3 layer (f o F 3 ) were observed after the occurrence of SSW. Large variations with quasi-16-day periodicity were observed in ∆H, f o F 2 and f o F 3. Both semidiurnal and quasi-16-day wave modulation observed after the 2009 SSW event are consistent with those reported earlier. We also noted quasi-6 day variations in ∆H and f o F 2 soon after the SSW commencement, not much reported before. During the counter-electrojet events linked with the SSW event, while equatorial E s (E sq ) disappeared as expected, there were no blanketing E s (E sb ), a finding not reported and discussed earlier. E sb was also not formed at the off-equatorial location, indicating the absence of required vertical wind shear, but E region plasma irregularities were observed by the ionosonde and radar with a close relationship between the two. Weak F region irregularities were observed in the post-midnight hours and case studies suggest the possible role of SSW related background electric field in the manifestation of post-midnight F region irregularities.

Beschreibung: The Defense Meteorological Satellite Program (DMSP) launches and maintains a network of satellites to monitor the meteorological, oceanographic, and solar-terrestrial physics environments. In the past decade, geomagnetic field modelers have focused much attention on magnetic measurements from missions such as CHAMP, Ørsted and SAC-C. With the completion of the CHAMP mission in 2010, there has been a multi-year gap in satellite-based vector magnetic field measurements available for main field modeling. In this study, we calibrate the Special Sensor Magnetometer (SSM) instrument onboard DMSP to create a dataset suitable for main field modeling. These vector field measurements are calibrated to compute instrument timing shifts, scale factors, offsets, and non-orthogonality angles of the fluxgate magnetometer cores. Euler angles are then computed to determine the orientation of the vector magnetometer with respect to a local coordinate system. We fit a degree 15 main field model to the dataset and compare with the World Magnetic Model (WMM) and Ørsted scalar measurements. We call this model DMSP-MAG-1 and its coefficients and software are available for download at http://geomag.org/models/dmsp.html . Our results indicate that the DMSP dataset will be a valuable source for main field modeling for the years between CHAMP and the recently launched Swarm mission.

Beschreibung: During eruptions onto low slopes, basaltic Pahoehoe lava can form thin lobes that progressively coalesce and inflate to many times their original thickness, due to a steady injection of magma beneath brittle and viscoelastic layers of cooled lava that develop sufficient strength to retain the flow. Inflated lava flows forming tumuli and pressure ridges have been reported in different kinds of environments, such as at contemporary subaerial Hawaiian-type volcanoes in Hawaii, La Réunion and Iceland, in continental environments (states of Oregon, Idaho, Washington), and in the deep sea at Juan de Fuca Ridge, the Galapagos spreading center, and at the East Pacific Rise (this study). These lava have all undergone inflation processes, yet they display highly contrasting morphologies that correlate with their depositional environment, the most striking difference being the presence of water. Lava that have inflated in subaerial environments display inflation structures with morphologies that significantly differ from subaqueous lava emplaced in the deep sea, lakes, and rivers. Their height is 2 to 3 times smaller and their length being 10 to 15 times shorter. Based on heat diffusion equation, we demonstrate that more efficient cooling of a lava flow in water leads to the rapid development of thicker (by 25%) cooled layer at the flow surface, which has greater yield strength to counteract its internal hydrostatic pressure than in subaerial environments, thus limiting lava breakouts to form new lobes, hence promoting inflation. Buoyancy also increases the ability of a lava to inflate by 60%. Together, these differences can account for the observed variations in the thickness and extent of subaerial and subaqueous inflated lava flows.

Beschreibung: [1]  Jupiter's large scale size and rapid planetary rotation period combine to produce the strong centrifugal force responsible for many unique properties of its magnetosphere. It was previously proposed that this centrifugal force and non-adiabatic field line stretching could cause the observed dawn-dusk asymmetry of Jupiter's plasma sheet, which is thickest near dusk. As flux tubes rotate and stretch between noon and dusk, particles bouncing along the field gain parallel energy and create pressure anisotropy. Because bounce times can be long compared with the outward expansion timescale, particles may respond non-adiabatically, and the resulting pressure anisotropy can drive the plasma sheet to instability. We used a large-scale kinetic (LSK) simulation to follow a collection of rotating particles as they move in a time-varying, rotating magnetic field designed to represent flux tube expansion in Jupiter's magnetosphere. The analysis quantifies the response of trapped particles by characterizing the pressure anisotropy and energy changes. We compare results of non-adiabatic and adiabatic outward expansion, and find that the non-adiabatic case leads to a large pitch angle anisotropy and higher total energy than for adiabatic expansion. Although the calculation was not handled fully self-consistently, the results support the proposition that plasma pressure changes lead to changes in the magnetic field structure with local time. Our findings are consistent with the idea that non-adiabatic effects in Jupiter's magnetosphere contribute to field dipolarization and the observed plasma sheet thickening between noon and dusk.

Beschreibung: [1]  The IMAGE network magnetic measurements during the period 1995-2009 (covering the whole of solar cycle 23) are used to characterize the annual variations in the westward electrojet. The results suggest that the annual variations in different local time sector are quite different due to the different sources. In the MLT-sector 2200-0100, the annual variations with maxima in winter suggest they are caused by the combined effects of the convective electric field and the conductivity associated with particle precipitation. Furthermore, the conductivity seems to play a more important role in the MLT-sector ~2200-2320, while the convective electric field appears to be more important in the MLT-sector ~2320-0100. In the MLT-sector 0300-0600, the annual variations with maxima in summer suggest they are caused by solar EUV conductivity effect and the equinoctial effect, which work in two fundamentally different ways. The solar EUV conductivity effect works by increasing ionospheric conductivity and enhancing the westward electrojet in summer, while the equinoctial effect works by decreasing solar wind-magnetosphere coupling efficiency and weakening the westward electrojet in winter. In the MLT-sector 0100-0300, the annual variations are relatively weak, and can be attributed to the combined effects of annual variations caused by all the previously mentioned effects. In addition, we find that a significant annual variation in substorm occurrence rate, mainly occurring in the premidnight region, is quite similar to that in the westward electrojet. We suggest that elevated solar wind driving during the winter months contributes to higher substorm occurrence in winter in the Northern Hemisphere.

Beschreibung: I analyze the seismic tomography around the Japanese Islands to elucidate the present slab morphology. A simple conversion from the high speed to the temperature anomaly is done based on the recent study of the conversion coefficient from the velocity to the thermal anomaly and on the study of the relation between the seismicity and the temperature in the slab. I find recognizable amount of cold temperature anomalies under the subducting slab. This suggests a rather continuous slab throughout the mantle consistent with the results of many recent numerical models. However, there still is a gap in the slab or diluted slab under the transition zone. To understand this, I construct a simple half-kinematic model of subduction zone in which a Byerlee's type yield stress and the depth independent yield stress are assumed. Taking into account the geologic history around the Japanese Islands, I find that the slab morphology similar to the tomographic image is obtained when the yield stress is O(100 MPa). The present study shows that the detailed studies of seismic tomography and the tectonic history of the surveyed area can provide the constraints on the slab dynamics.

Beschreibung: No abstract is available for this article.

Beschreibung: Chorus waves play an important role in energetic electron dynamics in the inner magnetosphere. In this work, we present a new hybrid code, DAWN, to simulate the generation of chorus waves. The DAWN code is unique in that it models cold electrons using linearized fluid equations and hot electrons using particle-in-cell techniques. The simplified fluid equations can be solved with robust and simple algorithms. We demonstrate that discrete chorus elements can be generated using the code. Waveforms of the generated element show amplitude modulation or “subpackets”, and the frequency sweep rate of the generated element is found to be consistent with that of observed chorus waves. Using the DAWN code, we then investigate the variation of wave intensity with respect to linear growth rates on the equatorial plane. Previous observations showed that the change in linear growth rates of whistler waves modulated by external processes such as density modulations is usually small , while the variation of the wave intensity is large . Using a chosen set of background plasma parameters, we demonstrate that a small change in linear growth rates can lead to significant variation of wave intensity only in the transition from the broadband whistler wave generation regime to the chorus wave generation regime. Our results demonstrate the importance of including nonlinear dynamics of chorus generation in understanding the whistler wave intensity modulation process in the inner magnetosphere.

Beschreibung: We report the observation of echo extreme horizontal drift speed (EEHS, ≥ 300 m s -1 ) during polar mesospheric (80-90 km) summer echoes (PMSEs) by the VHF (52-MHz) radar at Esrange, Sweden, in years of 2006 and 2008. The EEHS occur in PMSEs as correlated with high-speed solar wind streams (HSSs), observed at least once in 12-17 % of all hours of observation for the two summers. The EEHS rate peaks occur either during high solar wind speed in the early part of the PMSE season or during the arrival of interplanetary corotating interaction regions (CIRs) followed by peaks in PMSE occurrence rate after 1-4 days, in the latter part of the 2006 summer. The cause of EEHS rate peaks is likely under the competition between the interval of the CIR and HSS passage over the magnetosphere. A candidate process in producing EEHS is suggested to be localized strong electric field, which is caused by solar wind energy transfer from the interaction of CIR and HSS with the magnetosphere in a sequential manner. We suggest that EEHS are created by strong electric field, estimated as 〉 10-30 V m -1 at 85 km altitude, exceeding the mesospheric breakdown threshold field.

Beschreibung: As the supercontinent Rodinia was assembling ca. 1.1 billion years ago, there was extensive magmatism on at least five Proterozoic continents including the development of the North American Midcontinent Rift. New paleomagnetic data from 84 lava flows of the Osler Volcanic Group of the Midcontinent Rift reveal that there was a significant and progressive decrease in inclination between the initiation of extrusive volcanism in the region (ca. 1110 Ma) and ca. 1105 ± 2 Ma (the “early stage” of rift development). Paleomagnetic poles can be calculated for the lower portion of the reversed Osler Volcanic Group (40.9°N, 218.6 °E, A 95 =4.8°, N=30) and the upper portion of the reversed Osler Volcanic Group (42.5°N, 201.6 °E, A 95 =3.7°, N=59; this pole can be assigned the age of ca. 1105 ± 2 Ma). This result is a positive test of the hypothesis that there was significant plate motion during the early stage of rift development. In addition to being a time of widespread volcanism on Laurentia and other continents, this interval of the late Mesoproterozoic was characterized by rapid paleogeographic change.

Beschreibung: The relationship between magnetic hysteresis parameters and the degree of oxidation of ultrafine magnetite particles is examined by both experimental measurements (distributed particle assemblage with median grain size of ~80 nm and standard deviation 0.43) and micromagnetic simulations (single particles from 40 nm to 140 nm). Experimental results show that both coercivity (B c ) and the ratio of saturation remanence to saturation magnetization (M rs /M s ) increase slowly, as the oxidation parameter z increases from 0 to ~0.9. Thereafter both parameters decrease sharply as magnetite becomes completely oxidized to maghemite. Numerical simulations of hysteresis loop and microstructure using a micromagnetic model with a core-shell geometry (a stoichiometric core surrounded by an oxidized shell) show three categories of behavior for magnetic grains during oxidation. Firstly, the coercivity of SD particles decreases as oxidation proceeds, but their remanence magnetization remains in a uniform state. Secondly, for PSD sized particles near the critical SD boundary (80 nm to 100 nm), the initial vortex domain structure changes to a SD as oxidation occurs and returns to a vortex state upon complete maghemitization, resulting in an initial rise and then fall of B c and M rs . Finally, larger PSD grains remain a vortex state throughout the maghemitization, with less variations of B c and M rs . The predicted magnetic properties exhibit good agreement with experimental observations and suggest that the domain arrangement is likely to be dominated by a core-shell structure with strong exchange coupling at their interface. Overall, the partially oxidized magnetite in SD-PSD range can reliably record palaeomagnetic signals.

Beschreibung: We have developed a technique for estimating the temporal evolution of the plasmaspheric helium ion density based on a sequence of extreme ultraviolet (EUV) data obtained from the IMAGE satellite. In the proposed technique, the estimation is obtained by incorporating EUV images from IMAGE into a two-dimensional fluid model of the plasmasphere using a data assimilation approach based on the ensemble transform Kalman filter. Since the motion and the spatial structure of the helium plasmasphere is strongly controlled by the electric field in the inner magnetosphere, the electric field around the plasmapause can also be estimated using the ensemble transform Kalman filter. We performed an experiment using synthetic images that were generated from the same numerical model under a certain condition. It was confirmed that the condition that generated the synthetic images was successfully reproduced. We also present some results obtained using real EUV imaging data. Finally, we discuss the possibility of estimating the density profile along a magnetic field line. Since each EUV image was taken from a different direction due to the motion of the IMAGE satellite, we could obtain the information on the density profile along a field line by combining multiple images.

Beschreibung: Small-scale structure of the plasma convection and electron content within the subauroral polarization stream (SAPS) is investigated. We present ionospheric observations during the main phase of the geomagnetic storm on March 17, 2013, during which a sequence of intense, highly localized, and fast-moving electric field (EF) structures within SAPS was observed by the SuperDARN Christmas Valley West radar (CVW). The CVW EF measurements at 60-s resolution are analyzed in context of coincident GPS measurements of the total electron content (TEC) at 30-s resolution. The strong and narrow feature of the subauroral ion drift (SAID) was observed poleward of the TEC trough, with a TEC enhancement (peak) seen in the SAPS (SAID) region. The SAPS wave activity commenced ~2 hours (15 min) after first appearance of SAPS (SAID). The SAPS structures appeared near the poleward edge of the trough, propagated westward, and merged with SAID near TEC peak. The propagation velocity was comparable with convection velocity within each EF structure. The SAPS TEC exhibited a general decrease towards the end of the period. On a smaller time scale, TEC exhibited a small but appreciable decrease within EF structures. The wavelet spectra of EF and TEC showed similar variations, with wave period of ~5-min period near onset and increasing to 8–10 min towards the end of the period with significant wave activity. A scenario is discussed, in which the SAPS wave activity may modify the ionospheric conductance and TEC at small scales, with large-scale magnetosphere-ionosphere feedback acting to continuously deplete TEC where/when such activity does not occur.

Beschreibung: The intense inner radiation belt at Jupiter (〉50 MeV at 1.5 R J ) is generally accepted to be created by radial diffusion of electrons from further away from the planet. However, this requires a source with energies that exceed1 MeV outside the orbit of the moon Io at 5.9 R J , which has never been explained satisfactorily. Here we test the hypothesis that this source population could be formed from a very soft energy spectrum, by particle injection processes and resonant electron acceleration via whistler mode chorus waves. We use the BAS Radiation Belt Model to calculate the change in the electron flux between 6.5 and 15 R J ; these are the first simulations at Jupiter combining wave particle interactions and radial diffusion. The resulting electron flux at 100 keV and 1 MeV lies very close to the Galileo Interim Radiation Electron model spectrum after 1 and 10 days respectively. The primary driver for the increase in the flux is cyclotron resonant acceleration by chorus waves. A peak in phase space density forms such that inside L ≈ 9 radial diffusion transports electrons towards Jupiter, but outside L ≈ 9 radial diffusion acts away from the planet. The results are insensitiveto the softness of the initial energy spectrum but do depend on the value of the flux at the minimum energy boundary. We conclude by suggesting that the source population for the inner radiation belt at Jupiter could indeed be formed by wave-particle interactions.

Beschreibung: We report 3 He/ 4 He for 150 mid-ocean ridge basalt (MORB) glasses from the Southeast Indian Ridge (SEIR). Between 81°-101°E 3 He/ 4 He varies from 7.5 to 10.2 R A , encompassing more than half the MORB range away from ocean island hotspots. Abrupt transitions are present and in one case the full range occurs over ~10 km. Melting of lithologically heterogeneous mantle containing a few percent garnet pyroxenite or eclogite leads to lower 3 He/ 4 He, while 3 He/ 4 He above ~9 R A likely indicates melting of pyroxenite-free or eclogite-free mantle. Patterns in the length scales of variability represent a description of helium isotopic texture. We utilize four complementary methods of spectral analysis to evaluate this texture, including Periodogram, Redfit, Multi-Taper Method and Continuous Wavelet Transform. Long-wavelength lobes with prominent power at 1000 km and 500 km are present in all treatments, similar to hotspot-type spectra in Atlantic periodograms. The densely sampled region of the SEIR considered separately shows significant power at ~100 km and ~30-40 km, the latter scale resembling heterogeneity in the bimodal distribution of Hf and Pb isotopes in the same sample suite. Wavelet transform coherence reveals that 3 He/ 4 He varies in-phase with axial depth along the SEIR at ~1000 km length scale, suggesting a coupling between melt production, 3 He/ 4 He and regional variations in mantle temperature. Collectively, our results show that the length scales of MORB 3 He/ 4 He variability are dominantly controlled by folding and stretching of heterogeneities during regional (~1000 km) and mesoscale (~100 km) mantle flow, and by sampling during the partial melting process (~30 km).

Beschreibung: Enceladus’ southern plume is one of the major discoveries of the Cassini mission. The water neutrals and ice crystals ejected by the cryovolcanic activity populate Saturn's E-ring and the neutral torus, and they interact with the plasma environment of Saturn's magnetosphere. The plasma neutrality inside Enceladus’ plume has been shown by the Langmuir probe measurement to be modified by the presence of the dust particles. We present an independent method of determining the electron density inside the plume. Sometimes after dust impacts, plasma oscillations (ringing) were detected by the Cassini Radio and Plasma Wave Science (RPWS) instrument. The frequencies of these oscillations have been shown to be consistent with the local plasma frequency, thus providing a measurement of the local electron density.

Beschreibung: The study of complexity in two aspects of the magnetic activity in the Sun-Earth system is presented. We compare the temporal evolution of the magnetic fluctuations in the Earth's magnetosphere, and the spatial distribution of the magnetic field in the solar photosphere, by calculating fractal dimensions from the data. It is found that the fractal dimension of the D st data decreases during magnetic storm states, and is well correlated with other indexes of solar activity, such as the solar flare and coronal indexes. This correlation holds for individual storms, full year data, and the complete 23rd solar cycle. The fractal dimension from solar magnetogram data also correlates well with both the D st index and solar flare index, although the correlation is much more clear at the larger temporal scale of the 23rd solar cycle, showing a clear increase around solar maximum.

Beschreibung: The variations of the 3-D coronal magnetic fields associated with the X3.4-class flare of active region 10930 are studied in this paper. The coronal magnetic field data are reconstructed from the photospheric vector magnetograms obtained by the Hinode satellite and using the nonlinear force-free field extrapolation method developed in our previous work (He et al., 2011). The 3-D force-free factor α , 3-D current density, and 3-D magnetic energy density are employed to analyze the coronal data. The distributions of α and current density reveal a prominent magnetic connectivity with strong negative α values and strong current density before the flare. This magnetic connectivity extends along the main polarity inversionline and, is found to be totally broken after the flare. The distribution variation of magnetic energy density reveals the re-distribution of magnetic energy before and after the flare. In the lower space of the modeling volume the increase of magnetic energy dominates, and in the higher space the decrease of energy dominates. The comparison with the flare onset imaging observation exhibits that the breaking site of the magnetic connectivity and site with the highest values of energy density increase coincide with the location of flare initial eruption. We conclude that a cramped positive α region appearing in the photosphere causes the breaking of the magnetic connectivity. A scenario for flare initial eruption is proposed in which the Lorentz force acting on the isolated electric current at the magnetic connectivity breaking site lifts the associated plasmas and causes the initial ejection.

Beschreibung: Using a recent magnetohydrodynamic simulation of magnetotail dynamics we further investigate the build-up and evolution of the substorm current wedge (SCW), resulting from flow bursts generated by near-tail reconnection. Each flow burst generates an individual current wedge, which includes the reduction of cross-tail current and the diversion to region 1 (R1) type field-aligned currents (earthward on the dawn and tailward on the dusk side), connecting the tail with the ionosphere. Multiple flow bursts generate initially multiple SCW patterns, which at later times combine to a wider single SCW pattern. The standard SCW model is modified by the addition of several current loops, related to particular magnetic field changes: the increase of B z in a local equatorial region (“dipolarizationÓ), the decrease of | B x | away from the equator (“current disruptionÓ), and increases in | B y | resulting from azimuthally deflected flows. The associated loop currents are found to be of similar magnitude, 0.1-0.3 MA. The combined effect requires the addition of region 2 (R2) type currents closing in the near tail through dawnward currents but also connecting radially with the R1 currents. The current closure at the inner boundary, taken as a crude proxy of an idealized ionosphere, demonstrates westward currents as postulated in the original SCW picture as well as North-South currents connecting R1 and R2 type currents, which were larger than the westward currents by a factor of almost 2. However, this result should be applied with caution to the ionosphere because of our neglect of finite resistance and Hall effects.

Beschreibung: We report ARTEMIS dual-probe observations of two events in the terrestrial magnetotail lobe, both characterized by upward-moving heavy ions of lunar origin at one of the probes that is magnetically connected with the dayside lunar surface. By treating magnetic measurements at the other probe as the unperturbed lobe fields, we obtain background-subtracted magnetic perturbations (most significantly in B z ) when the first probe moved in the dawn-dusk direction across flux tubes magnetically connected to the Moon. These magnetic perturbations indicate the presence of field-aligned current above the lunar surface. By examining possible carriers of field-aligned current, we find that lunar heavy ions and accompanying electrons both contribute considerably to the current. Observations of the field-aligned current also suggest that the charging process at the dayside lunar surface and the associated lobe plasma environment, which have traditionally been viewed as a one-dimensional current balance problem, are actually more complicated. These observations give the first insights into how heavy ions affect the lunar dayside environment in terms of multi-species plasma dynamics.

Beschreibung: The solution of electric fields and currents in a height–resolved ionosphere is traditionally solved as an elliptic equation with Dirichlet or Neumann boundary condition in which the magnetosphere is represented as an unresponsive (prescribed) voltage generator or current source. In this paper we derive an alternative boundary condition based upon Alfvén waves in which only the Alfvén wave from the magnetosphere that is incident upon the ionosphere ( E xi ) is prescribed. For a uniform magnetosphere the new boundary condition reduces to ∂ φ /∂ z  =  (∂ 2 φ /∂ x 2  + 2∂ E xi /∂ x )/( μ 0 V A σ || ) and is evaluated at the magnetosphere–ionosphere interface. The resulting solution is interpreted as a responsive magnetosphere and establishes a key stage in the full, self-consistent and non-linear coupling of the magnetosphere and ionosphere.

Beschreibung: Reconstructing convective flow in the Earth's mantle is a crucial issue for a diversity of disciplines, from seismology to sedimentology. The common and fundamental limitation of these reconstructions based on geodynamic modelling is the unknown initial conditions. Because of the chaotic nature of convection in the Earth's mantle, errors in initial conditions grow exponentially with time and limit forecasting and hindcasting abilities. In this work we estimate for the first time the limit of predictability of Earth's mantle convection. Following the twin experiment method, we compute the Lyapunov time (i.e. e-folding time) for state-of-the art 3D spherical convection models, varying rheology and Rayleigh number. Our most Earth-like and optimistic solution gives a Lyapunov time of 136±13 My. Rough estimates of the uncertainties in best guessed initial conditions are around 5%, leading to a limit of predictability for mantle convection of 95 My. Our results suggest that error growth could produce unrealistic convective structures over timescales shorter than that of Pangea dispersal.

Beschreibung: . Sub-ionospheric radio-wave data from an AARDDVARK receiver located in Churchill, Canada, is analysed to determine the characteristics of electron precipitation into the atmosphere over the range 3 〈  L  〈 7. The study advances previous work by combining signals from two US transmitters from 20 July – 20 August 2010, allowing error estimates of derived electron precipitation fluxes to be calculated, including the application of time-varying electron energy spectral gradients. Electron precipitation observations from the NOAA POES satellites, and a ground-based riometer provide inter-comparison, and context, for the AARDDVARK measurements. AARDDVARK radiowave propagation data showed responses suggesting energetic electron precipitation from the outer radiation belt starting 27 July 2010, and lasting ~20 days. The uncertainly in 〉30 keV precipitation flux determined by the AARDDVARK technique was found to be ±10%. Peak 〉30 keV precipitation fluxes of AARDDVARK-derived precipitation flux during the main- and recovery-phase of the largest geomagnetic storm, that started on 04 August 2010, were 〉10 5 el. cm -2  s -1  sr -1 . The largest fluxes observed by AARDDVARK occurred on the dayside, and were delayed by several days from the start of the geomagnetic disturbance. During the main phase of the disturbances nightside fluxes were dominant. Significant differences in flux estimates between POES, AARDDVARK and the riometer were found after the main phase of the largest disturbance, with evidence provided to suggest that 〉700 keV electron precipitation was occurring. Currently the presence of such relativistic electron precipitation introduces some uncertainty in the analysis of AARDDVARK data, given the assumption of a power-law electron precipitation spectrum.