ALBERT

Description: [1]  Neutron monitors have recorded the flux of high energy Galactic cosmic rays for more than half a century. During the recent prolonged, deep minimum in solar activity many sources indicate that modulated Galactic cosmic rays have attained new space-age highs. However reported neutron monitor rates are ambiguous; some record new highs while others do not. This work examines the record of 15 long-running neutron monitors to evaluate cosmic ray fluxes during the recent extraordinary solar minimum in a long-term context. We show that ground-level neutron rates did reach a historic high during the recent solar minimum, and we present a new analysis of the cosmic ray energy spectrum in the year 2009 versus year 1987. To do this we define a reference as the average of eight high-latitude neutron monitors, four in the northern hemisphere (Apatity, Inuvik, Oulu, Thule) and four in the southern hemisphere (Kerguelen, McMurdo, Sanae, Terre Adelie). Most stations display changes in sensitivity, which we characterize by a simple linear trend. After correcting for the change in sensitivity, a consistent picture emerges. With our correction all stations considered display new highs at the recent solar minimum, approximately 3% above the previous record high. These increases are shown to be consistent with spacecraft observations.

Description: [1]  We present the occurrence frequency of downgoing auroral electron beams in magnetic local time and invariant latitude, and the dependence on solar cycle, as indicated by F10.7, on whether the ionospheric footpoint of the satellite is illuminated or dark, and on the energy flux carried by the electrons. As previously reported, we find that the occurrence of electron beams peaks in the pre-midnight local time sector and that solar illumination at the footpoint (solar zenith angle) reduces both the occurrence and energy of the electron beams. The effect of solar maximum conditions (indicated by F10.7) is almost as large as the effect of the solar zenith angle. The characteristic energy of the electron beams is dependent on the energy flux carried, in addition to both solar zenith angle and F10.7. The beam energy (and therefore the parallel potential drop) is ~1.6 times higher for during solar minimum than during solar maximum for both dark and illuminated footpoints. The beam energy during dark solar minimum conditions is a factor of ~3 more than during sunlit minimum conditions. The ‘area’ covered by intense aurora is also reduced during solar maximum, for both sunlit and dark conditions. There is no evidence that the statistical results are due to the fact that acceleration via parallel electric fields moves to lower latitudes during solar maximum.

Description: [1]  In this paper we describe and quantify the energy transfer, flow and distribution. Our high-resolution data-set covers 13 years of OMNI, SuperMAG and Kyoto data. We employ what we consider to be the best estimates for energy sinks, and relate these to SuperMAG indices for better coverage and spatial resolution. For the energy input we have used the method of dimensional analysis [ Vasyliunas et al ., 1982] that is presented in unit power and makes it appropriate for energy analysis. A cross-correlation analysis parameterizes the magnetospheric response on the solar wind parameters during a wide range of conditions, ranging from substorms and storms up to a decade. The determined functional form is then evaluated and scaled using superposed epoch analysis of geomagnetic storms, revealing that the effective area of interaction can not be considered static. Instead we present a dynamic area which depends to the first order on the cube of the IMF B z component. Also, we find that for longer time periods this area must be increased compared to the area used for geomagnetic storms. We argue that some of the terms in the energy coupling function are contributory to describing magnetosheath conditions, and discuss how our coupling function can be related to Maxwell stress components. Also, we quantify the relative importance of the different energy sinks during substorms, geomagnetic storms and long time series, and present the coupling efficiency of the solar wind. Our energy coupling functions is compared with the ɛ parameter [ Akasofu and Perreault , 1978] and performs better for almost any event.

Description: [1]  Plasma sheet flow bursts have been suggested to correspond to different types of auroral activity, such as poleward boundary intensifications (PBIs), ensuing auroral streamers, and substorms. The flow-aurora association leads to the important question of identifying the magnetotail source region for the flow bursts and how this region depends on magnetic activity. The present study uses the ARTEMIS spacecraft coordinated with conjugate ground-based auroral imager observations to identify flow bursts beyond 45 R E downtail and corresponding auroral forms. We find that quiet-time flows are directed dominantly earthward with a one-to-one correspondence with PBIs. Flow bursts during the substorm recovery phase, and during steady magnetospheric convection (SMC) periods are also directed earthward, and these flows are associated with a series of PBIs/streamers lasting for tens of minutes with similar durations to that of the series of earthward flows. Pre-substorm onset flows are also earthward and associated with PBIs/streamers. The earthward flows during those magnetic conditions suggest that the flow bursts, which lead to PBIs and streamers, originate from further downtail of ARTEMIS, possibly from the distant tail neutral line (DNL) or tailward-retreated near-Earth neutral line (NENL) rather than from the nominal NENL location in the mid-tail. We find that tailward flows are limited primarily to the substorm expansion phase. They continue throughout the period of auroral poleward expansion, indicating that the expansion-phase flows originate from the NENL and that NENL activity is closely related to the auroral expansion of the substorm expansion phase.

Description: [1]  The solar eclipse on 15 January 2010 traversed Asia and completed its travel on the Shandong Peninsula in China at sunset. Two vertical-incidence ionosondes at Wuhan and Beijing and the oblique-incidence ionosonde network in North China were implemented to record the ionospheric response to the solar eclipse. Following the initial electron density decrease caused by the eclipse, the ionosphere was characterized by a strong pre-midnight enhancement, and a subsequent ionospheric decay, and a ~10 hour later post-midnight enhancement. Neither geomagnetic disturbance occurred during the eclipse day, nor did obvious nighttime peak appear for the ten-day mean of the F2-layer critical frequency ( fo F2). The electron density profilogram of the Beijing ionosonde indicates that the two enhancements were the result of the plasma flux downward from the top ionosphere, possibly due to the steep decrease of the ionospheric electron density and plasma temperature during the solar eclipse. The two-dimensional differential fo F2 maps present the regional variations of the nighttime electron density peaks and decay. Both the pre- and post-midnight enhancements initially appeared in a belt almost in parallel with the eclipse track and then drifted southward. The different magnitudes of greatest eclipse in the umbra and outside tend to account for the different occurrence times of the plasma flux. The ionospheric decay following the pre-midnight enhancement is also considered as a consequence of the eclipse shade.

Description: [1]  We study a statistics of ∂  B z /∂  x in a thin stretching current sheet (substorm growth phases) observed by Cluster between 8 and 18 R E downtail. After 2005 spacecraft separation allowed to measure directly this derivative of B z along the tail axis. The near-tail events (within 14 R E ) exhibited a straight decrease of an initially large positive ∂  B z /∂  x to ∼ 1–2 nT/ R E . In the more stretched middle tail, usually the small | ∂  B z /∂  x | 〈 0.5 − 1 nT/ R E had no clear trend and fluctuated around zero with time scales 5–15 min. In general, negative ∂  B z /∂  x were ubiquitous. At some onsets larger negative ∂  B z /∂  x  〈 − 1 nT/ R E were associated with transient dipolarizations, propagating Earthward. There was no clear association of local plasma sheet activity onset with any value of ∂  B z /∂  x . We discuss relation of observations and recent modeling results.

Description: [1]  We investigated the evolution of seismicity and deformation in the unstable eastern flank of Etna volcano over a thirty-year period (from 1980 to 2012). A significant temporal correlation has been revealed between periods of flank acceleration and intensified seismic activity by comparing seismicity along the northern border (Pernicana fault system) of the sliding flank and the deformation of the eastern flank. Two marked phases have been observed in 1984-1986 and in the years following 2002. These two phases are separated by an intermediate phase from 1987 to 2001, in which the flank sliding slowed down and the seismicity dropped drastically. This common temporal evolution in the deformation rate and seismic release supports the hypothesis that the seismicity in the northern border can be viewed as a marker of the response to accommodate the stress exerted by the traction of the eastern flank sliding. This interplay has also been corroborated by Finite Element Method (FEM) numerical computations that highlight a good correlation between the seismicity pattern and areas of positive stress changes induced by the sliding surface. The two intense phases of flank acceleration are representative of two main different sources: volcano flank instability stretching the eastern sector in the first 1984-1986 phase and magmatic intrusions pushing the eastern flank seaward since the 2002-03 eruption. Establishing the relationship between flank acceleration and seismic activation, therefore, contributes to understanding Etna's mechanical behavior, and provides insights into the processes regulating the unstable flank response.

Description: [1]  Numerous studies in the Central Pyrenees have provided evidence for a rapid phase of exhumation of this mountain belt during the Late Eocene (37–30 Ma). Simultaneously, the closure of the Ebro foreland basin allowed the accumulation of sediments at the southern Piedmont, which partially covered the fold-and-thrust belt from Late Eocene ( e . g . when it was still actively deforming) to Miocene times. We aim here at understanding the consequences of such syn-tectonic sedimentation on the Southern Pyrenean fold-and-thrust belt by using a 2-D numerical model that reproduces the development of a thin-skinned wedge subject to different modes of sedimentation and erosion. The results show contrasting fold-and-thrust belt behavior when applying aggrading or prograding sedimentation, which we link to the critical state of the wedge. When the sediments are sourced from the hinterland (progradation), the thrusting propagates toward the foreland; whereas when the sediments aggrade from the basin, the thrusting sequence migrates backward. This latter mode shows patterns of deformation that compare favorably to the Pyrenean thrusting sequence observed during Eocene-Miocene times.

Description: [1]  ELF/VLF radio waves are difficult to generate with conventional antennas. Ionospheric HF heating facilities generate ELF/VLF waves via modulated heating of the lower ionosphere. HF heating of the ionosphere changes the lower ionospheric conductivity, which in the presence of natural currents such as the auroral electrojet, creates an antenna in the sky when heating is modulated at ELF/VLF frequencies. We present a summary of nearly 100 days of ELF/VLF wave generation experiments at the 3.6 MW HAARP facility near Gakona, Alaska, and provide a baseline reference of ELF/VLF generation capabilities with HF heating. Between February 2007 and August 2008, HAARP was operated on close to 100 days for ELF/VLF wave generation experiments, at a variety of ELF/VLF frequencies, seasons and times of day. We present comprehensive statistics of generated ELF/VLF magnetic fields observed at a nearby site, in the 500-3500 Hz band. Transmissions with a specific HF beam configuration (3.25 MHz, vertical beam, amplitude modulation) are isolated so the data comparison is self-consistent, across nearly 5 million individual measurements of either a tone or a piece of a frequency-time ramp. There is a minimum in the average generation close to local midnight. It is found that generation during local nighttime is on average weaker, but more highly variable, with a small number of very strong generation periods. Signal amplitudes from day to day may vary by as much as 20-30 dB. Generation strengthens by ~5 dB during the first ~30 minutes of transmission, which may be a signature of slow electron density changes from sustained HF heating. Theoretical calculations are made to relate the amplitude observed to the power injected into the waveguide and reaching250 km. The median power generated by HAARP and injected into the waveguide is ~0.05-0.1 W in this base-line configuration (vertical beam, 3.25 MHz, amplitude modulation), but may have generated hundreds of Watts for brief durations. Several efficiency improvements have improved the ELF/VLF wave generation efficiency further.

Description: [1]  In southern California, fault slip rate estimates along the San Andreas fault (SAF) and Garlock fault from geodetically-constrained kinematic models are systematically at the low end or lower than geologic slip rate estimates. The sum of geodetic model slip rates across the Eastern California Shear Zone is higher than the geologic sum. However, the ranges of reported model and geologic slip rate estimates in the literature are sufficiently large that it remains unclear whether these apparent discrepancies are real, or attributable to epistemic uncertainties in the two types of estimates. We further examine uncertainties in geodetically-derived slip rate estimates on major faults in southern California by conducting a suite of inversions with four kinematic models. Long-term-rigid elastic block models constrained by the geologic slip rates cannot fit the present-day GPS-derived velocity field. Deforming (permanent off-fault strain) elastic block models and viscoelastic earthquake cycle block models constrained by geologic slip rates can fit the present-day GPS-derived velocity field with 28-33% of the total geodetic moment rate occurring as distributed deformation off of the major faults. Models incorporating viscoelastic mantle flow predict systematically higher slip rates than purely elastic models on many of the the major southern California faults with ranges of (elastic/viscoelastic) 29-34/30-37 mm/yr for the Carrizo SAF segment, 20-24/20-32 mm/yr for the Mojave SAF segment, 14-17/18-22 mm/yr for the Coachella SAF segment, 13-19/14-22 mm/yr for the San Jacinto fault, and 5-11/5-11 mm/yr for the western Garlock fault.

Description: [1]  An inversion technique for estimating the properties of the magnetospheric plasma from the harmonic frequencies of the toroidal standing Alfvén waves has been used to derive the global equatorial mass density covering radial distances from 4 to 9 Earth radii ( R E ), within the local time sector spanning from 0300 to 1900 hours. This broad range of L shell extending to the outer magnetosphere allows us to examine the local time and radial dependence of the quiet-time equatorial mass density during solar minimum and thereby construct a global distribution of the equatorial mass density. The toroidal Alfvén waves were detected with magnetometers on the Active Magnetospheric Particle Tracer Explorers (AMPTE)/Charge Composition Explorer (CCE) during the nearly 5 year interval from August 1984 to January 1989 and on the Geostationary Operational Environmental Satellites (GOES) (10, 11 and 12) for 2 years from 2007 to 2008, both of which were operating during solar minimum years. The derived equatorial mass density, ρ eq , at geosynchronous orbit (GEO) monotonically increases with increasing magnetic local time (MLT) from the nightside towards the dusk sector. At other radial distances, ρ eq has the same MLT variation as that of GEO, while the magnitude logarithmically decreases with increasing L value. An investigation of the Dst and Kp dependence shows that the median value of ρ eq varies little in the daytime sector during moderately disturbed times, which agrees with previous studies. ρ eq calculated from the F 10.7 dependent empirical model shows good agreement with that of CCE but overestimates that of GOES probably due to the extreme solar cycle minimum in years 2007–2008.

Description: [1]  In this paper we compare observations of the high latitude cusp from DMSP data to simulations conducted using the Lyon-Fedder-Mobarry (LFM) global magnetosphere simulation. The LFM simulation is run for the 31 Aug 2005 to 02 Sep 2005 moderate storm, from which the solar wind data exhibits a wide range of conditions that enable a statistical representation of the cusp to be obtained. The location of the cusp is identified using traditional magnetic depression and plasma density enhancement at high altitude. A new diagnostic using the parallel ion number flux is also tested for cusp identification. The correlation of the cusp latitude and various solar wind IMF coupling functions is explored using the three different cusp identification methods. The analysis shows 1) the three methods give approximately the same location and size of the simulated cusp at high altitude; 2) the variations of the simulated cusp are remarkably consistent with the observed statistical variations of the low-altitude cusp. In agreement with observations a higher correlation is obtained using other solar wind coupling functions such as the Kan-Lee electric field. The MLT position of the simulated cusp is found to depend upon the IMF By component, with a lower linear correlation. The width of the simulated cusp in both latitude and MLT is also examined. The size of the cusp is found to increase with the solar wind dynamic pressure with saturation seen when the dynamic pressure is greater than 3 nPa.

Description: [1]  Daily profiles of phase measurements as observed on fixed VLF-paths generally show a transient phase advance, followed by a phase delay, for about 90 minutes after sunrise hours. This is indicative of a reflecting ionospheric C-region developing along the terminator line at an altitude below the normal D-region. The suggested occurrence of a C-region is consistent with rocket measurements made in the 1960's, showing a maximum of the electron density between 64 and 68 km, and by radio sounding in the 1980's. In order to correctly describe the properties of the phase effect associated with the presence of a C-region, it is important to understand the subionospheric propagation characteristics of the VLF-paths. In this paper, we analyze the variations presented by the temporal properties of the VLF narrow-band phase effect, and determined a parameter associated with the appearance of the C-region at sunrise hours observed by receivers from the South America VLF Network (SAVNET). Periodic patterns emerge from the parameter curves. Two distinct temporal behavior regimes can be identified: one exhibiting slow variations between March and October, and another one exhibiting faster variations between October and March. Solar illumination conditions and the geometrical configuration of the VLF paths relative to the sunrise terminator partly explain the slow variation regime. During periods of faster variations, we have observed good association with atmospheric temperature variability found in the measurements of the TIMED-SABER satellite instrument, which we assume to be related to the Winter Anomaly atmospheric phenomenon. However, when comparing the parameter time series with temperature curves, no direct one-to-one correspondence was found for transient events.

