Feed aggregator

Abstract

In March and April 2021, buoys were deployed in the Beaufort Sea, Arctic Ocean, to measure sea-ice horizontal deformation over spatial scales that had not been previously achieved. Geodetic-quality position measurements allowed measurements of strain-rate over lengths from about 200 m to 2 km. Conventional ice-drifters extended spatial coverage up to about 100 km. Past studies find there is multi-fractal behavior for horizontal sea-ice deformation from 10 to 1,000 km. Our results demonstrate that such behavior does not hold when including spatial scales below 10 km. We find that sea-ice deformation is not scale invariant between the scale of individual sea-ice floes and aggregates of floes. Therefore, we cannot expect the same physical laws or forcing to describe sea-ice kinematics over these regimes, nor can we assume log-log linear behavior for mean deformation. Using this scaling behavior as a metric to validate models that resolve sea ice floes and their interactions is hence not recommended.

Atmospheric odd nitrogen response to electron forcing from a 6D magnetospheric hybrid-kinetic simulation
Tuomas Häkkilä, Maxime Grandin, Markus Battarbee, Monika E. Szeląg, Markku Alho, Leo Kotipalo, Niilo Kalakoski, Pekka T. Verronen, and Minna Palmroth
Ann. Geophys. Discuss., https//doi.org/10.5194/angeo-2024-7,2024
Preprint under review for ANGEO (discussion: final response, 3 comments)
We study the atmospheric impact of auroral electron precipitation, by the novel combination of both magnetospheric and atmospheric modelling. We first simulate fluxes of auroral electrons, and then use these fluxes to model their atmospheric impact. We find an increase of up to 200 % in thermospheric odd nitrogen, and a corresponding decrease in stratospheric ozone of around 0.7 %. The produced auroral electron precipitation is realistic, and shows the potential for future studies.

Author(s): R. Liotard, B. Canaud, A. Pineau, A. Sollier, E. Lescoute, A. Colaïtis, and G. Duchateau

Laser direct drive (LDD) inertial confinement fusion (ICF) involves irradiating a spherical target of thermonuclear fuel coated with an ablator, usually made of polystyrene. Laser energy absorption near the target surface leads to matter ablation, hydrodynamic shocks, and ultimately capsule implosio…

[Phys. Rev. E 109, 065207] Published Wed Jun 12, 2024

Author(s): Ying Chen, Xi Liu, Zhenhua Chai, and Baochang Shi

In this paper we first present the general propagation multiple-relaxation-time lattice Boltzmann (GPMRT-LB) model and obtain the corresponding macroscopic finite-difference (GPMFD) scheme on conservative moments. Then based on the Maxwell iteration method, we conduct the analysis on the truncation …

[Phys. Rev. E 109, 065305] Published Wed Jun 12, 2024

Abstract

In the present study, we use broadband seismic data recorded by 190 stations of the EarthScope program's Transportable Array to construct a 3-D shear wave velocity model for the upper 180 km using a non-linear Bayesian Monte-Carlo joint inversion of receiver functions (RFs) and Rayleigh wave dispersion curves. Ambient noise and teleseismic data are used for obtaining Rayleigh wave phase velocity dispersion curves. A resonance removal filtering technique is applied to the RFs contaminated by reverberations from the thick sedimentary layers that cover most of the region. Our observations of the higher crustal shear velocities (∼3.40 km/s) beneath the Sabine Block (SB), along with the estimated relatively thicker crust (∼34.0 km) and lower crustal V p/V s estimates (∼1.80) in comparison with the rest of the Gulf Coastal Plain (GCP) (∼3.10 km/s for crustal shear velocities, ∼29.0 km for crustal thickness, and ∼1.90 for crustal V p/V s estimates), indicating that this crustal block has different crustal properties from the surrounding coastal plain regions. The southern Ouachita Mountains have a thin crust (∼30.0 km) and low mean crustal V p/V s value (∼1.73), suggesting that lower crustal delamination has occurred in this region. Low velocities in the upper mantle beneath most of the GCP are interpreted as a combined result of thin lithosphere, higher-than-normal temperatures, and possibly compositional variations.

Abstract

