Geophysics (physics.geo-ph)

We present a fully implicit formulation of coupled fluid flow and geomechanics for fluid injection/withdrawal in fractured reservoirs in the context of CO2storage. Utilizing a Galerkin finite-element approach, both flow and poroelasticity equations are discretized on a shared three-dimensional mesh. The fluid flow is assumed to be single-phase. The hydraulic behaviour of fractures is represented through a double-nodes flow element, which allows to efficiently model longitudinal and transversal fracture permeabilities. In addressing the mechanical subproblem, fractures are explicitly modelled using cohesive elements to account for contact, friction and opening phenomena. The nonlinear set of equations is solved implicitly through an iterative partitioned conjugate gradient procedure, extending its traditional application to continuous problems to those involving explicit discontinuities such as faults and fractures. The model's accuracy is verified against analytical solutions for different geomechanical problems, notably for the growth of a frictional slip rupture along a fault due to fluid injection. Such a particularly challenging benchmark for a critically stressed fault is here reproduced for the first time by a finite-element based scheme. The capabilities of the developed parallel solver are then illustrated through a scenario involving CO2 injection into a faulted aquifer. The original solver code, tutorials, and data visualization routines are publicly accessible.

Thermal tides are atmospheric tides caused by variations in day-night insolation, similar to gravitational tides but with key differences. While both result in delayed mass redistribution, energy dissipation, and angular momentum exchanges between the planet and its host star, thermal tides can drive a planet's dynamics away from the rotational equilibrium states predicted by classical tidal theory. In this work, we present a novel closed-form solution for the thermotidal response of rocky planets. This general solution is derived from first principles, assuming either dry or moist adiabatic temperature profiles for the planet's atmosphere, and can be readily applied to study the long-term evolution of exoplanets in the habitable zones of their host stars. Despite relying on a small number of parameters, the model successfully captures the key features of the thermotidal torque predicted by General Circulation Models (GCMs). It also accurately predicts Earth's current semidiurnal thermotidal response and provides new insights into the evolution of the length of day during the Precambrian era.

A new analytical model for the generation of axisymmetric tornado-type vortices has been developed. A solution to the nonlinear equation for the stream function in an unstable stratified atmosphere is obtained and analyzed within the framework of ideal hydrodynamics. The solution is sought by smooth connecting continuous solutions for the internal region ("eye"), the central region ("wall" with maximum velocities) and the external region of the tornado. Expressions describing radial dependences for the radial and vertical velocity components include combinations of Bessel functions. The vortex is spatially localized by radius and height. Convective instability of a stratified atmosphere leads to an increase in the radial and vertical components of velocities according to the hyperbolic sine law. A downward flow is observed near the tornado axis. The maximum speed of the upward flow is achieved at a certain radial distance at a certain height. Below this height, radial flows converge towards the central part of the tornado, and above this height there is an outflow from the "wall" to the axis and to the periphery. The radial structure of the azimuthal velocity is determined by the structure of the initial disturbance and can change with height. Maximum rotation is achieved in the tornado "wall" at a certain height. The increase in azimuthal velocity can occur according to a superexponential law. Possible structures of movements, scenarios for the development of a tornado and its dynamics are discussed.

Seismic phase picking, which aims to determine the arrival time of P- and S-waves according to seismic waveforms, is fundamental to earthquake monitoring. Generally, manual phase picking is trustworthy, but with the increasing number of worldwide stations and seismic monitors, it becomes more challenging for human to complete the task comprehensively. In this work, we explore multiple ways to do automatic phase picking, including traditional and learning-based methods.

We introduce the Seismic Language Model (SeisLM), a foundational model designed to analyze seismic waveforms -- signals generated by Earth's vibrations such as the ones originating from earthquakes. SeisLM is pretrained on a large collection of open-source seismic datasets using a self-supervised contrastive loss, akin to BERT in language modeling. This approach allows the model to learn general seismic waveform patterns from unlabeled data without being tied to specific downstream tasks. When fine-tuned, SeisLM excels in seismological tasks like event detection, phase-picking, onset time regression, and foreshock-aftershock classification. The code has been made publicly available on https://github.com/liutianlin0121/seisLM.

Data assimilation will be essential for the management and expansion of geological carbon storage operations. In traditional data assimilation approaches a fixed set of geological hyperparameters, such as mean and standard deviation of log-permeability, is often assumed. Such hyperparameters, however, may be highly uncertain in practical CO2 storage applications. In this study, we develop a hierarchical data assimilation framework for carbon storage that treats hyperparameters as uncertain variables characterized by hyperprior distributions. To deal with the computationally intractable likelihood function in hyperparameter estimation, we apply a likelihood-free (or simulation-based) inference algorithm, specifically sequential Monte Carlo-based approximate Bayesian computation (SMC-ABC), to draw independent posterior samples of hyperparameters given dynamic monitoring-well data. In the second step we use an ensemble smoother with multiple data assimilation (ESMDA) procedure to provide posterior realizations of grid-block permeability. To reduce computational costs, a 3D recurrent R-U-Net deep-learning surrogate model is applied for forward function evaluations. The accuracy of the surrogate model is established through comparisons to high-fidelity simulation results. A rejection sampling (RS) procedure for data assimilation is applied to provide reference posterior results. Detailed data assimilation results from SMC-ABC-ESMDA are compared to those from the reference RS method. These include marginal posterior distributions of hyperparameters, pairwise posterior samples, and history matching results for pressure and saturation at the monitoring location. Close agreement is achieved with 'converged' RS results, for two synthetic true models, in all quantities considered. Importantly, the SMC-ABC-ESMDA procedure provides speedup of 1-2 orders of magnitude relative to RS for the two cases.

The global demand for unconventional energy sources such as geothermal energy and white hydrogen requires new exploration techniques for precise subsurface structure characterization and potential reservoir identification. Magnetotelluric (MT) inversion is crucial for these tasks, providing critical information on the distribution of subsurface electrical resistivity at depths ranging from hundreds to thousands of meters. However, traditional iterative algorithm-based inversion methods require the adjustment of multiple parameters, demanding time-consuming and exhaustive tuning processes to achieve proper cost function minimization. Although recent advances have incorporated deep learning algorithms for MT inversion, these have been primarily based on supervised learning, which needs large labeled datasets for training. Therefore, it causes issues in generalization and model characteristics that are restricted to the neural network's features. This work utilizes TensorFlow operations to create a differentiable forward MT operator, leveraging its automatic differentiation capability. Moreover, instead of solving for the subsurface model directly, as classical algorithms perform, this paper presents a new deep unsupervised inversion algorithm guided by physics to estimate 1D MT models. Instead of using datasets with the observed data and their respective model as labels during training, our method employs a differentiable modeling operator that physically guides the cost function minimization, making the proposed method solely dependent on observed data. Therefore, the optimization problem is updating the network weights to minimize the data misfit. We test the proposed method with field and synthetic data at different acquisition frequencies, demonstrating that the resistivity models are more accurate than other results using state-of-the-art techniques.