Description: [1]  On 21 January 2005, one of the fastest interplanetary coronal mass ejections (ICME) of solar cycle 23, containing exceptionally dense plasma directly behind the sheath, hit the magnetosphere. We show from charge-state analysis that this material was a piece of the erupting solar filament, and further, based on comparisons to the simulation of a fast CME, that the unusual location of the filament material was a consequence of three processes. As the ICME decelerated, the momentum of the dense filament material caused it to push through the flux rope towards the nose. Diverging non-radial flows in front of the filament moved magnetic flux to the sides of the ICME. At the same time reconnection between the leading edge of the ICME and the sheath magnetic fields worked to peel away the outer layers of the flux rope creating a remnant flux rope and a trailing region of newly opened magnetic field lines. These processes combined to move the filament material into direct contact with the ICME sheath region. Within one hour after impact and under northward IMF conditions, a cold dense plasma sheet formed within the magnetosphere from the filament material. Dense plasma sheet material continued to move through the magnetosphere for more than 6 hours as the filament passed by the Earth. Densities were high enough to produce strong diamagnetic stretching of the magnetotail despite the northward IMF conditions and low levels of magnetic activity. The disruptions from the filament collision are linked to an array of unusual features throughout the magnetosphere, ionosphere and atmosphere. These results raise questions about whether rare collisions with solar filaments may, under the right conditions, be a factor in producing even more extreme events.

Description: [1]  We present observations from the Falkland Islands SuperDARN radar of the propagation of HF radio waves via the Weddell Sea ionospheric Anomaly (WSA), a region of enhanced austral summer nighttime ionospheric electron densities covering the southern Pacific and South Americas region. This anomaly is thought to be produced by uplift of the ionosphere by prevailing equatorwards thermospheric winds. Of particular interest are perturbations of the WSA-supported propagation, which suggest that during periods of geomagnetic disturbance the ionospheric layer can be lowered by several 10s of km and subsequently recover over a period of 1 to 2 hours. Perturbations can appear singly or as a train of 2 to 3 events. We discuss possible causes of the perturbations, and conclude that they are associated with equatorward-propagating large-scale atmospheric waves produced by magnetospheric energy deposition in the auroral or sub-auroral ionosphere. Changes in high/mid-latitude electrodynamics during geomagnetic storms may also account for the perturbations, but further modeling is required to fully understand their cause.

Description: [1]  Inner magnetosphere magnetic field and plasma flow data are examined during 228 steady magnetospheric convection events. We find that the B Z component of the magnetic field around geostationary orbit is weaker than during average conditions and the plasma flow speeds are higher than average in the dusk sector just beyond geostationary orbit. The SMC periods include more enhanced Earthward and tailward flow intervals than during average conditions. The steady convection period magnetic field is not steady: The near-geostationary nightside field grows increasingly taillike throughout the steady convection period. In the midtail, Earthward flows are enhanced in a wide region around the midnight sector, which leads to enhanced magnetic flux transport toward the Earth during the steady convection periods. Compared to well-known characteristics during magnetospheric substorms, the inner tail evolution resembles that during the substorm growth phase, while the midtail flow characteristics duringsteady convection periods are similar to those found during substorm recovery phases.

Description: [1]  A global total electron content (TEC) model response to geomagnetic activity described by the K p -index is built by using the Center for Orbit Determination of Europe (CODE) TEC data for full 13 years, January 1999 - December 2011. The model describes the most probable spatial distribution and temporal variability of the geomagnetically forced TEC anomalies assuming that these anomalies at a given modip latitude depend mainly on the K p -index, local time (LT) and longitude. The geomagnetic anomalies are expressed by the relative deviation of TEC from its 15-day median and are denoted as rTEC. The rTEC response to the geomagnetic activity is presented by a sum of two responses with different time delay constants and different sign of the cross-correlation function. It has been found that the mean dependence of rTEC on K p -index can be expressed by a cubic function. The LT dependence of rTEC is described by Fourier time series which includes the contribution of four diurnal components with periods 24, 12, 8 and 6 hours. The rTEC dependence on longitude is presented by Fourier series which includes the contribution of zonal waves with zonal wavenumbers up to 6. In order to demonstrate how the model is able to reproduce the rTEC response to geomagnetic activity three geomagnetic storms at different seasons and solar activity conditions are presented. The model residuals clearly reveal two types of the model deviation from the data: some underestimation of the largest TEC response to the geomagnetic activity and randomly distributed errors which are the data noise or anomalies generated by other sources. The presented TEC model fits to the CODE TEC input data with small negative bias of -0.204, root mean squares error RMSE  = 4.592 and standard deviation error STDE  = 4.588. The model offers TEC maps which depend on geographic coordinates (5 o x5 o in latitude and longitude), and universal time (UT) at given geomagnetic activity and day of the year. It could be used for both science and possible service (nowcasting and short-term prediction); for the latter a detailed validation of the model at different geophysical conditions has to be performed in order to clarify the model predicting quality.

Description: [1]  Magnetic reconnection (MR), a fundamental process in space plasmas that changes magnetic topology and converts magnetic energy into kinetic and thermal energies, is an ultimate driver of space weather. There exist two models of MR in the literatures, anti-parallel and component, associated with intensive studies of their generation and applications. In this paper we report an MR event observed by Cluster constellation in the geo-magnetotail where both types of MR were detected. By reconstructing the three-dimensional (3D) MR configuration we find that a pair of A-B nulls existed in both types of MR cases with two fan surfaces intersecting each other to form a separator line connecting the nulls. A weak or sizable magnetic field exists along the separator in the anti-parallel or component case, respectively. In the latter case, field strength is finite away from the two nulls and vanishes close to the nulls. Therefore, at least in the two cases observed, both anti-parallel and component MR geometries are local presentation of the separator MR configuration. This result supports the expectation that 3D nulls often occur as a crucial element of MR at least in the magnetotail and separator MR may play an important role in dynamics and reconfiguration of magnetic field in 3D MR processes.

Description: [1]  The development of equatorial plasma irregularity plumes can be well recorded by steerable backscatter radars operated at and off the magnetic equator due to the fact that the vertically extended plume structures are tracers of magnetically north-south aligned larger scale structures. From observations during March 2012, using two low latitude steerable backscatter radars in Southeast Asia, the Equatorial Atmosphere Radar (EAR) (0.2ºS, 100.3ºE; dip lat 10.4ºS) and the Sanya VHF radar (18.4ºN, 109.6ºE; dip lat 12.8ºN), the characteristics of backscatter plumes over the two sites separated in longitude by ~1000 km were simultaneously investigated. The beam steering measurements reveal frequent occurrences of multiple plumes over both radar sites, of which two cases are analyzed here. The observations on 30 March 2012 show plume structures initiated within the radar scanned area, followed by others drifting from the west of the radar beam over both stations. A tracing analysis on the onset locations of plasma plumes reveals spatially well-separated backscatter plumes, with a maximum east-west wavelength of about 1000 km, periodically generated in longitudes between 85ºE and 110ºE. The post-sunset backscatter plumes seen by the Sanya VHF radar are found to be due to the passage of sunset plumes initiated around the longitude of EAR. Most interestingly, the EAR measurements on the night of 21 March 2012 show multiple plume structures that developed successively in the radar scanned area with east-west separation of ~50 km, with however, no sunset plasma plume over Sanya. Co-located ionogram measurements show that spread F irregularities occurred mainly in the bottomside F-region at Sanya, whereas satellite traces in ionograms that are indications of large-scale wave structures, were observed on that night at both stations. Possible causes for the longitudinal difference in the characteristics of radar backscatter plumes are discussed.

Description: [1]  Determining the scale-length, magnitude, and distribution of heterogeneity in the lowermost mantle is crucial to understanding whole mantle dynamics, and yet it remains a much debated and ongoing challenge in geophysics. Common shortcomings of current seismically-derived lowermost mantle models are incomplete raypath coverage, arbitrary model parameterization, inaccurate uncertainty estimates, and an ad hoc definition of the misfit function in the optimization framework. In response, we present a new approach to global tomography. Apart from improving the existing raypath coverage using only high quality cross-correlated waveforms, the problem is addressed within a Bayesian framework where explicit regularization of model parameters is notrequired. We obtain high resolution images, complete with uncertainty estimates, of the lowermost mantle P-wave velocity structure using a hand-picked dataset of PKPab-df, PKPbc-df, and PcP-P differential traveltimes. Most importantly, our results demonstrate that the root mean square of the P-wave velocity variations in the lowermost mantle is approximately 0.87%, which is three times larger than previous global-scale estimates.

Description: [1]  The solar wind impacting the Earth varies over a wide range of time scales, driving a corresponding range of geomagnetic activity. Past work has strongly indicated the rate of merging on the frontside magnetosphere is the most important predictor for magnetospheric activity, especially over a few hours. However the magnetosphere exhibits variations on other time scales, including UT, seasonal, and solar cycle variations. Much of this geomagnetic variation cannot be reasonably attributed to changes in the solar wind driving – that is, it is not created by the original Russell-McPherron effect or any generalization thereof. In this paper we examine the solar cycle, seasonal, and diurnal effects based upon the frequency of substorm onsets, using a data set of 53,000 substorm onsets. These were identified through the SuperMAG collaboration and span three decades with continuous coverage. Solar cycle variations include a profound minima in 2009 (448 substorms) and peak in 2003 (3727). The magnitude of this variation (a factor of 8.3) is not explained through variations in estimators of the frontside merging rate (such as d Φ MP / dt ), even when the more detailed probability distribution functions are examined. Instead, v , or better, n 1/2 v 2 seems to be implicated in the dramatic difference between active and quiet years, even beyond the role of velocity in modulating merging. Moreover, we find that although most substorms are preceded by flux loading (78.5% are above the mean and 83.8% above median solar wind driving) a high solar wind v is almost as important (68.3% above mean, 74.8% above median). This and other evidence suggests either v or n 1/2 v 2 (but probably not p ) plays a strong secondary role in substorm onset. As for the seasonal and diurnal effects, the elliptical nature of the Earth's orbit, which is closest to the Sun in January, leads to a larger solar wind driving (measured by Bs , vBs , or d Φ MP / dt ) in November, as is confirmed by 22 years of solar wind observations. However substorms peak in October and March, and have a UT dependence best explained by whether a conducting path established by solar illumination exists in at least one hemisphere in the region where substorm onsets typically occur.

Description: [1]  The primary sources of energetic electron precipitation (EEP) which affect altitudes 〈100 km (〉30 keV) are expected to be from the radiation belts, and during substorms. EEP from the radiation belts should be restricted to locations between L  = 1.5-8, while substorm produced EEP is expected to range from L  = 4-9.5 during quiet geomagnetic conditions. Therefore, one would not expect any significant D-region impact due to electron precipitation at geomagnetic latitudes beyond about L  = 10. In this study we report on large unexpectedly high latitude D-region ionization enhancements, detected by an incoherent scatter radar at L  ≈ 16, which appear to be caused by electron precipitation from substorms. We go on to reexamine the latitudinal limits of substorm produced EEP using data from multiple low-Earth orbiting spacecraft, and demonstrate that the precipitation stretches many hundreds of kilometers polewards of the previously suggested limits. We find that a typical substorm will produce significant EEP over the IGRF L -shell range L  = 4.6 ± 0.2-14.5 ± 1.2, peaking at L  = 6-7. However, there is significant variability from event to event; in contrast to the median case, the strongest 25% of substorms have significant EEP in the range spaning L  = 4.1 ± 0.1-20.7 ± 2.2, while the weakest 25% of substorms have significant EEP in the range spaning L  = 5.5 ± 0.1-10.1 ± 0.7. We also examine the occurrence probability of very large substorms, focusing on those events which appear to be able to disable geostationary satellites when they are located near midnight MLT. On average these large substorms occur approximately 1-6 times per year, a significant rate given the potential impact on satellites.

Description: [1]  This paper presents new observations of the behaviour of simulated dust particles in space plasma based on a 3D particle in cell code. Multistep Monte Carlo collision is employed to simulate the dust charging process which is validated for the cases of charging of isolated dust particle and ensemble dust particles, where results indicate good agreement between simulation and theories. The code is then used to investigate plasma properties near a charged surface in a vicinity of a cloud of dust particles. The simulation reveals that a cloud of dust particle close to a spacecraft surface affects plasma densities around the spacecraft as well the spacecraft's surface potential. It is suggested that dust cloud causes the surface to charge to higher negative potential. The simulation also suggests that the combination of surface potential and dust cloud potential produces a region of trapped low energy electrons.

Description: [1]  We here document with magnetic field observations a passage of the MESSENGER spacecraft through Mercury's magnetosphere under conditions of a quasi-parallel bow shock, i.e., when the direction of the upstream interplanetary magnetic field (IMF) was within 45° of the bow shock normal. The spacecraft's fast transition of the magnetosheath and the steady solar wind conditions during the period analyzed allow both spatial and temporal properties of the shock crossing to be investigated. The observations show that the shock reformation process can be nearly periodic under stable solar wind conditions. Throughout the 25-min-long observation period, the pulsation duration deviated by at most ~10% from the average 10 s period measured. This quasi-periodicity allows us to study all aspects of the shock reconfiguration, including ultra-low-frequency waves in the upstream region and large-amplitude magnetic structures observed in the vicinity of the magnetosheath − solar wind transition region and inside the magnetosheath. We also show that bow shock reformation can be a substantial source of wave activity in the magnetosphere, on this occasion having given rise to oscillations in the magnetic field with peak-to-peak amplitudes of 40–50 nT over large parts of the dayside magnetosphere. The clean and cyclic behavior observed throughout the magnetosphere, the magnetosheath, and the upstream region indicates that the subsolar region was primarily influenced by a cyclic reformation of the shock front, rather than by a spatial and temporal patchwork of short large-amplitude magnetic structures (SLAMS), as is generally the case at the terrestrial bow shock under quasi-parallel conditions.

Description: [1]  The entry of solar wind into the magnetosphere is strongly influenced by kinetic-scale boundary layers where the rapid variation in the magnetic field and/or velocity can drive transport. In current layers with strong Alfvénic velocity shear, the generation of vortices from the Kelvin-Helmholtz instability can drive magnetic reconnection even in broader current sheets by locally compressing these layers as the vortices develop. Previous two-dimensional (2D) fully kinetic simulations of this vortex-induced reconnection process have demonstrated the copious formation of magnetic islands in regions of strongly compressed current between the vortices. Here we describe the first three-dimensional (3D) fully kinetic simulations of this process and demonstrate that the compressed current sheets give rise to magnetic flux ropes over a range of oblique angles and along the entire extent of the compressed current layer around the periphery of the vortex. These flux ropes propagate with the shear flow and eventually merge with the vortex. Over longer time scales, this basic scenario is repeated as the vortices drive new compressed current sheets. In the final stage, the vortices undergo a merging process that drives new compressed current sheets and flux ropes. Based on these simulations, a simple model is proposed that predicts the size of these flux ropes relative to their parent vortex. Both the relative sizes as well as the structure of the profiles across the vortex are in reasonable agreement with THEMIS observations at the Earth's low-latitude magnetopause.