Dichlorodifluoromethane (CFC-12) is an ozone-depleting substance and potent greenhouse gas, which was required to be phased out after 2010 under the Montreal Protocol. CFC-12 emissions need to be quantitatively traced. However, estimates of CFC-12 emissions in China based on atmospheric inversions are unavailable after 2010. Here, using atmospheric observations at nine sites across China and inversion techniques, we quantify CFC-12 emissions in China during 2011–2020 (on average 11.0 ± 0.6 Gg yr−1). The emissions derived from observations are 8.5 times larger than the previously reported inventories. Apart from emissions from eastern China revealed in previous studies, this study reveals that 71% of national total emissions were from other parts of China. Moreover, this study reconciled the global CFC-12 emissions during 2011–2020: 28% were traced to China by this study, 9% of emissions were traced in previous studies, while 63% remain untraced, indicating the need for more regional emission inversion studies.

Abstract

The subaerial exposure of the modern continental crust through time remains intensely debated, with estimates of the first exposure ranging from the late Archean to the Neoproterozoic. To constrain when and how much of the continental crust was exposed subaerially during Earth's history, we trained a supervised machine learning model on the compositions of modern subaerial and submerged basalts. Then, we applied this well-trained model to a refined worldwide data set of basaltic compositions and calculated the mean proportions of basalts erupted subaerially since 3.8 billion years ago (Ga). Our results suggest that ∼20% of the basalts were exposed subaerially in the early Archean, which may have driven the synthesis of biopolymers crucial to the origins and evolution of life. The proportion of subaerial basalts increased markedly during two stages between the late Archean and the Paleoproterozoic before reaching the present-day level no earlier than ∼1.8 Ga.

Abstract

Ground-level enhancements (GLEs) are sporadic events that signal the arrival of high fluxes of solar energetic particles (SEPs) that have been produced by solar eruptions. Ground-level enhancement events are characterized by a significant increase in the count rate of ground-based neutron monitors (NMs). The arrival of high-energy SEPs in the atmosphere leads to an enhancement of the radiation environment, with the enhancement at aviation altitudes being particularly hazardous to human health as pilots, crew, and airline passengers can be subjected to dangerous levels of radiation during a GLE. Through the use of a currently expanding library of analyzed GLEs and the application of a newly developed atmospheric radiation model, both of which have been created in-house, we found a strong statistically significant relationship between real-time NM data during GLE events and the radiation doses at aviation altitudes. This result provides a strong scientific basis for the use of real-time NM data as a proxy for radiation dose estimates during GLE events and aids in the development of future nowcasting models to help mitigate the dangerous impacts of future GLEs.

Abstract

The adiabatic invariants (M, J, Φ) and the related invariants (M, K, L*) have been established as effective coordinate systems for describing radiation belt dynamics at a theoretical level, and through numerical techniques, can be paired with in situ observations to order phase-space density. To date, methods for numerical techniques to calculate adiabatic invariants have focused on empirical models such the Tsyganenko models TS05, T96, and T89. In this work, we develop methods based on numerical integration and variable step size iteration for the calculation of adiabatic invariants, applying the method to the Lyon-Fedder-Mobarry (LFM) global magnetohydrodynamics (MHD) simulation code, with optional coupling to the Rice Convection Model (RCM). By opening the door to adiabatic invariant modeling with MHD magnetic fields, the opportunity for exploratory modeling work of radiation belt dynamics is enabled. Calculations performed using LFM are cross-referenced with the same code applied to the T96 and TS05 Tsyganenko models evaluated on the LFM grid. Important aspects of the adiabatic invariant calculation are reviewed and discussed, including (a) sensitivity to magnetic field model used, (b) differences in the problem between quiet and disturbed geomagnetic states, and (c) the selection of key parameters, such as the magnetic local time step size for drift shell determination. The rigorous development and documentation of this algorithm additionally acts as preliminary step for future thorough reassessment of in situ phase-space density results using alternative magnetic field models.

Abstract

