Feed aggregator

Abstract

We investigated velocity and anisotropic structure of the crust beneath the eastern margin of the Tibetan Plateau to better understand its deformation and evolution mechanism. We performed H-κ and Pms anisotropy analyses to obtain crustal thickness, Vp/Vs ratio, fast polarization direction, and splitting time from 711 stations, and further conducted quality control using slowness, harmonic and statistical analyses. The Songpan-Ganzi Block has a large splitting time and a fast polarization direction roughly parallel to the GPS motion and SKS fast direction. It also shows an overall high but complex distribution of Vp/Vs ratio, and large variations in crustal thickness, indicating that crustal deformation is likely caused by crustal shortening and lower crustal flow. The northern Sichuan-Yunnan Rhombic Block (SYRB) is featured by a thick crust and high Vp/Vs ratio, suggesting that the crust is likely inflated by partial melting lower crustal rocks. The subblock also exhibits a strong azimuthal anisotropy with a splitting time greater than 0.6 s. The fast polarization direction aligns with the nearly N-S extended direction and rotates clockwise in front of the Emeishan Large Igneous Province (ELIP). The observed anisotropy agrees with aligned amphibole minerals under a simple shear condition, supporting a southward lower crust flow being diverted by the ELIP. Anisotropy measurements on the southern SYRB are less robust and widely scattered, suggesting a deformation mechanism different from the northern SYRB. In addition, the southeastern margin of the Sichuan Basin shows a systematic pattern of crustal anisotropy consistent with a pure shear deformation mechanism.

Abstract

Shallow creep events provide opportunities to understand the mechanical properties and behavior of faults. However, due to physical limitations observing creep events, the precise spatio-temporal evolution of slip during creep events is not well understood. In 2023, the Superstition Hills and Imperial faults in California each experienced centimeter-scale slip events that were captured in unprecedented detail by satellite radar, sub-daily Global Navigation Satellite Systems, and creepmeters. In both cases, the slip propagated along the fault over 2–3 weeks. The Superstition Hills event propagated bilaterally away from its initiation point at average velocities of ∼9 km/day, but propagation velocities were locally much higher. The ruptures were consistent with slip from tens of meters to ∼2 km depths. These slowly propagating events reveal that the shallow crust of the Imperial Valley does not obey purely velocity-strengthening or velocity-weakening rate-and-state friction, but instead requires the consideration of fault heterogeneity or fault-frictional behaviors such as dilatant strengthening.

Abstract

The Superstition Hills Fault (SHF) exhibits a rich spectrum of slip modes, including M 6+ earthquakes, afterslip, quasi-steady creep, and both triggered and spontaneous slow slip events (SSEs). Following 13 years of quiescence, creepmeters recorded 25 mm of slip during 16–19 May 2023. Additional sub-events brought the total slip to 41 mm. The event nucleated on the northern SHF in early-May and propagated bi-laterally at rates on the order of kilometers per day. Surface offsets reveal a bi-modal slip distribution, with slip on the northern section of the fault being less localized and lower amplitude compared to the southern section. Kinematic slip models confirm systematic variations in the slip distribution along-strike and with depth and suggest that slip is largely confined to the shallow sedimentary layer. Observations and models of the 2023 SSE bear a strong similarity to previous slip episodes in 1999, 2006, and 2010, suggesting a characteristic behavior.

Abstract

The significant risk associated with fault reactivation often necessitates slip tendency analyses for effective risk assessment. However, such analyses are challenging, particularly in large areas with limited or absent reliable stress measurements and where the cost of extensive geomechanical analyses or simulations is prohibitive. In this paper, we propose a novel approach using a physics-informed neural network that integrates stress orientation and satellite displacement observations in a top-down multi-scale framework to estimate two-dimensional slip tendency analyses even in regions lacking comprehensive stress data. Our study demonstrates that velocities derived from a continental scale analysis, combined with reliable stress orientation averages, can effectively guide models at smaller scales to generate qualitative slip tendency maps. By offering customizable data selection and stress resolution options, this method presents a robust solution to address data scarcity issues, as exemplified through a case study of the South Australian Eyre Peninsula.

Abstract