Description: [1]  We survey the properties of electron pitch angle distributions in the magnetotail plasma sheet at a distance between 15 and 19 R E from the Earth, using data from the Cluster PEACE instrument. We limit our survey to those pitch angle distributions measured when the IMF had been steadily northward or steadily southward for the previous three hours. We find that, at sub-keV energies the plasma sheet electron pitch angle distribution has an anisotropy such that there is a higher differential energyflux of electrons in the (anti-) field-aligned directions. Fitting the measured pitch angle distributions with both a single and two component kappa distribution reveals that this anisotropy is the result of the presence of a second, cold, component of electrons that is observed more often than not, and occurs during both the northward and southward IMF intervals. We present evidence that suggests the cold electron component has an ionospheric, rather than magnetosheath, source and is linked tothe large scale field aligned current systems that couple the magnetosphere and ionosphere.

Description: [1]  The Dzhungarian strike-slip fault of Central Asia is one of a series of long, NW-SE right-lateral strike-slip faults that are characteristic of the northern Tien Shan region, and extends over 300 km from the high mountains into the Kazakh Platform. Our field-based and satellite observations reveal that the Dzhungarian fault can be characterised by three 100 km long sections based on variation in strike direction. Through morphological analysis of offset streams and alluvial fans, and through OSLdating, we find that the Dzhungarian fault has a minimum average late Quaternary slip rate of 2.2 ± 0.8 mm/yr and accommodates N-S shortening related to the India-Eurasia collision. This shortening may also be partly accommodated by counter-clockwise rotation about a vertical axis. Evidence for a possible paleo-earthquake rupture indicates that earthquakes up to at least Mw 7 can be associated with just the partitioned component of reverse slip on segments of the central section of the fault up to 30 km long. An event rupturing longer sections of the Dzhungarian fault has the potential to generate greater magnitude earthquakes ( Mw 8), however long time periods (e.g. thousands of years) are expected in order to accumulate enough strain to generate such earthquakes.

Description: [1]  Driver functions for the Earth's magnetosphere-ionosphere system are derived from physical principles. Two processes act simultaneously: a reconnection-coupled MHD generator G and a viscous interaction. G accounts for the dayside reconnection rate, the length of the reconnection X-line, and current-saturation limits for the solar-wind generator. Two viscous drivers are derived: Bohm viscosity B and the freestream-turbulence effect F. A problematic proxy effect is uncovered wherein the viscous driver functions also describe the strength of reconnection. Two magnetospheric-driver functions written in terms of upstream-solar-wind parameters are constructed: G+B and G+F. The driver functions are tested against 7 geomagnetic indices. The reaction of the geomagnetic indices to G+B and G+F is nonlinear: nonlinear versions of the driver functions are supplied. Applying the driver functions at multiple timesteps yields correlation coefficients of ~85% with the AE and Kp indices; it is argued that multiple timestepping removes high-frequency uncorrelated signal from the drivers. Autocorrelation-function analysis shows strong 1-d and 1-yr periodicities in the AE index that are not in the solar-wind driver functions; correspondingly, highpass and lowpass filtering finds uncorrelated signal at 1-d and 1-yr timescales. Residuals (unpredicted variance) between the geomagnetic indices and the driver functions are analyzed: the residuals are anti-correlated with the solar-wind velocity, the solar F 10.7 radio flux, and the solar-wind current-saturation parameter. Removing diurnal, semiannual, and annual trends from the indices improves their correlation with the solar-wind driver functions. Simplified versions of the driver functions are constructed: the simplified drivers perform approximately as well as the full drivers.

Description: [1]  Plasma properties of Saturn's pre-midnight tail region are surveyed using Cassini/CAPS ion observations from 2010. Only low-latitude (|lat| 〈 6°) intervals in which the CAPS viewing was roughly symmetric inward and outward around the corotation direction are used. Our numerical moments algorithm returns nonzero ion density for 70% (999) of the intervals selected. Of these, 642 had detectable water-group ion densities, and the remainder were dominantly, if not entirely, light ions. The derived plasma parameters are similar to those found in an earlier study for the post-midnight tail region, except that we find little evidence for the systematic outflows identified in that study, and we do find numerous significant inflow events. One such inflow is identified as a dipolarization event, the first reported plasma properties of such a structure at Saturn. A second, long-lasting event may be evidence for the existence at times of a quasi-steady reconnection region in the pre-midnight tail. The large majority of the plasma flows are found to be within 20° of the corotation direction, though with flow speeds significantly lower than full corotation. While the inflow events represent plausible evidence for internally-driven mass loss in the pre-midnight region, the absence of significant outflow events suggests that in the region surveyed here, tail reconnection has not yet proceeded to involve lobe field lines, so the disconnected plasma continues its general motion in the corotation direction.

Description: [1]  Following the first-time ionospheric imaging of a seismic fault, here we perform a case-study on retrieval of parameters of the extended seismic source ruptured during the great M9.0 Tohoku-oki earthquake. Using 1Hz ionospheric GPS data from the Japanese network of GPS receivers (GEONET) and several GPS satellites, we analyze spatio-temporal characteristics of co-seismic ionospheric perturbations and we obtain information on the dimensions and location of the sea surface uplift (seismic source). We further assess the criterion for the successful determination of seismic parameters from the ionosphere: the detection is possible when the line-of-sights from satellites to receivers cross the ionosphere above the seismic fault region. Besides, we demonstrate that the multi-segment structure of the seismic fault of the Tohoku-oki earthquake can be seen in high-rate ionospheric GPS-data. Overall, our results show that, under certain conditions, ionospheric GPS-derived TEC measurements could complement the currently working systems, or independent ionospherically-based system might be developed in the future.

Description: [1]  The non-dipolar portions of the Earth's main magnetic field constitute substantial differences between the two hemispheres. Beside the magnetic flux densities and patterns being different in the Northern Hemisphere (NH) and Southern Hemisphere (SH), also the offset between the invariant magnetic and the geographic poles is larger in the SH than in the NH. We investigated the effects of this magnetic field asymmetry on the high-latitude thermosphere and ionosphere using global numerical simulations, and compared our results with recent observations. While the effects on the high-latitude plasma convection are small, the consequences for the neutral wind circulation are substantial. The cross-polar neutral wind and ion drift velocities are generally larger in the NH than the SH, and the hemispheric difference shows a semidiurnal variation. The neutral wind vorticity is likewise larger in the NH than in the SH, with the difference probably becoming larger for higher solar activity. In contrast, the spatial variance of the neutral wind is considerably larger in the SH polar region, with the hemispheric difference showing a strong semidiurnal variation. Its phase is similar to the phase of the semidiurnal variation of the hemispheric magnitude differences. Hemispheric differences in ion drift and neutral wind magnitude are most likely caused partly by the larger magnetic flux densities in the near-polar regions of the SH, and partly by the larger offset between the invariant and geographic pole in the SH, while differences in spatial variance are probably just caused by the latter. We conclude that the asymmetry of the magnetic field, both in strength and in orientation, establishes substantial hemispheric differences in the neutral wind and plasma drift in the high-latitude upper atmosphere, which can help to explain observed hemispheric differences with EDI/Cluster and CHAMP.

Description: [1]  We have identified 8 events with double-belt structure in the outer radiation belt from 110 CME-driven magnetic storms and 223 CIR-driven storms during 1994 to 2003 based on the SAMPEX data sets. Among them, 3 cases are related to CME-driven magnetic storms and 5 cases are related CIR-driven storms. All double-belt structure events in the outer radiation belt are found during the recovery phase of a magnetic storm for both CME- and CIR-related events—they usually start to form within 3–4 days after the onset of the magnetic storm. The pre-conditions needed to form a double-belt structure, for all the CME-related events, are found to be high solar wind dynamic pressure ( P dy ) and southward IMF Bz; Nevertheless, for the CIR-related events, they are found to be associated with high speed stream with southward interplanetary magnetic field, which is enhanced by a suitably orientated B y component.It is further found that the flux distributions of the double-belt structure can be fitted well with a simply exponential decay function of L ∗ . Based on the RBC index, the proportion of the total number of 1.5-6.0 MeV electrons inside the position of maximum fluxes to that outside the maximum fluxes keeps rising during the double-belt period, which implies that the acceleration mainly occurs at regions inside the location of maximum fluxes. We suggest that the plasmapause and the strong wave-particle interactions with VLF and ULF waves near it, play an important role in the development of the double-belt structures.

Description: [1]  The solar minimum period during 2008–2009 was characterized by lower thermospheric density than the previous solar minimum, and lower than any previously measured. Recent work [ Solomon et al ., 2010; 2011] used the NCAR Thermosphere-Ionosphere-Electrodynamics General Circulation Model to show that the primary cause of density changes from 1996 to 2008 was a small reduction in solar extreme-ultraviolet irradiance (EUV), causing a decrease in thermospheric temperature and hence a contracted thermosphere. There are similar effects in the ionosphere, with most measurements showing an F -region ionosphere that is unusually low in density, and in peak altitude. This paper addresses the question of whether model simulations previously conducted, and their solar, geomagnetic, and anthropogenic inputs, produce ionospheric changes commensurate with observations. We conducted a 15-year model run and obtained good agreement with observations of the global mean thermospheric density at 400 km throughout the solar cycle, with a reduction of ~30% from the 1996 solar minimum to 2008–2009. We then compared ionosonde measurements of the mid-day peak density of the ionospheric F -region ( N m F 2 ) to the model simulations at various locations. Reasonable agreement was obtained between measurements and the model, supporting the validity of the neutral density comparisons. The global average N m F 2 was estimated to have declined between the two solar minima by ~15%. In these simulations, a 10% reduction of solar EUV plays the largest role in causing the ionospheric change, with a minor contribution from lower geomagnetic activity, and a very small additional effect from anthropogenic increase in CO 2 .

Description: [1]  We report the statistical features of sporadic sodium layers (SSLs) and the thermospheric enhanced sodium layers (TeSLs) observed by a lidar chain located at Beijing (40.2°N, 116.2°E), Hefei (31.8°N, 117.3°E), Wuhan (30.5°N, 114.4°E), and Haikou (19.5°N, 109.1°E). The average SSL occurrence rate was approximately 46.0, 12.3, 13.8, and 15.0 hr per SSL at Beijing, Hefei, Wuhan, and Haikou, respectively. However, the TeSLs occurred relatively infrequently and were more likely to appear at low and high latitudinal sites. Both the SSLs and TeSLs at four lidar sites showed evident summer enhancements and correlated well with Es ( f o E s  〉  4  MHz ). The co-observations of SSLs at three lidar site pairs, i.e., Hefei – Beijing, Hefei – Wuhan and Hefei – Beijing, indicated that a large-scale SSL extended horizontally for at least a few hundred kilometers and exhibited a tidal-induced modulation. Moreover, the SSLs were better correlated for the Hefei – Wuhan and Hefei – Haikou pairs than the Hefei – Beijing pair, which suggested a difference in the dynamical/chemical process in mesosphere and lower thermosphere (MLT) between the Beijing site and the other sites.

Description: [1]  A highly sensitive all-sky EMCCD airglow imager has been operative in Longyearbyen, Norway (78.1°N, 15.5°E) since October 2011. The imager obtains the 630.0 nm all-sky images with an exposure time of 4 sec, which is about 10 times shorter than the conventional cooled CCD airglow imagers. This new equipment allows us to image the on-going structuring of polar cap patches in 2D fashion. Here, we report a case in which faint undulations appeared along the trailing edge of patches propagating in the central polar cap. The separation between the fingers in the undulations was about 50–100 km and the e -folding time of their growth was ~5 min. We suggest that the gradient-drift instability (GDI) is one of the possible generation mechanisms of the undulating structures. The reasons for this interpretation are 1) the asymmetry in the preference of structuring between the leading and trailing edges is qualitatively consistent with the GDI mechanism, and 2) the linear growth rate of GDI calculated by using electron density estimates from simultaneous EISCAT Svalbard radar observations is roughly consistent with the observed growth time of the fingers. Such “unstable polar cap patches" could be important sources of seed irregularities, which would eventually be broken down to smaller-scale density perturbations affecting the trans-ionospheric satellite communications in the central polar cap.

Description: [1]  As the Polar spacecraft apogee precessed through the magnetic equator in 2001, Polar encountered numerous substorm events in the region between geosynchronous orbit and 10 Re geocentric distance; most of them in the plasma sheet boundary layers. Of these, a small number were recorded near the neutral sheet in the evening sector. Polar/TIDE provides a unique perspective on the lowest energy ion plasma, showing that these events exhibited a damped wavelike character, initiated by a burst of radially outward flow transverse to the local magnetic field at ~ 80 km/sec. They then exhibit strongly damped cycles of inward/outward flow with a period of several minutes. After one or two cycles, they culminated in a hot plasma electron and ion injection, quite similar to those observed at geosynchronous orbit. Cold plasmaspheric plasmas comprise the outward flow cycles, while the inward flow cycles contain counter-streaming field-parallel polar wind like flows. The observed wave-like structure, preceding the arrival of an Earthward moving substorm injection front, suggests an outward displacement driven by the inward motion at local times closer to midnight, that is, a “snowplow” effect. The damped in/out flows are consistent with interchange oscillations driven by the arrival at the observed local time by an injection originating at greater radius and local time.

Description: [1]  We performed shock compression experiments on preheated forsterite liquid (Mg 2 SiO 4 ) at an initial temperature of 2273 K and have revised the equation of state (EOS) that was previously determined by shock melting of initially solid Mg 2 SiO 4 (300 K). The linear Hugoniot, U S  = 2.674 ± 0.188 + 1.64 ± 0.06 u p km/s, constrains the bulk sound speed within a temperature and composition space as yet unexplored by 1-bar ultrasonic experiments. We have also revised the EOS for enstatite liquid (MgSiO 3 ) to exclude experiments that may have been only partially melted upon shock compression and also the EOS for anorthite liquid, which now excludes potentially un-relaxed experiments at low pressure. The revised fits and the previously determined EOS of fayalite and diopside were used to produce isentropes in the multicomponent CaO-MgO-Al 2 O 3 -SiO 2 -FeO system at elevated temperatures and pressures. Our results are similar to those previously presented for peridotite and simplified “chondrite” liquids such that regardless of where crystallization first occurs, the liquidus solid sinks upon formation. This process is not conducive to the formation of a basal magma ocean. We also examined the chemical and physical plausibility of the partial melt hypothesis to explain the occurrence and characteristics of ultralow velocity zones. We determined that the ambient mantle cannot produce an equilibrium partial melt and residue that is sufficiently dense to be a ULVZ mush. The partial melt would need to be segregated from its equilibrium residue and combined with a denser solid component to achieve a sufficiently large aggregate density.

Description: [1]  High resolution sparker and crustal-scale airgun seismic reflection data, coupled with repeat bathymetric surveys, document a region of repeated coseismic uplift on the portion of the Alaska subduction zone that ruptured in 1964. This area defines the western limit of Prince William Sound. Differencing of vintage and modern bathymetric surveys shows that the region of greatest uplift related to the 1964 Great Alaska earthquake was focused along a series of sub-parallel faults beneath Prince William Sound and the adjacent Gulf of Alaska shelf. Bathymetric differencing indicates that 12 m of coseismic uplift occurred along two faults that reached the sea floor as submarine terraces on the Cape Cleare bank southwest of Montague Island. Sparker seismic reflection data provide cumulative Holocene slip estimates as high as 9 mm/yr along a series of splay thrust faults within both the inner wedge and transition zone of the accretionary prism. Crustal seismic data show that these megathrust splay faults root separately into the subduction zone décollement. Splay fault divergence from this megathrust correlates with changes in mid-crustal seismic velocity and magnetic susceptibility values, best explained by duplexing of the subducted Yakutat terrane rocks above Pacific plate rocks along the trailing edge of the Yakutat terrane. Although each splay fault is capable of independent motion, we conclude that the identified splay faults rupture in a similar pattern during successive megathrust earthquakes and that the region of greatest seismic coupling has remained consistent throughout the Holocene.

