GRL

On 18 November 2023, SpaceX launched the Starship, the tallest and the most powerful rocket ever built. The Super Heavy engine separated from the Starship spacecraft and exploded at 90 km of altitude, while the main core Starship continued to rise up to 149 km and exploded after ∼8 min of flight. In this work, we used data from ground-based GNSS receivers and we analyzed total electron content (TEC) response to the Starship flight and the two explosions. For the first time, we observed large-distance northward propagation of intensive 2,000 km V-shaped ionospheric disturbances from the rocket trajectory. The observed perturbations, most likely, represent shock waves propagating with the cone angle of ∼14° on the North and ∼7° on the South against the flight track that corresponds to the Mach angle of the shock waves in the lower atmosphere. The Starship explosion also produced a non-chemical depletion in the ionospheric TEC.

The Arctic rapidly warms and sea ice retreats, a large fraction of organic carbon (OC), currently stored in coastal permafrost will be released into the marine system. Once reintroduced into the active carbon cycle, this material will either be decomposed or buried on the shelf depending on its hydrodynamic and chemical properties. Currently, carbon estimates are based on bulk measurements, which does not take the hydrodynamic pathway of different fractions into account. Therefore, eight coastal permafrost locations have been sampled along the Canadian Beaufort Sea Coast, hydrodynamically fractionated and analyzed for their C, N, 13C and 14C content. We found that the matrix-free fraction (low density <1.8 g/cm3, and high-density >1.8 g/cm3; <38 μm) account for 77%–98% of the OC. By using a coastal classification combined with field data, our results showed that short coastal segments can become key players in delivering matrix-free, easily degradable OC to the marine system.

Carbon and nitrogen are considered as candidate light elements present in planetary cores. However, there is limited understanding regarding the structure and physical properties of Fe-C-N alloys under extreme conditions. Here diamond anvil cell experiments were conducted, revealing the stability of hexagonal-structured Fe7(N0.75C0.25)3 up to 120 GPa and 2100 K, without undergoing any structural transformation or dissociation. Notably, the thermal expansion coefficient and Grüneisen parameter of the alloy exhibit a collapse at 55–70 GPa. First-principles calculations suggest that such anomaly is associated with the spin transition of iron within Fe7(N0.75C0.25)3. Our modeling indicates that the presence of ∼1.0 wt% carbon and nitrogen in liquid iron contributes to 9–12% of the density deficit of the Earth's outer core. The thermoelastic anomaly of the Fe-C-N alloy across the spin transition is likely to affect the density and seismic velocity profiles of (C,N)-rich planetary cores, thereby influencing the dynamics of such cores.

Tropical wetlands account for ∼20% of the global total methane (CH4) emissions, but uncertainties remain in emission estimation due to the inaccurate representation of wetland spatiotemporal variations. Based on the latest satellite observational inundation data, we constructed a model to map the long-term time series of wetland extents over the Sudd floodplain, which has recently been identified as an important source of wetland CH4 emissions. Our analysis reveals an annual, total wetland extent of 5.73 ± 2.05 × 104 km2 for 2003–2022, with a notable accelerated expansion rate of 1.19 × 104 km2 yr−1 during 2019–2022 driven by anomalous upstream precipitation patterns. We found that current wetland products generally report smaller wetland areas, resulting in a systematic underestimation of wetland CH4 emissions from the Sudd wetland. Our study highlights the pivotal role of comprehensively characterizing the seasonal and interannual dynamics of wetland extent to accurately estimate CH4 emissions from tropical floodplains.

Tectonic landforms and surficial lithologic age are essential data for producing quality late Quaternary fault maps and predicting coseismic fault rupture location before an earthquake. However, we lack a clear understanding of the relationship between tectonic landforms and shallow earthquake processes and how lithologic age relates to landform preservation. We assess how fault location error (rupture-to-fault separation distance) and coseismic displacement residual (difference between observed and predicted coseismic displacement) vary with tectonic landform and lithologic age for four historical earthquakes. Certain tectonic landforms identified before these earthquakes correlate with lower fault location errors and median displacements below model predictions. Faults cutting Holocene units exhibit the largest location errors, reflecting surface processes that erode or bury fault evidence. This study shows that tectonic landforms and lithologic age have a significant impact on fault location uncertainty and coseismic displacement, which should be considered in fault mapping and fault displacement assessment.