Mountain waves are known sources of fluctuations in the upper atmosphere. However, their effects over the Continental United States (CONUS) are considered modest as compared to hot spots such as the Southern Andes. Here, we present an observation-guided case study examining the dynamics of gravity waves (GWs) and their impacts on the ionosphere over the CONUS prior to the cold air outbreak in December 2022, which resulted from a significant distortion of the tropospheric polar vortex. The investigation relies on MERRA-2 and ERA5 reanalysis data sets for the climatological contextualization, analysis of GWs based on National Aeronautics and Space Administration Aqua satellite's Atmospheric Infrared Sounder, 557.7 and 630.0 nm airglow emission observations, and the measurements of ionospheric disturbances retrieved from Global Navigation Satellite System signal-based total electron content (TEC) and Super Dual Auroral Radar Network observations. We demonstrate that the tropospheric polar jet stream shifted toward the Rocky Mountains, generated large amplitude GWs (up to 11 K of brightness temperature), which, aided by winter-time winds over mid-latitudes, could propagate to mesospheric heights. The breaking of GWs plausibly led to the generation of a plethora of secondary acoustic and GWs that eventually emerged as the sources of extensive ionospheric fluctuations of ∼3–30 min periods and up to 0.7 TECu, observed across the entire CONUS for several days. This case offers a valuable demonstration of the interplay between tropospheric circulation and the ionosphere over CONUS, pointing to the need for a better understanding of wave-driven deep-atmosphere coupled dynamics.

Abstract

The regional ionospheric corrections have become one of the critical parts of Precise Point Positioning (PPP)-based Real-Time-Kinematic (RTK) services to achieve fast positioning convergence. Several methods for regional ionospheric corrections supporting PPP-RTK have been developed and implemented over recent years, but little attention is given to the theoretical foundation of existing ionospheric correction methods and their performance comparison to find an optimal method in some sense. The optimality criterion of such methods should not only be based on the precision of the ionospheric correction itself, but also on its broadcasting strategies, and implementation aspects. This contribution studies ionospheric correction generation methods within the best linear unbiased predictor (BLUP) framework. Comparing the accuracy performances of the methods, we demonstrate that the Kriging method with trend, as a special case of BLUP, is the most appropriate method for large-scale networks (above 500 km). A strategy for the evaluation of the uncertainty of the grid-interpolated ionospheric corrections is also developed. In contrast to other empirical methods, this new method is rigorous in the sense that it avoids the underestimation of the uncertainty of predicted ionospheric corrections, especially when reference stations are close to a grid point.

Nature Geoscience, Published online: 12 June 2024; doi:10.1038/s41561-024-01475-5

Artificial intelligence tools have the potential to revolutionize how scientists work and publish. We share our ground rules for managing the inherent risks.

Nature Geoscience, Published online: 12 June 2024; doi:10.1038/s41561-024-01463-9

A field experiment in Uganda shows how potassium and phosphorus keep leaves functioning during times of water scarcity, highlighting the need to consider ecosystem-scale processes in studying the response of forests to nutrient limitation.

SummaryDynamic topography is defined as the deflection of Earth's surface due to the convecting mantle. ASPECT (Advanced Solver for Planetary Evolution, Convection, and Tectonics) is a continually evolving, finite element code that uses modern numerical methods to investigate problems in mantle convection. With ASPECT version 2.0.0 a consistent boundary flux (CBF) algorithm, used to calculate radial stresses at the model boundaries, was implemented into the release version of ASPECT. It has been shown that the CBF algorithm improves the accuracy of dynamic topography calculations by approximately one order of magnitude. We aim to evaluate the influence of the CBF algorithm and explore the geophysical implications of these improved estimates of dynamic topography changes along the East Coast of the United States. We constrain our initial temperature conditions using the tomography models SAVANI, S40RTS, and TX2008 and combine them with a corresponding radial viscosity profile (2 for TX2008), and 2 different boundary conditions for a total of 8 experiments. We perform simulations with and without the CBF method, which takes place during post-processing and does not affect the velocity solution. Our dynamic topography calculations are spatially consistent in both approaches, but generally indicate an increase in magnitude using the CBF method (an average ∼15 per cent and ∼76 per cent absolute change in present-day instantaneous and rate of change of dynamic topography, respectively). This enhanced accuracy in dynamic topography calculations can be used to better evaluate the effects of mantle convection on surface processes including vertical land motions, sea-level changes, and sedimentation and erosion. We explore results along the US East Coast, where a Pliocene shoreline has been deformed by dynamic topography change. An increased accuracy in estimates of dynamic topography can improve Pleistocene and Pliocene sea-level reconstructions, which allow for a better understanding of past sea-level changes and ice sheet stability.