Description: [1]  We perform a time-lapse analysis of Rayleigh and Love wave anisotropy above an underground gas storage facility in the Paris Basin. The data were acquired with a three-component seismic array deployed during several days in April and November 2010. Phase velocity and back azimuth of Rayleigh and Love waves are measured in the frequency range 0.2-1.1 Hz using a three-component beamforming algorithm. In both snapshots, higher surface wave modes start dominating the signal above 0.4 Hz with a concurrent increase in back azimuth ranges. We fit anisotropy parameters to the array detections above 0.4 Hz using a bootstrap approach which also provides estimation uncertainty and enables significance testing. The isotropic phase velocity dispersion for Love and Rayleigh waves match for both snapshots. We also observe a stable fast direction of NNW-SSE for Love and Rayleigh waves which is aligned with the preferred orientation of known shallow (〈300 m) and deeper (~1000 m) fault systems in the area, as well as the maximum horizontal stress orientation. At lower frequencies corresponding to deeper parts of the basin, the anisotropic parameters exhibit higher magnitude in the November data. This may perhaps be caused by the higher pore-pressure changes in the gas reservoir in that depth range.

Description: [1]  Eruptive activity at the summit of Kilauea Volcano, Hawaii beginning in 2010 and continuing to the present time is characterized by transient outgassing bursts accompanied by very long period (VLP) seismic signals triggered by rockfalls from the vent walls impacting a lava lake in a pit within the Halemaumau pit crater. We use raw data recorded with a 11-station broadband network to model the source mechanism of signals accompanying two large rockfalls on August 29, 2012 and two smaller average rockfalls obtained by stacking over all events with similar waveforms to improve the signal-to-noise ratio. To determine the source centroid location and source mechanism, we minimize the residual error between data and synthetics calculated by the finite difference method for a point source embedded in a homogeneous medium that takes topography into account. We apply a new waveform inversion method that accounts for the contributions from both translation and tilt in horizontal seismograms through the use of Green's functions representing the seismometer response to translation and tilt ground motions. This method enables a robust description of the source mechanism over the period range 1-1000 s. The VLP signals associated with the rockfalls originate in a source region ~1 km below the eastern perimeter of the Halemaumau pit crater. The observed waveforms are well explained by a simple volumetric source with geometry composed of two intersecting cracks including an east striking crack (dike) dipping 80 ∘ to the north, intersecting a north striking crack (another dike) dipping 65 ∘ to the east. Each rockfall is marked by a similar step-like inflation trailed by decaying oscillations of the volumetric source, attributed to the efficient coupling at the source centroid location of the pressure and momentum changes induced by the rock mass impacting the top of the lava column. Assuming a simple lumped parameter representation of the shallow magmatic system, the observed pressure and volume variations can be modeled with the following attributes: rockfall volume (200 − 4500 m 3 ), length of magma column (120-210 m), diameter of pipe connecting the Halemaumau pit crater to the subjacent dike system (6 m), average thickness of the two underlying dikes (3 – 6 m), and effective magma viscosity (30–210 Pa s). Most rockfalls occur during episodes of sustained deflation of the Kilauea summit. The mass loss rate in the shallow magmatic system is estimated to be 1400  −  15, 000 kg s − 1 based on measurements of the temporal variation of VLP period in the two large rockfalls that occurred on August 29, 2012.

Description: [1]  Accurate knowledge of the global distribution of magnetospheric chorus waves is essential for radiation belt modeling because it provides a direct link to understanding radiation belt losses and acceleration processes. In this paper, we report on newly developed models of the global distribution of chorus amplitudes based on in-situ measurements of IMF and solar wind parameters as well as geomagnetic indices using an artificial neural network technique. We find that solar wind speed and IMF B Z are the most influential parameters that affect the evolution of the magnetospheric chorus. The variations of chorus amplitudes in the outer ( L  ≥ 7) and in the inner (5 ≤  L  〈 7) regions, respectively, are well correlated with the variations of solar wind speed and IMF B Z . In addition, the solar wind parameter-based chorus model generally results in a slightly higher correlation between measured and modeled chorus amplitudes than any other models including geomagnetic indices AE, Kp, and Dst. The developed model shows that the chorus is amplified near the pre-noon sector during the geomagnetically disturbed conditions. With increasing southward IMF B Z the location of peak chorus amplitude moves from the pre-noon sector to the midnight sector, which is due to the enhanced electron injection near midnight. We also present a comparison of diffusive transport simulations for radiation belt electrons interacting with two newly developed chorus models, solar wind parameter-based and geomagnetic index-based chorus models. The comparison between two models shows that the modeling outside the plasmapause can affect the dynamic even inside the plasmasphere because the populations outside the plasmapause can act as seed population to radiation belt particles inside the plasmapause. One weakness of our chorus modeling is that it is trained during the early phase of solar cycle 24 where very few strong storms occurred. Therefore, our model might not be valid in reproducing the chorus activity under extremely disturbed conditions, which should be updated in the future once chorus measurements for such conditions become available.

Description: [1]  We investigate the storm-scale morphology of the magnetospheric magnetic field as well as underlying distributions of electric currents, equatorial plasma pressure and entropy for four Steady Magnetospheric Convection (SMC) events that occurred during the May 2000 and October 2011 magnetic storms. The analysis is made using the empirical geomagnetic field model TS07D, in which the structure of equatorial currents is not predefined and it is dictated by data. The model also combines the strengths of statistical and event-oriented approaches in mining data for the reconstruction of the magnetic field. The formation of a near-Earth minimum of the equatorial magnetic field in the midnight sector is inferred from data without ad hoc assumptions of a special current system postulated in earlier empirical reconstructions. In addition, a new SMC class is discovered where the minimum equatorial field is substantially larger and located closer to Earth. The magnetic field tailward of the minimumis also much larger, and the corresponding area of accumulated magnetic flux may occupy a very short tail region. The equatorial current and plasma pressure are found to be strongly enhanced far beyond geosynchronous orbit and in a broad local time interval covering the whole nightside region. This picture is consistent with independent recent statistical studies of the SMC pressure distributions, global MHD and kinetic RCM-E simulations. Distributions of the flux tube volume and entropy inferred from data reveal different mechanisms of the magnetotail convection crisis resolution for two classes of SMC events.

Description: [1]  Drifts in the Parker spiral interplanetary magnetic field are known to be an important component in the propagation of galactic cosmic rays, while they are thought to be negligible for Solar Energetic Particles (SEPs). As a result they have so far been ignored in SEP propagation modelling and data analysis. We examine drift velocities in the Parker spiral within single particle first-order adiabatic theory, in a local coordinate system with an axis parallel to the magnetic field. We show that, in the presence of scattering in interplanetary space, protons at the high end of the SEP energy range experience significant gradient and curvature drift. In the scatter-free case, drift due to magnetic field curvature is present. The magnitude of drift velocity increases by more than an order of magnitude at high heliographic latitudes compared to near the ecliptic; it has a strong dependence on radial distance r from the Sun, reaching a maximum at r ~1 AU at low heliolatitudes and r ~10 AU at high heliolatitudes. Due to the mass over charge dependence of drift velocities, the effect of drift for partially ionised SEP heavy ions is stronger than for protons. Drift is therefore likely to be a considerable source of cross field transport for high energy SEPs.

Description: [1]  We present a novel method for the automatic retrieval of local plasma density measurements from the Mars advanced radar for sub-surface and ionospheric sounding (MARSIS) active ionospheric sounder (AIS) instrument. The resulting large data set is then used to study the configuration of the Martian ionosphere at altitudes above ~300 km. An empirical calibration routine is used, which relates the local plasma density to the measured intensity of multiple harmonics of the local plasma-frequency oscillation, excited in the plasma surrounding the antenna in response to the transmission of ionospheric sounding pulses. Enhanced accuracy is achieved in higher-density ( n e  〉 150 cm − 3 ) plasmas, when MARSIS AIS is able to directly measure the fundamental frequency of the local plasma oscillation. To demonstrate the usefulness of this data set, the derived plasma densities are binned by altitude and solar zenith angle (SZA) in regions over weak ( nT) and strong ( nT) crustal magnetic fields, and we find clear and consistent evidence for a significant asymmetry between these two regions. We show that within the ~300-1200 km altitude range sampled, the median plasma density is substantially higher on the dayside in regions of relatively stronger crustal fields than under equivalent illuminations in regions of relatively weaker crustal fields. Conversely, on the nightside, median plasma densities are found to be higher in regions of relatively weaker crustal fields. We suggest that the observed asymmetry arises as a result of the modulation of the efficiency of plasma transport processes by the irregular crustal fields, and the generally horizontal draped interplanetary magnetic field (IMF).

Description: [1]  Global models of the Van Allen radiation belts usually include resonant wave-particle interactions as a diffusion process, but there is a large uncertainty over the diffusion rates. Here we present a new diffusion matrix for whistler mode chorus waves that can be used in such models. Data from seven satellites are used to construct 3,536 power spectra for upper and lower band chorus for 1.5 ≤  L ∗  ≤ 10, MLT, magnetic latitude 0 o  ≤ | λ m | ≤ 60 o and five levels of K p . Five density models are also constructed from the data. Gaussian functions are fitted to the spectra and capture typically 90% of the wave power. The frequency maxima of the power spectra vary with L ∗ and are typically lower than that used previously. Lower band chorus diffusion increases with geomagnetic activity and is largest between 21:00 and 09:00 MLT. Energy diffusion extends to a few MeV at large pitch angles 〉 60 o and at high energies exceeds pitch angle diffusion at the loss cone. Most electron diffusion occurs close to the geomagnetic equator (〈 12 o ). Pitch angle diffusion rates for lower band chorus increase with L ∗ and are significant at L ∗  = 8 even for low levels of geomagnetic activitywhile upper band chorus is restricted to mainly L ∗  〈 6. The combined drift and bounce averaged diffusion rates for upper and lower band chorus extend from a few keV near the loss cone up to several MeV at large pitch angles indicating loss at low energies and net acceleration at high energies.

Description: [1]  We report that radio science (RS) experiment onboard Mars Express (MEX) have observed three plasma layers in the nighttime ionosphere of Mars at altitudes ~80-100 km, ~120 km and ~160 km, which are reproduced by model calculation due to impact of meteoroid, solar wind proton and electron respectively. The densities of 21 ions (Mg + , Fe + , Si + , MgO + , MgCO 2 + , MgO 2 + , MgN 2 + , FeO + , FeO 2 + , FeN 2 + , FeCO 2 + , SiO + , SiCO 2 + , SiN 2 + , SiO 2 + , CO 2 + , N 2 + , O + , O 2 + , CO + , and NO + ) have been computed between altitude 50 km and 200 km. The model shows that all atmospheric ions (CO 2 + , N 2 + , O + , CO + , O 2 + and NO + ) are produced above 100 km due to solar wind electron and proton impact ionizations. The metallic ions are formed between 50 km and 100 km due to ablation of micrometeoroids. It is found that mass ~3.0 × 10 -4  g of incoming meteoroid is sufficient for meteor ablation and its characteristic flux ~4.0 × 10 -15  cm -2  s -1 could produce the nighttime metallic layer observed by MEX. The calculated electron densities are also compared with the occultation measurements made by Mars 4/5 in the nighttime ionosphere of Mars.

Description: [1]  Ultra Low Frequency (ULF) waves transfer energy in the Earth's magnetosphere through a variety of mechanisms that impact the Earth's ionosphere, radiation belts, and other plasma populations. Measurements of the electromagnetic portion of the energy transfer rate are an important source of information for assessing the importance of ULF waves relative to other energy transfer mechanisms as well as a diagnostic for studying the behavior of ULF waves. Using THEMIS satellite data, we examine the time averaged electromagnetic energy transfer rate, or Poynting vector, as a function of frequency and region of the magnetosphere; for this study, we focus on the direction and rate of energy transfer relative to the background magnetic field, comparing perpendicular and parallel transfer rates. This study extends earlier studies of the ULF wave Poynting vector that focused on narrower frequency ranges or specific regions of the magnetosphere; here, we consider the 3–50 mHz frequency range, all local time sectors, radial distances from 3 to 13 Re, and magnetic latitudes close to the equatorial plane. We measure time averaged Poynting vectors that range from  10 − 11 to 10 − 5 W / m 2 , with larger Poynting vector magnitudes occurring at largerradial distances and smaller frequencies. In every spatial region and frequency we examined, we found a large degree of scatter in both the Poynting vector magnitude and direction. The Poynting vector tends to be anisotropic at all frequencies,with more energy transferred along rather than across the background magnetic field. This preference for parallel energy transfer near the magnetic equator suggests that Joule dissipation in the ionosphere and the acceleration of auroral electrons are the largest sinks of ULF wave energy in the magnetosphere.

Description: [1]  Previous statistical studies have found a close relationship between high speed streams and high latitude geomagnetic activity. The speed by itself would increase the geoeffectivity of the solar wind. But it is also believed that pure Alfvénic fluctuations, often found in the trailing part of the streams, play a role in the solar wind driving of geomagnetic activity by amplifying the north-south component of the magnetic field ( B Z ), and thereby the dayside reconnection electric field. By automatically identifying slow and fast solar wind streams and by analysing them for more than one solar cycle, we aimed to study the relation between speed, Alfvénicity and B Z in the solar wind. It was found out that streams whose trailing parts are dominated by pure Alfvénic fluctuations, are the most geoeffective streams, on average. However, it is not the pure Alfvénic fluctuations themselves which cause the streams to be more geoeffective. There is only a variation of about 10 % in B Z due to the Alfvénicity of the fluctuations. Instead the streams are more geoeffective because the pure Alfvénic fluctuations tend to occur during high solar wind speed and strong IMF. There is a substantial variation within the solar cycle of how Alfvénic the solar wind streams are, and years with many extremely Alfvénic streams tends to have more days with moderately large geoeffectivity. The list of solar wind streams is included as extra material to this paper.

Description: [1]  Data Interpolating Empirical Orthogonal Functions (DINEOFs) are a data-based method for determining a few orthogonal basis functions that optimally reproduce a given data set. This technique is applied to meridional drift measurements performed by the Coupled Ion Neutral Dynamics Investigation (CINDI) onboard the Communication/Navigation Outage Forecasting System (C/NOFS) as well as electron density profiles derived from GPS Radio Occulations (RO) performed by COSMIC. The low densities of the equatorial ionosphere spanning 2009 - 2010 restricted quality drift measurements by CINDI to altitudes near perigee, limiting the local time coverage of measurements. Full local time descriptions may be obtained as perigee moves through all local times though this requires a minimum 67 day season. To increase the data coverage of the ionosphere CINDI data is supplemented with COSMIC GPS RO data. DINEOFs are applied to median meridional drift measurements as well as COSMIC measurements spanning 2009-10 and are used to make a best estimate of the equatorial ionosphere at locations not observed. The scattered distribution of COSMIC profiles as well as the physical relationship between meridional ion drifts and the distribution of density with altitude improve thequality of the reconstructions compared to using CINDI alone. The DINEOF reconstructions demonstrate that the annual anomomly of reduced ionospheric densities in June compared to December measured by COSMIC is coincident with a change in the meridional ion drifts at the geomagnetic equator measured by CINDI.