Basement relief gravimetry is crucial in geophysics, especially for oil exploration and mineral prospecting. It involves solving an inverse problem to infer geological model parameters from observed data. The model represents basement relief with constant-density prisms, and the data reflect gravitational anomalies from these prisms. Inverse problems are often ill-posed, meaning small data changes can lead to large solution variations. To mitigate this, regularization techniques like Tikhonov's are used to stabilize solutions. This study compares regularization methods applied to gravimetric inversion, including Smoothness Constraints, Total Variation, Discrete Cosine Transform (DCT), and Discrete Wavelet Transform (DWT) using Daubechies D4 wavelets. Optimization, particularly with Genetic Algorithms (GA), is used to find prism depths that best match observed anomalies. GA, inspired by natural selection, selects the best solutions to minimize the objective function. The results, evaluated through fit metrics and error analysis, show the effectiveness of all regularization methods and GA, with the Smoothness constraint performing best in synthetic models. For the real data model, all methods performed similarly.

We have been developing a novel interferometer formed directly between a Very Long Baseline Interferometry (VLBI) radio telescope and a Global Navigation Satellite Systems (GNSS) antenna/receiver. This interferometer is enabled by the High Rate Tracking Receiver (HRTR), a high-performance GNSS software-defined receiver that records baseband data similar to a VLBI receiving system. Of particular interest is the potential of the technique to produce precise local tie vectors directly between the geodetic reference points of VLBI and GNSS antennas for use in the determination of combination terrestrial reference frames. In this paper, we will describe the unique setup of this interferometer before discussing current developmental work in estimating local tie vectors from experimental data collected using this technique. In particular, much recent work has been in exploring the possibility of direct comparison of the GNSS and VLBI processing techniques. A significant advantage of GNSS-VLBI co-observation of GNSS satellites comes from the ability to cross-check position and clock estimates derived from differential measurements in VLBI-style processing with differential GNSS processing and absolute GNSS positioning with Precise Point Positioning (PPP). We show a preliminary VLBI analysis of interferometric data including a differential positioning solution using phase delays.

Artificial intelligence (AI) has become a buzz word since Google's AlphaGo beat a world champion in 2017. In the past five years, machine learning as a subset of the broader category of AI has obtained considerable attention in the research community of granular materials. This work offers a detailed review of the recent advances in machine learning-aided studies of granular materials from the particle-particle interaction at the grain level to the macroscopic simulations of granular flow. This work will start with the application of machine learning in the microscopic particle-particle interaction and associated contact models. Then, different neural networks for learning the constitutive behaviour of granular materials will be reviewed and compared. Finally, the macroscopic simulations of practical engineering or boundary value problems based on the combination of neural networks and numerical methods are discussed. We hope readers will have a clear idea of the development of machine learning-aided modelling of granular materials via this comprehensive review work.

Ground-penetrating radar (GPR) is a mature geophysical method that has gained increasing popularity in planetary science over the past decade. GPR has been utilised both for Lunar and Martian missions providing pivotal information regarding the near surface geology of Terrestrial planets. Within that context, numerous processing pipelines have been suggested to address the unique challenges present in planetary setups. These processing pipelines often require manual tuning resulting to ambiguous outputs open to non-unique interpretations. These pitfalls combined with the large number of planetary GPR data (kilometers in magnitude), highlight the necessity for automatic, objective and advanced processing and interpretation schemes. The current paper investigates the potential of deep learning for interpreting and processing GPR data. The one-shot multi-offset configuration is investigated via a coherent numerical case study, showcasing the potential of deep learning for A) reconstructing the dielectric distribution of the the near surface of Terrestrial planets, and B) filling missing or bad-quality traces. Special care was taken for the numerical data to be both realistic and challenging. Moreover, the generated synthetic data are properly labelled and made publicly available for training future data-driven pipelines and contributing towards developing pre-trained foundation models for GPR.

The methodology developed by McNamara and Bulland (2004) for computing Power Spectral Densities (PSDs) has gained popularity due to its low computational cost and reduction of spectral variance. This methodology is widely used in seismic noise studies and station performance evaluations and is implemented in tools like ISPAQ, MUSTANG, and PQLX. However, concerns have been raised about its appropriateness in certain contexts, particularly when high-resolution spectral detail is required. This study evaluates McNamara and Bulland's methodology by comparing it with Welch's method across three Alaskan stations with differing microseism conditions. When calculating seismic power across a band of frequencies--for example, the 5-10s secondary microseism--we find that both methodologies produce time series with nearly identical trends, albeit with slight differences in absolute power values. Our results demonstrate that McNamara and Bulland's methodology is fully appropriate for certain applications, specifically ones that rely on averaged seismic energy over a frequency band as opposed to a single discrete frequency.

The increased brightness temperature of young rocky protoplanets during their magma ocean epoch makes them potentially amenable to atmospheric characterization to distances from the solar system far greater than thermally equilibrated terrestrial exoplanets, offering observational opportunities for unique insights into the origin of secondary atmospheres and the near surface conditions of prebiotic environments. The Large Interferometer For Exoplanets (LIFE) mission will employ a space-based mid-infrared nulling interferometer to directly measure the thermal emission of terrestrial exoplanets. Here, we seek to assess the capabilities of various instrumental design choices of the LIFE mission concept for the detection of cooling protoplanets with transient high-temperature magma ocean atmospheres, in young stellar associations in particular. Using the LIFE mission instrument simulator (LIFEsim) we assess how specific instrumental parameters and design choices, such as wavelength coverage, aperture diameter, and photon throughput, facilitate or disadvantage the detection of protoplanets. We focus on the observational sensitivities of distance to the observed planetary system, protoplanet brightness temperature using a blackbody assumption, and orbital distance of the potential protoplanets around both G- and M-dwarf stars. Our simulations suggest that LIFE will be able to detect (S/N $\geq$ 7) hot protoplanets in young stellar associations up to distances of $\approx$100 pc from the solar system for reasonable integration times (up to $\sim$hours). Detection of an Earth-sized protoplanet orbiting a solar-sized host star at 1 AU requires less than 30 minutes of integration time. M-dwarfs generally need shorter integration times. The contribution from wavelength regions $<$6 $\mu$m is important for decreasing the detection threshold and discriminating emission temperatures.