Nickel stable isotopes (δ60Ni) provide insight to Ni biogeochemistry in the modern and past oceans. Here, we present the first Pacific Ocean high-resolution dissolved Ni concentration and δ60Ni data, from the US GEOTRACES GP15 cruise. As in other ocean basins, increases in δ60Ni toward the surface ocean are observed across the entire transect, reflecting preferential biological uptake of light Ni isotopes, however the observed magnitude of fractionation is larger in the tropical Pacific than the North Pacific Subtropical Gyre. Such surface ocean fractionation by phytoplankton should accumulate isotopically lighter Ni in the deep Pacific, yet we find that North Pacific deep ocean δ60Ni is similar to previously reported values from the deep Atlantic. Finally, we find that seawater dissolved δ60Ni in regions with hydrothermal input can be either higher or lower than background deep ocean δ60Ni, depending on vent geochemistry and proximity.

The last deglaciation experienced the retreat of massive ice sheets and a transition from the cold Last Glacial Maximum to the warmer Holocene. Key simulation challenges for this period include the timing and extent of ice sheet decay and meltwater input into the oceans. Here, major uncertainties and forcing factors for the last deglaciation are evaluated. Two sets of transient simulations are performed based on the novel ice-sheet reconstruction PaleoMist and the more established GLAC1D. The simulations reveal that the proximity of the Atlantic meridional overturning circulation (AMOC) to a bifurcation point, where it can switch between on- and off-modes, is primarily determined by the interplay of greenhouse gas concentrations, orbital forcing and freshwater forcing. The PaleoMist simulation qualitatively replicates the Bølling-Allerød (BA)/Younger Dryas (YD) sequence: a warming in Greenland and Antarctica during the BA, followed by a cooling northern North Atlantic and an Antarctic warming during the YD.

The 2024 Mw 7.5 Noto Peninsula, Japan, earthquake was initiated within the source region of intense swarm activity. To reveal the mainshock early process, we relocated the earthquake hypocenters and found that many key phenomena, including the mainshock initiation, foreshocks, swarm earthquakes, and deep aseismic slip, occurred at parts of a previously unrecognized fault in intricate fault network. This fault is subparallel (several kilometers deeper) to a known active fault, and the mainshock initiation and foreshocks occurred at the front of a 2-year westward swarm migration. The initiation location coincides with the destination of the upward migration of a deeper earthquake cluster via several smaller faults. Fluid supply, small earthquakes, and aseismic slip on the fault likely triggered the mainshock, leading to the first major rupture at the western region, propagating further to the west and east sides, resulting in an Mw7.5 event, exceeding 100 km in length.

As a practical reflection of the opportunity window of Madden-Julian Oscillation (MJO), there are intermittent periods of relatively high forecasting skills, namely the forecast skill windows. Robust forecast skill windows are identified based on the subseasonal-seasonal reforecast database, during which the majority of models show high forecast skills. A total of 15 MJO forecast skill windows during 1993–2020 have been identified. Most of the forecast skill windows are closely associated with active MJO events with high amplitude. Whether a high-skill forecast window appears significantly depends on the magnitude of MJO intensity during the same period. The maintenance of active strong MJO events is potentially related with the warmer surface sea temperature anomalies in the western Pacific. Further research into such processes may unveil the MJO development mechanism and improve the MJO forecast skill.

This study analyzes hydrometeor evolution during rapid intensification (RI) and tropical cyclone (TC) intensity dependence using satellite data. Previous studies have suggested ice cloud water or non-convective precipitation as a predictor of RI from different perspectives. However, few studies have focused on the impact of TC intensity or comprehensive comparisons to identify better indicators. During RI, hydrometeor contents in weak TCs increase over the entire region, whereas they increase mainly in the inner-core region and decrease in advance in the outer-core region for strong TCs. Hydrometeor contents in the inner-core are higher in RI than in slow intensification, and their maxima location is related to TC intensity and intensification rate. Cloud water path (CWP) in the inner-core region is most correlated with the intensification rate, especially in weak TCs. Therefore, the CWP can serve as a predictor of RI and can be applied to all TC intensities.