The absorption of atmospheric carbon dioxide (CO2) in the Southern Ocean represents a critical component of the global oceanic carbon budget. Previous assessments of air-sea carbon flux variations and long-term trends in polar regions during winter have faced limitations due to scarce field data and the lack of ocean color satellite imagery, causing uncertainties in estimating CO2 flux estimation. This study utilized the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation satellite to construct a continuous 16-year (2006–2021) time series of sea surface partial pressure of CO2 (pCO2) in the Southern Ocean. Our findings revealed that the polar region in South Ocean acts as a carbon sink in winter, with CO2 flux of ∼30 TgC in high-latitude areas (South of 50°S). This work highlights the efficacy of active remote sensing for monitoring sea surface pCO2 and contributes insights into the dynamic carbonate systems of the Southern Ocean.

Abstract

We explore the links between elevation variability of the Antarctic Ice Sheet (AIS) and large-scale climate modes. Using multiple linear regression, we quantify the time-cumulative effects of El Niño Southern Oscillation (ENSO) and the Southern Annular Mode (SAM) on gridded AIS elevations. Cumulative ENSO and SAM explain a median of 29% of the partial variance and up to 85% in some coastal areas. After spatial smoothing, these signals have high spatial correlation with those from GRACE gravimetry (r∼ = 0.65 each). Much of the signal is removed by a firn densification model but inter-model differences exist especially for ENSO. At the lower parts of the Thwaites and Pine Island glaciers, near their grounding line, we find the Amundsen Sea Low (ASL) explains ∼90% of the observed elevation variability. There, modeled firn effects explain only a small fraction of the variability, suggesting significant height changes could be a response to climatological ice-dynamics.

Abstract

We construct a high-resolution shear-wave velocity (VS) model for the uppermost 100 m using ambient noise tomography near the West Antarctic Ice Sheet Divide camp. This is achieved via joint inversion of Rayleigh wave phase velocity and H/V ratio, whose signal-to-noise ratios are boosted by three-station interferometry and phase-matched filtering, respectively. The VS shows a steep increase (0.04–0.9 km/s) in the top 5 m, with sharp interfaces at ∼8–12 m, followed by a gradual increase (1.2–1.8 km/s) between 10 and 45 m depth, and to 2 km/s at ∼65 m. The compressional-wave velocity and empirically-obtained density profile compares well with the results from Herglotz–Wiechert inversion of diving waves in active-source shot experiments and ice core analysis. Our approach offers a tool to characterize high-resolution properties of the firn and shallow ice column, which helps to infer the physical properties of deeper ice sheets, thereby contributes to improved understanding of Earth's cryosphere.

Abstract

Using high-resolution data from the Magnetospheric Multiscale mission, an electron-only reconnection current sheet is found between two successive dipolarization fronts (DFs). The electron-only reconnection occurs between the northward component of the magnetic field of the flux pileup region (FPR) of the first DF (DF1) and the southward component of the magnetic dip of the second DF (DF2). The faster DF2 compresses the FPR of DF1, which constitutes an anti-parallel topology and reduces the thickness of the current sheet. Further analysis shows that the current sheet is unstable to the electron tearing instability, which may power the onset of the reconnection. Our results suggest that these two DFs may merge into one by the reconnection, which sheds light on the evolution of DFs during their earthward propagation.

Abstract

The magnetic fields of terrestrial planets are created by core convection. Molten silicate mantles could also generate magnetic fields through their convective motion, known as a silicate dynamo. Recent computational studies have suggested that silicate melts may exhibit high electrical conductivity (EC) at temperatures above 4000 K due to strong electronic conduction, which could activate a silicate dynamo. We determined the EC of dense molten MgSiO3 up to 71 GPa and 4490 K by static compression experiments. It jumped by one order of magnitude upon melting, but 57(27) S/m at 4490 K is much lower than previous predictions, suggesting that molten MgSiO3 carries charge via ions rather than predicted electronic conduction. Nevertheless, the strong temperature dependence of the ionic conductivity found in this study suggests that super-Earths’ hotter magma ocean with larger-scale convection could power a dynamo that drives magnetic fields, which plays key roles in sustaining planetary surface environments.

Abstract

An unseasonal equatorial plasma bubble (EPB) event occurred in the East/Southeast Asian sector during the geomagnetic storm on 1 December 2023, causing strong amplitude scintillations from equatorial to middle latitudes. Based on the observations from multiple instruments over a large latitudinal and longitudinal region, the spatial features of the super EPB were investigated. The EPB developed vertically at a fast rising speed ∼470 m/s over the magnetic equator and extended to a very high middle latitude more than 40°N, despite that the storm intensity was not very strong with the minimum SYM-H index −132 nT. In the zonal direction, the super EPB covered over a specific region ∼95–140°E, where the local sunset roughly coincided with southward turning of interplanetary magnetic field (IMF) Bz component. Before the onset of the super EPB, significant upward plasma drift up to ∼110 m/s was observed over the magnetic equator, which could amplify the growth rate of Rayleigh-Taylor instability and lead to the generation of the super EPB. The significant drift was likely caused by eastward penetration electric field (PEF) due to sharp southward turning of IMF Bz. The local time of storm onset and duration of IMF Bz southward turning during the storm main phase may partly determine the onset region and zonal coverage of the EPB.

