Feed aggregator

A mini-neutron monitor (MNM) was installed at the German Antarctic Neumayer III station, measuring the variation of galactic cosmic rays and searching for Forbush Decreases (FDs) caused by solar activities. Running continuously from 2014 until the end of 2017, the long-term stability of the detector could be investigated. After correcting the air pressure and normalization to the 27 days running mean averages of the SANAE and TERA Neutron Monitors (NMs), the daily running mean count rates are compared with the SANAE and TERA NMs also installed in Antarctica. For most of the 14 FDs with magnitudes greater than 3, taken from the list compiled by the IZMIRAN group (http://spaceweather.izmiran.ru/eng/dbs.html), the three detectors show consistent particle flux variation, although the average rate of the MNM is more than a hundred times smaller. The light and low-cost MNM is an ideal alternative to heavy and old NMs, especially at high altitudes and remote environments.

In this paper, we study Titan's magnetotail using Cassini data from the T122-T126 flybys. These consecutive flybys had a similar flyby geometry and occurred at similar Saturn magnetospheric conditions, enabling an analysis of the magnetotail's structure. Using measurements from Cassini's magnetometer (MAG) and Radio and Plasma Wave System/Langmuir probe (RPWS/LP) we identify several features consistent with reported findings from earlier flybys, for example, T9, T63 and T75. We find that the so-called ’split’ signature of the magnetotail becomes more prominent at distances of at least 3,260 km (1.3 R T ) downstream of Titan. We also identify a specific signature of the sub-alfvenic interaction of Titan with Saturn, the Alfvén wings, which are observed during the T123 and T124 flyby. A coordinate transformation is applied to mitigate variations in the upstream magnetic field, and all the flybys are projected into a new reference frame—aligned to the background magnetic field reference frame (BFA). We show that Titan's magnetotail is confined to a narrow region of around ∼4 R T Y BFA . Finally, we analyze the general draping pattern in Titan's magnetotail throughout the TA to T126 flybys.

The SEM-2 (Space Environment Monitor-2) instrument embedded on the NOAA-15 Low Earth Orbit satellite provides measurements of trapped protons in the Van Allen inner belt from 1998 to nowadays. This continuous set of measurements enables us to study the dynamics of the South Atlantic Anomaly (SAA) over more than two solar cycles, particularly, its temporal evolution. We observe that the area of the SAA is anti-correlated with the solar activity. Two physical processes explain this anticorrelation. First, the more the Sun is active the more it disables the cosmic rays to reach the Earth's magnetosphere and fill in the inner radiation belt with protons. Then, when the Sun is more active, the upper atmosphere is warmer and therefore absorbs more protons from the radiation belt. Then, we investigate the protons flux centroid of the SAA. The temporal evolution of its position, latitude and, longitude is studied over the same time interval (1998–2022). We notice that the latitude of the centroid is also anti-correlated with the solar activity whereas the longitude seems absolutely independent. The temporal evolution of the position of the centroid shows a drift of the SAA. Indeed from 1998 to 2022 the SAA drifted of about 7° West. The SEM-2 instrument measures flux for protons of different energies (16, 36, 70, and, 140 MeV). For each energy, the SAA dynamic has a similar trend but with different values. These differences are investigated and the results are discussed.

Perchlorate in the environment originates from both natural and anthropogenic sources. A previous study of a 300-year Greenland ice core perchlorate record found that anthropogenic impact on environmental perchlorate became significant starting around 1980, while natural formation is the only significant source of environmental perchlorate prior to that. The study also found increased perchlorate deposition in the Arctic following certain volcanic eruptions and suggested that at least some volcanic eruptions could enhance natural perchlorate production. Here we compare the perchlorate record with the volcanic record from sulfate in the same Greenland ice core and find that only stratospheric eruptions—large eruptions injecting volcanic substances directly into the stratosphere—enhance perchlorate production. No contribution to naturally formed perchlorate is detected from non-stratospheric eruptions. The high-resolution ice core perchlorate data are used to quantify contributions from volcanic eruptions, non-volcanic natural processes, as well as from human activities during different periods. For the location in the Arctic in the perchlorate Pre-Anthropogenic Era (1701–1979), the magnitude (0.26 μg m−2 yr−1 on average) of perchlorate produced during sporadic stratospheric eruptions is comparable to that (0.23 μg m−2 yr−1) produced by non-volcanic natural processes. In the Anthropogenic Era (1980–2006), the magnitude of both the volcanic and non-volcanic natural perchlorate production is similar to the enhancement (0.29 μg m−2 yr−1) by human activities.

Previous studies have suggested possible connections between the decreasing Arctic sea-ice and long-duration (>5 days, LD) cold weather events in Eurasia and North America. Here we document the occurrences of weather regimes in winter by their durations, based on the empirical orthogonal function analyses of the daily geopotential height fields at 500 hPa (z500) for the months of November–March 1979–2019. Significant changes in the occurrence frequency and persistence of Ural ridge (UR) and weak stratospheric polar vortex (PV) were found between winters following high and low autumn sea-ice covers (SIC) in the Barents and Kara seas. It is shown that a strengthening of the UR is accompanied with a weakening of the PV, and a weak PV favors Greenland ridge (GR). Cold spells in East Asia persist for 5 more days after an LDUR. Cold spells from Canada to the U.S. occur 2–5 days after an LD Ural trough (UT) and are associated with a z500 anomaly dipole centered over Alaska (+) and Hudson Bay (−). Cold spells in the eastern U.S. occur 1–4 days after an LDGR due to circulations resembling the Pacific-North America pattern. Increased occurrences of UR in winter are associated with a decreased eastward propagation of synoptic waves from the North Atlantic to Japan and the North Pacific.

As the number of satellites on orbit grows it is increasingly important to understand their operating environment. Physics-based models can simulate the behavior of the Earth's radiation belts by solving a Fokker-Planck equation. Three-dimensional models use diffusion coefficients to represent the interactions between electromagnetic waves and the electrons. One-dimensional radial diffusion models neglect the effects of energy diffusion and represent the losses due to the waves with a loss timescale. Both approaches may use pitch angle distributions (PADs) to create boundary conditions, to map observations from low to high equatorial pitch angles and to calculate phase-space density from observations. We present a comprehensive set of consistent PADs and loss timescales for 2 ≤ L* ≤ 7, 100 keV ≤ E ≤ 5 MeV and all levels of geomagnetic activity determined by the Kp index. These are calculated from drift-averaged diffusion coefficients that represent all the VLF waves that typically interact with radiation belt electrons and show good agreement with data. The contribution of individual waves is demonstrated; magnetosonic waves have little effect on loss timescales when lightning-generated whistlers are present, and chorus waves contribute to loss even in low levels of geomagnetic activity. The PADs vary in shape depending on the dominant waves. When chorus is dominant the distributions have little activity dependence, unlike the corresponding loss timescales. Distributions peaked near 90° are formed by plasmaspheric hiss for L* ≤ 3 and E < 1 MeV, and by EMIC waves for L* > 3 and E > 1 MeV. When hiss dominates, increasing activity broadens the distribution but when EMIC waves dominate increasing activity narrows the distribution.

Positive leaders branch less frequently than negative counterpart, and the physical processes and properties of positive leader branching remain a mystery. We investigated 10 m laboratory discharges under four positive voltages using a high-speed video camera. Positive leaders differ from negative leaders by either directly splitting or connecting with floating bidirectional leaders to form branching, and the number of leader branches shows a positive correlation with the applied voltage, that is, the branched channels increased from 1 to 4 when the voltage increased by a factor of 1.5. Grounding points are positioned beneath the electrode and are more concentrated with lower voltage. During the stable progression of the leader, there is a slight increase in its development speed as the applied voltage rises. When the voltage is increased by 70%, the average breakdown time decreases by 40%. These characteristics provide insights into the branching mechanism of positive leaders.

We report results from a global simulation of the September 2017 geomagnetic storm. The global model comprises the ionospheric code SAMI3 and the atmosphere/thermosphere code WACCM-X. We show that a train of large-scale EPBs form in the Pacific sector during the storm recovery phase on 8 September 2017. The EPBs are associated with storm-induced modification of the zonal and meridional winds. These changes lead to an eastward electric field which in turn causes an upward E × B drift in the post-midnight sector. A large decrease in the Pedersen conductance caused by meridional equatorward winds leads to an increase in the growth rate of the generalized Rayleigh-Taylor instability that causes EPBs to develop. Interestingly, several EPBs reach altitudes above 3,000 km.

Damage zones are common around faults, but their effects on earthquake mechanics are still incompletely understood. Here, we investigate how damage affects rupture patterns, source time functions (STF) and ground motions in 2D fully-dynamic cycle models. We find that back-propagating rupture fronts emerge in large faults and can be triggered by residual stresses left by previous ruptures or by damage-induced pulse-to-crack transitions. Damage-induced back-propagating fronts are modulated by slip rate oscillations, amplify high-frequency radiation, and sharpen the multiple peaks in STF even in the absence of frictional heterogeneity or fault segmentation. Near-field ground motion is predominantly controlled by stress heterogeneity left by prior seismicity, and further amplified within the damage zone by trapped waves and outside it by secondary rupture fronts. This study refines our knowledge on damage zone effects on earthquake rupture and identifies their potentially observable signatures in the near and far field.

A devastating earthquake with moment magnitude 7.5 occurred in the Noto Peninsula in central Japan on 1 January 2024. We estimate the rupture evolution of this earthquake from teleseismic P-wave data using the potency-density tensor inversion method, which provides information on the spatiotemporal slip distribution including fault orientations. The results show a long and quiet initial rupture phase that overlaps with regions of preceding earthquake swarms and associated aseismic deformation. The following three major rupture episodes evolve on segmented, differently oriented faults bounded by the initial rupture region. The irregular initial rupture process followed by the multi-scale rupture growth is considered to be controlled by the preceding seismic and aseismic processes and the geometric complexity of the fault system. Such a discrete rupture scenario, including the triggering of an isolated fault rupture, adds critical inputs on the assessment of strong ground motion and associated damages for future earthquakes.

Using a two-way coupled magnetohydrodynamics with embedded kinetic physics model, we perform a substorm event simulation to study electron velocity distribution functions (VDFs) evolution associated with Bursty Bulk Flows (BBFs). The substorm was observed by Magnetospheric Multiscale satellite on 16 May 2017. The simulated BBF macroscopic characteristics and electron VDFs agree well with observations. The VDFs from the BBF tail to its dipolarization front (DF) during its earthward propagation are revealed and they show clear energization and heating. The electron pitch angle distributions (PADs) at the DF are also tracked, which show interesting energy dependent features. Lower energy electrons develop a “two-hump” PAD while the higher energy ones show persist “pancake” distribution. Our study reveals for the first time the evolution of electron VDFs as a BBF moves earthward using a two-way coupled global and kinetic model, and provides valuable contextual understanding for the interpretation of satellite observations.

Whakaari/White Island has been the most active New Zealand volcano in the 21st century, producing small phreatic and phreatomagmatic eruptions, which are hard to predict. The most recent eruption occurred in 2019, tragically claiming the lives of 22 individuals and causing numerous injuries. We employed shear-wave splitting analyses to investigate variations in anisotropy between 2018 and 2020, during quiescence, unrest, and the eruption. We examined spatial and temporal variations in 3,499 shear-wave splitting and 2,656 V p /V s ratio measurements. Comparing shear-wave splitting parameters from similar earthquake paths across different times indicates that the observed temporal changes are unlikely to result from variations in earthquake paths through media with spatial variability. Instead, these changes may stem from variations in anisotropy over time, likely caused by changes in crack alignment due to stress or varying fluid content.

Mixed-phase stratiform clouds contain numerous liquid, mixed-phase, and ice clusters, quantifying the cluster scales is potentially helpful to improve the parameterizations of microphysics and radiation models. However, the scales of hydrometeor clusters at different levels of stratiform clouds are not well understood. In this study, using airborne measurements and a large eddy simulation, we show that turbulence plays an important role in controlling the clusters with length of a few hundred meters, while the scales of larger clusters have stronger vertical variations from cloud base to top. The liquid clusters are the largest near the cloud top, while the lengths of ice clusters decrease from cloud base to top. The lengths of mixed-phase clusters depend on the glaciation process, a faster glaciation results in smaller mixed-phase clusters. The results improve our understanding on how the liquid and ice are mixed at different levels in stratiform clouds.

The impact of an episodic river flood is intimately linked to its duration. Yet it is still unclear how often should a river be observed to accurately determine the occurrence and duration of extreme events. Here we assess flow statistics along with peak flow event detection and duration as a function of the discharge sampling period for large tributaries of the Mississippi basin using hourly gages over 2010–2022. Median event durations above high quantiles spatially vary from around 2 days upstream to 30 days downstream. Discharge mean, standard deviation, and quantiles can all be estimated within 2.5% error for sampling periods up to 8 days. A minimum temporal sampling 4× (2×) finer than peak flow event median duration is required to detect 95 ± 3% (85 ± 5%) of events and to estimate their duration within 90 ± 5% (75 ± 10%) median accuracy. Our findings have direct implications for future satellite missions concerned with capturing flood events.

Seismic anisotropy is a powerful tool to map deformation processes in the deep Earth. Below 660 km, however, observations are scarce and conflicting. In addition, the underlying crystal scale mechanisms, leading to microstructures and crystal orientations, remain poorly constrained. Here, we use multigrain X-ray diffraction in the laser-heated diamond anvil cell to investigate the orientations of hundreds of grains in pyrolite, a model composition of the Earth's mantle, at in situ pressure and temperature. Bridgmanite in pyrolite exhibits three regimes of microstructures, due to transformation and deformation at low and high pressure. These microstructures result in predictions of 1.5%–2% shear wave splitting between 660 and 2,000 km with reversals in fast S-wave polarization direction at about 1,300 km depth. Anisotropy can develop in pyrolite at lower mantle conditions, but pressure has a significant impact on the plastic behavior of bridgmanite, and hence seismic observations, which may explain conflicting anisotropy observations.

The processing of hundreds of Synthetic Aperture Radar (SAR) images acquired by two satellite systems: Sentinel-1 and COSMO-SkyMed reveals a decade of ground deformation for a ∼0.5 km diameter area around the summit crater of the only active carbonatitic volcano on Earth: Ol Doinyo Lengai in Tanzania. Further decomposing ascending and descending orbits when the appropriate SAR data sets overlap allow us to interpret the imaged deformation as ground subsidence with a significant rate of ∼3.6 cm/yr for the pixels located just north of the summit crater. Using geodetic modeling and inverting the highest spatial resolution COSMO-SkyMed data set, we show that the mechanism explaining this subsidence is most likely a deflating very shallow (≤1 km depth below the summit crater at the 95% confidence level) magma reservoir, consistent with geochemical-petrological and seismo-acoustic studies.

This paper presents the first ever observations on aspect-sensitive characteristics of 205 MHz stratosphere–troposphere (ST) radar located at a tropical station Cochin (10.04°N, 76.3°E) using volume scanning. The most significant and new observation is that the signal-to-noise ratio in zenith and off-zenith beams are nearly equal in some height region, indicating the presence of isotropic turbulence. Signal strength decreases by 0.75 dB per degree from 0 to 10 degree off-zenith, 0.9 dB per degree from 10 to 20 degree off-zenith and 0.3 dB per degree beyond 20 degree off-zenith. Different causative mechanisms are discussed on the basis of various estimated parameters associated with aspect sensitivity. Maximum aspect sensitivity is observed between 12 and 17 km, indicating the presence of dynamic instability arising due to strong wind shear and atmospheric stability. When both the square of wind shear and stability parameters are above 0.25 × 10−3 s−2, the scatterers become mostly isotropic. The study also shows a power difference in the symmetric beams as well as azimuth angle dependency. Analysis suggests that this asymmetry is due to the tilting of layers by the action of atmospheric gravity waves generated through Kelvin-Helmholtz instability. The present configuration of radar can provide a better understanding of three-dimensional structures of turbulence and instabilities.

Recent observations beneath Kanto, Japan have shown that seismic activity and seismic attenuation within the overlying continental plate change with time due to drainage caused by slow-slip events (SSEs) along the upper boundary of the Philippine Sea plate. However, associated changes in rock properties have not been investigated. In this study, we estimate frequency-dependent shear-wave anisotropy to provide a detailed insight into the structural change associated with drainage. We perform shear-wave splitting analysis in frequency ranges of 1–4, 2–6, and 4–8 Hz for 306 earthquakes that occur during September 2009–August 2021 and recorded at the Metropolitan Seismic Observation network. Obtained time differences between fast and slow S waves (delay time) range from almost zero to 0.16–0.18 s, exhibiting spatio-temporal variation and frequency dependence. The fast S-wave polarization directions are generally consistent with the direction of the maximum horizontal compressional axis in the study region, which suggests that the observed anisotropy is probably caused by the NE–SW-oriented fractures developed under the regional stress field. The temporal variation in delay times is correlated with SSEs activity with a lag time of 0.0–0.1 year. Furthermore, comparisons between observed frequency-dependent delay times and numerical calculation of fracture-induced anisotropy suggest that the average fracture radius is almost constant (0.30–0.35 m) over time but fracture density temporally varies from 0.025 to 0.035. We infer that the fracture density is probably enhanced by opening of the NE–SW-oriented fractures during the upward migration of fluids that are expelled from the plate interface.

To gain a deeper understanding of the extensive and varied lithospheric deformations beneath northern Oman, we examine seismic anisotropy in this region using splitting analysis of teleseismic shear wave data. Our study utilizes data from a dense network consisting of 13 permanent and 45 temporary seismic stations, which were operational for approximately 2.5 years starting from October 2013. By examining the azimuthal distribution of shear wave splitting (SWS) parameters, we were able to divide the study area into three sub-regions. The stations located to the west of the Hawasina window exhibit relatively azimuthally invariant SWS parameters suggesting a single anisotropic layer. On the other hand, most of the stations located in the central and eastern regions display variations versus back-azimuth, indicating the potential presence of depth-dependent anisotropy. The General NW-SE trend of the Fast Polarization Directions (FPDs) of the one-layer anisotropy in the west and FPDs of the upper layers in the east is concordant with the strike of the structures resulting from the collision between the continental and oceanic plates. A clear contrast in SWS parameters is observed in the Semail Gap Fault Zone (SGFZ), suggesting that the SGFZ can be a lithospheric-scale structure that hampers the intrusion of mafic magma from the southeast. Furthermore, the FPDs of the lower layer in the east exhibit an NE-SW trend, which may be indicative of the large-scale mantle flow resulting from the present-day plate motion.

We developed an enhanced Kalman-based approach to quantify abrupt changes and significant non-linearity in vertical land motion (VLM) along the coast of Chile and the Antarctic Peninsula using a combination of multi-mission satellite altimetry (ALT), tide gauge (TG), and GPS data starting from the early 1990s. The data reveal the spatial variability of co-seismic and post-seismic subsidence at TGs along the Chilean subduction zone in response to the Mw8.8 Maule 2010, Mw8.1 Iquique 2014, and Mw8.3 Illapel 2015 earthquakes that are not retrievable from the interpolation of sparse GPS observations across space and time. In the Antarctic Peninsula, where continuous GPS data do not commence until ∼1998, the approach provides new insight into the ∼2002 change in VLM at the TGs of +5.3 ± 2.2 mm/yr (Palmer) and +3.5 ± 2.8 mm/yr (Vernadsky) due to the onset of ice-mass loss following the Larsen-B Ice Shelf breakup. We used these data to constrain viscoelastic Earth model parameters for the northern Antarctic Peninsula, obtaining a preferred lithosphere thickness of 115 km and upper mantle viscosity of 0.9 × 1018 Pa s. Our estimates of regionally-correlated ALT systematic errors are small, typically between ∼±0.5–2.5 mm/yr over single-mission time scales. These are consistent with competing orbit differences and the relative errors apparent in ALT crossovers. This study demonstrates that, with careful tuning, the ALT-TG technique can provide improved temporal and spatial sampling of VLM, yielding new constraints on geodynamic models and assisting sea-level change studies in otherwise data sparse regions and periods.