Inversion of Meteor Rayleigh Waves on Earth and Modeling of Air Coupled Rayleigh Waves on Mars

Meteor impacts and/or meteor events generate body and surface seismic waves on the surface of a planet. When meteoroids burst in the atmosphere, they generate shock waves that subsequently convert into acoustic waves in the atmosphere and seismic waves in the ground. This effect can be modeled as the amplitude of Rayleigh and other Spheroidal modes excitation, due to atmospheric/ground coupling effects.First, an inversion of the seismic source of Chelyabinsk superbolide is performed. We develop an approach in order to model a line source in the atmosphere, corresponding to the consecutive generation of shock waves by the interaction with the atmosphere. The model is based on the known trajectory. We calculate the synthetic seismograms of Rayleigh waves associated with the event by the summation of normal modes of a model of the solid part and the atmosphere of the planet. Through an inversion technique based on singular value decomposition, we perform a full Rayleigh wave inversion and we provide solutions for the moment magnitude.SEIS will likely detect seismic waves generated by impacts and the later might be further located by remote sensing differential processing. In the case of Mars, we use the same method to obtain waveforms associated with impacts on the planetary surface or in low altitudes in the Martian atmosphere. We show that the contribution of the fundamental spheroidal solid mode is dominating the waveforms, compared to that of the first two overtones. We perform an amplitude comparison and we show that small impactors (diameter of 0.5 to 2 m), can be detected by the SEIS VBB seismometer of InSight mission, even in short epicentral distances, in the higher frequencies of the Rayleigh waves. We perform an analysis based on impact rate estimations and we calculate the number of detectable events of 1 meter diameter meteor impacts to be 6.7 to 13.4 per 1 Martian year for a \(Q=500\).     

Evolution of the Sunspot Number and Solar Wind $$ Time Series

The past two decades have witnessed significant changes in our knowledge of long-term solar and solar wind activity. The sunspot number time series (1700-present) developed by Rudolf Wolf during the second half of the 19th century was revised and extended by the group sunspot number series (1610–1995) of Hoyt and Schatten during the 1990s. The group sunspot number is significantly lower than the Wolf series before ～1885. An effort from 2011–2015 to understand and remove differences between these two series via a series of workshops had the unintended consequence of prompting several alternative constructions of the sunspot number. Thus it has been necessary to expand and extend the sunspot number reconciliation process. On the solar wind side, after a decade of controversy, an ISSI International Team used geomagnetic and sunspot data to obtain a high-confidence time series of the solar wind magnetic field strength (\(B\)) from 1750-present that can be compared with two independent long-term (> ～600 year) series of annual \(B\)-values based on cosmogenic nuclides. In this paper, we trace the twists and turns leading to our current understanding of long-term solar and solar wind activity.     

Single Degenerate Models for Type Ia Supernovae: Progenitor’s Evolution and Nucleosynthesis Yields

We review how the single degenerate models for Type Ia supernovae (SNe Ia) works. In the binary star system of a white dwarf (WD) and its non-degenerate companion star, the WD accretes either hydrogen-rich matter or helium and undergoes hydrogen and helium shell-burning. We summarize how the stability and non-linear behavior of such shell-burning depend on the accretion rate and the WD mass and how the WD blows strong wind. We identify the following evolutionary routes for the accreting WD to trigger a thermonuclear explosion. Typically, the accretion rate is quite high in the early stage and gradually decreases as a result of mass transfer. With decreasing rate, the WD evolves as follows: (1) At a rapid accretion phase, the WD increase its mass by stable H burning and blows a strong wind to keep its moderate radius. The wind is strong enough to strip a part of the companion star’s envelope to control the accretion rate and forms circumstellar matter (CSM). If the WD explodes within CSM, it is observed as an “SN Ia-CSM”. (X-rays emitted by the WD are absorbed by CSM.) (2) If the WD continues to accrete at a lower rate, the wind stops and an SN Ia is triggered under steady-stable H shell-burning, which is observed as a super-soft X-ray source: “SN Ia-SSXS”. (3) If the accretion continues at a still lower rate, H shell-burning becomes unstable and many flashes recur. The WD undergoes recurrent nova (RN) whose mass ejection is smaller than the accreted matter. Then the WD evolves to an “SN Ia-RN”. (4) If the companion is a He star (or a He WD), the accretion of He can trigger He and C double detonations at the sub-Chandrasekhar mass or the WD grows to the Chandrasekhar mass while producing a He-wind: “SN Ia-He CSM”. (5) If the accreting WD rotates quite rapidly, the WD mass can exceed the Chandrasekhar mass of the spherical WD, which delays the trigger of an SN Ia. After angular momentum is lost from the WD, the (super-Chandra) WD contracts to become a delayed SN Ia. The companion star has become a He WD and CSM has disappeared: “SN Ia-He WD”. We update nucleosynthesis yields of the carbon deflagration model W7, delayed detonation model WDD2, and the sub-Chandrasekhar mass model to provide some constraints on the yields (such as Mn) from the comparison with the observations. We note the important metallicity effects on ⁵⁸Ni and ⁵⁵Mn.     