Description: [1]  On 5 September 2012, a large thrust earthquake (M w 7.6) ruptured a densely-instrumented seismic gap on the shallow-dipping plate boundary beneath the Nicoya Peninsula, Costa Rica. Ground motion recordings directly above the rupture zone provide a unique opportunity to study the detailed source process of a large shallow megathrust earthquake using very nearby land observations. Hypocenter relocation using local seismic network data indicates that the event initiated with small emergent seismic waves from a hypocenter ~10 km offshore, 13 km deep on the megathrust. A joint finite-fault inversion using high-rate GPS, strong-motion ground velocity recordings, GPS static offsets, and teleseismic P waves reveals that the primary slip zone (slip 〉 1 m) is located beneath the peninsula. The rupture propagated down-dip from the hypocenter with a rupture velocity of ~3.0 km/s. The primary slip zone extends ~70 km along strike and ~30 km along dip, with an average slip of ~2 m. The associated static stress drop is ~3 MPa. The seismic moment is 3.5 x 10 20  Nm, giving M w  = 7.6. The co-seismic large-slip patch directly overlaps an onshore inter-seismic locked region indicated by geodetic observations, and extends down-dip to the intersection with the upper plate Moho. At deeper depths, below the upper plate Moho, seismic tremor and low frequency earthquakes have been observed. Most tremor locates in adjacent areas of the megathrust that have little co-seismic slip; a region of prior slow slip deformation to the southeast also has no significant co-seismic slip or aftershocks. An offshore locked patch indicated by geodetic observations does not appear to have experienced co-seismic slip, and aftershocks do not overlap this region, allowing the potential for a comparable size rupture offshore in the future.

Description: [1]  Seismic shear waves emitted by earthquakes can be modeled as plane (transverse) waves. When entering an anisotropic medium they can be split into two orthogonal components moving at different speeds. This splitting occurs along an axis, the fast polarization, that is determined by geologic conditions. We present here a comprehensive analysis of the Silver and Chan (1991) method, used to obtain shear wave splitting parameters, comprising theoretical derivations and statistical tests of the assumptions used to construct the standard errors. We find discrepancies in the derivations of equations in their article, with the most important being a mistake in how the standard errors are calculated. Our simulations suggest that the degrees of freedom are being overestimated by this method and consequently the standard errors are too small. Using a set of S waveforms from very similar shallow earthquakes on Reunion Island, we perform a statistical analysis on the noise of these replicates and find that the assumption of Gaussian noise does not hold. Further, the properties of background noise differ substantially from the noise obtained from the shear wave splitting analysis. However, we find that the standard errors for the fast polarization are comparable to the spread in the fast polarization parameters. Delay time errors appear to be comparable to delay time estimates once cycle skipping is accounted for. Future work using synthetic seismograms with simulated noise should be conducted to confirm this is the case for earthquakes in general.

Description: [1]  The sliver strike-slip Great Sumatra Fault (GSF) traverses mainland Sumatra from the Sunda Strait in the southeast to Banda Aceh in the northwest, and defines the present day plate boundary between the Sunda Plate in the north and the Burmese Sliver Plate in the south. It has been well studied on mainland Sumatra but poorly north of Banda Aceh in the Andaman Sea. Here we present deep seismic reflection images along the northward extension of the GSF over 700 km until it joins the Andaman Sea Spreading Centre and we interpret these images in the light of earthquake, gravity, and bathymetry data. We find that the GSF has two strands between Banda Aceh and Nicobar Island: a transpression in the south and a deep narrow active rift system in the north dotted with volcanoes in the center, suggesting that the volcanic arc is coincident with rifting. Further north of Nicobar Island, an active strike-slip fault, the Andaman-Nicobar Fault, cuts through a rifted deep basin until its intersection with the Andaman Sea Spreading Centre. The volcanic arc lies just east of the rift basin. The western margin of this basin seems to be a rifted continental margin, tilted westward, and flooring the Andaman-Nicobar forearc basin. The Andaman-Nicobar forearc basin is bounded in the west by backthrusts similar to the West Andaman and Mentawai faults. The cluster of seismicity after the 2004 great Andaman-Sumatra earthquake just north of Nicobar Island coincides with the intersection of two strike-slip fault systems.

Description: [1]  The Japan Tohoku-Oki earthquake (9.0 Mw) of 11 March 2011 has left signatures in the Earth's gravity field that are detectable by data of the GRACE mission. Because ESA's satellite gravity mission GOCE – launched in 2009 – aims at high spatial resolution, its measurements could complement the GRACE information on coseismic gravity changes, although time-variable gravity was not foreseen as goal of the GOCE mission. We modeled the coseismic earthquake geoid signal and converted this signal to vertical gravity gradients at GOCE satellite altitude. We combined the single gradient observations in a novel way reducing the noise level, required to detect the coseismic gravity change, subtracted a global gravity model, and applied tailored outlier detection to the resulting gradient residuals. Furthermore, the measured gradients were along-track filtered using different gradient bandwidths where in the space domain Gaussian smoothing has been applied. One year periods before and after earthquake occurrence have been compared with the modeled gradients. The comparison reveals that the earthquake signal is well above the accuracy of the vertical gravity gradients at orbital height. Moreover, the obtained signal from GOCE shows a 1.3 times higher amplitude compared with the modeled signal. Besides the statistical significance of the obtained signal, it has a high spatial correlation of ~0.7 with the forward modeled signal. We conclude therefore that the coseismic gravity change of the Japan Tohoku-Oki earthquake left a statistically significant signal in the GOCE measured gravity gradients.

Description: [1]  The Ninetyeast Ridge (NER), one of the longest linear volcanic features on the Earth, offers an excellent opportunity of understanding the isostatic response to the interactions of mantle plume with the migrating mid-ocean ridge. Bathymetry, geoid and gravity (ship-borne and satellite) data along 72 closely spaced transects and 17 overlapping grids on the NER are analyzed and modeled to determine the effective elastic thickness ( Te ) beneath the entire ridge. The results of 2-D and 3-D flexural modeling of the NER show large spatial variations in Te values ranging from 4 to 35 km, suggesting that the ridge was compensated along its length by different isostatic mechanisms. The southern (south of 22°S latitude) and northern (north of 2°N latitude) parts of the NER have Te values of 〉10 and 〉23 km, respectively, revealing that the southern part was emplaced on a lithosphere of intermediate strength possibly on flank of the Indian plate, whereas the northern part was emplaced in an intraplate setting. In contrast, in the central part of the NER (between latitudes 22°S and 2°N), highly variable Te values (4–22 km) are estimated. The scattered Te values in the central NER suggest that this part may have evolved due to the occurrence of more frequent ridge jumps caused by the interaction of Kerguelen hot spot with rapid northward migration of the Wharton spreading ridge. Residual Mantle Bouguer Anomaly (RMBA) map of the NER and adjacent basins reveals that the entire length of the NER is associated with a significant negative anomaly up to 200 mGal, indicating the presence of thickened crust or less dense mantle beneath the ridge. 3-D crustal thickness map of the NER, generated by inversion of the RMBA data, shows a thick crust ranging from 15 to 19 km. The present study clearly shows that NER possesses a highly segmented isostatic pattern with the occurrence of sub-crustal underplating or sub-surface loading.

Description: [1]  We conducted deep-sea magnetic measurements using autonomous underwater vehicles in the Bayonnaise knoll caldera, the Izu-Ogasawara island arc, which hosts the large Hakurei hydrothermal field. We improved the conventional correction method applied for removing the effect of vehicle magnetization, thus greatly enhancing the precision of the resulting vector anomalies. The magnetization distribution obtained from the vector anomaly data shows an ∼ 2-km-wide belt of high magnetization, trending NNW–SSE going through the caldera, and a low magnetization zone ∼ 300 m by ∼ 500 m in area, extending over the Hakurei site. Comparison between the results obtained using the vector anomaly and the total intensity anomaly shows that the magnetic field is determined more accurately, especially in areas of sparse data distribution, when the vector anomaly rather than the total intensity anomaly is used. We suggest a geologically motivated model that basaltic volcanism associated with the backarc rifting occurred after the formation of the caldera, resulting in the formation of the high magnetization belt underneath the silicic caldera. The Hakurei hydrothermal field lies in the intersection of the basaltic volcanism belt and the caldera wall fault, suggesting a mechanism that hot water generated by the heat of the volcanic activity has been spouting out through the caldera wall fault. The deposit apparently extends beyond the low magnetization zone, climbing up the caldera wall. This may indicate that hot water rising from the deep through the alteration zone is transported laterally when it comes near the seafloor along fissures and fractures in the caldera wall.

Description: [1]  Using multiple ScS reverberations we examine mantle reflectivity structure beneath northeast China and the northwest Pacific. We find several upper mantle discontinuities, including a melt layer with a mean thickness of 64 km atop the 410-km discontinuity, present on both sides of the subducting slab near the Nankai trench. The transition zone contains a split 520-km discontinuity in several paths, and tomographic images show stagnant slabs at this depth. We believe this may be slab-related based on experimental work (Saikia, A., Frost, D. J., Rubie, D. C., 2008. Splitting of the 520-kilometer seismic discontinuity and chemical heterogeneity in the mantle. Science 319 (5869), 1515–1518). A negative reflector is found in one path beneath the northeast China craton at a depth of 598 km. Mid-mantle reflectors are found in all of our paths and are present throughout a wide depth range (~750 – 1600 km).

Description: [1]  Three-dimensional P - and S -wave velocity (V P , V S ) models and high-resolution earthquake relocations are determined for the New Madrid Seismic Zone using double-difference local earthquake tomography. The data set consists of arrival times and differential times recorded by the Cooperative New Madrid Seismic Network (CNMSN) from 2000-2007 and the 1989-1992 Portable Array Network and Data Acquisition deployment. Waveform cross-correlation derived differential times for the CNMSN data are also incorporated. The velocity solutions are compatible with previous solutions centered on the active arms of seismicity and cover a broader area including mafic intrusions along the margin of the Reelfoot Rift. Major features include elevated V P and V S associated with the mafic plutons and reduced V P and V S along and southeast of the Axial fault (AF), a major arm of seismicity trending along the rift axis. Low V P extends to a depth of at least 20 km along the portion of the AF that extends south of the Missouri bootheel. A locally high V P /V S anomaly imaged along the central portion of the Reelfoot fault is spatially correlated with a significant change in fault trend and is interpreted as a region containing high pore pressure and/or water-filled microcracks.

Description: [1]  We present a catalog of InSAR constraints on deformation that occurred during earthquake sequences in southern Iran between 1992-2011, and explore the implications on the accommodation of large-scale continental convergence between Saudi Arabia and Eurasia within the Zagros Mountains. The Zagros Mountains, a salt-laden fold-and-thrust-belt involving ~10 km of sedimentary rocks overlying Precambrian basement rocks, have formed as a result of ongoing continental collision since 10-20 Ma that is currently occurring at a rate of ~3 cm/yr. We first demonstrate that there is a biased misfit in earthquake locations in global catalogs that likely results from neglect of 3D velocity structure. Previous work involving two M ~ 6 earthquakes with well-recorded aftershocks has shown that the deformation observed with InSAR may represent triggered slip on faults much shallower than the primary earthquake, which likely occurred within the basement rocks (〉10 km depth). We explore the hypothesis that most of the deformation observed with InSAR spanning earthquake sequences is also due to shallow, triggered slip above a deeper earthquake, effectively doubling the moment release for each event. We quantify the effects that this extra moment release would have on the discrepancy between seismically and geodetically constrained moment rates in the region, finding that even with the extra triggered fault slip, significant aseismic deformation during the interseismic period is necessary to fully explain the convergence between Eurasia and Saudi Arabia.

Description: [1]  It is well accepted that the propagation of electromagnetic ion cyclotron (EMIC) waves are bidirectional near their source regions and unidirectional when away from these regions. The generally believed source region for EMIC waves is around the magnetic equatorial plane. Here we describe a series of EMIC waves in the Pc1 (0.2-5 Hz) frequency band above the local He + cyclotron frequency observed in situ by all four Cluster spacecraft on 9 April 2005 at mid-magnetic latitudes (MLAT = ~33°-49°) with L = 10.7-11.5 on the dayside (MLT = 10.3-10.4). A Poynting vector spectrum shows that the wave packets consist of multiple groups of packets propagating bidirectionally, rather than unidirectionally, away from the equator, while the local plasma conditions indicate that the spacecraft are entering into a region sufficient for local wave excitation. One possible interpretation is that, while part of the observed waves are inside their source region, the others are either close enough to the source region, or mixed with the wave packets from multiple source regions at different latitudes.

Description: [1]  The identification and evaluation of trigger mechanisms for volcano flank instabilities and/or collapse represent a key issue for risk assessment in densely populated volcanic areas, as well as in long distance settings, particularly in case of island or coastal volcanoes. Here, we address quantitatively the effects of external (seismic) and inner (magmatic) forcing on the stress-strain state associated to flank instabilities at Mt. Etna (Sicily, southern Italy) by means of a 2-D Finite-Difference-Method numerical modelling. Modelled seismic actions include strong near-field, strong far-field and low-magnitude near-field earthquakes. Magmatic actions consider the inner presssure changes induced by energetic lava fountains in the summit crater area, sub-vertical and oblique dyke ascent below the summit area. Model results are validated in light of available monitoring data and recent eruptive activity. Numerical results show that the main strain effects are produced by high-magnitude near-field earthquakes (expected return time of ~10 3  yrs), and by vertical uprise of a magma dyke below the volcano summit area. Maximum displacements in the order of tens of centimetres may involve the summit area, up to some 10 6  m 3 /m over some km laterally. Stress releases up to 10 7  Pa may affect a limited portion of the magmatic conduit, thus favouring major effusive flank eruptions. Major catastrophic events, such as volcano flank collapse, should not be expected by applying, either individually or combined, the aforementioned actions.

Description: [1]  Seismicity closely related to hydrological impacts has been observed in several locations worldwide; particularly in intraplate areas where tectonic stressing rates are small. The triggering mechanism is usually explained by a poroelastic response of the seismogenic crust to surface water flux, leading to pore pressure changes at depth. To explain the earthquake triggering in response of those small stress changes, however, the crust has to be near a critical state in which other transient processes might be significant. One of the prominent examples is the Mt. Hochstaufen in SW Germany, where seismicity is known to vary seasonally. A previous analysis showed that the seismicity in 2002 was highly correlated with model forecasts based on fluid diffusion and rate- and state-dependent frictional nucleation. Here we revisit this case by accounting additionally for poroelastic effects, as well as for thermoelastic and tidal stresses. We also test whether the model can explain the observations of the subsequent eight years between 2003 and 2010. Our analysis confirms that rainfall is the dominant driving force in this region. The model not only fits the year 2002 activity very well, but provides with the same parameters a reasonable fit to the subsequent period, with a probability gain of about 4 per event in comparison to a time-independent Poisson model.

Description: ABSTRACT [1]  The absolute magnitude of stress in the crust and the shear strength of faults are poorly known, yet fundamental quantities, in lithospheric dynamics. While stress magnitude cannot be measured directly, deviatoric stress state can be inferred indirectly from focal mechanism solutions collected before and after an earthquake. We extend a standard stress inversion for normalized stresses to invert for the 3D spatial distribution of absolute deviatoric stress and variation of fault strength with depth using focal mechanism solutions and coseismic stress changes produced by large earthquakes. We apply the method to the 2011 M9 Tohoku-oki, Japan earthquake. The northern Japan forearc crust between 5 and 15 km depth appears to be weak with fault strength of 40–90 MPa, consistent with a coefficient of friction of 0.2-0.5. The M9 Tohoku-oki coseismic stress change was large enough, relative to the ambient stress, to rotate the principal stress directions typically ~20° in the upper 20 km of the crust. The data from Japan require a heterogeneous ambient deviatoric stress field with short wavelength (~20-50 km) fluctions in principal stress orientations.

