Search

On 10 January 2001, Cassini briefly entered into the magnetosphere of Jupiter, en route to Saturn. During this excursion into the Jovian magnetosphere, the Cassini Magnetosphere Imaging Instrument/Charge‐Energy‐Mass Spectrometer detected oxygen and sulfur ions. While Charge‐Energy‐Mass Spectrometer can distinguish between oxygen and sulfur charge states directly, only 95.9 ± 2.9 keV/e ions were sampled during this interval, allowing for a long time integration of the tenuous outer magnetospheric (~200 RJ) plasma at one energy. For this brief interval for the 95.9 keV/e ions, 96% of oxygen ions were O+, with the other 4% as O2+, while 25% of the energetic sulfur ions were S+, 42% S2+, and 33% S3+. The S2+/O+ flux ratio was observed to be 0.35 (±0.06 Poisson error)., Key Points Cassini measured the relative charge state abundances of 95.9 keV/e oxygen and sulfur in the outer magnetosphere (~200 RJ) of JupiterThe flux of 95.9 keV/e O+ was higher than that of S2+, with S2+ being the most abundant charge state among sulfur ionsThe relative abundances of 95.9 keV/e heavy ions are compared to thermal ions in the inner to middle magnetosphere

The past developments in the predictability of weather and climate are discussed from the point of view of nonlinear dynamical systems. The problems ahead for long‐range predictability extending into the climate time scale are also presented. The sensitive dependence of chaos on initial conditions and the imperfections in the models limit reliable predictability of the instantaneous state of the weather to less than 10 days in present‐day operational forecasts. The existence of slowly varying components such as the sea surface temperature, soil moisture, snow cover, and sea ice may provide basis for predicting certain aspects of climate at long range. The regularly varying nonlinear oscillations, such as the Madden‐Julian Oscillation, monsoon intraseasonal oscillations, and El Niño‐Southern Oscillation, are also possible sources of extended‐range predictability at the climate time scale. A prediction model based on phase space reconstruction has demonstrated that monsoon intraseasonal oscillation can be better predicted at long leads., Key Points The predictability of weather forecast models is limited to less than 10 days because of the limit imposed by chaos and model imperfectionsThe existence of slowly varying components of the climate system and regularly varying phenomena provide basis for climate predictabilityA prediction model of monsoon intraseasonal oscillation using phase space reconstruction shows that some aspects of climate can be predicted

Extended periods of northward interplanetary magnetic field (IMF) lead to the formation of a cold, dense plasma sheet due to the entry of solar wind plasma into the magnetosphere. Identifying the paths that the solar wind takes to enter the magnetosphere, and their relative importance has remained elusive. Any theoretical model of entry must satisfy observational constraints, such as the overall entry rate and the dawn‐dusk asymmetry observed in the cold, dense plasma sheet. We model, using a combination of global magnetohydrodynamic and test particle simulations, solar wind ion entry into the magnetosphere during northward IMF and compare entry facilitated by the Kelvin‐Helmholtz instability to cusp reconnection. For Kelvin‐Helmholtz entry we reproduce transport rates inferred from observation and kinetic modeling and find that intravortex reconnection creates buoyant flux tubes, which provides, through interchange instability, a mechanism of filling the central plasma sheet with cold magnetosheath plasma. For cusp entry we show that an intrinsic dawn‐dusk asymmetry is created during entry that is the result of alignment of the westward ion drift with the dawnward electric field typically observed during northward IMF. We show that both entry mechanisms provide comparable mass but affect entering plasma differently. The flank‐entering plasma is cold and dawn‐dusk symmetric, whereas the cusp‐entering plasma is accelerated and preferentially deflected toward dawn. The combined effect of these entry mechanisms results in a plasma sheet population that exhibits dawn‐dusk asymmetry in the manner that is seen in nature: a two‐component (hot and cold) dusk flank and hotter, broadly peaked dawn population., Key Points Global MHD model with kinetic test particles shows comparable solar wind entry from KH vortices and cusp reconnection under northward IMFSolar wind entry exhibits intrinsic dawn‐dusk asymmetry due to the preferential energization and dawnward deflection of cusp‐entering ionsCombined flank‐ and cusp‐entering solar wind ions recreates observed two‐component (hot and cold) dusk and broad‐peaked dawn population

Vb cyclones typically emerge in the Western Mediterranean and propagate to the Northeast into Central Europe. This paper explores the temporal characteristics of Vb cyclone occurrence based on cyclone tracks identified at the atmospheric levels of Z700 and sea level pressure, using JRA‐55 reanalysis data for the period 1959–2015. The risk of Vb occurrence was significantly high in the 1960s and has remained at a lower level since then. Vb cyclones do not occur fully randomly according to a Poisson point process. Eleven well‐separated and distinct clusters as well as 11 hiatus periods are identified, with average occurrence rates of 21.5 and 5.2 yrea−1, respectively. During the event of Vb, the large‐scale atmospheric circulation is changed into a state favoring the development of successive Vb cyclones. Clustering is very prominent in the case of Genoan Vb cyclones in summer as well as those Vb cyclones developing over the Iberian Peninsula or the North African Coast in winter. Superposition of the polar and the subtropical jet stream over the Western Mediterranean is identified as a main feature at the onset of Vb cyclones. Vb cyclone occurrence appears to be synchronized with the Northern Atlantic Oscillation (NAO; at Z500) and Arctic Oscillation (AO; at Z1000). Clusters have occurred when both NAO and AO were negative. This relation applies to Western Mediterranean cyclones not following a Vb track as well, however to a much weaker extent. In contrast, Vb cyclone frequency was particularly low from 1988 to 1997 during a sustained positive phase of both NAO and AO., Key Points In the past 50 years, the frequency of Vb cyclones was high in the 1960s and lower since thenVb cyclones were strongly synchronized with NAO and AO; high frequency clusters occurred when NAO and AO were negativeThe coupling of the polar and the subtropical jet stream over the Western Mediterranean is a driving mechanism of the onset of Vb cyclones

The COVID-19 pandemic led to widespread reductions in mobility and induced observable changes in atmospheric emissions. Recent work has employed novel mobility data sets as a proxy for trace gas emissions from traffic by scaling CO2 emissions linearly with those near-real-time mobility data. Yet, there has been little work evaluating these emission numbers. Here, we systematically compare these mobility data sets to traffic data from local governments in seven diverse urban and national/state regions to characterize the magnitude of errors that result from using the mobility data. We observe differences in excess of 60% between these mobility data sets and local traffic data. We could not find a general functional relationship between the mobility data and traffic flow over all the regions and observe higher deviations from using such general relationships than the original data. Finally, we give an overview of the potential errors that come from estimating CO2 emissions using (mobility or traffic) activity data. Future work should be cautious while using these mobility metrics for emission estimates., Journal of Geophysical Research: Atmospheres, 126 (11), ISSN:0148-0227, ISSN:2169-897X

Greenland’s outlet glaciers have been a leading source of mass loss and accompanying sea‐level rise from the Greenland Ice Sheet (GrIS) over the last 25 years. The dynamic component of outlet glacier mass loss depends on both the ice flux through the terminus and the inland extent of glacier thinning, initiated at the ice‐ocean interface. Here, we find limits to the inland spread of thinning that initiates at glacier termini for 141 ocean‐terminating outlet glaciers around the GrIS. Inland diffusion of thinning is limited by steep reaches of bed topography that we call “knickpoints.” We show that knickpoints exist beneath the majority of outlet glaciers but they are less steep in regions of gentle bed topography, giving glaciers in gentle bed topography the potential to contribute to ongoing and future mass loss from the GrIS by allowing the diffusion of thinning far into the ice sheet interior., Key Points Greenland glaciers with high potential to contribute to mass loss identified using novel metricThe majority of thinning limits are set by a knickpoint, a steep drop in bed topographyGlaciers that transmit thinning far inland are preferentially in regions with gentle bed topography