In this study, prominent dust source areas are identified along with their plume extent using high temporal frequency satellite observations. Hourly dust plume observations of the Dust Belt from geostationary-orbit satellites are analyzed for the 2017-12–2022-11 period. To identify dust source areas and their extents, we back-track plumes to their source, assessing source areas in terms of emission frequency, contribution, and plume extent patterns. This method advances over traditional source allocation techniques that rely on polar-orbiting satellites based on a few daily passes and meteorological wind fields for backtracking. Our findings indicate that Boreal summer is the most intense season for most sources, except in the Southern Sahara, which experiences winterly winds. Our analysis also reveals significant contributions from regions within the Sahara that experience expansive but infrequent dust storms, highlighting the importance of considering both frequency and magnitude in understanding dust emissions.

Radio-echo sounding (RES) shows large-scale englacial stratigraphic folds are ubiquitous in Greenland's ice sheet. However, there is no consensus yet on how these folds form. Here, we use the full-Stokes code Underworld2 to simulate ice movements in three-dimensional convergent flow, mainly considering ice anisotropy due to a crystallographic preferred orientation, vertical viscosity and density gradients in ice layers, and bedrock topography. Our simulated folds show complex patterns and are classified into: large-scale folds (>100 m amplitude), small-scale folds (<<100 m) and basal-shear folds. The amplitudes of large-scale folds tend to be at their maximum in the middle of the ice column or just below, in accordance with observations in RES data. We conclude that ice anisotropy amplifies the perturbations in ice layers (mainly due to bedrock topography) into large-scale folds during flow. Density differences between the warm deep ice and cold ice above may enhance fold amplification.

Ion-acoustic waves (IAWs) commonly occur near interplanetary (IP) shocks. These waves are important because of their potential role in the dissipation required for collisionless shocks to exist. We study IAW occurrence statistically at different heliocentric distances using Solar Orbiter to identify the processes responsible for IAW generation near IP shocks. We show that close to IP shocks the occurrence rate of IAW increases and peaks at the ramp. In the upstream region, the IAW activity is highly variable among different shocks and increases with decreasing distance from the Sun. We show that the observed currents near IP shocks are insufficient to reach the threshold for the current-driven instability. We argue that two-stream proton distributions and suprathermal electrons are likely sources of the waves.

The Central Asian Orogenic Belt (CAOB) was formed by the aggregation and collage of numerous Paleozoic subduction-accretion assemblages and Precambrian microcontinental blocks. However, the tectonic nature of the southeastern CAOB remains controversial, which complicates the reconstruction of the Paleo-Asian Ocean. To address this issue, a deep seismic reflection survey was initiated across the southeastern CAOB and reveals broad gentle sub-horizontal reflectors in the middle-lower crust and a relatively transparent zone in the upper crust. Combining with the Precambrian geological outcrops and other geophysical features, we support a microcontinental block, the Xilinhot Block, existed in the Paleo-Asian domain. Thus, the Paleo-Asian Ocean was separated into two branches that underwent north-dipping and double-dipping oceanic plate subduction, respectively, to form the Hegenshan-Heihe and Solonker sutures. Multiple relics beneath Hegenshan-Heihe Suture indicate that multiple sets of unidirectional oceanic subduction-accretion and magmatism were important mechanisms of continental growth.

Severe storms produce ocean waves with periods of 18–26 s, corresponding to wavelengths 500–1,055 m. These waves radiate globally as swell, generating microseisms and affecting coastal areas. Despite their significance, long waves often elude detection by existing remote sensing systems when their height is below 0.2 m. The new Surface Water Ocean Topography (SWOT) satellite offers a breakthrough by resolving these waves in global sea level measurements. Here we show that SWOT can detect 25-s waves with heights as low as 3 cm, and resolves period and direction better than in situ buoys. SWOT provides detailed maps of wave height, wavelength, and direction across ocean basins. These measurements unveil intricate spatial patterns, shedding light on wave generation in storms, currents that influence propagation, and refraction, diffraction and reflection in shallow regions. Notably, the magnitude of reflections exceeds previous expectations, illustrating SWOT's transformative impact.