SummaryThermal convection in planetary solid (rocky or icy) mantles sometimes occurs adjacent to liquid layers with a phase equilibrium at the boundary. The possibility of a solid-liquid phase change at the boundary has been shown to greatly help convection in the solid layer in spheres and plane layers and a similar study is performed here for a spherical shell with a radius-independent central gravity subject to a destabilising temperature difference. The solid-liquid phase change is considered as a mechanical boundary condition and applies at either or both horizontal boundaries. The boundary condition is controlled by a phase change number, Φ, that compares the time-scale for latent heat exchange in the liquid side to that necessary to build a topography at the boundary. We introduce a numerical tool, available at https://github.com/amorison/stablinrb, to carry out the linear stability analysis of the studied setup as well as other similar situations (cartesian geometry, arbitrary temperature and viscosity depth-dependent profiles). Decreasing Φ makes the phase change more efficient, which reduces the importance of viscous resistance associated to the boundary and makes the critical Rayleigh number for the onset of convection smaller and the wavelength of the critical mode larger, for all values of the radii ratio, γ. In particular, for a phase change boundary condition at the top or at both boundaries, the mode with a spherical harmonics degree of 1 is always favoured for Φ ≲ 10−1. Such a mode is also favoured for a phase change at the bottom boundary for small (γ ≲ 0.45) or large (γ ≳ 0.75) radii ratio. Such dynamics could help explaining the hemispherical dichotomy observed in the structure of many planetary objects.

SummaryWe conduct a thorough analysis of seismic and acoustic data purported to be from the so-called ‘Interstellar Meteor’ which entered the Earth’s atmosphere off the coast of Papua New Guinea on 2014-01-08. Previous work had suggested that this meteor may have been caused by an alien spacecraft burning up in the atmosphere. We conclude that both previously-reported seismic signals are spurious - one has characteristics suggesting a local vehicular-traffic based origin; whilst the other is statistically indistinguishable from the background noise. As such, previously-reported localisations based on this data are unreliable. Analysis of acoustic data provides a best fit location estimate which is very far (∼170 km) from the reported fireball location. Accordingly, we conclude that material recovered from the seafloor and purported to be from this event is almost certainly unrelated to it, and is likely of more mundane (non-interstellar) origin.

SummaryIn this study, we propose a systematic and effective method, that is, an extended version of the generalized reflection/transmission (R/T) coefficient method, for computing the phase-velocity (${c}_r$) dispersion curves, attenuation coefficient ($\alpha $) curves, and eigenfunctions of both Rayleigh and Love waves as well as the ellipticity of Rayleigh waves in layered viscoelastic-vertical transversely isotropic (VTI) media. The numerical scheme of combining the root-searching method with the local optimization method is designed for determining the complex-valued modal solutions (i.e., complex wavenumber $k = {\omega / {{c_{r}} - i\alpha }}$) of surface waves. The near-surface sedimentary geological environment is taken as the model example because it is typical viscoelastic-VTI media. Besides the anisotropic-viscoelastic (AV) media, our algorithm can also compute surface waves in isotropic-elastic (IE), isotropic-viscoelastic (IV), and anisotropic-elastic (AE) media by resetting the corresponding parameters. Using the six-layer half-space models and in these four media, we verify the correctness of our algorithm by benchmarking the modal solutions against those from other methods. In the four-layer half-space model, by comparing the results of IE, IV, AE, and AV media, we analyze the effects of velocity anisotropy, viscoelasticity and attenuation anisotropy on the dispersion and attenuation characteristics of both Rayleigh and Love waves in detail. Our study can provide a theoretical basis and useful tool for surface wave imaging considering the anisotropy and/or viscoelasticity of the medium, which has the potential to better investigate the solid Earth's internal structure.

The Periadriatic Fault system (PAF) ranks among the largest post-collisional structures of the European Alps. Recent geodetic data suggest that a fraction of the Adria–Europe convergence is still being accommo...

Emissions of nitrous oxide—a greenhouse gas more potent than carbon dioxide or methane—continued unabated between 1980 and 2020, a year when more than 10-million metric tons were released into the atmosphere primarily through farming practices, according to a new report by the Global Carbon Project.

Climate models are a key technology in predicting the impacts of climate change. By running simulations of the Earth's climate, scientists and policymakers can estimate conditions like sea level rise, flooding, and rising temperatures, and make decisions about how to appropriately respond. But current climate models struggle to provide this information quickly or affordably enough to be useful on smaller scales, such as the size of a city.