Description: [1]  Using mesospheric temperature information obtained by Rayleigh lidar observations over Gadanki (13.5°N, 79.2°E) and temperature, chemical constituents and their exothermic chemical reaction heating rates obtained by Sounding of Atmosphere by Broadband Emission Radiometry (SABER) instrument onboard Thermosphere Ionosphere Mesosphere Energetics and Dynamics (TIMED) satellite over and near to Gadanki location, it is demonstrated that, the dominant role of chemical heating for the occurrence of a few large mesospheric inversion layer (MIL) events during January-February 2011. The heating rates due to exothermic chemical reactions among atomic oxygen (O), atomic hydrogen (H), molecular oxygen (O 2 ), ozone (O 3 ) and hydroxyl (OH) species are the key factors involved in the MIL formation. As relatively larger (lower) ozone mixing ratios are observed at the top (bottom) of the inversion layers, large heating rate (~ 15 K/day) is found to occur between 80 and 85 km due to the exothermic reaction between atomic hydrogen and ozone H + O 3  → OH + O 2 , which dominates all other important exothermic chemical reactions during these MIL events and the large volume emission rates of energy of O 2 , OH are found to be mostly due to the reaction between H and O 3. It is also demonstrated that dynamics plays a negligible role in the formation of these MILs as conditions for the convective and dynamic instability are found to be not satisfied during these MIL events.

Description: [1]  Bursty bulk flows (BBFs) play an important role for the mass, energy and magnetic flux transport in the plasma sheet, and the flow pattern in and around a BBF has important consequences for the localized energy conversion between the electromagnetic and plasma mechanical energy forms. The plasma flow signature in and around BBFs is often rather complicated. Return flows and plasma vortices are expected to exist at the flanks of the main flow channel, especially near the inner plasma sheet boundary, but also further down-tail. A dipolarization front (DF) is often observed at the leading edge of a BBF, and a flux pileup region (FPR) behind the DF. Here we present Cluster data of three FPRs associated with vortex flows observed in the mid-tail plasma sheet on August 15, 2001. According to the principles of Fu et al . [2011] and Fu et al . [2012c], two of the FPRs are considered to be in an early stage of evolution (growing FPRs). The third FPR is in a later stage of evolution (decaying FPR). For the first time, the detailed energy conversion properties during various stages of the FPR evolution have been measured. We show that the later stage FPR has a more complex vortex pattern than the two earlier stage FPRs. The two early stage FPR correspond to generators, E  ⋅  J  〈 0, while the later stage FPR only shows weak generator characteristics and is instead dominated by load signatures at the DF, E  ⋅  J  〉 0. Moreover, to our knowledge this is one of the first times BBF related plasma vortices have been observed to propagate over the spacecraft in the mid-tail plasma sheet at geocentric distances of about 18 R E . Our observations are compared to recent simulation results [ Wiltberger et al ., 2000; Birn et al ., 2011] and previous observations [ Fu et al ., 2011, 2012c; Pang et al ., 2012].

Description: [1]  We review marine heat flow data along the Nankai Trough and show that observations 〉 30 km seaward of the deformation front are 20% below conductive predictions (129–94 mW m -2 ) but consistent with the global heat flow average for oceanic crust of the same age (16-28 Ma). Heat flow values  〈  30 km seaward of the deformation front are generally 20% higher than conductive predictions. This heat flow pattern is consistent with the advection of heat by fluid flow in the subducting oceanic crust and explains both the high heat flux in the vicinity of the trench, 〉 200 and 〉 140 mW m -2 , and steep landward declines to values of approximately 60 mW m -2 over distances of 65 and 50 km along the Muroto and Kumano transects, respectively. Along the Ashizuri transect the lack of heat flow data preclude a definitive interpretation. We conclude that fluid flow in the subducting oceanic crust leads to temperatures that are generally 25 ° C higher near the toe of the margin wedge and 50 - 100 ° C lower near the downdip limit of the seismogenic zone than estimated by purely conductive models.

Description: [1]  We use the measurements by the Vector Electric Field Investigation (VEFI) on board the C/NOFS satellite to study the local time, seasonal, and longitudinal dependent climatology of equatorial F region zonal plasma drifts during 2008-2011. These drifts are essentially westward during the day and eastward at night. Over Jicamarca, the satellite measurements are in good agreement with incoherent scatter radar drifts. Our data show strong longitudinal variations, particularly in the South American sector during the solstices. The equinoctial data exhibit short-lived and largely enhanced westward drifts near sunrise, and wave-4 structures from theearly afternoon to late night equinoctial periods. The nighttime eastward drifts are largest near the western American sector at all seasons. The June solstice post-midnight eastward drifts decrease sharply at longitudes of about 310˚ and have much smaller values in the entire eastern hemisphere. We also briefly discuss the relationship between the longitude dependent vertical and zonal plasma drifts.

Description: [1]  We investigate whether predictions of mantle structure from tectonic reconstructions are in agreement with a detailed tomographic image of seismic P-wave velocity structure under the Caribbean region. In the upper mantle, positive seismic anomalies are imaged under the Lesser Antilles and Puerto Rico. These anomalies are interpreted as remnants of Atlantic lithosphere subduction and confirm tectonic reconstructions that suggest at least 1100 km of convergence at the Lesser Antilles island arc during the past ~45 Myr. The imaged Lesser-Antilles slab consists of a northern and southern anomaly, separated by a low velocity anomaly across most of the upper mantle, which we interpret as the subducted North America-South America plate boundary. The southern edge of the imaged Lesser Antilles slab agrees with vertical tearing of South America lithosphere. The northern Lesser Antilles slab is continuous with the Puerto Rico slab along the northeastern plate boundary. This results in an amphitheatre-shaped slab and it is interpreted as westward subducting North America lithosphere that remained attached to the surface along the northeastern boundary of the Caribbean plate. At the Muertos Trough, however, material is imaged until a depth of only 100 km, suggesting a small amount of subduction. The location and length of the imaged South Caribbean slab agrees with proposed subduction of Caribbean lithosphere under the northern South America plate. An anomaly related to proposed Oligocene subduction at the Nicaragua rise is absent in the tomographic model. Beneath Panama, a subduction window exists across the upper mantle, which is related to the cessation of subduction of the Nazca plate under Panama since 9.5 Ma and possibly the preceding subduction of the extinct Cocos-Nazca spreading center. In the lower mantle two large anomaly patterns are imaged. The westernmost anomaly agrees with the subduction of Farallon lithosphere. The second lower mantle anomaly is found east of the Farallon anomaly and is interpreted as a remnant of the late Mesozoic subduction of North and South America oceanic lithosphere at the Greater Antilles, Aves ridge and Leeward Antilles. The imaged mantle structure does not allow us to discriminate between an ‘Intra-Americas’ origin and a ‘Pacific origin’ of the Caribbean plate.

Description: [1]  Long-term variation in the geomagnetic Sq field and the cause of the variation were examined. 　The amplitude of the geomagnetic Y component (Sq(Y)) in equinox was averaged for each year, and adopted as a proxy of the Sq field. Sq(Y) was combined with the ionospheric conductivity estimated by the IRI model to determine the dynamo electric field and neutral wind velocity by using the geomagnetic main field strength. It was found that the solar activity dependence of the Sq field could be almost completely attributed to the conductivity variation, and neutral winds tend to decrease when the solar activity increases Although the long-term variation in the dynamo field differed among observatories, these differences were mostly attributed to the locality of the geomagnetic secular variation, whereas the variations in neutral wind amplitude were nearly the same in all regions. On the other hand, no clear long-term variation in neutral wind was detected other than that by solar activity.

Description: [1]  We have measured interseismic deformation across the Ashkabad strike-slip fault using 13 Envisat interferograms covering a total effective timespan of ~30 years. Atmospheric contributions to phase delay are significant and variable due to the close proximity of the Caspian Sea. In order to retrieve the pattern of strain accumulation, we show it is necessary to use data from Envisat's Medium Resolution Imaging Spectrometer (MERIS) instrument, as well numerical weather model outputs from the European Centre for Medium-Range Weather Forecasting (ECMWF), to correct interferograms for differences in water vapour and atmospheric pressure respectively. This has enabled us to robustly estimate the slip rate and locking depth for the Ashkabad fault using a simple elastic dislocation model. Our data are consistent with a slip rate of 5–12 mm/yr below a locking depth of 5.5–17 km for the Ashkabad fault, and synthetic tests support the magnitude of the uncertainties on these estimates. Our estimate of slip rate is 1.25–6 times higher than some previous geodetic estimates, with implications for both seismic hazard and regional tectonics, in particular supporting fast relative motion between the South Caspian Block and Eurasia. This result reinforces the importance of correcting for atmospheric contributions to interferometric phase for small strain measurements. We also attempt to validate a recent method for atmospheric correction based on ECMWF ERA-Interim model outputs alone and find that this technique does not work satisfactorily for this region when compared to the independent MERIS estimates.

Description: [1]  The seafloor within the Perth Abyssal Plain (PAP), offshore Western Australia, is the only section of crust that directly records the early spreading history between India and Australia during the Mesozoic breakup of Gondwana. However, this early spreading has been poorly constrained due to an absence of data, including marine magnetic anomalies and data constraining the crustal nature of key tectonic features. Here, we present new magnetic anomaly data from the PAP that shows that the crust in the western part of the basin was part of the Indian Plate – the conjugate flank to the oceanic crust immediately offshore the Perth margin, Australia. We identify a sequence of M2 and older anomalies in the west PAP within crust that initially moved with the Indian Plate, formed at intermediate half-spreading rates (35 mm/yr) consistent with the conjugate sequence on the Australian Plate. More speculatively, we reinterpret the youngest anomalies in the east PAP, finding that the M0-age crust initially formed on the Indian Plate was transferred to the Australian Plate by a westward jump or propagation of the spreading ridge shortly after M0 time. Samples dredged from the Gulden Draak and Batavia Knolls (at the western edge of the PAP) reveal that these bathymetric features are continental fragments rather than igneous plateaus related to Broken Ridge. These microcontinents rifted away from Australia with Greater India during initial breakup at ~130 Ma, then rifted from India following the cessation of spreading in the PAP (~101-103 Ma).

Description: [1]  We develop a three-step Maximum-A-Posteriori probability (MAP) method for coseismic rupture inversion, which aims at maximizing the a posterior probability density function (PDF) of elastic deformation solutions of earthquake rupture. The method originates from the Fully Bayesian Inversion (FBI) and Mixed linear-nonlinear Bayesian inversion (MBI) methods, shares the same posterior PDF with them, while overcoming difficulties with convergence when large numbers of low-quality data are used and greatly improving the convergence rate using optimization procedures. A highly-efficient global optimization algorithm, Adaptive Simulated Annealing (ASA), is used to search for the maximum of a posterior PDF (" mode " in statistics) in the first step. The second step inversion approaches the " true" solution further using the Monte Carlo Inversion (MCI) technique with positivity constraints, with all parameters obtained from step one as the initial solution. Then slip artifacts are eliminated from slip models in the third step using the same procedure of the second step, with fixed fault geometry parameters. [2]  We first design a fault model with 45°-dip angle and oblique slip, and produce corresponding synthetic InSAR datasets to validate the reliability and efficiency of the new method. We then apply this method to InSAR data inversion for the coseismic slip-distribution of the April 14, 2010 Mw 6.9 Yushu, China earthquake. Our preferred slip model is composed of three segments with most of the slip occurring within 15 km depth and the maximum slip reaches 1.38 m at the surface. The seismic moment released is estimated to be 2.32e + 19 Nm, consistent with the seismic estimate of 2.50e + 19 Nm.

Description: [1]  We analyze data recorded from October 2010 to September 2011, during the ascending phase of the 24 th solar cycle, from an AIS-INGV (Advanced Ionospheric Sounder - Istituto Nazionale di Geofisica e Vulcanologia) ionosonde and a GISTM (GPS Ionospheric Scintillation and TEC Monitor) scintillation receiver, co-located at low-latitude in the Southern American longitudinal sector (Tucumán, 26.9°S, 294.6°E, mag. Lat. 15.5°S, Argentina). The site offers the opportunity to perform spread-F and GPS scintillation statistics of occurrence under the southern crest of the equatorial ionospheric anomaly. Spread-F signatures, classified into four types (strong range spread-F (SSF), range spread-F (RSF), frequency spread-F (FSF) and mixed spread-F (MSF)), the phase and amplitude scintillation index ( σ Φ and S 4 , respectively), the total electron content ( TEC ) and the rate of TEC ( ROT ) parameter, marker of the TEC gradients, that can cause scintillations, are considered. The seasonal behavior results as follows: the occurrence of all four types of spread-F is higher in summer and lower in winter, while the occurrence of scintillations peaks at equinoxes in the post-sunset sector and shows a minimum in winter. The correspondence between SSF and scintillations seems to be systematic, and a possible correlation between S 4 and FSF peaks is envisaged at the terminator. The investigation focused also on two particular periods, from 12 to 16 March 2011, and from 23 to 29 September 2011, both characterized by the simultaneous presence of SSF signatures and scintillation phenomena allowing to discuss the role of Travelling Ionospheric Disturbances as a strong candidate causing ionospheric irregularities.

Description: ABSTRACT [1]  Neutrons contribute a significant radiation dose at commercial passenger airplane altitudes. With cosmic ray energies 〉 1  GeV , these effects could in principle be propagated to ground level. Under current conditions, the cosmic ray spectrum incident on the Earth is dominated by particles with energies 〈 1  GeV . Astrophysical shocks from events such as supernovae accelerate high energy cosmic rays (HECRs) well above this range. The Earth is likely episodically exposed to a greatly increased HECR flux from such events. Solar events of smaller energies are much more common and short lived, but still remain a topic of interest due to the ground level enhancements (GLEs) they produce. The air showers produced by cosmic rays (CRs) ionize the atmosphere and produce harmful secondary particles such as muons and neutrons. Although the secondary spectra from current day terrestrial cosmic ray flux are well known, this is not true for spectra produced by many astrophysical events. This work shows the results of Monte Carlo simulations quantifying the neutron flux due to CRs at various primary energies and altitudes. We provide here lookup tables (described later) that can be used to determine neutron fluxes from proton primaries with total energies 1  MeV – 1 PeV . By convolution, one can compute the neutron flux for any arbitrary CR spectrum. This contrasts with all other similar work, which is spectrum dependent. Our results demonstrate the difficulty in deducing the nature of primaries from the spectrum of ground level neutron enhancements.

Description: [1]  Rayleigh lidar temperature observations over Gadanki show a few mesospheric inversion layer (MIL) events during 20-25 January 2007. The zonal mean removed SABER temperature shows warm anomalies around 50°E and 275°E indicating the presence of planetary wave of zonal wavenumber 2. The MIL amplitudes in SABER temperature averaged for 10°N-15°N and 70°E-90°E show a clear 2-day wave modulation during 20-28 January 2007. Prior to 20 January 2007, a strong 2-day wave (zonal wavenumber 2) is observed in the height region of 80-90 km and it gets largely suppressed during 20-26 January 2007 as the condition for vertical propagation is not favourable, though it prevails at lower heights. The 10-day mean zonal wind over Tiruneveli shows deceleration of eastward winds indicating the westward drag due to wave dissipation. The nightly mean MF radar observed zonal winds show the presence of alternating eastward and westward winds during the period of 20-26 January 2007. The two dimensional spectrum of Rayleigh lidar temperature observations available for the nights of 20, 22 and 24 January 2007 shows the presence of gravity wave activity with periods 18 min, 38 min, 38 min, and vertical wavelengths 6.4 km, 4.0 km, 6.4 km respectively. From the dispersion relation of gravity waves, it is inferred that, these waves are internal gravity waves rather than inertia gravity waves with the horizontal phase speeds of ~40 m/s, ~37 m/s and ~50 m/s respectively. Assuming the gravity waves are eastward propagating waves, they get absorbed only in the eastward local wind fields of the planetary wave thereby causing turbulence and eddy diffusion which can be inferred from the estimation of large drag force due to the breaking of gravity wave leading to the formation of large amplitude inversion events in alternate nights. The present study shows that, the mesospheric temperature inversion is caused mainly due to the gravity wave breaking and the inversion amplitude may get modulated by the interaction between gravity waves and planetary waves. The eddy diffusion associated with gravity wave drag may also cause suppression in the planetary wave activity.