Electromagnetic ion cyclotron (EMIC) waves play important roles in particle loss processes in the magnetosphere. Determining the evolution of EMIC waves as they propagate and how this evolution affects wave‐particle interactions requires accurate knowledge of the wave vector, k. We present a technique using the curl of the wave magnetic field to determine k observationally, enabled by the unique configuration and instrumentation of the Magnetospheric MultiScale (MMS) spacecraft. The wave curl analysis is demonstrated for synthetic arbitrary electromagnetic waves with varying properties typical of observed EMIC waves. The method is also applied to an EMIC wave interval observed by MMS on October 28, 2015. The derived wave properties and k from the wave curl analysis for the observed EMIC wave are compared with the Waves in Homogenous, Anisotropic, Multi‐component Plasma (WHAMP) wave dispersion solution and with results from other single‐ and multi‐spacecraft techniques. We find good agreement between k from the wave curl analysis, k determined from other observational techniques, and k determined from WHAMP. Additionally, the variation of k due to the time and frequency intervals used in the wave curl analysis is explored. This exploration demonstrates that the method is robust when applied to a wave containing at least 3–4 wave periods and over a rather wide frequency range encompassing the peak wave emission. These results provide confidence that we are able to directly determine the wave vector properties using this multi‐spacecraft method implementation, enabling systematic studies of EMIC wave k properties with MMS., Key Points The wave curl analysis is a new implementation of determining k using observed wave magnetic field and associated current densityThe wave curl analysis reliably determines k for both synthetic waves and Magnetospheric MultiScale observations of electromagnetic ion cyclotron wavesThe calculated k is robust relative to time and frequency ranges used in the analysis, and agrees well with linear dispersion theory

Lakes are often defined by seasonal cycles. The seasonal timing, or phenology, of many lake processes are changing in response to human activities. However, long‐term records exist for few lakes, and extrapolating patterns observed in these lakes to entire landscapes is exceedingly difficult using the limited number of available in situ observations. Limited landscape‐level observations mean we do not know how common shifts in lake phenology are at macroscales. Here, we use a new remote sensing data set, LimnoSat‐US, to analyze U.S. summer lake color phenology between 1984 and 2020 across more than 26,000 lakes. Our results show that summer lake color seasonality can be generalized into five distinct phenology groups that follow well‐known patterns of phytoplankton succession. The frequency with which lakes transition from one phenology group to another is tied to lake and landscape level characteristics. Lakes with high inflows and low variation in their seasonal surface area are generally more stable, while lakes in areas with high interannual variations in climate and catchment population density show less stability. Our results reveal previously unexamined spatiotemporal patterns in lake seasonality and demonstrate the utility of LimnoSat‐US, which, with over 22 million remote sensing observations of lakes, creates novel opportunities to examine changing lake ecosystems at a national scale., Key Points Summer lake color phenology can be generalized into five distinct seasonal patterns of greening and blueing eventsSince the mid‐1990s, the number of lakes with color patterns corresponding to eutrophic waterbodies has been increasingWe observe these patterns using a new U.S. lake remote sensing data set that contains over 22 million lake observations

Explosive magnetotail activity has long been understood in the context of its auroral manifestations. While global models have been used to interpret and understand many magnetospheric processes, the temporal and spatial scales of some auroral forms have been inaccessible to global modeling creating a gulf between observational and theoretical studies of these phenomena. We present here an important step toward bridging this gulf using a newly developed global magnetosphere‐ionosphere model with resolution capturing ≲ 30 km azimuthal scales in the auroral zone. In a global magnetohydrodynamic (MHD) simulation of the growth phase of a synthetic substorm, we find the self‐consistent formation and destabilization of localized magnetic field minima in the near‐Earth magnetotail. We demonstrate that this destabilization is due to ballooning‐interchange instability which drives earthward entropy bubbles with embedded magnetic fronts. Finally, we show that these bubbles create localized field‐aligned current structures that manifest in the ionosphere with properties matching observed auroral beads., Key Points We present the first global magnetosphere simulation to reveal ballooning‐interchange instability of a narrow Bz minimum in the near‐Earth magnetotailThe instability is prominent during the substorm growth phase and generates earthward entropy bubbles with embedded magnetic frontsThe bubbles drive mesoscale ionospheric field‐aligned currents and auroral structures (beads) with properties matching to those observed

In situ measurements of the solar wind have been available for almost 60 years, and in that time plasma physics simulation capabilities have commenced and ground‐based solar observations have expanded into space‐based solar observations. These observations and simulations have yielded an increasingly improved knowledge of fundamental physics and have delivered a remarkable understanding of the solar wind and its complexity. Yet there are longstanding major unsolved questions. Synthesizing inputs from the solar wind research community, nine outstanding questions of solar wind physics are developed and discussed in this commentary. These involve questions about the formation of the solar wind, about the inherent properties of the solar wind (and what the properties say about its formation), and about the evolution of the solar wind. The questions focus on (1) origin locations on the Sun, (2) plasma release, (3) acceleration, (4) heavy‐ion abundances and charge states, (5) magnetic structure, (6) Alfven waves, (7) turbulence, (8) distribution‐function evolution, and (9) energetic‐particle transport. On these nine questions we offer suggestions for future progress, forward looking on what is likely to be accomplished in near future with data from Parker Solar Probe, from Solar Orbiter, from the Daniel K. Inouye Solar Telescope (DKIST), and from Polarimeter to Unify the Corona and Heliosphere (PUNCH). Calls are made for improved measurements, for higher‐resolution simulations, and for advances in plasma physics theory., Key Points Nine outstanding questions of solar wind physics are synthesized from inputs from the heliospheric research communityNew ways of viewing these questions are put forth, and suggestions for future progress are offeredCalls are made for improved measurements, simulations, and plasma physics theory

The Vlasiator hybrid‐Vlasov code was developed to investigate global magnetospheric dynamics at ion‐kinetic scales. Here we focus on the role of magnetic reconnection in the formation and evolution of magnetic islands at the low‐latitude magnetopause, under southward interplanetary magnetic field conditions. The simulation results indicate that (1) the magnetic reconnection ion kinetics, including the Earthward pointing Larmor electric field on the magnetospheric side of an X‐point and anisotropic ion distributions, are well‐captured by Vlasiator, thus enabling the study of reconnection‐driven magnetic island evolution processes, (2) magnetic islands evolve due to continuous reconnection at adjacent X‐points, “coalescence” which refers to the merging of neighboring islands to create a larger island, “erosion” during which an island loses magnetic flux due to reconnection, and “division” which involves the splitting of an island into smaller islands, and (3) continuous reconnection at adjacent X‐points is the dominant source of magnetic flux and plasma to the outer layers of magnetic islands resulting in cross‐sectional growth rates up to + 0.3 RE 2/min. The simulation results are compared to the Magnetospheric Multiscale (MMS) measurements of a chain of ion‐scale flux transfer events (FTEs) sandwiched between two dominant X‐lines. The MMS measurements similarly reveal (1) anisotropic ion populations and (2) normalized reconnection rate ~0.18, in agreement with theory and the Vlasiator predictions. Based on the simulation results and the MMS measurements, it is estimated that the observed ion‐scale FTEs may grow Earth‐sized within ~10 min, which is comparable to the average transport time for FTEs formed in the subsolar region to the high‐latitude magnetopause. Future simulations shall revisit reconnection‐driven island evolution processes with improved spatial resolutions., Key Points Anisotropic ion distributions are reported in Vlasiator simulations and MMS observations of reconnection inflow regionsThe 2‐D simulations suggest magnetic islands grow mainly via continuous reconnection. Island coalescence, erosion, and division are also presentBased on simulation results, ion‐scale FTEs are estimated to grow at