The elementary theory of relaxation of the source cooling down after the main shock of an earthquake is presented axiomatically. The names of the objects under study are given and the relationships between them are determined. A new basic concept of earthquake source deactivation is introduced and a procedure for calculating the deactivation coefficient from aftershock frequency measurement data is indicated. An important property of the system is that the axioms do not contain statements regarding the actual process of source relaxation. From two simple axioms a number of meaningful statements (theorems) logically follow. A two-stage mode of source relaxation was discovered. The sharp transition between stages has the character of a bifurcation. It is shown that the classical Omori law has limited applicability. It describes the evolution of aftershocks only at the first stage of relaxation. The well-known Hirano-Utsu law is not applicable to describe aftershocks either at the first or second stages of relaxation. The conclusions of the elementary theory allow for interesting generalizations that expand the possibilities of experimental and theoretical study of the source. Keywords: earthquake source, foreshocks, main shock, aftershocks, source state, deactivation coefficient, bifurcation, Omori's law, Hirano-Utsu's law.

In recent years, Full-Waveform Inversion (FWI) has been extensively used to derive high-resolution subsurface velocity models from seismic data. However, due to the nonlinearity and ill-posed nature of the problem, FWI requires a good starting model to avoid producing non-physical solutions. Moreover, conventional optimization methods fail to quantify the uncertainty associated with the recovered solution, which is critical for decision-making processes. Bayesian inference offers an alternative approach as it directly or indirectly evaluates the posterior probability density function. For example, Markov Chain Monte Carlo (MCMC) methods generate multiple sample chains to characterize the solution's uncertainty. Despite their ability to theoretically handle any form of distribution, MCMC methods require many sampling steps; this limits their usage in high-dimensional problems with computationally intensive forward modeling, as is the FWI case. Variational Inference (VI), on the other hand, provides an approximate solution to the posterior distribution in the form of a parametric or non-parametric proposal distribution. Among the various algorithms used in VI, Stein Variational Gradient Descent (SVGD) is recognized for its ability to iteratively refine a set of samples to approximate the target distribution. However, mode and variance-collapse issues affect SVGD in high-dimensional inverse problems. This study aims to improve the performance of SVGD within the context of FWI by utilizing, for the first time, an annealed variant of SVGD and combining it with a multi-scale strategy. Additionally, we demonstrate that Principal Component Analysis (PCA) can be used to evaluate the performance of the optimization process. Clustering techniques are also employed to provide more rigorous and meaningful statistical analysis of the particles in the presence of multi-modal distributions.

Astronomical solutions provide calculated orbital and rotational parameters of solar system bodies based on the dynamics and physics of the solar system. Application of astronomical solutions in the Earth sciences has revolutionized our understanding in at least two areas of active research. (i) The Astronomical (or Milankovic) forcing of climate on time scales > ~10 kyr and (ii) the dating of geologic archives. The latter has permitted the development of the astronomical time scale, widely used today to reconstruct highly accurate geological dates and chronologies. The tasks of computing vs. applying astronomical solutions are usually performed by investigators from different backgrounds, which has led to confusion and recent inaccurate results on the side of the applications. Here we review astronomical solutions and Milankovic forcing in the Earth sciences, primarily aiming at clarifying the astronomical basis, applicability, and limitations of the solutions. We provide a summary of current up-to-date and outdated astronomical solutions and their valid time span. We discuss the fundamental limits imposed by dynamical solar system chaos on astronomical calculations and geological/astrochronological applications. We illustrate basic features of chaotic behavior using a simple mechanical system, i.e., the driven pendulum. Regarding so-called astronomical "metronomes", we point out that the current evidence does not support the notion of generally stable and prominent metronomes for universal use in astrochronology and cyclostratigraphy. We also describe amplitude and frequency modulation of astronomical forcing signals and the relation to their expression in cyclostratigraphic sequences. Furthermore, the various quantities and terminology associated with Earth's axial precession are discussed in detail. Finally, we provide some suggestions regarding practical considerations.

Single borehole radar detection suffers from azimuthal ambiguity, while borehole optical tests only provide information about the borehole wall. These limitations prevent either detection method from revealing the complete spatial patterns of rock fractures on their own. In this paper, we address these challenges by proposing a joint imaging method that combines the advantages of both borehole detection methods. Geological azimuthal parameters are extracted from optical images by fitting the fracture curves to sinusoidal functions. A 2D Kirchhoff time migration is then implemented using radar common offset gather. Up-dip and down-dip events are separated by the f-k transform or z-s transform, depending on their geometric relation. The complete fracture planes, including trend, dip angle, gap width, and extension length, are finally reconstructed in 3D space by mapping the migration profile using azimuthal information from optical images. The method is proven reliable and high-resolution through both numerical tests and real field data.

In subsurface imaging, learning the mapping from velocity maps to seismic waveforms (forward problem) and waveforms to velocity (inverse problem) is important for several applications. While traditional techniques for solving forward and inverse problems are computationally prohibitive, there is a growing interest in leveraging recent advances in deep learning to learn the mapping between velocity maps and seismic waveform images directly from data. Despite the variety of architectures explored in previous works, several open questions still remain unanswered such as the effect of latent space sizes, the importance of manifold learning, the complexity of translation models, and the value of jointly solving forward and inverse problems. We propose a unified framework to systematically characterize prior research in this area termed the Generalized Forward-Inverse (GFI) framework, building on the assumption of manifolds and latent space translations. We show that GFI encompasses previous works in deep learning for subsurface imaging, which can be viewed as specific instantiations of GFI. We also propose two new model architectures within the framework of GFI: Latent U-Net and Invertible X-Net, leveraging the power of U-Nets for domain translation and the ability of IU-Nets to simultaneously learn forward and inverse translations, respectively. We show that our proposed models achieve state-of-the-art (SOTA) performance for forward and inverse problems on a wide range of synthetic datasets, and also investigate their zero-shot effectiveness on two real-world-like datasets.

Modern climate change presents unprecedented challenges, posing critical crises that threaten sustainable development, human well-being, and planetary health. A significant concern is the potential for global warming to cause irreversible disruptions to the water-carbon cycle, a topic that remains underexplored. This study seeks to address a crucial knowledge gap by examining how increasing wet extremes impact ecosystem productivity. The research agenda focuses on three primary questions: 1) How do the intensity and duration of various wet extremes affect evapotranspiration across different watersheds and terrestrial biomes? 2) How do immediate and lagged responses to wet extremes vary across different biomes, and what insights do these temporal patterns provide about the causal and predictive relationships between wet extreme and evapotranspiration? 3) To what extent do watershed characteristics (such as soil properties, hydrological conditions, and vegetation factors) modulate the relationship between wet extremes and ecosystem productivity? As climate change alters precipitation patterns, understanding these complex ecosystem responses becomes crucial for developing adaptive strategies and improving food and water resource management.