Abstract

How the solar wind influences the magnetospheres of the outer planets is a fundamentally important question, but is difficult to answer in the absence of consistent, simultaneous monitoring of the upstream solar wind and the large-scale dynamics internal to the magnetosphere. To compensate for the relative lack of in-situ solar wind data, propagation models are often used to estimate the ambient solar wind conditions at the outer planets for comparison to remote observations or in-situ measurements. This introduces another complication: the propagation of near-Earth solar wind measurements introduces difficult-to-assess uncertainties. Here, we present the Multi-Model Ensemble System for the outer Heliosphere (MMESH) to begin to address these issues, along with the resultant multi-model ensemble (MME) of the solar wind conditions near Jupiter. MMESH accepts as input any number of solar wind models together with contemporaneous in-situ spacecraft data. From these, the system characterizes typical uncertainties in model timing, quantifies how these uncertainties vary under different conditions, attempts to correct for systematic biases in the input model timing, and composes a MME with uncertainties from the results. For the Juno-era (04/07/2016–04/07/2023) MME hindcast for Jupiter presented here, three solar wind propagation models were compared to in-situ measurements from the near-Jupiter spacecraft Ulysses and Juno spanning diverse geometries and phases of the solar cycle across >14,000 hr of data covering 2.5 decades. The MME gives the most-probable near-Jupiter solar wind conditions for times within the tested epoch, outperforming the input models and returning quantified estimates of uncertainty.

Science, Volume 384, Issue 6703, Page 1389-1389, June 2024.

A near-global multiyear climate data record of the fine-mode and coarse-mode components of atmospheric pure dust
Emmanouil Proestakis, Antonis Gkikas, Thanasis Georgiou, Anna Kampouri, Eleni Drakaki, Claire L. Ryder, Franco Marenco, Eleni Marinou, and Vassilis Amiridis
Atmos. Meas. Tech., 17, 3625–3667, https://doi.org/10.5194/amt-17-3625-2024, 2024
A new four-dimensional, multiyear, and near-global climate data record of the fine-mode (submicrometer diameter) and coarse-mode (supermicrometer diameter) components of atmospheric pure dust is presented. The dataset is considered unique with respect to a wide range of potential applications, including climatological, time series, and trend analysis over extensive geographical domains and temporal periods, validation of atmospheric dust models and datasets, and air quality.

An empirical characterization of the aerosol Ångström exponent interpolation bias using SAGE III/ISS data
Robert P. Damadeo, Viktoria F. Sofieva, Alexei Rozanov, and Larry W. Thomason
Atmos. Meas. Tech., 17, 3669–3678, https://doi.org/10.5194/amt-17-3669-2024, 2024
Comparing different aerosol data sets for scientific studies often requires converting aerosol extinction data between different wavelengths. A common approximation for the spectral behavior of aerosol is the Ångström formula; however, this introduces biases. Using measurements across many different wavelengths from a single instrument, we derive an empirical relationship to both characterize this bias and offer a correction for other studies that may employ this analysis approach.

Calibrating and validating the Integrated Valuation of Ecosystem Services and Tradeoffs (InVEST) urban cooling model: case studies in France and the United States
Perrine Hamel, Martí Bosch, Léa Tardieu, Aude Lemonsu, Cécile de Munck, Chris Nootenboom, Vincent Viguié, Eric Lonsdorf, James A. Douglass, and Richard P. Sharp
Geosci. Model Dev., 17, 4755–4771, https://doi.org/10.5194/gmd-17-4755-2024, 2024
The InVEST Urban Cooling model estimates the cooling effect of vegetation in cities. We further developed an algorithm to facilitate model calibration and evaluation. Applying the algorithm to case studies in France and in the United States, we found that nighttime air temperature estimates compare well with reference datasets. Estimated change in temperature from a land cover scenario compares well with an alternative model estimate, supporting the use of the model for urban planning decisions.