In this study, the ion composition of flux transfer events (FTEs) observed within the magnetosheath proper is examined. These FTEs were observed just upstream of the Earth's postnoon magnetopause by the National Aeronautics and Space Administration (NASA) Magnetospheric Multiscale (MMS) spacecraft constellation. The minor ion characteristics are described using energy spectrograms, flux distributions, and ion moments as the constellation encountered each FTE. In conjunction with electron data and magnetic field observations, such observations provide important contextual information on the formation, topologies, and evolution of FTEs. In particular, minor ions, when combined with the field‐aligned streaming of electrons, are reliable indicators of FTE topology. The observations are also placed (i) in context of the solar wind magnetic field configuration, (ii) the connection of the sampled flux tube to the ionosphere, and (iii) the location relative to the modeled reconnection line at the magnetopause. While protons and alpha particles were often depleted within the FTEs relative to the surrounding magnetosheath plasma, the He+ and O+ populations showed clear enhancements either near the center or near the edges of the FTE, and the bulk plasma flow directions are consistent with magnetic reconnection northward of the spacecraft and convection from the dayside toward the flank magnetopause., Key Points Long‐duration magnetosheath FTEs observed under at postnoon local times and under similar solar wind conditions are compared and contrastedMinor ions from the magnetosphere are observed at higher energies at the edges than at the center of the FTEsThe cores of the magnetosheath FTEs are on closed field lines

We discuss a family of Vlasov-Maxwell equilibrium distribution functions for current sheet equilibria that are intermediate cases between the Harris sheet and the force-free (or modified) Harris sheet. These equilibrium distribution functions have potential applications to space and astrophysical plasmas. The existence of these distribution function had been briefly discussed in by Harrison and Neukirch 2009, but here it is shown that their approach runs into problems in the limit where the guide field goes to zero. The nature of this problem will be discussed and an alternative approach will be suggested that avoids the problem. This is achieved by considering a slight variation of the magnetic field profile, which allows a smooth transition between the Harris and force-free Harris sheet cases., 13 pages, 6 figures, accepted by Journal of Plasma Physics

At the Earth's low‐latitude magnetopause, clear signatures of the Kelvin‐Helmholtz (KH) waves have been frequently observed during periods of the northward interplanetary magnetic field (IMF), whereas these signatures have been much less frequently observed during the southward IMF. Here, we performed the first 3‐D fully kinetic simulation of the magnetopause KH instability under the southward IMF condition. The simulation demonstrates that fast magnetic reconnection is induced at multiple locations along the vortex edge in an early nonlinear growth phase of the instability. The reconnection outflow jets significantly disrupt the flow of the nonlinear KH vortex, while the disrupted turbulent flow strongly bends and twists the reconnected field lines. The resulting coupling of the complex field and flow patterns within the magnetopause boundary layer leads to a quick decay of the vortex structure, which may explain the difference in the observation probability of KH waves between northward and southward IMF conditions., Key Points Three‐dimensional fully kinetic simulation of Kelvin‐Helmholtz instability at the Earth's magnetopause under the southward IMF condition is performedFast reconnection causes a rapid decay of the nonlinear vortex structure in the early nonlinear growth phase of the instabilityThe vortex decay can lead to a lower probability of observing magnetopause Kelvin‐Helmholtz waves/vortices during southward IMF periods

On 5 May 2017, MMS observed a crater‐type flux rope on the dawnside tailward magnetopause with fluctuations. The boundary‐normal analysis shows that the fluctuations can be attributed to nonlinear Kelvin‐Helmholtz (KH) waves. Reconnection signatures such as flow reversals and Joule dissipation were identified at the leading and trailing edges of the flux rope. In particular, strong northward electron jets observed at the trailing edge indicated midlatitude reconnection associated with the 3‐D structure of the KH vortex. The scale size of the flux rope, together with reconnection signatures, strongly supports the interpretation that the flux rope was generated locally by KH vortex‐induced reconnection. The center of the flux rope also displayed signatures of guide‐field reconnection (out‐of‐plane electron jets, parallel electron heating, and Joule dissipation). These signatures indicate that an interface between two interlinked flux tubes was undergoing interaction, causing a local magnetic depression, resulting in an M‐shaped crater flux rope, as supported by reconstruction., Key Points MMS observed a magnetic flux rope formed on the boundary of a nonlinear Kelvin‐Helmholtz (KH) waveBoth in‐plane and midlatitude reconnection associated with the 3‐D structure of the KH vortex‐inducedflux rope were identifiedA current sheet at the flux rope center supported by reconstruction indicates two flux tubes interlinked to form a crater‐type flux rope

We have designed and built a versatile modularized software library—ODYN—that wraps a comprehensive set of advanced data analysis methods meant to facilitate the study of turbulence, nonlinear dynamics, and complexity in space plasmas. The Python programming language is used for the algorithmic implementation of models and methods devised to understand fundamental phenomena of space plasma physics like elements of spectral analysis, probability distribution functions and their moments, multifractal analysis, or information theory. ODYN is an open‐source software analysis tool and freely available to any user interested in turbulence and nonlinear dynamics analysis and provides a tool to perform automatic analysis on large collections of space measurements, in situ or simulations, a feature that distinguishes ODYN from other similar software. A user‐friendly configurator is provided, which allows customization of key parameters of the analysis methods, most useful for nonprogrammers., Key Points We describe an open‐source software data analysis tool based on PythonThe software includes a large portfolio of methods to analyze turbulence and nonlinear dynamics and visualize the resultsThe software tool is adapted to ingest and process large collections of spacecraft data as well as results of numerical simulations

During the September 2015 to March 2016 duration (sometimes referred to as Phase 1A) of the Magnetospheric Multiscale Mission, the Dual Electron Spectrometers (DES) were configured to generously utilize lossy compression. While this maximized the number of velocity distribution functions downlinked, it came at the expense of lost information content for a fraction of the frames. Following this period of lossy compression, the DES was reconfigured in a way that allowed for 95% of the frames to arrive to the ground without loss. Using this high‐quality set of frames from on‐orbit observations, we compressed and decompressed the frames on the ground to create a side‐by‐side record of the compression effect. This record was used to drive an optimization method that (a) derived basis functions capable of approximating the lossless sample space and with nonnegative coefficients and (b) fitted a function which maps the lossy frames to basis weights that recreate the frame without compression artifacts. This method is introduced and evaluated in this paper. Data users should expect a higher level of confidence in the absolute scale of density/temperature measurements and notice less sinusoidal bias in the velocity X and Y components (GSE), Key Points 57% of the burst resolution frames in the September 2015 to March 2016 duration of the MMS mission were affected by some form of lossy compression errorThis data can be corrected by optimizing a model with side‐by‐side examples of the before/after compression effectData users should expect a higher level of confidence in the absolute scale of density/temperature measurements and notice less sinusoidal bias in the velocity X and Y components (GSE)