Ionospheric irregularities and associated scintillations under geomagnetically active/quiet conditions have detrimental effects on the reliability and performance of space- and ground-based navigation satellite systems, especially over the low-latitude region. The current work investigates the low-latitude ionospheric irregularities using the phase screen theory and the corresponding temporal Power Spectral Density (PSD) analysis to present an estimate of the outer irregularity scale sizes over these locations. The study uses simultaneous L5 signal C/N$_o$ observations of NavIC (a set of GEO and GSO navigation satellite systems) near the northern crest of EIA (Indore: 22.52$^\circ$N, 75.92$^\circ$E, dip: 32.23$^\circ$N) and in between the crest and the dip equator (Hyderabad: 17.42$^\circ$N, 78.55$^\circ$E, dip: 21.69$^\circ$N). The study period (2017-2018) covers disturbed and quiet-time conditions in the declining phase of the solar cycle 24. The PSD analysis brings forward the presence of irregularities, of the order of a few hundred meters during weak-to-moderate and quiet-time conditions and up to a few km during the strong event, over both locations. The ROTI values validate the presence of such structures in the Indian region. Furthermore, only for the strong event, a time delay of scintillation occurrence over Indore, with values of 36 minutes and 50 minutes for NavIC satellites (PRNs) 5 and 6, respectively, from scintillation occurrence at Hyderabad is observed, suggesting a poleward evolution of irregularity structures. Further observations show a westward propagation of these structures on this day. This study brings forward the advantage of utilizing continuous data from the GEO and GSO satellite systems in understanding the evolution and propagation of the ionospheric irregularities over the low-latitude region.

This paper presents a proof of concept for spatial prediction of rock saturation probability using classifier ensemble methods on the example of the giant Groningen gas field. The stages of generating 1481 seismic field attributes and selecting 63 significant attributes are described. The effectiveness of the proposed method of augmentation of well and seismic data is shown, which increased the training sample by 9 times. On a test sample of 42 wells (blind well test), the results demonstrate good accuracy in predicting the ensemble of classifiers: the Matthews correlation coefficient is 0.7689, and the F1-score for the "gas reservoir" class is 0.7949. Prediction of gas reservoir thicknesses within the field and adjacent areas is made.

Life is a planetary feature that depends on its environment, but it has also strongly shaped the physical conditions on Earth, having created conditions highly suitable for a productive biosphere. Clearly, the second law of thermodynamics must apply to these dynamics as well, but how? What insights can we gain by placing life and its effects on planetary functioning in the context of the second law? In Kleidon (2010), I described a thermodynamic Earth system perspective by placing the functioning of the Earth system in terms of the second law. The Earth system is represented by a planetary hierarchy of energy transformations that are driven predominantly by incoming solar radiation, these transformations are constrained by the second law, but they are also modified by the feedbacks from various dissipative activities. It was then hypothesised that life evolves its dissipative activity to the limits imposed by this hierarchy and evolves feedbacks aimed at pushing these limits to higher levels of dissipative activity. Here I provide an update of this perspective. I first review applications to climate and global climate change to demonstrate its success in predicting magnitudes of physical processes, particularly regarding temperatures, heat redistribution and hydrological cycling. I then focus on the limits to dissipative activity of the biosphere. It would seem that the limitations by thermodynamics act indirectly by imposing limitations associated with transport and material exchange. I substantiate this interpretation and discuss the broader implications for habitability, the emergence and evolution of life, and the contemporary biosphere.

Due to the saturation of the body ($m_b$) and surface ($M_S$) earthquake magnitudes, the moment magnitude ($M_W$) is a more convenient parameter for representing earthquake energies. We use the HRVD data, including 18,569 earthquakes recorded from 1976 January 1 to 2010 December 31, for five regions on the Earth (\textbfAsia-Oceania, Europe, Africa, North America, and South America). Applying the error propagation technique on the moment tensor, we estimate the standard error for $M_W$ showing a distribution that deviates from Gaussian errors which recommends using the Monte Carlo Markov chain (MCMC) to obtain model parameters instead of the least square approach. We investigate the linear relationship between the body and surface waves in small ($\le$ 6.1) and large ($>$6.1) magnitudes via MCMC. The slope varies from 0.811 to 1.565 ($m_b - M_W$) and 0.566 to 0.971 ($M_S - M_W$) across five regions. For 18,569 global earthquakes, we obtained the slope for $m_b \le 6.1$ ($m_b>$6.1) 0.863 (1.374), while we find the slope for $M_S\le$6.1 ($M_S >$ 6.1) 0.581 (0.921). These transition relationships of the magnitude moments are useful for increasing earthquake monitoring capacity.

We present X-ray observations of the upper atmospheric density disturbance caused by the explosive eruption of the Hunga Tonga-Hunga Ha'apai (HTHH) volcano on 15 January 2022. From 14 January to 16 January, the Chinese X-ray astronomy satellite, Insight-HXMT, was observing the supernova remnant Cassiopeia A. The X-ray data obtained during Earth's atmospheric occultations allowed us to measure neutral densities in the altitude range of ~90-150 km. The density profiles above 110 km altitude obtained before the major eruption are in reasonable agreement with expectations by both GAIA and NRLMSIS 2.0 models. In contrast, after the HTHH eruption, a severe density depletion was found up to ~1,000 km away from the epicenter, and a relatively weak depletion extending up to ~7,000 km for over 8 hr after the eruption. In addition, density profiles showed wavy structures with a typical length scale of either ~20 km (vertical) or ~1,000 km (horizontal). This may be caused by Lamb waves or gravity waves triggered by the volcanic eruption.

Magnetotelluric deep learning (DL) inversion methods based on joint data-driven and physics-driven have become a hot topic in recent years. When mapping observation data (or forward modeling data) to the resistivity model using neural networks (NNs), incorporating the error (loss) term of the inversion resistivity's forward modeling response--which introduces physical information about electromagnetic field propagation--can significantly enhance the inversion accuracy. To efficiently achieve data-physical dual-driven MT deep learning inversion for large-scale 3-D MT data, we propose using DL forward modeling networks to compute this portion of the loss. This approach introduces pseudo-physical information through the forward modeling of NN simulation, further guiding the inversion network fitting. Specifically, we first pre-train the forward modeling networks as fixed forward modeling operators, then transfer and integrate them into the inversion network training, and finally optimize the inversion network by minimizing the multinomial loss. Theoretical experimental results indicate that despite some simulation errors in DL forward modeling, the introduced pseudo-physical information still enhances inversion accuracy and significantly mitigates the overfitting problem during training. Additionally, we propose a new input mode that involves masking and adding noise to the data, simulating the field data environment of 3-D MT inversion, thereby making the method more flexible and effective for practical applications.