Description: [1]  To evaluate whether ground magnetic responses traditionally viewed as signatures of substorm onset may be related to post-onset auroral streamers rather than to substorm onset seen in aurora, we have analyzed ground magnetic and all-sky image responses for 14 substorm onsets selected solely on imager coverage. We find evidence supporting this possibility, the response being approximately coincident for mid-latitude positive bay and Pi2 pulsation enhancements and for abrupt auroral zone H decreases at stations near a streamer. These are all signatures of current wedge formation, indicating that substorm current wedge formation may be more of a response to plasma sheet flow channels associated with post-onset streamers than to the process leading to auroral onset, and that the current wedge may develop via a series of distinctly narrower wedge-like structures. We also find evidence that periods of streamers can occur more than once during a substorm expansion phase, and can give magnetic signatures of multiple onsets even if there are no additional auroral onsets. Furthermore, the peak auroral zone H decreases are seen in association with streamers and at times varying from just a few minutes to well over a half hour after substorm auroral onset if there is a prolonged period of streamers. This indicates peak substorm AE and AL values are likely related to features of post-onset streamers and their locations relative to observing stations, and are not a measure of the strength of the substorm auroral onset processes. We discuss possible additional implications on substorm expansion phase development.

Description: [1]  We conducted total magnetic field and Bouguer gravity measurements to investigate the shallow structure beneath the summit caldera of Kīlauea Volcano, Hawai‘i. Two significant and distinctive magnetic anomalies were identified within the caldera. One is interpreted to be associated with a long-lived pre-historic eruptive centre, the Observatory vent, located ~1 km east of the Hawaiian Volcano Observatory. The second magnetic anomaly corresponds to a set of eruptive fissures that strike northeast from Halema‘uma‘u Crater, suggesting this is an important transport pathway for magma. The Bouguer gravity data were inverted to produce 3D models of density contrasts in the upper 2 km beneath Kīlauea. The models detect 3.0 km 3 of material, denser than 2800 kg m -3 , beneath the caldera that may represent an intrusive complex centred northeast of Halema‘uma‘u. Recent temporal gravity studies indicate continual addition of mass beneath the caldera during 1975–2008 centred west of Halema‘uma‘u and suggest this is due to filling of void space. The growth of a large intrusive complex, apparent cyclical caldera formation, and continual mass addition without inflation, however, can also be explained by extensional rifting caused by the continual southward movement of Kīlauea's unstable south flank.

Description: [1]  The 2006-2007 doublet of M W  〉 8 earthquakes in the Kuril subduction zone caused postseismic transient motion in the asthenosphere, which we observed on the Kuril GPS Array in 2007–2011. Here we show that the Maxwell asthenospheric viscosity that best fits the geodetic data increased by nearly an order of magnitude over the interval of four years, from 2 × 10 17 to 1 × 10 18  Pa s. These effective values of viscosity can be explained by a power-law rheology for which strain rate is proportional to stress raised to a power n  〉 1. The apparent change in viscosity can also be caused by other factors such as coupling between afterslip and viscoelastic flow. The open and intriguing question in connection with postseismic data after the Kuril earthquake doublet is the magnitude of the long-term asthenospheric viscosity, which shall be revealed by continued observations. An asthenosphere with viscosity of about 1 × 10 19  Pa s is favored by the postseismic deformation still observed several decades after the 1960 Chile and 1964 Alaska M W ~9 earthquakes. However, postseismic deformation associated with the 1952 southern Kamchatka M W ~9 earthquake currently is not observed in the northern Kurils, an indication that the long-term asthenospheric viscosity in the Kurils is lower than in Chile and Alaska.

Description: [1]  The region of central Chile offers a unique opportunity to study the links between the subducting Juan Fernandez Ridge, the flat slab, the Double Seismic Zone (DSZ) and the absence of modern volcanism. Here, we report the presence and characteristics of the first observed DSZ within the intermediate-depth Nazca slab using two temporary seismic catalogues (OVA99 and CHARSME). The lower plane of seismicity (LP) is located 20–25 km below the upper plane (UP), begins at 50 km depth and merges with the lower plane at 120 km depth, where the slab becomes horizontal. Focal mechanism analysis and stress tensor calculations indicate that the slab's state of stress is dominantly controlled by plate convergence and overriding crust thickness: Above 60–70 km depth, the slab is in horizontal compression, and below, it is in horizontal extension, parallel to plate convergence, which can be accounted for by vertical loading of the overriding lithosphere. Focal mechanisms below 60–70 km depth are strongly correlated with offshore outer rise bend faults, suggesting the reactivation of pre-existing faults below this depth. The large interplane distances for all Nazca DSZs can be related to the slab's unusually cold thermal structure with respect to its age. Since LPs globally seem to mimic mantle mineral dehydration paths, we suggest that fluid migration and dehydration embrittlement provide the mechanism necessary to weaken the rock and that the stress field determines the direction of rupture.

Description: [1]  We statistically examine the occurrence rate of electromagnetic ion cyclotron (EMIC) waves observed by AMPTE/CCE. We use the 8-Hz magnetic field dataset that covers the whole CCE mission period of nearly 4.5 years from Augutst 1984 to Janauary 1989, which is more than three times the period studied by Anderson et al. [1992] (~452 days). The large data volume allows us to evaluate the storm phase dependence of the spatial occurence pattern of EMIC waves. The major results of this study are summarized as follows. (a) The occurrence rate is below 5% on the nightside at all L. On the dayside the rate is 〈5% in the inner magnetosphere (L 〈 6), while it is higher than 5% in the outer magnetosphere (L ≥ 6), up to 25%. The highest rate appears in the afternoon sector. (b) The overall occurrence rate is higher for H-band events than He-band events, except for the opposite feature seen in the inner magnetosphere on the early afternoon-to-post midnight side (L 〈 6, 14 h 〈 MLT 〈 22 h). (c) H-band events occur frequently in the outer magnetosphere (L ≥ 7) in the afternoon sector, regardless of geomagnetic activity. Under quiet conditions, H-band events also occur in the outer magnetosphere on the morning side (4 h ≤ MLT 〈 8 h). (d) He-band events frequently occur in the inner magnetosphere (L 〈 7) on the pre-noon to dusk side (10 h ≤ MLT 〈 19 h) under disturbed conditions (Dst ≤ −50 nT). (e) The storm-time He-band waves are generated more frequently during the main phase than the recovery phase, with the main-phase wave excitation seen toward the afternoon-side outer magnetosphere (L 〉 7). The results indicate two independent major processes that cause EMIC wave excitation in the Earth's magnetosphere: one externally triggers H-band waves on the dayside, and the other internally excites He-band waves on the dusk to afternoon side. We suggest that the former is due to solar wind compression which leads to perpendicular adiabatic ion heating and in turn an increase in temperature anisotropy, and that the latter is caused by injections of new, highly energetic ion population from the plasma sheet, with its velocity distributions becoming pancake-like on the dusk-to-afternoon side. The frequent occurrence seen on the afternoon side, at a wide L range, and during the main/development phase, strongly suggests the significant role of the sunward surge of the plasmasphere and plasma plumes in the injection-associated (i.e., storm-time) EMIC wave generation.

Description: [1]  We performed a statistical study of downward moving protons and alpha particles of ~ keV energy (assumed to be of solar wind origin), inside the Martian induced magnetosphere from July 2006 to July 2010. Ion and electron data are from the Analyzer of Space Plasma and Energetic Atoms (ASPERA-3) package onboard Mars Express. We investigated the solar wind ion entry into the ionosphere, excluding intervals of low altitude magnetosheath encounters. The study compares periods of quiet solar wind conditions and periods of solar wind pressure pulses, including interplanetary coronal mass ejections and corotating interaction regions. The solar wind ion precipitation appears localized and/or intermittent, consistent with previous measurements. Precipitation events are less frequent and the precipitating fluxes do not increase during pressure pulse encounters. During pressure pulses the occurrence frequency of observed proton precipitation events is reduced by a factor ~3 and for He 2+ events the occurrence frequency is reduced by a factor ~2. One explanation is that during pressure pulse periods, the mass loading of the solar wind plasma increases due to a deeper penetration of the interplanetary magnetic flux tubes into the ionosphere. The associated decrease of the solar wind speed thus increases the pile-up of the interplanetary magnetic field on the dayside of the planet. The magnetic barrier becomes thicker in terms of solar wind ion gyroradii, causing the observed reduction of H + /He 2+ precipitations.

Description: [1]  We present THEMIS observations of a " Spontaneous Hot Flow Anomaly" (SHFA) upstream from the pre-noon bow shock at 0431 UT on August 12, 2007. Although the SHFA exhibited the greatly heated and deflected solar wind plasmas used to identify hot flow anomalies (HFAs), it did not result from the standard mechanism invoked for the formation of HFAs, namely the interaction of an interplanetary magnetic field (IMF) discontinuity with the bow shock. We employ THEMIS A, B, C, and D observations to describe the evolution of the event from a proto-SHFA exhibiting regions of depressed magnetic field strength and density but little evidence for plasma heating or flow deflection, to a well-developed SHFA further downstream. These observations show that SHFA can be generated without the presence of an IMF discontinuity and are therefore a new category of HFAs.

Description: [1]  Among the different types of waves embedded in seismic noise, body waves present appealing properties but are still challenging to extract. Here we first validate recent improvements in numerical modeling of microseismic compressional ( P ) body waves and then show how this tool allows fast detection and location of their sources. We compute sources at ~ 0.2 Hz within typical P teleseismic distances (30-90 degrees) from the South California Seismic Network (SCSN) and analyze the most significant discrete sources. The locations and relative strengths of the computed sources are validated by the good agreement with beam-forming analysis. These ~75 noise sources exhibit a highly heterogeneous distribution, and cluster along the usual storm tracks in the Pacific and Atlantic oceans. They are mostly induced in the open ocean, at or near water depths of 2800 and 5600 km, most likely within storms or where ocean waves propagating as swell meet another swell or wind sea. We then emphasize two particularly strong storms to describe how they generate noise sources in their wake. We also use these two specific noise bursts to illustrate the differences between microseismic body- and surface-waves in terms of source distribution and resulting recordable ground motion. The different patterns between body- and surface-waves result from distinctive amplification of ocean wave-induced pressure perturbation and different seismic attenuation. Our study demonstrates the potential of numerical modeling to provide fast and accurate constraints on where and when to expect microseismic body waves, with implications for seismic imaging and climate studies.

Description: [1]  Earthquakes that rupture across steps between faults can be larger than those predicted from individual fault lengths, making understanding multifault events critical to assessing earthquake hazard. Empirical data from earthquake surface ruptures suggest that the distances between faults that rupture together can range from 〈1 km to 5 km. Dynamic and quasi-static models of planar faults determine similar distances. However, studies of interactions between realistic, 3D non-planar faults are few. A general comparison of quasi-static stress perturbations and triggering potentials with mechanical models incorporating either planar or non-planar faults highlight the sensitivity of planar fault models to model parameters and reveal no clear relationship between mean fault slip and triggering potential. More specifically, planar fault models predict triggering across a 3 km extensional step, while models incorporating non-planar faults indicate that a connecting fault is necessary to transfer slip through a 3 km step along the 1992 Landers, California earthquake rupture. The mechanical approach taken captures the stress changes as well as the total stress following fault slip, improving the criterion used to determine triggered failure potential. This underscores the need for additional constraint on fault strength and cohesion. The focus on complex fault geometry restricts analyses to the quasi-static realm, limiting the application of results to fault interactions over the short distances and slow rupture velocities for which the quasi-static stress field is relevant or approximates the dynamic stress field.

Description: [1]  The two primary methods responsible for solar wind magnetosphere coupling are magnetic reconnection and the viscous interaction. The viscous interaction is generated due to the antisunward dragging of plasma inside the magnetopause by the plasma flowing in the magnetosheath, creating a return flow deeper inside the magnetosphere and producing a circulation pattern. This viscous circulation pattern is mapped into the ionosphere via magnetic field lines, which results in ionospheric electric field in the non-rotating Earth's frame. We measure this interaction in terms of an electric potential, the viscous potential. In this paper, we use the results obtained from the LFM simulation model during periods of purely northward IMF for different solar wind velocity and ionospheric conductivity, showing a reduction of the viscous potential with increasing magnitude of northward IMF. The viscous potential is found to settle around 5-10 kV for large + B z values. The decrease in viscous potential was found to be associated with a weak or non-existent sunward plasma flow in the nightside plasmasheet. Instead, the return flow to the dayside occurs at high latitudes and is associated with the reconnection topology and dynamics that occur during northward IMF periods. We also show that the magnetosphere remains closed during purely northward IMF, except for two small regions-one on each hemisphere, where the magnetic reconnection occur. We argue that the reduction of the viscous potential is due to a reduction of the velocity shear across the magnetopause and the lack of sunward convection in the equatorial tail.

Description: [1]  We present an analytical model of the Poynting flux that is generated by the interaction between the plume of Enceladus and Saturn's magnetospheric plasma. Our purpose is to analyze the influence of two key elements of Enceladus’ magnetospheric interaction on the electromagnetic energy radiated away in the moon's Alfvén wings. First, the north-south asymmetry of the obstacle generates a system of hemisphere coupling currents which allows a transport of electromagnetic energy into Saturn's northern hemisphere, even if the field-aligned currents connecting to the plume are completely blocked at the non-conducting icy crust of Enceladus. Second, the presence of electron-absorbing dust grains within the plume was recently found to drastically modify the electromagnetic field configuration within Enceladus’ Alfvén wings (Anti-Hall effect), thereby also altering the energy flux radiated away from the interaction region. By systematically studying the impact of varying strengths of the hemisphere coupling currents and varying electron absorption fractions on the energy flux, we come to the following conclusions: (1) The integrated Poynting flux into Saturn's southern hemisphere always exceeds the integrated flux into the northern hemisphere. In particular, the power transmitted towards the south may become several orders of magnitude larger than the power transmitted towards the north. (2) The search for Enceladus’ auroral footprint has so far mainly focused on Saturn's northern hemisphere. However, based on the Poynting fluxes radiated away by the interaction, detections of the footprint should occur more likely in the giant planet's southern hemisphere, if no other far-field effects play a role. (3) Electron absorption by the dust grains within the plume makes a measurable contribution to the energy flux. Thus, changes in the brightness of the observed footprints cannot exclusively be ascribed to variations in the incident magnetospheric flow conditions or the neutral gas content of the plume, but variations in the properties of the dust component need to be taken into account as well.

Description: [1]  Suprathermal singly-charged molecular ions, O 2 + (at ~32 Da/e) and the Mass-28 ion group 28 M + (ions at ~28 Da/e, with possible contributions from C 2 H 5 + , HCNH + , N 2 + , and/or CO + ), are present throughout Saturn's ~4-20 Rs (1 Saturn radius, Rs = 60,268 km) near-equatorial magnetosphere from mid-2004 until mid-2012. These ~83-167 keV/e heavy ions measured by Cassini's CHarge-Energy-Mass Spectrometer, CHEMS, have long-term temporal profiles that differ from each other and differ relative to the dominant water group ions, W + (O + , OH + , H 2 O + , and H 3 O + ). O 2 + /W + , initially ~0.05, declined steadily until equinox in mid-2009 by a factor of ~6, and 28 M + /W + , initially ~0.007, declined similarly until early-2007 by a factor of ~2. The O 2 + /W + decline is consistent with Cassini's in-situ ring-ionosphere thermal ion measurements, and with proposed and modeled seasonal photolysis of Saturn's rings for thermal O 2 and O 2 + . The water ice dominated main rings and Enceladus plume depositions thereon are the two most likely O 2 + sources. Enceladus' dynamic plumes, though, have no known long-term dependence. After declining, O 2 + /W + and 28 M + /W + levels remained low until late-2011 when O 2 + /W + increased, but 28 M + /W + did not. The O 2 + /W + increase was steady and became statistically significant by mid-2012, indicating a clear increase after a decline, that is, a possibly delayed O 2 + “seasonal” recovery. Ring insolation is driven by solar UV flux which itself varies with the sun's 11-year activity cycle. The O 2 + /W + and 28 M + /W + declines are consistent with seasonal ring insolation. No O 2 + /W + response to the late-2008 solar-cycle UV minimum and recovery is evident. However, the O 2 + /W + recovery from the post-equinox baseline levels in late-2011 coincided with a strong solar UV enhancement. We suggest a scenario/framework in which the O 2 + observations can be understood.