The evolution of a linearly-polarized, long-wavelength Alfv\'en wave --propagating in a collisionless magnetized plasma with a sheared parallel-directed velocity flow-- is here studied by means of two-dimensional hybrid Vlasov-Maxwell (HVM) simulations. The unperturbed sheared flow has been represented by an exact solution of the HVM set of equations (Malara {\it et al.}, Phys. Rev. E 97, 053212), this avoiding spurious oscillations that would arise from the non-stationary initial state and inevitably affect the dynamics of the system. We have considered the evolution of both a small and a moderate amplitude Alfv\'en wave, in order to separate linear wave-shear flow couplings from kinetic effects, the latter being more relevant for larger wave amplitudes. The phase-mixing generated by the shear flow modifies the initial perturbation, leading to the formation of small-scale transverse fluctuations at scales comparable with the proton inertial length. By analyzing both the polarization and group velocity of perturbations in the shear regions, we identify them as Kinetic Alfv\'en Waves (KAWs). In the moderate amplitude run, kinetic effects distort the proton distribution function in the shear region. This leads to the formation of a proton beam, at the Alfv\'en speed and parallel to the magnetic field. Such a feature, due to the parallel electric field associated with KAWs, positively compares with solar-wind observations of suprathermal ions' populations, suggesting that it may be related to the presence of ion-scales KAW-like fluctuations., Comment: 20 pages, 8 figures

Microbursts are an impulsive increase of electrons from the radiation belts into the atmosphere and have been directly observed in low Earth orbit and the upper atmosphere. Prior work has estimated that microbursts are capable of rapidly depleting the radiation belt electrons on the order of a day; hence, their role to radiation belt electron losses must be considered. Losses due to microbursts are not well constrained, and more work is necessary to accurately quantify their contribution as a loss process. To address this question, we present a statistical study of >35 keV microburst sizes using the pair of AeroCube‐6 CubeSats. The microburst size distribution in low Earth orbit and the magnetic equator was derived using both spacecraft. In low Earth orbit, the majority of microbursts were observed, while the AeroCube‐6 separation was less than a few tens of kilometers, mostly in latitude. To account for the statistical effects of random microburst locations and sizes, Monte Carlo and analytic models were developed to test hypothesized microburst size distributions. A family of microburst size distributions were tested, and a Markov chain Monte Carlo sampler was used to estimate the optimal distribution of model parameters. Finally, a majority of observed microbursts map to sizes less than 200 km at the magnetic equator. Since microbursts are widely believed to be generated by scattering of radiation belt electrons by whistler mode waves, the observed microburst size distribution was compared to whistler mode chorus size distributions derived in prior literature., Key Points The dual AeroCube‐6 CubeSats simultaneously observed >35 keV microbursts at a variety of spatial separations ranging from 2 to 100 kmIn low Earth orbit the majority of microbursts have a size on the order of a few tens of kmMapped to the magnetic equator, the majority of microbursts are less than 200 km in size, corresponding to the size of chorus wave packets

Inward radial diffusion driven by ULF waves has long been known to be capable of accelerating radiation belt electrons to very high energies within the heart of the belts, but more recent work has shown that radial diffusion values can be highly event‐specific, and mean values or empirical models may not capture the full significance of radial diffusion to acceleration events. Here we present an event of fast inward radial diffusion, occurring during a period following the geomagnetic storm of 17 March 2015. Ultrarelativistic electrons up to ∼8 MeV are accelerated in the absence of intense higher‐frequency plasma waves, indicating an acceleration event in the core of the outer belt driven primarily or entirely by ULF wave‐driven diffusion. We examine this fast diffusion rate along with derived radial diffusion coefficients using particle and fields instruments on the Van Allen Probes spacecraft mission., Key Points Fast radial diffusion of ultrarelativistic electrons is observed days after storm main phaseEvent‐specific radial diffusion can be orders of magnitude higher than statistical valuesULF‐wave driven acceleration can account for intense particle enhancement observed in inner magnetosphere

The goal of this paper is to understand the processes by which solar wind electrons are energized in the Martian magnetosphere and how this compares to processes at Venus and Earth. Each is unique in the source of its magnetic field topology and how this influences electron energization. To achieve this goal, 24 million spectra spanning 13 years have been examined using the electron spectrometer from the Mars Express spacecraft between about 12,000 km and about 250 km altitude, and from all latitudes and local times. The top 10 largest differential energy flux at energies above the differential energy flux peak have been found: seven spectra from the magnetosheath near noon, three from the dark tail (the largest two from the middle and ionospheric edge of the magnetosheath). Spectral comparisons show a decade range in the peak of the electron distributions; however, all distributions show a similar energy maximum dictated by solar wind/planet interaction. Similarly derived, the largest Venus spectrum occurred near the magnetosheath bow shock and had the same shape as the most intense Mars inner magnetosheath spectrum. The Mars and Venus dayside spectra compared to the Mars nightside spectrum that included an enhanced optical signal attributed to discrete “auroral” precipitation show a similar shape. These spectra are also compared to a selected auroral zone electron spectra from the Earth. The Mars and Venus results suggest that there is no more energy needed to generate electrons forming the nightside precipitation than is gained during the solar wind/planet interaction., Key Points Top 10 orbit‐selected electron energy spectra from Mars are identified from MEx ASPERA‐3 ELS between 2004 and 2016The Mars and Venus dayside electron spectra are similar, and a discrete “auroral” spectrum from the nightside tail of Mars is also similarNo more energy than that supplied by the solar wind‐planet interaction is required to produce discrete auroral spectrum from the nightside tail of Mars

Cold pools influence convective initiation and organization, dust lofting, and boundary layer properties, but little is known about their interactions with the land surface, particularly in dry continental environments. In this study, two‐way cold pool‐land surface interactions are investigated using high‐resolution idealized simulations of an isolated, transient cold pool evolving in a dry convective boundary layer. Results using a fully interactive land surface demonstrate that sensible heat fluxes are suppressed at the center of the cold pool but enhanced at the edge due to the spatial patterns of land surface cooling and the air temperature and wind speed perturbations. This leads to cold pool dissipation from the edge inward. Latent heat fluxes are primarily suppressed within the cold pool, and the magnitude of this suppression is controlled by competition between atmospheric and land surface effects. By comparing the fully interactive land surface simulation to a simulation with imposed surface fluxes, the land surface‐cold pool feedbacks are shown to reduce the cold pool lifetime, extent, and intensity by up to 50% and influence the pattern of boundary layer turbulent kinetic energy recovery, which have significant implications for cold pool‐induced convective initiation., Key Points An isolated cold pool evolving in a dry convective boundary layer is simulated at high resolutionSensible heat fluxes are enhanced at the cold pool's edge and suppressed in the center, leading to dissipation from the edge inwardLand surface interactions reduce the cold pool lifetime, area, and intensity by up to 50% and affect the boundary layer evolution

We report global observations of high‐m poloidal waves during the recovery phase of the 22 June 2015 magnetic storm from a constellation of widely spaced satellites of five missions including Magnetospheric Multiscale (MMS), Van Allen Probes, Time History of Events and Macroscale Interactions during Substorm (THEMIS), Cluster, and Geostationary Operational Environmental Satellites (GOES). The combined observations demonstrate the global spatial extent of storm time poloidal waves. MMS observations confirm high azimuthal wave numbers (m ~ 100). Mode identification indicates the waves are associated with the second harmonic of field line resonances. The wave frequencies exhibit a decreasing trend as L increases, distinguishing them from the single‐frequency global poloidal modes normally observed during quiet times. Detailed examination of the instantaneous frequency reveals discrete spatial structures with step‐like frequency changes along L. Each discrete L shell has a steady wave frequency and spans about 1 R E, suggesting that there exist a discrete number of drift‐bounce resonance regions across L shells during storm times., Key Points Observed long‐lasting high‐m poloidal waves associated with second harmonics of field line resonances during a major magnetic stormDemonstrated global spatial extent of storm time poloidal FLR region using observations from a constellation of widely spaced satellitesRevealed discrete spatial structures of resonant L shells with step‐like frequency changes