Accurate flood modeling is crucial for effective analysis and forecasting. Full momentum hydrodynamic models often require extensive computational time, sometimes exceeding the forecast horizon. In contrast, low-complexity models, like local-inertial approximations, provide accurate results in subcritical flows but may have limited skillfulness in supercritical conditions. This paper explores two main aspects: (i) the impact of urban infrastructure on 2D hydrodynamic modeling without detailed sewer and drainage data, and (ii) the accuracy of 2D local-inertial modeling using three numerical schemes (original formulation, s-centered, and s-upwind) in a dam-break scenario on complex, flat terrain. The HydroPol2D model is benchmarked against HEC-RAS 2D full momentum solver. We present one numerical case study and three real-world scenarios in São Paulo, Brazil: a detention pond with a $1$ in $100$-year inflow, a highly urbanized catchment with a $1$ in $50$-year hyetograph, and a dam-break scenario threatening a coastal city of nearly 200,000 residents. Results show that the model accurately simulates internal boundary conditions, achieving peak errors under 5\% compared to HEC-RAS 2D. However, neglecting urban infrastructure can lead to a 17.5\% difference in peak discharges at the outlet and significant mismatches in hydrographs, with computational times nearly doubling. The dam-break scenario demonstrates good predictive performance for maximum flood depths (CSI = $0.95$ for the original model, $0.92$ for s-centered, and $0.89$ for s-upwind), though the model's lack of convective inertia results in faster flood wave propagation than the full momentum solver. Notably, HydroPol2D is 23 times faster than HEC-RAS 2D, making it well-suited for simulating dam collapses in forecasting systems and capable of modeling urban drainage infrastructure such as orifices, weirs, and pumps.

Full Waveform Inversion (FWI) is a vital technique for reconstructing high-resolution subsurface velocity maps from seismic waveform data, governed by partial differential equations (PDEs) that model wave propagation. Traditional machine learning approaches typically map seismic data to velocity maps by encoding seismic waveforms into latent embeddings and decoding them into velocity maps. In this paper, we introduce a novel framework that reframes FWI as a joint diffusion process in a shared latent space, bridging seismic waveform data and velocity maps. Our approach has two key components: first, we merge the bottlenecks of two separate autoencoders-one for seismic data and one for velocity maps-into a unified latent space using vector quantization to establish a shared codebook. Second, we train a diffusion model in this latent space, enabling the simultaneous generation of seismic and velocity map pairs by sampling and denoising the latent representations, followed by decoding each modality with its respective decoder. Remarkably, our jointly generated seismic-velocity pairs approximately satisfy the governing PDE without any additional constraint, offering a new geometric interpretation of FWI. The diffusion process learns to score the latent space according to its deviation from the PDE, with higher scores representing smaller deviations from the true solutions. By following this diffusion process, the model traces a path from random initialization to a valid solution of the governing PDE. Our experiments on the OpenFWI dataset demonstrate that the generated seismic and velocity map pairs not only exhibit high fidelity and diversity but also adhere to the physical constraints imposed by the governing PDE.

The ubiquitous phenomena of crystallization and melting occur in various geophysical contexts across many spatial and temporal scales. In particular, they take place in the iron core of terrestrial planets and moons, profoundly influencing their dynamics and magnetic field generation. Crystallization and melting entail intricate multiphase flows, buoyancy effects, and out-of-equilibrium thermodynamics, posing challenges for theoretical modeling and numerical simulations. Besides, due to the inaccessible nature of the planetary deep interior, our understanding relies on indirect data from seismology, mineral physics, geochemistry, and magnetism. Consequently, phase-change-driven flows in planetary cores constitute a compelling yet challenging area of research. This paper provides an overview of the role of laboratory fluid dynamics experiments in elucidating the solid-liquid phase change phenomena occurring thousands of kilometers beneath our feet and within other planetary depths, along with their dynamic repercussions. Drawing parallel with metallurgy, it navigates through all scales of phase change dynamics, from microscopic processes (nucleation and crystal growth) to macroscopic consequences (solid-liquid segregation and large-scale flows). The review delves into the two primary planetary solidification regimes, top-down and bottom-up, and elucidates the formation of mushy and/or slurry layers in the various relevant configurations. Additionally, it outlines remaining challenges, including insights from ongoing space missions poised to unveil the diverse planetary regimes.

GPR full-waveform inversion optimizes the subsurface property model iteratively to match the entire waveform information. However, the model gradients derived from wavefield continuation often contain errors, such as ghost values and excessively large values at transmitter and receiver points. Furthermore, models updated based on these gradients frequently exhibit unclear characterization of anomalous bodies or false anomalies, making it challenging to obtain accurate inversion results. To address these issues, we introduced a novel full-waveform inversion (FWI) framework that incorporates an embedded convolutional neural network (CNN) to adaptively filter model parameters and gradients. Specifically, we embedded the CNN module before the forward modeling process and ensured the entire FWI process remains differentiable. This design leverages the auto-grad tool of the deep learning library, allowing model values to pass through the CNN module during forward computation and model gradients to pass through the CNN module during backpropagation. Experiments have shown that filtering the model parameters during forward computation and the model gradients during backpropagation can ultimately yield high-quality inversion results.

In this brief paper, we analyze space weather events that occurred on May 11 and 12, 2024, from the perspective of an operational space weather center that provides advisories for civil aviation. One of the key metrics monitored by the center is the radiation dose rate at operational flight altitudes. A model implemented by the center provides the dose rate in real time. The model showed that dangerous levels were momentarily exceeded just above the usual 30,000 feet level during the events. This paper highlights differences in models used by various space weather centers, emphasizing the need for harmonization.

The recent development of deep learning (DL) methods for computer vision has been driven by the creation of open benchmark datasets on which new algorithms can be tested and compared with reproducible results. Although DL methods have many applications in geophysics, few real seismic datasets are available for benchmarking DL models, especially for denoising real data, which is one of the main problems in seismic data processing scenarios in the oil and gas industry. This article presents a benchmark dataset composed of synthetic seismic data corrupted with noise extracted from a filtering process implemented on real data. In this work, a comparison between two well-known DL-based denoising models is conducted on this dataset, which is proposed as a benchmark for accelerating the development of new solutions for seismic data denoising. This work also introduces a new evaluation metric that can capture small variations in model results. The results show that DL models are effective at denoising seismic data, but some issues remain to be solved.

Classical point process models, such as the epidemic-type aftershock sequence (ETAS) model, have been widely used for forecasting the event times and locations of earthquakes for decades. Recent advances have led to Neural Point Processes (NPPs), which promise greater flexibility and improvements over classical models. However, the currently-used benchmark dataset for NPPs does not represent an up-to-date challenge in the seismological community since it lacks a key earthquake sequence from the region and improperly splits training and testing data. Furthermore, initial earthquake forecast benchmarking lacks a comparison to state-of-the-art earthquake forecasting models typically used by the seismological community. To address these gaps, we introduce EarthquakeNPP: a collection of benchmark datasets to facilitate testing of NPPs on earthquake data, accompanied by a credible implementation of the ETAS model. The datasets cover a range of small to large target regions within California, dating from 1971 to 2021, and include different methodologies for dataset generation. In a benchmarking experiment, we compare three spatio-temporal NPPs against ETAS and find that none outperform ETAS in either spatial or temporal log-likelihood. These results indicate that current NPP implementations are not yet suitable for practical earthquake forecasting. However, EarthquakeNPP will serve as a platform for collaboration between the seismology and machine learning communities with the goal of improving earthquake predictability.