Description: [1]  We analyze ion temperature data near 350 km altitude over the years 1966–2012 to seek explanations for three outstanding questions concerning the long-term cooling observed in the upper thermosphere: (1) why is the cooling so much larger than expected, (2) why has the cooling lasted so long, and (3) why is the thermospheric density response to the cooling so small? We speculate that gravity waves may cause this cooling and provide answers to these questions. Recent simulations have shown that gravity waves are expected to cool the upper thermosphere by an amount comparable to that observed over our data timeline. A gravity-wave proxy formed from the non-tidal fluctuations in temperature shows a positive long-term trend throughout its timeline, consistent with the increasing cooling observed. The time scales of the long-term trend and its decadal fluctuations are characteristic of the ocean, not the atmosphere. We suggest that the following scenario may explain these behaviors: (a) the climate regime shift of 1976–1977 launched slow Rossby waves across the oceans which continue to propagate to this day, (b) winds over this increasingly corrugated ocean have launched increasing fluxes of gravity waves into the atmosphere, (c) these increasing fluxes of gravity waves have propagated to the thermosphere to produced increasing amounts of cooling. The strong thermospheric cooling seen would be expected to produce thermospheric density declines much larger than those observed via satellite drag. These temperature and density results would be compatible if the turbopause were lowered 4 km over the time span of observations.

Description: [1]  Batch and flow-through experiments were performed on quartz-feldspar granular aggregates and sandstone samples to investigate time-dependent effects of fluid-rock interactions on fluid and rock conductivity, respectively. The experiments were conducted at temperatures up to 164, at confining and pore pressures up to 10 and 5 MPa, respectively, and for up to 136 days. It showed that changes in rock conductivity were unequivocally related to changes in pore fluid conductivity. It is inferred that these changes were dependent on kinetically controlled dissolution reactions between the mineral grains and the fluid. The time-dependent signature of rock conductivity implied a detectable transition from initial dissolution towards some state of equilibrium. The response of rock conductivity to temperature changes followed an Arrhenius-type behavior. An exploratory kinetic evaluation of the conductivity data for sandstone samples yielded an apparent activation energy of approximately 32 kJ/mol. A concurrent chemical fluid analysis showed that this is an integrated value over all reactions occurring in parallel within a sample. These reactions namely concern silica and silicate dissolutionbut also the dissolution of accessory salt minerals. It is concluded that measuring the evolution of rock conductivity in combination with chemical pore fluid analysis constitutes a powerful and quantitative tool for monitoring time-dependent changesin pore fluid chemistry and thus fluid-rock interactions in real time.

Description: [1]  We describe a multi-parameter experiment at Erebus volcano, Antarctica, employing Doppler radar, video, acoustic, and seismic observations to estimate the detailed energy budget of large (up to 40-m-diameter) bubble bursts from a persistent phonolite lava lake. These explosions are readily studied from the crater rim at ranges of less than 500 m, and present an ideal opportunity to constrain the dynamics and mechanism of magmatic bubble bursts that can drive Strombolian and Hawaiian eruptions. We estimate the energy budget of the first second of a typical Erebus explosion as a function of time and energy type, and constrain gas pressures and forces using an analytic model for the expansion of a gas bubble above a conduit that incorporates conduit geometry and magma and gas parameters. The model, consistent with video and radar observations, invokes a spherical bulging surface with a base diameter equal to that of the lava lake. The model has no ad hoc free parameters, and geometrical calculations predict zenith height, velocity and acceleration during shell expansion. During explosions, the energy contained in hot over-pressured gas bubbles is freed and partitioned into other energy types, where by far the greatest non-thermal energy component is the kinetic and gravitational potential energy of the accelerated magma shell (〉10 9 J). Seismic source energy created by explosions is estimated from radar measurements and is consistent with source energy determined from seismic observations. For the generation of the infrasonic signal, a dual mechanism incorporating a terminally disrupted slug is proposed, which clarifies previous models and provides good fits to observed infrasonic pressures. A new and straightforward method is presented for determining gas volumes from slug explosions at volcanoes from remote infrasound recordings.

Description: [1]  We explore the application of GPS data to earthquake early warning and investigate whether the co-seismic ground deformation can be used to provide fast and reliable magnitude estimations and ground shaking predictions. We use an algorithm to extract the permanent static offset from GPS displacement time series and invert for the slip distribution on the fault plane, which is discretized into a small number of rectangular patches. We developed a completely “self-adapting” strategy in which the initial fault plane model is built based on a quick, approximate magnitude estimation, and is then allowed to increase in size based on the evolutionary magnitude estimation resulting from the slip inversion. Two main early warning outputs are delivered in real-time: magnitude and the along-strike extent of the rupture area. These are finally used to predict the expected ground shaking due to the finite source. We tested the proposed strategy by simulating real-time environments for three earthquakes. For the Mw 9.0, 2011 Tohoku-Oki earthquake our algorithm provides the first magnitude estimate of 8.2 at 39 sec after the origin time, and then gradually increases to 8.9 at 120 sec. The estimated rupture length remains constant from the outset at ~360 km. For the Mw 8.3, 2003 Tokachi-Oki earthquake the initial magnitude estimate is 8.5 at 24 sec and drops to 8.2 at 40 sec with a rupture length of 290 km. Finally, for the Mw 7.2, 2010 El Mayor-Cucapah earthquake the magnitude estimate is 7.0 from the outset with a rupture length of 140 km. The accuracy of the ground shaking prediction using the GPS-based magnitude and finite extent is significantly better than existing seismology-based point source approaches. This approach would also facilitate more rapid tsunami warnings

Description: ABSTRACT [1]  Analysis of Lake Bonneville shorelines using LIDAR digital elevation data challenges accepted models of Wasatch fault deformation since the late Pleistocene. While footwall deformation of the Weber segment of the Wasatch fault is consistent with back-rotation of the footwall block and greatest displacement rate towards the center of the segment, shorelines along the footwall of the Salt Lake City segment decrease in elevation towards the interior and are highest at the segment boundaries, an opposite pattern of footwall deformation than predicted for boundaries arresting or strongly inhibiting displacement during earthquakes. The spatial pattern of footwall rebound implies that some of the proposed persistent fault segment boundaries do not stop earthquake ruptures that originate on adjacent fault segments, nor constrain ruptures initiated within the Salt Lake City segment. Net vertical fault displacement at the boundary between the Salt Lake and Provo segments is 16—20 m over the past 16.3—18.5 ka, corresponding to a vertical displacement rate of 0.8—1.2 mm/yr, a net fault slip rate of 2.0—2.8 mm/yr and horizontal extension rate of 1.8—2.6 mm/yr on the 25 o west-southwest dipping fault that forms the southern Salt Lake City segment boundary. Shoreline analysis suggests isostatic rebound caused by a drop in lake level was concentrated during a relatively short (~2000 yr) time period following the Bonneville flood at ~16 ka. LIDAR-derived topography in conjunction with robust geomorphic datums improves our ability to map deformation associated with lithospheric flexure and faulting while demonstrating the limitation of lacustrine shorelines in this type of analysis.

Description: [1]  We use recent results on statistical analysis of seismicity to present a robust method for comprehensive detection and analysis of earthquake clusters. The method is based on nearest-neighbor distances of events in space-time-energy domain. The method is applied to a 1981–2011 relocated seismicity catalog of southern California having 111,981 events with magnitudes m  ≥ 2, and corresponding synthetic catalogs produced by the Epidemic Type Aftershock Sequence (ETAS) model. Analysis of the ETAS model demonstrates that the cluster detection results are accurate and stable with respect to (i) three numerical parameters of the method, (ii) variations of the minimal reported magnitude, (iii) catalog incompleteness, and (iv) location errors. Application of the method to the observed catalog separates the 111,981 examined earthquakes into 41,393 statistically significant clusters comprised of foreshocks , mainshocks and aftershocks . The results reproduce the essential known statistical properties of earthquake clusters, which provide overall support for the proposed technique. In addition, systematic analysis with our method allows us to detect several new features of seismicity that include (i) existence of a significant population of single-event clusters ; (ii) existence of foreshock activity in natural seismicity that exceeds expectation based on the ETAS model; and (iii) dependence of all cluster properties, except area, on the magnitude difference of events from mainshocks but not on their absolute values. The classification of detected clusters into several major types, generally corresponding to singles, burst-like and swarm-like sequences, and correlations between different cluster types and geographic locations is addressed in a companion paper.

Description: [1]  For the first time, we report the amplitude variation with angle (AVA) pattern of bottom simulating reflectors (BSRs) beneath fracture-filled gas hydrate deposits when the effective medium is anisotropic. The common depth point (CDP) gathers of two mutually perpendicular multi-channel seismic profiles, located in the vicinity of Site NGHP-01-10, are appropriately processed such that they are fit for AVA analysis. AVA analysis of the BSR shows normal-incidence reflection coefficients of -0.04 to -0.11 with positive gradients of 0.04 to 0.31 indicating class IV pattern. The acoustic properties from isotropic rock physics model predict class III AVA pattern which cannot explain the observed class IV AVA pattern in Krishna-Godavari basin due to the anisotropic nature of fracture-filled gas hydrate deposits. [2]  We modeled the observed class IV AVA of the BSR by assuming that the gas hydrate bearing sediment can be represented by horizontally transversely isotropic (HTI) medium after accounting for anisotropic wave propagation effects on BSR amplitudes. The effective medium properties are estimated using Backus averaging technique and the AVA pattern of BSRs is modeled using the properties of overlying HTI and underlying isotropy/HTI media with or without free gas. Anisotropic AVA analysis of the BSR from the inline seismic profile shows 5–30 % gas hydrate concentration (equivalent to fracture density) and the azimuth of fracture system (fracture orientation) with respect to the seismic profile is close to 45°. Free gas below the base of gas hydrate stability zone is interpreted in the vicinity of fault system (F1).

Description: [1]  We introduce a new model to study Jupiter's magnetosphere and how it interacts with the solar wind. We first derive a set of one-fluid MHD equations to consistently include the ion-neutral collisions in Jupiter's ionosphere, and the mass-loading in the Io torus. The mass-loading and the subsequent radial mass transport in Jupiter's magnetosphere leads to a deviation from full corotation of the magnetospheric plasma. Ion-neutral collisions in Jupiter's ionosphere and subsequent transport of angular momentum out into the magnetosphere acts to spin up the magnetosphere's plasma. Our model explicitly includes mass-loading in the Io plasma torus and an inner boundary region, which represents the effects of Jupiter's ionosphere. We present the results of five model runs where different mass-loading rates and ionospheric conductances are used. For these model runs, we consider an anti-parallel interplanetary magnetic field and a strong solar wind dynamic pressure, resulting in a compressed magnetosphere. The results are compared with analytical models, in situ measurements, and remote-sensing observations. Our azimuthal velocity profiles and the position of the corotation break-down are in quantitative agreement with theoretical predictions by Hill [1979, 2001] and Saur et al . [2004a], and Voyager observations [ McNutt et al ., 1981]. The total current flowing into and out of the ionosphere is 48.7 MA, which is in agreement with estimates from measurements and analytical models. Using the field aligned electric current j ║ to determine the position of the aurorae, we find that our main auroral oval is associated, as expected, with the position of the corotation break-down (between 20.6 R J and 30.1 R J for the different model runs). The discontinuity in the main oval observed by Radioti et al . [2008] is also present in our results, where it is caused by an asymmetry in the pressure distribution, due to the interaction between the rotating plasma and the magnetopause.

Description: [1]  In this study we quantify the contribution of individual large-scale waves to ionospheric electrodynamics, and examine the dependence of the ionospheric perturbations on solar activity. We focus on migrating diurnal tide (DW1) plus mean winds, migrating semidiurnal tide (SW2), quasi-stationary planetary wave 1 (QSPW1), and nonmigrating semidiurnal westward wave 1 (SW1) under northern winter conditions, when QSPW1 and SW1 are climatologically strong. From TIME-GCM simulations under solar minimum conditions, it is found that the mean winds and DW1 produce a wave 2 pattern in equatorial vertical E  ×  B drift that is upward in the morning and around dusk. The modeled SW2 also produces a wave 2 pattern in the ionospheric vertical drift that is nearly a half wave cycle out of phase with that due to mean winds and DW1. SW1 can cause large vertical drifts around dawn, while QSPW1 does not have any direct impact on the vertical drift. Wind components of both SW2 and SW1 become large at mid to high latitudes in the E-region, and kernel functions obtained from numerical experiments reveal that they can significantly affect the equatorial ion drift, likely through modulating the E-region wind dynamo. The most evident changes of total ionospheric vertical drift when solar activity is increased are seen around dawn and dusk, reflecting the more dominant role of large F-region Pedersen conductivity and of the F-region dynamo under high solar activity. Therefore, the lower atmosphere driving of the ionospheric variability is more evident under solar minimum conditions, not only because variability is more identifiable in a quieter background, but also because the E-region wind dynamo is more significant. These numerical experiments also demonstrate that the amplitudes, phases and latitudinal and vertical structures of large-scale waves are important in quantifying the ionospheric responses.

Description: [1]  An improved dispersion relation, with thermal corrections retained, for parallel propagating electromagnetic waves in a warm plasma is developed for both left-hand (L-) and right-hand (R-) polarized modes. Compared with the cold plasma dispersion relation, the derived dispersion relation is in much better agreement with the full hot plasma dispersion relation (including the wave growth rate). The pitch angle scattering rates of energetic electrons are compared between using cold and full dispersion relation. Significant differences are found when evaluating pitch angle scattering rate of MeV electron caused by He + band waves in multiple ion plasmas. Due to He + ion cyclotron absorption, the He + band EMIC waves, which are able to resonate with MeV electrons (at large wave number), tend to be strongly suppressed. Less significant differences in scattering rate of electrons between using cold and hot dispersion relation are found in the case of L-mode waves in single H + ion plasma, and in the case of R-mode waves.

Description: [1]  The fundamental features of ~0.1–0.2 s duration ~0.5 s spaced ionospheric electron precipitation “microbursts”, ~5 to 15 s microburst “trains”, and 5–15 s electron precipitation pulsations are reviewed in light of similar temporal structures of electromagnetic whistler mode “chorus” waves detected in the outer magnetosphere. Past observations of microbursts point to extremely rapid (ms timescale) wave-particle interactions, probably between lower-band chorus subelements (durations of ~10 to 100 ms) and energetic ~10 to 100 keV electrons. A recent theory explaining such rapid interaction rates observed in microbursts is briefly reviewed. Arguments are given why ~5–15 s x-ray (and optical) pulsations are also associated with chorus scattering of energetic electrons. Comments about relativistic (E 〉 1 MeV) microbursts are also provided. There are, however, many other unsolved problems of outer zone energetic electron precipitation. The authors will attempt to indicate several of these for the interested reader. Finally an appendix is provided for a brief review of two-frequency chorus and some current problems with that topic.