We use measurements from NASA's Van Allen Probes to calculate the decay time constants for electrons over a wide range of energies (30 keV to 4 MeV) and L values ( L = 1.3–6.0) in the Earth's radiation belts. Using an automated routine to identify flux decay events, we construct a large database of lifetimes for near‐equatorially mirroring electrons over a 5‐year interval. We provide the first accurate estimates of the long decay timescales in the inner zone ( ∼100 days), which are highly resolved in energy and free from proton contamination. In the slot region and outer zone, we compare our lifetime calculations with prior empirical estimates and find good quantitative agreement (lifetimes ∼1–20 days). The comparisons suggest that some prior estimates may overestimate electron lifetimes between L≈ 2.5–4.5 due to instrumental effects and/or background contamination. Previously reported two‐stage decays are explicitly demonstrated to be a consequence of using integral fluxes., Key Points A large database of radiation belt electron decay timescales is calculated from Van Allen Probes MagEIS measurementsWe provide the first accurate estimates of these timescales over a wide range of energies in the inner zone, free from proton contaminationOuter zone decay timescales generally agree well with prior estimates; some differences exist and may be due to instrumental effects

We compute quasilinear diffusion rates due to pitch angle scattering by various mechanisms in the Earth's electron radiation belts. The calculated theoretical lifetimes are compared with observed decay rates, and we find excellent qualitative agreement between the two. The overall structure of the observed lifetime profiles as a function of energy and L is largely due to plasmaspheric hiss and Coulomb scattering. The results also reveal a local minimum in lifetimes in the inner zone at lower energy ( ∼50 keV), attributed to enhanced scattering via ground‐based very low frequency transmitters, and a reduction in lifetimes at higher L and energy ( >1 MeV), attributed to enhanced electromagnetic ion cyclotron wave scattering. In addition, we find significant quantitative disagreement at L, Key Points Electron scattering by Coulomb collisions, hiss, EMIC, and VLF transmitter waves is demonstrated in observed loss timescalesThe energy‐ L structure of the observed electron lifetimes is in excellent qualitative agreement with quasilinear theoryTheoretical lifetimes in the inner region ( L

The solar wind is slowed, deflected, and heated as it encounters Venus's induced magnetosphere. The importance of kinetic plasma processes to these interactions has not been examined in detail, due to a lack of constraining observations. In this study, kinetic‐scale electric field structures are identified in the Venusian magnetosheath, including plasma double layers. The double layers may be driven by currents or mixing of inhomogeneous plasmas near the edge of the magnetosheath. Estimated double‐layer spatial scales are consistent with those reported at Earth. Estimated potential drops are similar to electron temperature gradients across the bow shock. Many double layers are found in few high cadence data captures, suggesting that their amplitudes are high relative to other magnetosheath plasma waves. These are the first direct observations of plasma double layers beyond near‐Earth space, supporting the idea that kinetic plasma processes are active in many space plasma environments., Key Points Plasma double layers are detected near the Venusian bow shockMultiple double layers are identified in a small amount of burst dataKinetic processes may help mediate interaction between the solar wind and induced magnetospheres

Energetic neutral atoms (ENAs) created by charge‐exchange of ions with the Earth's hydrogen exosphere near the subsolar magnetopause yield information on the distribution of plasma in the outer magnetosphere and magnetosheath. ENA observations from the Interstellar Boundary Explorer (IBEX) are used to image magnetosheath plasma and, for the first time, low‐energy magnetospheric plasma near the magnetopause. These images show that magnetosheath plasma is distributed fairly evenly near the subsolar magnetopause; however, low‐energy magnetospheric plasma is not distributed evenly in the outer magnetosphere. Simultaneous images and in situ observations from the Magnetospheric Multiscale (MMS) spacecraft from November 2015 (during the solar cycle declining phase) are used to derive the exospheric density. The ~11–17 cm−3 density at 10 RE is similar to that obtained previously for solar minimum. Thus, these combined results indicate that the exospheric density 10 RE from the Earth may have a weak dependence on solar cycle., Key Points ENA cameras image both magnetosheath and magnetospheric plasmas in the vicinity of the subsolar magnetopauseMagnetospheric plasma is not distributed evenly across the dayside near the magnetopauseThe exospheric hydrogen density near the magnetopause may have a weak dependence on solar F10.7

We develop an optimization approach to model the magnetic field configuration of magnetic clouds, based on a linear-force free formulation in three dimensions. Such a solution, dubbed the Freidberg solution, is kin to the axi-symmetric Lundquist solution, but with more general "helical symmetry". The merit of our approach is demonstrated via its application to two case studies of in-situ measured magnetic clouds. Both yield results of reduced $\chi^2\approx 1$. Case 1 shows a winding flux rope configuration with one major polarity. Case 2 exhibits a double-helix configuration with two flux bundles winding around each other and rooted on regions of mixed polarities. This study demonstrates the three-dimensional complexity of the magnetic cloud structures., Comment: Submitted to ApJL on Aug 24, 2020; rejected by Editor without review... Resubmitted to GRL, in press (see DOI; Open Access)

Stochastic schemes, designed to represent unresolved subgrid‐scale variability, are frequently used in short and medium‐range weather forecasts, where they are found to improve several aspects of the model. In recent years, the impact of stochastic physics has also been found to be beneficial for the model's long‐term climate. In this paper, we demonstrate for the first time that the inclusion of a stochastic physics scheme can notably affect a model's projection of global warming, as well as its historical climatological global temperature. Specifically, we find that when including the “stochastically perturbed parametrization tendencies” (SPPT) scheme in the fully coupled climate model EC‐Earth v3.1, the predicted level of global warming between 1850 and 2100 is reduced by 10% under an RCP8.5 forcing scenario. We link this reduction in climate sensitivity to a change in the cloud feedbacks with SPPT. In particular, the scheme appears to reduce the positive low cloud cover feedback and increase the negative cloud optical feedback. A key role is played by a robust, rapid increase in cloud liquid water with SPPT, which we speculate is due to the scheme's nonlinear interaction with condensation., Key Points The inclusion of a stochastic scheme reduces climate sensitivity in a general circulation modelThis reduction, of around 10%, is linked to changes in cloud cover and cloud optical depth feedbacksWell‐calibrated stochastic schemes may give more accurate global warming projections

著者人数: 20名, Accepted: 2019-05-28, 資料番号: SA1190061000

To achieve breakthroughs in the areas of heliospheric and magnetospheric energetic neutral atom (ENA) imaging, a new class of instruments is required. We present a high angular resolution ENA imager concept aimed at the suprathermal plasma populations with energies between 0.5 and 20 keV. This instrument is intended for understanding the spatial and temporal structure of the heliospheric boundary recently revealed by Interstellar Boundary Explorer instrumentation and the Cassini Ion and Neutral Camera. The instrument is also well suited to characterize magnetospheric ENA emissions from low‐altitude ENA emissions produced by precipitation of magnetospheric ions into the terrestrial upper atmosphere, or from the magnetosheath where solar wind protons are neutralized by charge exchange, or from portions of the ring current region. We present a new technique utilizing ultrathin carbon foils, 2‐D collimation, and a novel electron optical design to produce high angular resolution (≤2°) and high‐sensitivity (≥10−3 cm2 sr/pixel) ENA imaging in the 0.5–20 keV energy range., Key Points The LENI instrument is capable of high angular resolution ENA imagingThe LENI instrument utilizes novel MCP collimation for ENA imagingThe LENI instrument is designed for heliospheric and magnetospheric imaging