Altered and gangue quartz in hydrothermal veins from the Kubi Gold deposit in Dunkwa on Offin in the central region of Ghana are investigated for possible Au-associated indicator minerals and to provide the understanding and increase the knowledge of the mineral hosting and alteration processes in quartz. X-ray diffraction, air annealing furnace, differential scanning calorimetry, energy dispersive X-ray spectroscopy, and transmission electron microscopy have been applied on different quartz types outcropping from surface and bedrocks at the Kubi Gold Mining to reveal the material properties at different temperatures. From the diffraction results of the fresh and annealed quartz samples, we find that the samples contain indicator and the impurity minerals iron disulfide, biotite, titanium oxide, and magnetite. These minerals, under oxidation process between 574-1400 \degC temperatures experienced hematite alterations and a transformation from \alpha-quartz to e̱ta-quartz and further to cristobalite as observed from the calorimetry scans for hydrothermally exposed materials. The energy dispersive spectroscopy revealed elemental components of Fe, S, Mg, K, Al, Ti, Na, Si, O, and Ca contained in the samples, and these are attributed to the impurity phase minerals observed in the diffraction. The findings also suggest that during the hydrothermal flow regime, impurity minerals and metals can be trapped by voids and faults. Under favorable temperature conditions, the trapped minerals can be altered to change color at different depositional stages by oxidation and reduction processes leading to hematite alteration which is a useful indicator mineral in mineral exploration.

Air-sea exchange processes have been identified as essential for both short- and long-term atmospheric and ocean forecasts. The two phases of the fluid layer covering our planet interact across a vast range of scales that we need to explore to achieve a better understanding of the exchange processes. While satellites provide a distributed large-scale view of the sea surface situation, highly detailed measurements, e.g., from oceanographic towers, are necessarily local. The required intermediate solution (i.e., data that are both accurate and distributed) can be provided by swarms of miniature surface buoys. As size, weight, and cost are reduced, these can be deployed in large numbers to investigate specific processes that are at present only crudely parameterized in our models as a result of scarcity of good measurements. Perhaps the most crucial process is white-capping in stormy conditions, where air-sea exchanges are enhanced by one or two orders of magnitude. Other applications include wave-current interactions, wave-ice interactions, and plunging breakers in the coastal zone. Following a dedicated workshop, we summarize here the main findings and possibilities derived from the different approaches, and summarize the state of the art for a selection of miniature buoys. We list the solutions presented at the workshop, as well as other similar buoys, with their main characteristics and range of application. We describe the various possibilities of practical use and the scientific and engineering problems to be solved. Looking to the future, we also point out where the present technological improvements are leading to.

Garnet has been widely used to decipher the pressure-temperature-time history of rocks, but its physical properties such as elasticity and diffusion are strongly affected by trace amounts of hydrogen. Experimental measurements of H diffusion in garnet are limited to room pressure. We use atomistic simulations to study H diffusion in perfect and defective garnet lattices, focusing on protonation defects at the Si and Mg sites, which are shown to be energetically favored. The ab-initio simulation of H diffusion is computationally challenging due to a transient trapping of H, which is overcome with machine learning techniques by training a deep neural network that encodes the interatomic potential. Our results show high mobility of hydrogen in defect free garnet lattices, whereas H diffusivity is significantly diminished in defective lattices. Tracer simulation focusing on H alone highlights the vital role of atomic vibrations of heavier atoms like Mg in the untrapping of H atoms. Two regimes of H diffusion are identified: a diffuser-dominated regime at high hydrogen content with low activation energies, due to saturation of vacancies by hydrogen, and a vacancy-dominated regime at low hydrogen content with high activation energies, due to trapping of H atoms at vacancy sites. These regimes account for experimental observations, such as a H-concentration dependent diffusivity and the discrepancy in activation energy between deprotonation and D-H exchange experiments. This study underpins the crucial role of vacancies in H diffusion and demonstrates the utility of machine-learned interatomic potentials in studying kinetic processes in the Earth's interior.

Monitoring networks increasingly aim to assimilate data from a large number of diverse sensors covering many sensing modalities. Bayesian optimal experimental design (OED) seeks to identify data, sensor configurations, or experiments which can optimally reduce uncertainty and hence increase the performance of a monitoring network. Information theory guides OED by formulating the choice of experiment or sensor placement as an optimization problem that maximizes the expected information gain (EIG) about quantities of interest given prior knowledge and models of expected observation data. Therefore, within the context of seismo-acoustic monitoring, we can use Bayesian OED to configure sensor networks by choosing sensor locations, types, and fidelity in order to improve our ability to identify and locate seismic sources. In this work, we develop the framework necessary to use Bayesian OED to optimize a sensor network's ability to locate seismic events from arrival time data of detected seismic phases at the regional-scale. Bayesian OED requires four elements: 1) A likelihood function that describes the distribution of detection and travel time data from the sensor network, 2) A Bayesian solver that uses a prior and likelihood to identify the posterior distribution of seismic events given the data, 3) An algorithm to compute EIG about seismic events over a dataset of hypothetical prior events, 4) An optimizer that finds a sensor network which maximizes EIG. Once we have developed this framework, we explore many relevant questions to monitoring such as: how to trade off sensor fidelity and earth model uncertainty; how sensor types, number, and locations influence uncertainty; and how prior models and constraints influence sensor placement.

Ground motion prediction equations (GMPEs) play a key role in seismic hazard assessment (SHA). Considering the seismo-tectonic, geophysical and geotectonic characteristics of a target region, all the GMPEs may not be suitable in predicting the observed ground motion effectively. With a fairly large number of published GMPEs, the selection and ranking of suitable GMPEs for the design of logic trees in SHA for a particular target region has become a necessity of late. This paper presents a detailed quantitative evaluation of performance of 16 GMPEs against recorded ground motion data in two target regions, characterized by distinct seismo-tectonic, geophysical and geotectonical nature. The data set comprises of 465 three-component spectral accelerograms corresponding to 122 earthquake events. The suitability of a GMPE is tested by two widely accepted data-driven statistical methods, namely, likelihood (LH) and log-likelihood (LLH) method. Different suites of GMPEs are shown suitable for different periods of interest. The results will be useful to scientists and engineers for microzonation and estimation of seismic design parameters for the design of earthquake-resistant structures in these regions.