Crustal S-Wave Velocity from Apparent Incidence Angles: A Case Study in Preparation for InSight

Retrieval of crustal structure and thickness of Mars is among the main goals of InSight. Here we investigate which constraints on the crust at the landing site can be provided by apparent P-wave incidence angles derived from P-receiver functions. We consider receiver functions for six different Mars models, calculated from synthetic seismograms generated via Instaseis from the Green’s function databases of the Marsquake Service, in detail. To allow for a larger range of crustal thicknesses and structures, we additionally analyze data from five broad-band stations across Central Europe. We find that the likely usable epicentral distance range for P-wave receiver functions on Mars lies between \(35^{\circ}\) and the core shadow, and can be extended to more than \(150^{\circ}\) by also using the PP-phase. Comparison to models for the spatial distribution of Martian seismicity indicates that sufficient seismicity should occur within the P-wave distance range around InSight within the nominal mission duration to allow for the application of our method. Apparent P-wave incidence angles are derived from the amplitudes of vertical and radial receiver functions at the P-wave onset within a range of period bands, up to 120 s. The apparent incidence angles are directly related to apparent S-wave velocities, which are inverted for the subsurface S-wave velocity structure via a grid search. The veracity of the forward calculated receiver functions and apparent S-wave velocities is ensured by benchmarking various algorithms against the Instaseis synthetics. Results indicate that apparent S-wave velocity curves provide valuable constraints on crustal thickness and structure, even without any additional constraints, and considering the location uncertainty and limited data quantity of InSight. S-wave velocities in the upper half of the crust are constrained best, but if reliable measurements at long periods are available, the curves also provide constraints down to the uppermost mantle. Besides, it is demonstrated that the apparent velocity curves can differentiate between crustal velocity models that are indistinguishable by other methods.     

The Origin,Early Evolution and Predictability of Solar Eruptions

Coronal mass ejections (CMEs) were discovered in the early 1970s when space-borne coronagraphs revealed that eruptions of plasma are ejected from the Sun. Today, it is known that the Sun produces eruptive flares, filament eruptions, coronal mass ejections and failed eruptions; all thought to be due to a release of energy stored in the coronal magnetic field during its drastic reconfiguration. This review discusses the observations and physical mechanisms behind this eruptive activity, with a view to making an assessment of the current capability of forecasting these events for space weather risk and impact mitigation. Whilst a wealth of observations exist, and detailed models have been developed, there still exists a need to draw these approaches together. In particular more realistic models are encouraged in order to asses the full range of complexity of the solar atmosphere and the criteria for which an eruption is formed. From the observational side, a more detailed understanding of the role of photospheric flows and reconnection is needed in order to identify the evolutionary path that ultimately means a magnetic structure will erupt.     

Venus: The Atmosphere,Climate, Surface,Interior and Near-Space Environment of an Earth-Like Planet