We present a statistical survey of ultralow‐frequency wave activity within the Hermean magnetosphere using the entire MErcury Surface, Space ENvironment, GEochemistry, and Ranging magnetometer data set. This study is focused upon wave activity with frequencies, Key Points Evidence that the Kelvin‐Helmholtz instability is driving Hermean ULF wave activityObservations suggest that Earth‐like field line resonance is possible at MercuryWave power mapped to Mercury's surface reveals the likely average polar cap boundary location

Energetic charged particle detectors characterize a portion of the plasma distribution function that plays critical roles in some physical processes, from carrying the currents in planetary ring currents to weathering the surfaces of planetary objects. For several low‐resource missions in the past, the need was recognized for a low‐resource but highly capable, mass‐species‐discriminating energetic particle sensor that could also obtain angular distributions without motors or mechanical articulation. This need led to the development of a compact Energetic Particle Detector (EPD), known as the “Puck” EPD (short for hockey puck), that is capable of determining the flux, angular distribution, and composition of incident ions between an energy range of ~10 keV to several MeV. This sensor makes simultaneous angular measurements of electron fluxes from the tens of keV to about 1 MeV. The same measurements can be extended down to approximately 1 keV/nucleon, with some composition ambiguity. These sensors have a proven flight heritage record that includes missions such as MErcury Surface, Space ENvironment, GEochemistry, and Ranging and New Horizons, with multiple sensors on each of Juno, Van Allen Probes, and Magnetospheric Multiscale. In this review paper we discuss the Puck EPD design, its heritage, unexpected results from these past missions and future advancements. We also discuss high‐voltage anomalies that are thought to be associated with the use of curved foils, which is a new foil manufacturing processes utilized on recent Puck EPD designs. Finally, we discuss the important role Puck EPDs can potentially play in upcoming missions., Key Points Review of the compact Energetic Particle Detector known as the PuckPuck EPD heritage includes five successful scientific missionsUnexpected results and potential advancements of the Puck are discussed