Stress changes in the Earth's crust induced by an earthquake typically trigger other earthquakes. Identifying such aftershocks is an important step in seismic hazard assessment but has remained challenging. Here we demonstrate that earthquake catalogs without triggering can lead to spurious Omori-Utsu and productivity relations for the commonly used nearest-neighbor declustering method. However, we show that the robustness of the Omori-Utsu exponent on newly introduced parameters of the method allows to identify true aftershocks. For a natural swarm catalog, which lacks a clear distinction between triggered and background events, we reveal that the catalog is dominated by aftershocks.

Downhole earthquake monitoring, without the complex effects from the near surface, can record more and better seismic data than monitoring on surface. The San Andreas Fault Observatory at Depth (SAFOD) is a borehole observatory equipped with different instruments inside to study the earthquake mechanism of the San Andreas fault at Parkfield, California. During April to May in 2005, Paulsson deployed an 80-level 3-component geophone array in the SAFOD main hole, and continuously recorded seismic data for about 13 days. We located 125 local earthquakes from the borehole earthquake monitoring data using a homogeneous velocity model and compared it with 35 earthquakes' locations from surface earthquake monitoring by the United State Geological Survey (USGS) during the same monitoring time. The borehole earthquake locating is assumably more accurate in the borehole's vicinity. We also compared the result with 1,074 earthquakes' locations from the surface earthquake monitoring in the last 9 years from 2015 to 2024. The hypocenters from our nearly 2 weeks' borehole earthquake monitoring form similar structures as that from the 9 years' surface earthquake monitoring by the USGS.

Induced seismicity has emerged as a source of a significant earthquake hazard associated with recent development of unconventional energy resources. Therefore, it is imperative to develop stochastic models that can accurately describe the observed seismicity rate and forecast its evolution. In this study, a mechanism suggested by linear response theory is incorporated into a stochastic earthquake model to account for changes in the seismicity rate. It is derived that the induced rate can be modelled as a convolution of the forcing, related to fluid injection operations, and a specific response kernel. The model is incorporated into a Bayesian framework to compute the probabilities for the occurrence of the largest expected events during future time intervals. The applicability of the model is illustrated by analyzing the injection and seismicity data at the Geysers geothermal field in California. The suggested approach provides further insight into the probabilistic assessment of earthquake hazard associated with fluid injection operations. It also can be used for probing the rheological properties of the subsurface by analysing the inherent characteristic time-scales associated with the subsurface response to external forcing.

Salt precipitation triggered by the evaporation of formation brine into injected supercritical CO2 can cause injectivity and containment issues in near-wellbore regions. Predicting the distribution of precipitated salts and their impact on near-wellbore properties remains challenging. This study investigates the influence of surface wettability on CO2-induced halite precipitation and growth within hydrophilic and hydrophobic microfluidic chips designed to mimic rock-structure porous geometries. A series of high-pressure brine-CO2 flow experiments, direct microscopic observations, and detailed image processing were conducted to explore how substrate wettability affects salt precipitation. The experiments show that wettability markedly controls residual brine relocation, film flow movement, solute consumption, and salt formation. Tracking halite precipitation dynamics revealed distinct crystal formation signatures: hydrophilic chips exhibited irregular, larger aggregation patches, while the hydrophobic network showed more numerous, smaller, and limited aggregations. Large individual crystals were observed in both chips, with a notable dominance in the hydrophobic one. Crystallization dynamics varied, with nucleation and growth occurring earlier, progressing faster, and forming bulkier aggregates in the hydrophilic chip. Despite these differences, the three identified temporal stages of brine evaporation and halite surface coverage were comparable. Spatial analysis along the chips indicated that crystal aggregation properties, such as size and distribution, were position-dependent, with hydrophilic chips exhibiting greater probabilistic variability. These observations underscore the impact of surface wettability on salt precipitation through brine accessibility and capillarity, with implications for mitigating and remediating salt issues in saline aquifers.

The striped patterns in rocks, characterized by thin dark bands and thick light bands, are a result of the natural processes that have shaped the rocks through in situ metamorphism. Notably, the dark and light bands are composed of similar minerals, with the dark bands containing additional trapped materials such as fluid pockets or extra mineral grains. Here we employed the Gray-Scott Reaction Diffusion model to investigate the interaction between two sets of virtual chemicals, denoted as 'u' and 'v'. The differing diffusion and reaction rates of 'u' and 'v' chemical orders lead to the formation of 180 degree out-of-phase chemical domains, resulting in striped patterns. Utilizing the Gray-Scott model in this manner, we gain valuable insights into the early microscopic stages of these geologically significant striped patterns in metamorphic rocks.

A theoretical model describing the one-dimensional large strain consolidation of the modified Cam Clay soil is presented in this paper. The model is based on the Lagrangian formulation, and is capable of featuring the variability of soil compressibility (inherently so due to the direct incorporation of the specific Cam Clay plasticity model) and permeability, as well as the impact of overconsolidation ratio. The derivation starts from the establishment of the incremental stress-strain relations for both purely elastic and elastoplastic deformations under one-dimensional compression condition, and thereafter the coefficients of compressibility/volume change that are essential to the consolidation analysis. The governing partial differential equation is then neatly deduced in conjunction with the continuity and equilibrium conditions for the soil, with the vertical effective stress being the privileged unknown to be solved for. Subsequently, semi-analytical solution to the developed rigorous poroelastoplastic large strain consolidation model is obtained and verified with the ABAQUS finite element numerical results. Parametric analyses are finally provided to investigate in detail the influences of the soil overconsolidation ratio, large strain configuration, and the variability of the soil permeability on the calculated one-dimensional consolidation response.

Accurate calibration of gravity adaptation zones is of great significance in fields such as underwater navigation, geophysical exploration, and marine engineering. With the increasing application of gravity field data in these areas, traditional calibration methods based on single features are becoming inadequate for capturing the complex characteristics of gravity fields and addressing the intricate interrelationships among multidimensional data. This paper proposes an attention-enhanced algorithm for gravity adaptation zone calibration. By introducing an attention mechanism, the algorithm adaptively fuses multidimensional gravity field features and dynamically assigns feature weights, effectively solving the problems of multicollinearity and redundancy inherent in traditional feature selection methods, significantly improving calibration accuracy and robustness.In addition, a large-scale gravity field dataset with over 10,000 sampling points was constructed, and Kriging interpolation was used to enhance the spatial resolution of the data, providing a reliable data foundation for model training and evaluation. We conducted both qualitative and quantitative experiments on several classical machine learning models (such as SVM, GBDT, and RF), and the results demonstrate that the proposed algorithm significantly improves performance across these models, outperforming other traditional feature selection methods. The method proposed in this paper provides a new solution for gravity adaptation zone calibration, showing strong generalization ability and potential for application in complex environments. The code is available at \hrefthis link https://github.com/hulnifox/RF-ATTN.