This is a review of current knowledge about Earth’s nearest planetary neighbour and near twin, Venus. Such knowledge has recently been extended by the European Venus Express and the Japanese Akatsuki spacecraft in orbit around the planet; these missions and their achievements are concisely described in the first part of the review, along with a summary of previous Venus observations. The scientific discussions which follow are divided into three main sections: on the surface and interior; the atmosphere and climate; and the thermosphere, exosphere and magnetosphere. These reports are intended to provide an overview for the general reader, and also an introduction to the more detailed topical surveys in the following articles in this issue, where full references to original material may be found.     

Towards a Global Unified Model of Europa’s Tenuous Atmosphere

Despite the numerous modeling efforts of the past, our knowledge on the radiation-induced physical and chemical processes in Europa’s tenuous atmosphere and on the exchange of material between the moon’s surface and Jupiter’s magnetosphere remains limited. In lack of an adequate number of in situ observations, the existence of a wide variety of models based on different scenarios and considerations has resulted in a fragmentary understanding of the interactions of the magnetospheric ion population with both the moon’s icy surface and neutral gas envelope. Models show large discrepancy in the source and loss rates of the different constituents as well as in the determination of the spatial distribution of the atmosphere and its variation with time. The existence of several models based on very different approaches highlights the need of a detailed comparison among them with the final goal of developing a unified model of Europa’s tenuous atmosphere. The availability to the science community of such a model could be of particular interest in view of the planning of the future mission observations (e.g., ESA’s JUpiter ICy moons Explorer (JUICE) mission, and NASA’s Europa Clipper mission). We review the existing models of Europa’s tenuous atmosphere and discuss each of their derived characteristics of the neutral environment. We also discuss discrepancies among different models and the assumptions of the plasma environment in the vicinity of Europa. A summary of the existing observations of both the neutral and the plasma environments at Europa is also presented. The characteristics of a global unified model of the tenuous atmosphere are, then, discussed. Finally, we identify needed future experimental work in laboratories and propose some suitable observation strategies for upcoming missions.     

A Hitch-hiker’s Guide to Stochastic Differential Equations

In this review, an overview of the recent history of stochastic differential equations (SDEs) in application to particle transport problems in space physics and astrophysics is given. The aim is to present a helpful working guide to the literature and at the same time introduce key principles of the SDE approach via “toy models”. Using these examples, we hope to provide an easy way for newcomers to the field to use such methods in their own research. Aspects covered are the solar modulation of cosmic rays, diffusive shock acceleration, galactic cosmic ray propagation and solar energetic particle transport. We believe that the SDE method, due to its simplicity and computational efficiency on modern computer architectures, will be of significant relevance in energetic particle studies in the years to come.     

Daytime Ionosphere Retrieval Algorithm for the Ionospheric Connection Explorer (ICON)

The NASA Ionospheric Connection Explorer Extreme Ultraviolet spectrograph, ICON EUV, will measure altitude profiles of the daytime extreme-ultraviolet (EUV) OII emission near 83.4 and 61.7 nm that are used to determine density profiles and state parameters of the ionosphere. This paper describes the algorithm concept and approach to inverting these measured OII emission profiles to derive the associated \(\mathrm{O}^{+}\) density profile from 150–450 km as a proxy for the electron content in the F-region of the ionosphere. The algorithm incorporates a bias evaluation and feedback step, developed at the U.S. Naval Research Laboratory using data from the Special Sensor Ultraviolet Limb Imager (SSULI) and the Remote Atmospheric and Ionospheric Detection System (RAIDS) missions, that is able to effectively mitigate the effects of systematic instrument calibration errors and inaccuracies in the original photon source within the forward model. Results are presented from end-to-end simulations that convolved simulated airglow profiles with the expected instrument measurement response to produce profiles that were inverted with the algorithm to return data products for comparison to truth. Simulations of measurements over a representative ICON orbit show the algorithm is able to reproduce hmF2 values to better than 5 km accuracy, and NmF2 to better than 12% accuracy over a 12-second integration, and demonstrate that the ICON EUV instrument and daytime ionosphere algorithm can meet the ICON science objectives which require 20 km vertical resolution in hmF2 and 18% precision in NmF2.