Assessment of a tiling energy budget approach in a land surface model, ORCHIDEE-MICT (r8205)
Yi Xi, Chunjing Qiu, Yuan Zhang, Dan Zhu, Shushi Peng, Gustaf Hugelius, Jinfeng Chang, Elodie Salmon, and Philippe Ciais
Geosci. Model Dev., 17, 4727–4754, https://doi.org/10.5194/gmd-17-4727-2024, 2024
The ORCHIDEE-MICT model can simulate the carbon cycle and hydrology at a sub-grid scale but energy budgets only at a grid scale. This paper assessed the implementation of a multi-tiling energy budget approach in ORCHIDEE-MICT and found warmer surface and soil temperatures, higher soil moisture, and more soil organic carbon across the Northern Hemisphere compared with the original version.

The ddeq Python library for point source quantification from remote sensing images (version 1.0)
Gerrit Kuhlmann, Erik Koene, Sandro Meier, Diego Santaren, Grégoire Broquet, Frédéric Chevallier, Janne Hakkarainen, Janne Nurmela, Laia Amorós, Johanna Tamminen, and Dominik Brunner
Geosci. Model Dev., 17, 4773–4789, https://doi.org/10.5194/gmd-17-4773-2024, 2024
We present a Python software library for data-driven emission quantification (ddeq). It can be used to determine the emissions of hot spots (cities, power plants and industry) from remote sensing images using different methods. ddeq can be extended for new datasets and methods, providing a powerful community tool for users and developers. The application of the methods is shown using Jupyter notebooks included in the library.

Abstract

Stratovolcanoes are common globally, with high-altitude summit regions that are often glacier-clad and intersect the seasonal and perennial snow line. During an eruption, interaction between snow/ice and hot, pyroclastic deposits will potentially lead to extensive melt and steam production. This is particularly pertinent when pyroclastic density currents (PDCs) are emplaced onto and propagate over glacierised substrates. Generated melt and steam are incorporated into the flow, which can cause a transformation from a hot, dry granular flow, to a water-saturated, sediment-laden flow, termed a lahar. Both PDCs and ice-melt lahars are highly hazardous due to their high energy during flow and long runout distances. Knowledge of the physics that underpin these interactions and the transformation to ice-melt lahar is extremely limited, preventing accurate descriptions within hazard models. To physically constrain the thermal interactions we conduct static melting experiments, where a hot granular layer was emplaced onto an ice substrate. The rate of heat transfer through the particle layer, melt and steam generation were quantified. Experiments revealed systematic increases in melt and steam with increasing particle layer thicknesses and temperatures. We also present a one-dimensional numerical model for heat transfer, calibrated against experimental data, capable of accurately predicting temperature and associated melting. Furthermore, similarity solutions are presented for early-time melting which are used to benchmark our numerical scheme, and to provide rapid estimates for meltwater flux hydrographs. These data are vital for predicting melt volume and incorporation into PDCs required to facilitate the transformation to and evolution of ice-melt lahars.

Abstract

The Mediterranean Sea releases approximately 1 Sv of water into the North Atlantic through the Gibraltar Straits, forming the saline Mediterranean Outflow Water (MOW). Its impact on large-scale flow and specifically its northbound Lagrangian pathways are widely debated, yet a comprehensive overview of MOW pathways over recent decades is lacking. We calculate and analyze synthetic Lagrangian trajectories in 1980–2020 reanalysis velocity data. Sixteen percent of the MOW follow a direct northbound path to the sub-polar gyre, reaching a 1,000 m depth crossing window at the southern tip of Rockall Ridge in about 10 years. Surprisingly, time-dependent chaotic advection, not steady currents, drives over half of the northbound transport. Our results suggest a potential 15–20 years predictability in the direct northbound transport. Additionally, monthly variability appears more significant than inter-annual variability in Lagrangian mixing and spreading the MOW.

Abstract

Large-scale seepage faces occur on tidal flats with gentle slope, which are widely distributed worldwide. Evaporation on these seepage faces, leading to salt retention and accumulation, may significantly impact the density-dependent groundwater flow beneath the tidal flats. However, due to nonlinear complexities of the groundwater flow and solute transport on seepage faces, explicit boundary conditions and numerical models to quantify these processes are lacking. In this study, we present both mathematical and numerical models to quantify these processes. Compared to the results of our previous study, this paper shows that seepage-face evaporation can (a) significantly increase the groundwater salinity in the upper intertidal zone, and form multiple groundwater circulation cells in the intertidal zone, (b) cause the disappearance of multiple seepage-faces and reduce the spatial extent of seepage faces notably, (c) and intensify the groundwater and salt exchange as well as the seawater-groundwater circulation through the intertidal zone.