Significant progress has been made in understanding planetary core dynamics using numerical models of rotating convection (RC) in spherical shell geometry. However, the behaviour of forces in these models within various dynamic regimes of RC remains largely unknown. Directional anisotropy, scale dependence, and the role of dynamically irrelevant gradient contributions in incompressible flows complicate the representation of dynamical balances. In this study, we systematically compare integrated and scale-dependent representations of mean and fluctuation forces and curled forces (which contain no gradient contributions) separately for the three components ($\hat{r},\hat{\theta},\hat{\phi}$). The analysis is performed with simulations in a range of convective supercriticality $Ra/Ra_c=1.2-1967$ and Ekman number $E=10^{-3}-10^{-6}$, with fixed Prandtl number $Pr=1$, no-slip and fixed flux boundaries. We exclude $10$ velocity boundary layers from the spherical shell boundaries, to get a consistent representation of the dynamics between the force and curled force balance. Radial, azimuthal and co-latitudinal components exhibit distinct balances. The total magnitudes of the mean forces and mean curled forces exhibit a primary thermal wind (TW) balance; the corresponding fluctuating forces are in a quasi-geostrophic (QG) primary balance, while the fluctuating curled forces transition from a Viscous-Archimedean-Coriolis balance to an Inertia-Viscous-Archimedean-Coriolis balance with increasing $Ra/Ra_c$. The curled force balances are weakly scale-dependent compared to the forces and do not show clear cross-over scales. The fluctuating force and curled force balances are consistent with three regimes of RC (weakly nonlinear, rapidly rotating, and weakly rotating), but do not exhibit sharp changes with $Ra/Ra_c$, which inhibits the identification of precise regime boundaries from these balances.

Ensemble data assimilation techniques form an indispensable part of numerical weather prediction. As the ensemble size grows and model resolution increases, the amount of required storage becomes a major issue. Data compression schemes may come to the rescue not only for operational weather prediction, but also for weather history archives. In this paper, we present the design and implementation of an easy-to-use framework for evaluating the impact of lossy data compression in large scale ensemble data assimilation. The framework leverages robust statistical qualifiers to determine which compression parameters can be safely applied to the climate variables. Furthermore, our proposal can be used to apply the best parameters during operation, while monitoring data integrity. We perform an exemplary study on the Lorenz96 model to identify viable compression parameters and achieve a 1/3 saving in storage space and an effective speedup of 6% per assimilation cycle, while monitoring the state integrity.

On 24 September 2023, the Origins, Spectral Interpretation, Resource Identification, and Security Regolith Explorer (OSIRIS-REx) Sample Return Capsule entered the Earth's atmosphere after successfully collecting samples from an asteroid. The known trajectory and timing of this return provided a rare opportunity to strategically instrument sites to record geophysical signals produced by the capsule as it traveled at hypersonic speeds through the atmosphere. We deployed two optical-fiber distributed acoustic sensing (DAS) interrogators to sample over 12 km of surface-draped, fiber-optic cables along with six co-located seismometer-infrasound sensor pairs, spread across two sites near Eureka, NV. This campaign-style rapid deployment is the first reported recording of a sample return capsule entry with any distributed fiber optic sensing technology. The DAS interrogators recorded an impulsive arrival with an extended coda which had features that were similar to recordings from both the seismometers and infrasound sensors. While the signal-to-noise of the DAS data was lower than the seismic-infrasound data, the extremely dense spacing of fiber-optic sensors allowed for more phases to be clearly distinguished and the continuous transformation of the wavefront as it impacted the ground could be visualized. Unexpectedly, the DAS recordings contain less low-frequency content than is present in both the seismic and infrasound data. The deployment conditions strongly affected the recorded DAS data, in particular, we observed that fiber selection and placement exert strong controls on data quality.

Convection is a ubiquitous process driving geophysical/astrophysical fluid flows, which are typically strongly constrained by planetary rotation on large scales. A celebrated model of such flows, rapidly rotating Rayleigh-Bénard convection, has been extensively studied in direct numerical simulations (DNS) and laboratory experiments, but the parameter values attainable by state-of-the-art methods are limited to moderately rapid rotation (Ekman numbers $Ek\gtrsim10^{-8}$), while realistic geophysical/astrophysical $Ek$ are significantly smaller. Asymptotically reduced equations of motion, the nonhydrostatic quasi-geostrophic equations (NHQGE), describing the flow evolution in the limit $Ek\to 0$, do not apply at finite rotation rates. The geophysical/astrophysical regime of small but finite $Ek$ therefore remains currently inaccessible. Here, we introduce a new, numerically advantageous formulation of the Navier-Stokes-Boussinesq equations informed by the scalings valid for $Ek\to0$, the \textitrescaled incompressible Navier-Stokes equations (RiNSE). We solve the RiNSE using a spectral quasi-inverse method resulting in a sparse, fast algorithm to perform efficient DNS in this previously unattainable parameter regime. We validate our results against the literature across a range of $Ek$, and demonstrate that the algorithmic approaches taken remain accurate and numerically stable at $Ek$ as low as $10^{-15}$. Like the NHQGE, the RiNSE derive their efficiency from adequate conditioning, eliminating spurious growing modes that otherwise induce numerical instabilities at small $Ek$. We show that the time derivative of the mean temperature is inconsequential for accurately determining the Nusselt number in the stationary state, significantly reducing the required simulation time, and demonstrate that full DNS using RiNSE agree with the NHQGE at very small $Ek$.

Geomagnetic indices derived from ground magnetic measurements characterize the intensity of solar-terrestrial interaction. The \textitKp index derived from multiple magnetic observatories at mid-latitude has commonly been used for space weather operations. Yet, its temporal cadence is low and its intensity scale is crude. To derive a new generation of geomagnetic indices, it is desirable to establish a geomagnetic `baseline' that defines the quiet-level of activity without solar-driven perturbations. We present a new approach for deriving a baseline that represents the time-dependent quiet variations focusing on data from Chambon-la-Forêt, France. Using a filtering technique, the measurements are first decomposed into the above-diurnal variation and the sum of 24h, 12h, 8h, and 6h filters, called the daily variation. Using correlation tools and SHapley Additive exPlanations, we identify parameters that dominantly correlate with the daily variation. Here, we predict the daily `quiet' variation using a long short-term memory neural network trained using at least 11 years of data at 1h cadence. This predicted daily quiet variation is combined with linear extrapolation of the secular trend associated with the intrinsic geomagnetic variability, which dominates the above-diurnal variation, to yield a new geomagnetic baseline. Unlike the existing baselines, our baseline is insensitive to geomagnetic storms. It is thus suitable for defining geomagnetic indices that accurately reflect the intensity of solar-driven perturbations. Our methodology is quick to implement and scalable, making it suitable for real-time operation. Strategies for operational forecasting of our geomagnetic baseline 1 day and 27 days in advance are presented.