MMS3 spacecraft passed the vicinity of the electron diffusion region of magnetotail reconnection on 3 July 2017, observing discrepancies between perpendicular electron bulk velocities and E→×B→ drift, and agyrotropic electron crescent distributions. Analyzing linear wave dispersions, Burch et al. (2019, https://doi.org/10.1029/2019GL082471) showed the electron crescent generates high‐frequency waves. We investigate harmonics of upper‐hybrid (UH) waves using both observation and particle‐in‐cell (PIC) simulation, and the generation of electromagnetic radiation from PIC simulation. Harmonics of UH are linearly polarized and propagate along the perpendicular direction to the ambient magnetic field. Compared with two‐dimensional PIC simulation and nonlinear kinetic theory, we show that the nonlinear beam‐plasma interaction between the agyrotropic electrons and the core electrons generates harmonics of UH. Moreover, PIC simulation shows that agyrotropic electron beam can lead to electromagnetic (EM) radiation at the plasma frequency and harmonics., Key Points Harmonics of upper‐hybrid waves are observed near an electron diffusion region of tail reconnectionNonlinear beam‐plasma interaction between the core and crescent electrons is responsible for harmonics generationParticle‐in‐cell simulation result suggests that radio emission is a result of the beam‐plasma interaction

A simple Monte Carlo model is presented that considers the effects of spacecraft orbital sampling on the inferred distribution of magnetic flux ropes, generated through magnetic reconnection in the magnetotail current sheet. When generalized, the model allows the determination of the number of orbits required to constrain the underlying population of structures: It is able to quantify this as a function of the physical parameters of the structures (e.g., azimuthal extent and probability of generation). The model is shown adapted to the Hermean magnetotail, where the outputs are compared to the results of a recent survey. This comparison suggests that the center of Mercury's neutral line is located dawnward of midnight by 0.37−1.02+1.21RM and that the flux ropes are most likely to be wide azimuthally (∼50% of the width of the Hermean tail). The downtail location of the neutral line is not self‐consistent or in agreement with previous (independent) studies unless dissipation terms are included planetward of the reconnection site; potential physical explanations are discussed. In the future the model could be adapted to other environments, for example, the dayside magnetopause or other planetary magnetotails., Key Points Monte Carlo model allows the estimation of X‐line location and reconnection frequency given sampling with a single spacecraftMercury's magnetotail reconnection site is consistent with a center offset ( 0.37‐1.02+1.21RM) dawnward of midnightMercury's downtail X‐line location is only self‐consistent if dissipation terms are included planetward of the X‐line

Satellite charging is one of the most important risks for satellites on orbit. Satellite charging can lead to an electrostatic discharge resulting in component damage, phantom commands, and loss of service and in exceptional cases total satellite loss. Here we construct a realistic worst case for a fast solar wind stream event lasting 5 days or more and use a physical model to calculate the maximum electron flux greater than 2 MeV for geostationary orbit. We find that the flux tends toward a value of 106 cm−2·s−1·sr−1 after 5 days and remains high for another 5 days. The resulting flux is comparable to a 1 in 150‐year event found from an independent statistical analysis of electron data. Approximately 2.5 mm of Al shielding would be required to reduce the internal charging current to below the National Aeronautics and Space Administration‐recommended guidelines, much more than is currently used. Thus, we would expect many satellites to report electrostatic discharge anomalies during such an event with a strong likelihood of service outage and total satellite loss. We conclude that satellites at geostationary orbit are more likely to be at risk from fast solar wind stream event than a Carrington‐type storm., Key Points We calculate the electron flux at geostationary orbit for a realistic worst‐case eventApproximately 2.4 mm of Al is required to shield electronic componentsSatellites at geostationary orbit are more at risk from a fast solar wind stream event than a Carrington‐type event

A Monte Carlo method of investigating the effects of placing selection criteria on the magnetic signature of in situ encounters with flux ropes is presented. The technique is applied to two recent flux rope surveys of MESSENGER data within the Hermean magnetotail. It is found that the different criteria placed upon the signatures will preferentially identify slightly different subsets of the underlying population. Quantifying the selection biases first allows the distributions of flux rope parameters to be corrected, allowing a more accurate estimation of the intrinsic distributions. This is shown with regard to the distribution of flux rope radii observed. When accounting for the selection criteria, the mean radius of Hermean magnetotail quasi‐force‐free flux ropes is found to be 589−269+273 km. Second, it is possible to weight the known identifications in order to determine a rate of recurrence that accounts for the presence of the structures that will not be identified. In the case of the Hermean magnetotail, the average rate of quasi‐force‐free flux ropes is found to 0.12 min−1 when selection effects are accounted for (up from 0.05 min−1 previously inferred from observations)., Key Points A Monte Carlo method is presented to estimate and correct for sampling biases in spacecraft surveys of magnetic flux ropesMethod allows the correction of observed distributions (e.g., flux rope radii) according to the selection criteria employedAccounting for unidentified flux ropes increases the average rate of flux ropes in Mercury's magnetotail from 0.05 to 0.12 min−1

The radiation belts and magnetospheres of Jupiter and Saturn show significant intensities of relativistic electrons with energies up to tens of megaelectronvolts (MeV). To date, the question on how the electrons reach such high energies is not fully answered. This is largely due to the lack of high‐quality electron spectra in the MeV energy range that models could be fit to. We reprocess data throughout the Galileo orbiter mission in order to derive Jupiter's electron spectra up to tens of MeV. In the case of Saturn, the spectra from the Cassini orbiter are readily available and we provide a systematic analysis aiming to study their acceleration mechanisms. Our analysis focuses on the magnetospheres of these planets, at distances of L > 20 and L > 4 for Jupiter and Saturn, respectively, where electron intensities are not yet at radiation belt levels. We find no support that MeV electrons are dominantly accelerated by wave‐particle interactions in the magnetospheres of both planets at these distances. Instead, electron acceleration is consistent with adiabatic transport. While this is a common assumption, confirmation of this fact is important since many studies on sources, losses, and transport of energetic particles rely on it. Adiabatic heating can be driven through various radial transport mechanisms, for example, injections driven by the interchange instability or radial diffusion. We cannot distinguish these processes at Saturn with our technique. For Jupiter, we suggest that the dominating acceleration process is radial diffusion because injections are never observed at MeV energies., Key Point Electrons in Saturn's and Jupiter's magnetosphere, outside the most intense radiation belts, are accelerated through adiabatic transport

Using the Helium Oxygen Proton Electron (HOPE) and Electric Field and Waves (EFW) instruments from the Van Allen Probes, we explored the relationship between electron energy fluxes in the eV and keV ranges and spacecraft surface charging. We present statistical results on spacecraft charging within geosynchronous orbit by L and MLT. An algorithm to extract the H+ charging line in the HOPE instrument data was developed to better explore intense charging events. Also, this study explored how spacecraft potential relates to electron number density, electron pressure, electron temperature, thermal electron current, and low‐energy ion density between 1 and 210 eV. It is demonstrated that it is imperative to use both EFW potential measurements and the HOPE instrument ion charging line for examining times of extreme spacecraft charging of the Van Allen Probes. The results of this study show that elevated electron energy fluxes and high‐electron pressures are present during times of spacecraft charging but these same conditions may also occur during noncharging times. We also show noneclipse significant negative charging events on the Van Allen Probes., Key Points Statistical analysis of Van Allen Probes spacecraft chargingIntense negative charging occurs in sunlightElectron Pressure, temperature, and fluxes affect charging occurrence and intensity

Intervals of enhanced magnetic fluctuations have been frequently observed in the solar wind. But it remains an open question as to whether these waves are generated at the Sun and then transported outward by the solar wind or generated locally in the interplanetary medium. Magnetic field and plasma measurements from the Wind spacecraft under slow solar wind conditions on 19 March 2005 demonstrate seven events of enhanced magnetic fluctuations at spacecraft‐frame frequencies somewhat above the proton cyclotron frequency and propagation approximately parallel or antiparallel to the background magnetic field B o. The proton velocity distributions during these events are characterized by two components: a more dense, slower core and a less dense, faster beam. Observed plasma parameters are used in a kinetic linear dispersion equation analysis for electromagnetic fluctuations at k x B o = 0; for two events the most unstable mode is the Alfvén‐cyclotron instability driven by a proton component temperature anisotropy T⊥/T|| > 1 (where the subscripts denote directions relative to B o), and for three events the most unstable mode is the right‐hand polarized magnetosonic instability driven primarily by ion component relative flows. Thus, both types of ion anisotropies and both types of instabilities are likely to be local sources of these enhanced fluctuation events in the solar wind., Key Points Ion temperature anisotropies and proton beam/core flows are sources of enhanced field observationsFor two events Alfven‐cyclotron modes are most unstableFor three events magnetosonic modes are most unstable

A polarizing heliospheric imager is a critical next generation tool for space weather monitoring and prediction. Heliospheric imagers can track coronal mass ejections (CMEs) as they cross the solar system, using sunlight scattered by electrons in the CME. This tracking has been demonstrated to improve the forecasting of impact probability and arrival time for Earth‐directed CMEs. Polarized imaging allows locating CMEs in three dimensions from a single vantage point. Recent advances in heliospheric imaging have demonstrated that a polarized imager is feasible with current component technology.Developing this technology to a high technology readiness level is critical for space weather relevant imaging from either a near‐Earth or deep‐space mission. In this primarily technical review, we developpreliminary hardware requirements for a space weather polarizing heliospheric imager system and outline possible ways to flight qualify and ultimately deploy the technology operationally on upcoming specific missions. We consider deployment as an instrument on NOAA's Deep Space Climate Observatory follow‐on near the Sun‐Earth L1 Lagrange point, as a stand‐alone constellation of smallsats in low Earth orbit, or as an instrument located at the Sun‐Earth L5 Lagrange point. The critical first step is the demonstration of the technology, in either a science or prototype operational mission context., Key Points Polarized heliospheric imaging is a necessary next step for improved space weather predictionPolarized heliospheric imaging is feasible with current technologyPolarization is a useful augmentation to current instruments, either on or off the Earth‐Sun line

We investigate the plasmaspheric influence on the resonant mode coupling of magnetospheric ultralow frequency (ULF) waves using the Lyon‐Fedder‐Mobarry (LFM) global magnetohydrodynamic (MHD) model. We present results from two different versions of the model, both driven by the same solar wind conditions: one version that contains a plasmasphere (the LFM coupled to the Rice Convection Model, where the Gallagher plasmasphere model is also included) and another that does not (the stand‐alone LFM). We find that the inclusion of a cold, dense plasmasphere has a significant impact on the nature of the simulated ULF waves. For example, the inclusion of a plasmasphere leads to a deeper (more earthward) penetration of the compressional (azimuthal) electric field fluctuations, due to a shift in the location of the wave turning points. Consequently, the locations where the compressional electric field oscillations resonantly couple their energy into local toroidal mode field line resonances also shift earthward. We also find, in both simulations, that higher‐frequency compressional (azimuthal) electric field oscillations penetrate deeper than lower frequency oscillations. In addition, the compressional wave mode structure in the simulations is consistent with a radial standing wave oscillation pattern, characteristic of a resonant waveguide. The incorporation of a plasmasphere into the LFM global MHD model represents an advance in the state of the art in regard to ULF wave modeling with such simulations. We offer a brief discussion of the implications for radiation belt modeling techniques that use the electric and magnetic field outputs from global MHD simulations to drive particle dynamics., Key Points Magnetosphere responds as a resonant waveguide to ULF fluctuations in solar wind dynamic pressureInclusion of a plasmasphere has a substantial impact on the nature of the simulated ULF wavesInclusion of a plasmasphere leads to a deeper penetration of azimuthal electric field oscillations

A significant fraction of the energy released by magnetotail reconnection appears to go into ion heating, but this heating is generally anisotropic. We examine ARTEMIS dual‐spacecraft observations of a long‐duration magnetotail exhaust generated by antiparallel reconnection in conjunction with particle‐in‐cell simulations, showing spatial variations in the anisotropy across the outflow far (>100d i) downstream of the X line. A consistent pattern is found in both the spacecraft data and the simulations: While the total temperature across the exhaust is rather constant, near the boundaries T i,|| dominates. The plasma is well above the firehose threshold within patchy spatial regions at |B X|∈[0.1,0.5]B 0, suggesting that the drive for the instability is strong and the instability is too weak to relax the anisotropy. At the midplane ( |BX|≲0.1B0), T i,⊥>T i,|| and ions undergo Speiser‐like motion despite the large distance from the X line., Key Points Parallel temperature peaks at the edges, perpendicular at the neutral planeFirehose instability threshold is greatly exceeded, indicating strong drivingSpeiser‐like motion persists far from X line despite large fluctuations

Magnetospheric Multiscale observations are used to probe the structure and temperature profile of a guide field reconnection exhaust ~100 ion inertial lengths downstream from the X‐line in the Earth's magnetosheath. Asymmetric Hall electric and magnetic field signatures were detected, together with a density cavity confined near 1 edge of the exhaust and containing electron flow toward the X‐line. Electron holes were also detected both on the cavity edge and at the Hall magnetic field reversal. Predominantly parallel ion and electron heating was observed in the main exhaust, but within the cavity, electron cooling and enhanced parallel ion heating were found. This is explained in terms of the parallel electric field, which inhibits electron mixing within the cavity on newly reconnected field lines but accelerates ions. Consequently, guide field reconnection causes inhomogeneous changes in ion and electron temperature across the exhaust., Key Points A guide field reconnection exhaust was encountered by MMS in the magnetosheath ~100 ion inertial lengths downstream from the X‐lineA density cavity forms on one edge of the exhaust with embedded electron jetting toward the X‐line and electron holes on the cavity edgeThe parallel electric field causes electron cooling and ion heating in the cavity and inhomogeneous temperature profiles across the exhaust

Sea ice models with the traditional viscous‐plastic (VP) rheology and very small horizontal grid spacing can resolve leads and deformation rates localized along Linear Kinematic Features (LKF). In a 1 km pan‐Arctic sea ice‐ocean simulation, the small‐scale sea ice deformations are evaluated with a scaling analysis in relation to satellite observations of the Envisat Geophysical Processor System (EGPS) in the Central Arctic. A new coupled scaling analysis for data on Eulerian grids is used to determine the spatial and temporal scaling and the coupling between temporal and spatial scales. The spatial scaling of the modeled sea ice deformation implies multifractality. It is also coupled to temporal scales and varies realistically by region and season. The agreement of the spatial scaling with satellite observations challenges previous results with VP models at coarser resolution, which did not reproduce the observed scaling. The temporal scaling analysis shows that the VP model, as configured in this 1 km simulation, does not fully resolve the intermittency of sea ice deformation that is observed in satellite data., Key Points Sea ice deformation localizes along lines in high‐resolution viscous‐plastic sea ice modelsThe model reproduces spatial scaling properties also observed in satellite dataThe model underestimates temporal scaling compared to satellite data

We report a 511‐keV photon flux enhancement that was observed inside a thundercloud and is a result of positron annihilation. The observation was made with the In‐flight Lightning Damage Assessment System (ILDAS) on board of an A340 test aircraft. The aircraft was intentionally flying through a thunderstorm at 12‐km altitude over Northern Australia in January 2016. Two gamma ray detectors showed a significant count rate increase synchronously with fast electromagnetic field variations registered by an on‐board antenna. A sequence of 10 gamma ray enhancements was detected, each lasted for about 1 s. Their spectrum mainly consists of 511‐keV photons and their Compton component. The local electric activity during the emission was identified as a series of static discharges of the aircraft. A full‐scale Geant4 model of the aircraft was created to estimate the emission area. Monte Carlo simulation indicated that the positrons annihilated in direct vicinity or in the aircraft body., Key Points A sequence of 10 X‐ray pulses was observed from an aircraft inside thundercloud at 12‐km altitudeThe pulses contain strong enhancement of the positron annihilation line (511 keV)Each pulse lasts for 1 s and is accompanied by a static discharge from the aircraft

We present a comprehensive statistical analysis of mirror mode waves and the properties of their plasma surroundings in sheath regions driven by interplanetary coronal mass ejection (ICME). We have constructed a semi-automated method to identify mirror modes from the magnetic field data. We analyze 91 ICME sheath regions from January 1997 to April 2015 using data from the Wind spacecraft. The results imply that similarly to planetary magnetosheaths, mirror modes are also common structures in ICME sheaths. However, they occur almost exclusively as dip-like structures and in mirror stable plasma. We observe mirror modes throughout the sheath, from the bow shock to the ICME leading edge, but their amplitudes are largest closest to the shock. We also find that the shock strength (measured by Alfvén Mach number) is the most important parameter in controlling the occurrence of mirror modes. Our findings suggest that in ICME sheaths the dominant source of free energy for mirror mode generation is the shock compression. We also suggest that mirror modes that are found deeper in the sheath are remnants from earlier times of the sheath evolution, generated also in the vicinity of the shock. Keywords. Interplanetary physics (plasma waves and turbulence; solar wind plasma) – space plasma physics (waves and instabilities)

Fast plasma flows produced as outflow jets from reconnection sites or X lines are a key feature of the dynamics in the Earth's magnetosphere. We have used a polar plane simulation of the hybrid-Vlasov model Vlasiator, driven by steady southward interplanetary magnetic field and fast solar wind, to study fast plasma sheet ion flows and related magnetic field structures in the Earth's magnetotail. In the simulation, lobe reconnection starts to produce fast flows after the increasing pressure in the lobes has caused the plasma sheet to thin sufficiently. The characteristics of the earthward and tailward fast flows and embedded magnetic field structures produced by multi-point tail reconnection are in general agreement with spacecraft measurements reported in the literature. The structuring of the flows is caused by internal processes: interactions between major X points determine the earthward or tailward direction of the flow, while interactions between minor X points, associated with leading edges of magnetic islands carried by the flow, induce local minima and maxima in the flow speed. Earthward moving flows are stopped and diverted duskward in an oscillatory (bouncing) manner at the transition region between tail-like and dipolar magnetic fields. Increasing and decreasing dynamic pressure of the flows causes the transition region to shift earthward and tailward, respectively. The leading edge of the train of earthward flow bursts is associated with an earthward propagating dipolarization front, while the leading edge of the train of tailward flow bursts is associated with a tailward propagating plasmoid. The impact of the dipolarization front with the dipole field causes magnetic field variations in the Pi2 range. Major X points can move either earthward or tailward, although tailward motion is more common. They are generally not advected by the ambient flow. Instead, their velocity is better described by local parameters, such that an X point moves in the direction of increasing reconnection electric field strength. Our results indicate that ion kinetics might be sufficient to describe the behavior of plasma sheet bulk ion flows produced by tail reconnection in global near-Earth simulations. Keywords. Magnetospheric physics (magnetospheric configuration and dynamics; plasma sheet) – space plasma physics (numerical simulation studies)

著者人数: 19名, Accepted: 2018-01-07, 資料番号: SA1170245000

Observations from Magnetospheric MultiScale (~8 Re) and Van Allen Probes (~5 and 4 Re) show that the initial dayside response to a small interplanetary shock is a double‐peaked dawnward electric field, which is distinctly different from the usual bipolar (dawnward and then duskward) signature reported for large shocks. The associated E × B flow is radially inward. The shock compressed the magnetopause to inside 8 Re, as observed by Magnetospheric MultiScale (MMS), with a speed that is comparable to the E × B flow. The magnetopause speed and the E × B speeds were significantly less than the propagation speed of the pulse from MMS to the Van Allen Probes and GOES‐13, which is consistent with the MHD fast mode. There were increased fluxes of energetic electrons up to several MeV. Signatures of drift echoes and response to ULF waves also were seen. These observations demonstrate that even very weak shocks can have significant impact on the radiation belts., Key Points Magnetospheric response to small shocks not typical bipolar electric field pulse seen for large shocksPropagation of compression on VAP, MMS, and GOES‐13 is consistent with fast mode speed; large‐scale structure is not fast modeElectric field pulse associated with small shock energizes electrons to >50 keV; flux enhancements seen up to several MeV