Megathrusts at convergent plate boundaries generate the largest and some of the most hazardous earthquakes on Earth. However, their physical properties, including those influencing fault slip accumulation and release and earthquake-related surface displacements, are still poorly constrained at critical depths. Here, we combine seismic imaging and geodetic modeling to investigate the structure and mechanical behavior of the Main Himalayan Thrust fault (MHT) in the center of the 2015 Mw 7.8 Gorkha rupture in Nepal. Our results from two independent observations consistently suggest the presence of a channel associated with the MHT with high compliance (shear modulus as low as ∼4 GPa) and strain anisotropy (stiffer in the vertical orientation than in the horizontal), likely arising from a weak subducting layer with north-dipping foliation. Such mechanical heterogeneity significantly influences the quantification of short-term fault kinematics and associated earthquake potential, with implications on across-scale dynamics of plate boundaries in Himalaya and elsewhere.

The recently-observed poleward shift in western North Pacific tropical cyclone (TC) genesis has increased the TC threat to East Asia. We find that the poleward shift of TC genesis since 1979 is linked to mega-ENSO. A downscaling analysis of TC genesis latitude given the constraint of mega-ENSO using 30 models from the Coupled Model Intercomparison Project Phase 6 (CMIP6) show a continued increasing poleward shift with additional warming. We use the dynamic genesis potential index as a TC proxy in future CMIP6 simulations. These simulations show enhanced TC formation in the subtropics and decreased TC formation in the tropics. Modeled TCs in CMIP6 high-resolution models that well represent mega-ENSO project future poleward shifts in TC genesis. Both observations and simulations show that extra-tropical North Pacific sea surface temperature warming associated with mega-ENSO are the primary driver of the TC genesis poleward shift. Our study provides new insights into climate change-driven TC migration.

Brown carbon (BrC) in China is of great interest to the regional and global climate due to its strong absorption of sunlight. However, the contribution of BrC to total carbonaceous aerosol light absorption and its direct radiative effects (DRE) in China remains largely uncertain. To better assess its climate impact in China, we develop an explicit BrC scheme and implement it in a global climate model, which includes optical parameters of primary BrC derived from local measurements, secondary BrC absorption, and a photobleaching parameterization of BrC. By comparing with multi-type observational data, we find that with the implementation of this scheme, the model captures the seasonal variations of BrC light absorption well in China. The model estimates that BrC contributes 19% and 12% to the total light absorption of carbonaceous aerosol in China in winter and summer, resulting in 0.110 and 0.205 W m−2 of DRE, respectively.

Mars' polar ice deposits are thought to preserve a record of climate throughout their evolution. In addition to the large north polar layered deposits (NPLD) at Mars' north pole, smaller ice deposits are preserved in craters nearby. These outlying deposits were potentially formed by the same mechanisms that drive NPLD formation, or may represent more local mechanisms. Distinguishing between these possibilities would help elucidate the spatial homogeneity of Martian climate processes. Here, we analyzed SHARAD radar depth profiles from 34 locations across the NPLD and 5 locations within the Korolev crater ice deposit using Fourier transform analysis and dynamic time warping to quantitatively assess the similarity between the internal layered stratigraphy of the two deposits. We identify broad stratigraphic similarities between the Korolev deposit and the NPLD, suggesting they likely formed due to the same climate forcing mechanism, with local variability also observed across the NPLD.

The evolution of charge structure involved in lightning discharge of a thunderstorm over the central Tibetan Plateau is investigated for the first time, based on the data from very high frequency interferometer, radar and sounding. During the developing-mature stage, the TP thunderstorm exhibited a tripolar charge structure evolved from an initial inverted dipole. At the mature stage, a bottom-heavy tripole charge structure is clearly presented, with a strong lower positive charge center (LPCC) at temperatures above −10°C, a middle negative charge region between −30°C and −15°C, and an upper positive charge region at T < −30°C. As the LPCC was depleted, the charge structure evolved into a normal tripole with a pocket LPCC. The merging between different convective cells resulted in the formation of two adjacent negative charge regions located directly and obliquely above the LPCC, and horizontally arranged different charge regions were simultaneously involved in the same lightning discharge.