Showing papers in "arXiv: Geophysics in 2017"

Machine Learning Predicts Laboratory Earthquakes

Abstract: Forecasting fault failure is a fundamental but elusive goal in earthquake science. Here we show that by listening to the acoustic signal emitted by a laboratory fault, machine learning can predict the time remaining before it fails with great accuracy. These predictions are based solely on the instantaneous physical characteristics of the acoustical signal, and do not make use of its history. Surprisingly, machine learning identifies a signal emitted from the fault zone previously thought to be low-amplitude noise that enables failure forecasting throughout the laboratory quake cycle. We hypothesize that applying this approach to continuous seismic data may lead to significant advances in identifying currently unknown signals, in providing new insights into fault physics, and in placing bounds on fault failure times.

Natural Time, Nowcasting and the Physics of Earthquakes: Estimation of Seismic Risk to Global Megacities

Abstract: This paper describes the use of the idea of natural time to propose a new method for characterizing the seismic risk to the world's major cities at risk of earthquakes. Rather than focus on forecasting, which is the computation of probabilities of future events, we define the term seismic nowcasting, which is the computation of the current state of seismic hazard in a defined geographic region.

Three-dimensional density structure of La Soufrie\`ere de Guadeloupe lava dome from simultaneous muon radiographies and gravity data

Abstract: Muon imaging has recently emerged as a powerful method to complement standard geophysical tools in the understanding of the Earth's subsurface. Muon measurements can yield a radiography of the average density along the muon path, allowing to image large volumes of a geological body from a single observation point. Here we jointly invert muon data from three simultaneous telescope acquisitions together with gravity data to estimate the three-dimensional density structure of the La Soufriere de Guadeloupe lava Dome. Our unique dataset allows us to achieve an unprecedented spatial resolution with this novel technique. The retrieved density model reveals an extensive, low-density anomaly where the most active part of the volcanic hydrothermal system is located, supporting previous studies that indicate this region as the most likely to be involved in a partial edifice collapse.

The Cosmic-Ray Neutron Rover - Mobile Surveys of Field Soil Moisture and the Influence of Roads

Abstract: Measurements of root-zone soil moisture across spatial scales of tens to thousands of meters have been a challenge for many decades. The mobile application of Cosmic-Ray Neutron Sensing (CRNS) is a promising approach to measure field soil moisture non-invasively by surveying large regions with a ground-based vehicle. Recently, concerns have been raised about a potentially biasing influence of local structures and roads. We employed neutron transport simulations and dedicated experiments to quantify the influence of different road types on the CRNS measurement. We found that the presence of roads introduces a bias in the CRNS estimation of field soil moisture compared to non-road scenarios. However, this effect becomes insignificant at distances beyond a few meters from the road. Measurements from the road could overestimate the field value by up to 40 % depending on road material, width, and the surrounding field water content. The bias could be successfully removed with an analytical correction function that accounts for these parameters. Additionally, an empirical approach is proposed that can be used on-the-fly without prior knowledge of field soil moisture. Tests at different study sites demonstrated good agreement between road-effect corrected measurements and field soil moisture observations. However, if knowledge about the road characteristics is missing, any measurements on the road could substantially reduce the accuracy of this method. Our results constitute a practical advancement of the mobile CRNS methodology, which is important for providing unbiased estimates of field-scale soil moisture to support applications in hydrology, remote sensing, and agriculture.

Revealing the cluster of slow transients behind a large slow slip event

Abstract: Capable of reaching similar magnitudes to large megathrust earthquakes ($M_w>7$), slow slip events play a major role in accommodating tectonic motion on plate boundaries. These slip transients are the slow release of built-up tectonic stress that are geodetically imaged as a predominantly aseismic rupture, which is smooth in both time and space. We demonstrate here that large slow slip events are in fact a cluster of short-duration slow transients. Using a dense catalog of low-frequency earthquakes as a guide, we investigate the $M_w7.5$ slow slip event that occurred in 2006 along the subduction interface 40~km beneath Guerrero, Mexico. We show that while the long-period surface displacement as recorded by GPS suggests a six month duration, motion in the direction of tectonic release only sporadically occurs over 55 days and its surface signature is attenuated by rapid relocking of the plate interface.These results demonstrate that our current conceptual model of slow and continuous rupture is an artifact of low-resolution geodetic observations of a superposition of small, clustered slip events. Our proposed description of slow slip as a cluster of slow transients implies that we systematically overestimate the duration $T$ and underestimate the moment magnitude $M$ of large slow slip events.

Characterization of Marcellus shale fracture properties through size effect tests and computations

Abstract: Mechanical characterization of shale-like rocks requires understanding the scaling of the measured properties to enable the extrapolation from small scale laboratory tests to field study. In this paper, the size effect of Marcellus shale was analyzed, and the fracture properties were obtained through size effect tests. A number of fracture tests were conducted on Three-Point-Bending (TPB) specimens with increasing size. Test results show that the nominal strength decreases with increasing specimen size, and can be fitted well by Bazant's Size Effect Law (SEL). It is shown that SEL accounts for the effects of both specimen size and geometry, allowing an accurate identification of the initial fracture energy of the material, Gf, and the effective Fracture Process Zone (FPZ) length, cf. The obtained fracture properties were verified by the numerical simulations of the investigated specimens using standard Finite Element technique with cohesive model. Significant anisotropy was observed in the fracture properties determined in three principal notch orientations: arrester, divider, and short-transverse. The size effect of the measured structural strength and apparent fracture toughness was discussed. Neither strength-based criterion which neglects size effect, nor classic LEFM which does not account for the finiteness of the FPZ can predict the reported size effect data, and nonlinear fracture mechanics of the quasibrittle type is instead applicable.

Estimations of the seismic pressure noise on Mars determined from Large Eddy Simulations and demonstration of pressure decorrelation techniques for the InSight mission

Abstract: The atmospheric pressure fluctuations on Mars induce an elastic response in the ground that creates a ground tilt, detectable as a seismic signal on the InSight seismometer SEIS. The seismic pressure noise is modeled using Large Eddy Simulations of the wind and surface pressure at the InSight landing site and a Green's function ground deformation approach that is subsequently validated via a detailed comparison with two other methods based on Sorrells' theory (Sorrels 1971; Sorrels et al. 1971). The horizontal acceleration as a result of the ground tilt due to the LES turbulence-induced pressure fluctuations are found to be typically ~2 - 40 nm/s^2 in amplitude, whereas the direct horizontal acceleration is two orders of magnitude smaller and is thus negligible in comparison. The vertical accelerations are found to be ~0.1 - 6 nm/s^2 in amplitude. We show that under calm conditions, a single-pressure measurement is representative of the large-scale pressure field (to a distance of several kilometers), particularly in the prevailing wind direction. However, during windy conditions, small-scale turbulence results in a reduced correlation between the pressure signals, and the single-pressure measurement becomes less representative of the pressure field. Nonetheless, the correlation between the seismic signal and the pressure signal is found to be higher for the windiest period because the seismic pressure noise reflects the atmospheric structure close to the seismometer. In the same way that we reduce the atmospheric seismic signal by making use of a pressure sensor that is part of the InSight APSS (Auxiliary Payload Sensor Suite), we also the use the synthetic noise data obtained from the LES pressure field to demonstrate a decorrelation strategy.

High Performance Clocks and Gravity Field Determination

Abstract: Time measured by an ideal clock crucially depends on the gravitational potential and velocity of the clock according to general relativity. Technological advances in manufacturing high-precision atomic clocks have rapidly improved their accuracy and stability over the last decade that approached the level of 10$^{-18}$. Based on a fully relativistic description of the background gravitational physics, we discuss the impact of those highly-precise clocks on the realization of reference frames and time scales used in geodesy. We discuss the current definitions of basic geodetic concepts and come to the conclusion that the advances in clocks and other metrological technologies will soon require the re-definition of time scales or, at least, clarification to ensure their continuity and consistent use in practice. The relative frequency shift between two clocks is directly related to the difference in the values of the gravity potential at the points of clock's localization. According to general relativity the relative accuracy of clocks in 10$^{-18}$ is equivalent to measuring the gravitational red shift effect between two clocks with the height difference amounting to 1 cm. We show how clock measurements can provide geopotential numbers for the realization of gravity-field-related height systems and can resolve discrepancies in classically-determined height systems as well as between national height systems. Another application of clocks is the direct use of observed potential differences for the improved recovery of regional gravity field solutions. Finally, clock measurements for space-borne gravimetry are analyzed along with closely-related deficiencies of this method like an extra-ordinary knowledge of the spacecraft velocity, etc. For all these applications besides the near-future prospects, we also discuss the challenges that are related to using those novel clock data in geodesy.

Joint inversion in hydrogeophysics and near-surface geophysics

Abstract: The near-surface environment is often too complex to enable inference of hydrological and environmental variables using one geophysical data type alone. Joint inversion and coupled inverse modeling involving numerical flow- and transport simulators have, in the last decade, played important roles in pushing applications towards increasingly challenging targets. Joint inversion of geophysical data that is based on structural constraints is often favored over model coupling based on explicit petrophysical relationships. More specifically, cross-gradient joint inversion has been applied to a wide range of near-surface applications and geophysical data types. To infer hydrological subsurface properties, the most appropriate approach is often to use temporal changes in geophysical data that can be related to hydrological state variables. This allows using geophysical data as indirect hydrological observables, while the coupling with a flow- and transport simulator ensures physical consistency. Future research avenues include investigating the validity of different coupling strategies at various scales, the spatial statistics of near-surface petrophysical relationships, the influence of the model conceptualization, fully probabilistic joint inversions, and how to include complex prior information in the joint inversion.

Machine learning for graph-based representations of three-dimensional discrete fracture networks

Abstract: Structural and topological information play a key role in modeling flow and transport through fractured rock in the subsurface. Discrete fracture network (DFN) computational suites such as dfnWorks are designed to simulate flow and transport in such porous media. Flow and transport calculations reveal that a small backbone of fractures exists, where most flow and transport occurs. Restricting the flowing fracture network to this backbone provides a significant reduction in the network's effective size. However, the particle tracking simulations needed to determine the reduction are computationally intensive. Such methods may be impractical for large systems or for robust uncertainty quantification of fracture networks, where thousands of forward simulations are needed to bound system behavior. In this paper, we develop an alternative network reduction approach to characterizing transport in DFNs, by combining graph theoretical and machine learning methods. We consider a graph representation where nodes signify fractures and edges denote their intersections. Using random forest and support vector machines, we rapidly identify a subnetwork that captures the flow patterns of the full DFN, based primarily on node centrality features in the graph. Our supervised learning techniques train on particle-tracking backbone paths found by dfnWorks, but run in negligible time compared to those simulations. We find that our predictions can reduce the network to approximately 20% of its original size, while still generating breakthrough curves consistent with those of the original network.

Objective Estimation of Spatially Variable Parameters of Epidemic Type Aftershock Sequence Model: Application to California

Abstract: The ETAS model is widely employed to model the spatio-temporal distribution of earthquakes, generally using spatially invariant parameters. We propose an efficient method for the estimation of spatially varying parameters, using the Expectation-Maximization (EM) algorithm and spatial Voronoi tessellation ensembles. We use the Bayesian Information Criterion (BIC) to rank inverted models given their likelihood and complexity and select the best models to finally compute an ensemble model at any location. Using a synthetic catalog, we also check that the proposed method correctly inverts the known parameters. We apply the proposed method to earthquakes included in the ANSS catalog that occurred within the time period 1981-2015 in a spatial polygon around California. The results indicate a significant spatial variation of the ETAS parameters. We find that the efficiency of earthquakes to trigger future ones (quantified by the branching ratio) positively correlates with surface heat flow. In contrast, the rate of earthquakes triggered by far-field tectonic loading or background seismicity rate shows no such correlation, suggesting the relevance of triggering possibly through fluid-induced activation. Furthermore, the branching ratio and background seismicity rate are found to be uncorrelated with hypocentral depths, indicating that the seismic coupling remains invariant of hypocentral depths in the study region. Additionally, triggering seems to be mostly dominated by small earthquakes. Consequently, the static stress change studies should not only focus on the Coulomb stress changes caused by specific moderate to large earthquakes but also account for the secondary static stress changes caused by smaller earthquakes.

The sinking of the El Faro: predicting real world rogue waves during Hurricane Joaquin

Abstract: We present a study on the prediction of rogue waves during the 1-hour sea state of Hurricane Joaquin when the Merchant Vessel El Faro sank east of the Bahamas on October 1, 2015. High-resolution hindcast of hurricane-generated sea states and wave simulations are combined with novel probabilistic models to quantify the likelihood of rogue wave conditions. The data suggests that the El Faro vessel was drifting at an average speed of approximately~$2.5$~m/s prior to its sinking. As a result, we estimated that the probability that El Faro encounters a rogue wave whose crest height exceeds 14 meters while drifting over a time interval of 10~(50) minutes is $\sim1/400$~$(1/130)$. The largest simulated rogue wave has similar generating mechanism and characteristics of the Andrea, Draupner and Killard rogue waves as the constructive interference of elementary waves enhanced by bound nonlinearities.

High spatial resolution trace element determination of geological samples by laser ablation quadrupole plasma mass spectrometry: implications for glass analysis in volcanic products

Abstract: Increasing the spatial resolution of Laser Ablation Inductively Coupled Plasma Mass Spectrometry (LA-ICP-MS) is a challenge in microanalysis of geological samples. Smaller sizes for the laser beam will allow for (1) high resolution determination of trace element compositions, (2) accurate estimation of crystal/melt partition coefficients, (3) detailed characterization of diffusion profiles, and (4) analysis of fine volcanic glasses. Here, we report about the figures of merit for LA-ICP Quadrupole MS down to a spatial resolution of 5 {\mu}m. This study highlights the possibility to achieve suitable limits of detection, accuracy and precision for geological samples even at spatial resolutions of the order of 5 {\mu}m. At a beam size of 15 {\mu}m precision (measured as one sigma) and accuracy (expressed as relative deviation from the reference value) are of the order of 10%. At a smaller beam size of 8um, precision decreases to 15% for concentration above 1.7 {\mu}g g-1. As the beam size is decreased to ~5 {\mu}m, precision declines to about 15% and 20% for concentrations above 10 {\mu}g g-1 using 42Ca and 29Si as internal standard, respectively. Accuracy is better or equal to 10% and 13 % at beam sizes of 15 and 10 {\mu}m respectively. When the spatial resolution is increased to 8 {\mu}m, accuracy remains better than 15% and 20% for 42Ca and 29Si as internal standard, respectively. We employed such high-resolution techniques to volcanic glasses in ash particles of the 2010 Eyjafjallajokull eruption. Our results are well consistent with the previously reported data obtained at lower spatial resolution, supporting the reliability of the method.

A simple hysteretic constitutive model for unsaturated flow

Abstract: In this paper we present a constitutive model to describe unsaturated flow that considers the hysteresis phenomena. This constitutive model provides simple mathematical expressions for both saturation and hydraulic conductivity curves, and a relationship between permeability and porosity. The model is based on the assumption that the porous media can be represented by a bundle of capillary tubes with throats or "ink-bottles" and a fractal pore size distribution. Under these hypotheses, hysteretic curves are obtained for saturation and relative hydraulic conductivity in terms of pressure head. However, a non-hysteretic relationship is obtained when relative hydraulic conductivity is expressed as a function of saturation. The proposed relationship between permeability and porosity is similar to the well-known Kozeny-Carman equation but depends on the fractal dimension. The performance of the constitutive model is tested against different sets of experimental data and previous models. In all of the cases the proposed expressions fit fairly well the experimental data and predicts values of permeability and hydraulic conductivity better than others models.

New Variable Compliance Method for Estimating In-Situ Stress and Leak-Off from DFIT Data

Abstract: It is shown that using Carter leak-off is an oversimplification that leads to significant errors in the interpretation of DFIT data. Most importantly, this article reveals that previous methods of estimating minimum in-situ stress often lead to significant over or underestimates. Based on our modeling and simulation results, we propose a much more accurate and reliable method to estimate the minimum in-situ stress and fracture pressure dependent leak-off rate.

Earthquake precursors in the light of peroxy defects theory: critical review of systematic observations

Abstract: The starting point of the present review is to acknowledge that there are innumerable reports of non-seismic types of earthquake precursory phenomena that are intermittent and seem not to occur systematically, while associated reports are not widely accepted by the geoscience community at large because no one could explain their origins. We review a unifying theory for a solid-state mechanism, based on decades of research bridging semi-conductor physics, chemistry and rock physics. A synthesis has emerged that all pre-earthquake phenomena could trace back to one fundamental physical process: the activation of electronic charges (electrons and positive holes) in rocks subjected to ever-increasing tectonic stresses prior to any major seismic activity, via the rupture of peroxy bonds. In the second part of the review, we critically examine satellite and ground station data, recorded before past large earthquakes, as they have been claimed to provide evidence that precursory signals tend to become measurable days, sometimes weeks before the disasters. We review some of the various phenomena that can be directly predicted by the peroxy defect theory , namely, radon gas emanations, corona discharges, thermal infrared emissions, air ionization, ion and electron content in the ionosphere, and electro-magnetic anomalies. Our analysis demonstrates the need for further systematic investigations, in particular with strong continuous statistical testing of the relevance and confidence of the precursors. Only then, the scientific community will be able to assess and improve the performance of earthquake forecasts.

Experimental constraints on the rheology, eruption and emplacement dynamics of analog lavas comparable to Mercury's northern volcanic plains

Abstract: We present new viscosity measurements of a synthetic silicate system considered an analogue for the lava erupted on the surface of Mercury. In particular, we focus on the northern volcanic plains (NVP), which correspond to the largest lava flows on Mercury and possibly in the Solar System. High-temperature viscosity measurements were performed at both superliquidus (up to 1736 K) and subliquidus conditions (1569-1502 K) to constrain the viscosity variations as a function of crystallinity (from 0 to 28\%) and shear rate (from 0.1 to 5 s 1). Melt viscosity shows moderate variations (4-16 Pa s) in the temperature range of 1736-1600 K. Experiments performed below the liquidus temperature show an increase in viscosity as shear rate decreases from 5 to 0.1 s 1, resulting in a shear thinning behavior, with a decrease in viscosity of 1 log unit. The low viscosity of the studied composition may explain the ability of NVP lavas to cover long distances, on the order of hundreds of kilometers in a turbulent flow regime. Using our experimental data we estimate that lava flows with thickness of 1, 5, and 10 m are likely to have velocities of 4.8, 6.5, and 7.2 m/s, respectively, on a 5 degree ground slope. Numerical modeling incorporating both the heat loss of the lavas and its possible crystallization during emplacement allows us to infer that high effusion rates (>10,000 m3/s) are necessary to cover the large distances indicated by satellite data from the MErcury Surface, Space ENvironment, GEochemistry, and Ranging spacecraft.

Revisiting slope influence in turbulent bedload transport: consequences for vertical flow structure and transport rate scaling

Abstract: Gravity-driven turbulent bedload transport has been extensively studied over the past century in regard to its importance for Earth surface processes such as natural riverbed morphological evolution. In the present contribution, the influence of the longitudinal channel inclination angle on gravity-driven turbulent bedload transport is studied in an idealised framework considering steady and uniform flow conditions. From an analytical analysis based on the two-phase continuous equations, it is shown that : (i) the classical slope correction of the critical Shields number is based on an erroneous formulation of the buoyancy force, (ii) the influence of the slope is not restricted to the critical Shields number but affects the whole transport formula and (iii) pressure-driven and gravity-driven turbulent bedload transport are not equivalent from the slope influence standpoint. Analysing further the granular flow driving mechanisms, the longitudinal slope is shown to not only influence the fluid bed shear stress and the resistance of the granular bed, but also to affect the fluid flow inside the granular bed - responsible for the transition from bedload transport to debris flow. The relative influence of these coupled mechanisms allows us to understand the evolution of the vertical structure of the granular flow and to predict the transport rate scaling law as a function of a rescaled Shields number. The theoretical analysis is validated with coupled fluid-discrete element simulations of idealised gravity-driven turbulent bedload transport, performed over a wide range of Shields number values, density ratios and channel inclination angles. In particular, all the data are shown to collapse onto a master curve when considering the sediment transport rate as a function of the proposed rescaled Shields number.

Hydrogeophysical characterization of transport processes in fractured rock by combining push-pull tracer tests and single-hole ground penetrating radar

Abstract: The in situ characterization of transport processes in fractured media is particularly challenging due to the considerable spatial uncertainty on tracer pathways and dominant controlling processes, such as dispersion, channeling, trapping, matrix diffusion, ambient and density driven flows. We attempted to reduce this uncertainty by coupling push-pull tracer experiments with single-hole ground penetrating radar (GPR) time-lapse imaging. The experiments involved different injection fractures, chaser volumes and resting times, and were performed at the fractured rock research site of Ploemeur in France (H+ network, this http URL). For the GPR acquisitions, we used both fixed and moving antenna setups in a borehole that was isolated with a flexible liner. During the fixed-antenna experiment, time-varying GPR reflections allowed us to track the spatial and temporal dynamics of the tracer during the push-pull experiment. During the moving antenna experiments, we clearly imaged the dominant fractures in which tracer transport took place, fractures in which the tracer was trapped for longer time periods, and the spatial extent of the tracer distribution (up to 8 m) at different times. This demonstrated the existence of strongly channelized flow in the first few meters and radial flow at greater distances. By varying the resting time of a given experiment, we identified regions affected by density-driven and ambient flow. These experiments open up new perspectives for coupled hydrogeophysical inversion aimed at understanding transport phenomena in fractured rock formations.

On the anomalous changes of seismicity and geomagnetic field prior to the 2011 $M_w$ 9.0 Tohoku earthquake

Abstract: Xu et al. [J. Asian Earth Sci. {\bf 77}, 59-65 (2013)] It has just been reported that approximately 2 months prior to the $M_w$9.0 Tohoku earthquake that occurred in Japan on 11 March 2011 anomalous variations of the geomagnetic field have been observed in the vertical component at a measuring station about 135 km from the epicenter for about 10 days (4 to 14 January 2011). Here, we show that this observation is in striking agreement with independent recent results obtained from natural time analysis of seismicity in Japan. In particular, this analysis has revealed that an unprecedented minimum of the order parameter fluctuations of seismicity was observed around 5 January 2011, thus pointing to the initiation at that date of a strong precursory Seismic Electric Signals activity accompanied by the anomalous geomagnetic field variations. Starting from this date, natural time analysis of the subsequent seismicity indicates that a strong mainshock was expected in a few days to one week after 08:40 LT on 10 March 2011.

Heat Flow and Boundary Heterogeneity in Rotating Convection

Abstract: Convection in the metallic cores of terrestrial planets is likely to be subjected to lateral variations in heat flux through the outer boundary imposed by creeping flow in the overlying silicate mantles. Boundary anomalies can significantly influence global diagnostics of core convection when the Rayleigh number, $Ra$, is weakly supercritical; however, little is known about the strongly supercritical regime appropriate for planets. We perform numerical simulations of rapidly rotating convection in a spherical shell geometry and impose two patterns of boundary heat flow heterogeneity: a hemispherical $Y_1^1$ spherical harmonic pattern; and one derived from seismic tomography of Earth's lower mantle. We consider Ekman numbers $10^{-4} \le E \le 10^{-6}$, flux-based Rayleigh numbers up to $\sim800$ times critical, and Prandtl number unity. The amplitude of the lateral variation in heat flux is characterised by $q_{\rm L}^*=0, 2.3, 5.0$, the peak-to-peak amplitude of the outer boundary heat flux divided by its mean. We find that the Nusselt number, $Nu$, can be increased by up to $\sim25$\% relative to the equivalent homogeneous case due to boundary-induced correlations between the radial velocity and temperature anomalies near the top of the shell. The $Nu$ enhancement tends to become greater as the amplitude and length scale of the boundary heterogeneity are increased and as the system becomes more supercritical. This $Ra$ dependence can steepen the $Nu \propto Ra^\gamma$ scaling in the rotationally dominated regime, with $\gamma$ for our most extreme case approximately $20\%$ greater than the equivalent homogeneous scaling. Therefore, it may be important to consider boundary heterogeneity when extrapolating numerical results to planetary conditions.

Three-Dimensional Numerical Modeling of Shear Stimulation of Naturally Fractured Reservoirs

Abstract: Shear dilation based hydraulic stimulations enable exploitation of geothermal energy from reservoirs with inadequate initial permeability. While contributing to enhancing the reservoir's permeability, hydraulic stimulation processes may lead to undesired seismic activity. Here, we present a three dimensional numerical model aiming to increase understanding of this mechanism and its consequences. The fractured reservoir is modeled as a network of explicitly represented large scale fractures immersed in a permeable rock matrix. The numerical formulation is constructed by coupling three physical processes: fluid flow, fracture deformation, and rock matrix deformation. For flow simulations, the discrete fracture matrix model is used, which allows the fluid transport from high permeable conductive fractures to the rock matrix and vice versa. The mechanical behavior of the fractures is modeled using a hyperbolic model with reversible and irreversible deformations. Linear elasticity is assumed for the mechanical deformation and stress alteration of the rock matrix. Fractures are modeled as lower dimensional surfaces embodied in the domain, subjected to specific governing equations for their deformation along the tangential and normal directions. Both the fluid flow and momentum balance equations are approximated by finite volume discretizations. The new numerical model is demonstrated considering a three dimensional fractured formation with a network of 20 explicitly represented fractures. The effects of fluid exchange between fractures and rock matrix on the permeability evolution and the generated seismicity are examined for test cases resembling realistic reservoir conditions.

Seismic wave propagation in icy ocean worlds

Abstract: Seismology was developed on Earth and shaped our model of the Earth's interior over the 20th century. With the exception of the Philae lander, all in situ extraterrestrial seismological effort to date was limited to other terrestrial planets. All have in common a rigid crust above a solid mantle. The coming years may see the installation of seismometers on Europa, Titan and Enceladus, so it is necessary to adapt seismological concepts to the setting of worlds with global oceans covered in ice. Here we use waveform analyses to identify and classify wave types, developing a lexicon for icy ocean world seismology intended to be useful to both seismologists and planetary scientists. We use results from spectral-element simulations of broadband seismic wavefields to adapt seismological concepts to icy ocean worlds. We present a concise naming scheme for seismic waves and an overview of the features of the seismic wavefield on Europa, Titan, Ganymede and Enceladus. In close connection with geophysical interior models, we analyze simulated seismic measurements of Europa and Titan that might be used to constrain geochemical parameters governing the habitability of a sub-ice ocean.

Comparing Topology of Engineered and Natural Drainage Networks

Abstract: We investigated the scaling and topology of engineered urban drainage networks (UDNs) in two cities, and further examined UDN evolution over decades. UDN scaling was analyzed using two power-law characteristics widely employed for river networks: (1) Hack's law of length ($L$)-area ($A$) scaling [$L \propto A^{h}$], and (2) exceedance probability distribution of upstream contributing area $(\delta)$ [$P(A\geq \delta) \sim a \delta^{-\epsilon}$]. For the smallest UDNs ($ \text{km}^2$), length-area scales linearly ($h\sim 1$), but power-law scaling emerges as the UDNs grow. While $P(A\geq \delta)$ plots for river networks are abruptly truncated, those for UDNs display exponential tempering [$P(A\geq \delta) \>\text{=}\> a \delta^{-\epsilon}\exp(-c\delta)$]. The tempering parameter $c$ decreases as the UDNs grow, implying that the distribution evolves in time to resemble those for river networks. However, the power-law exponent $\epsilon$ for large UDNs tends to be slightly larger than the range reported for river networks. Differences in generative processes and engineering design constraints contribute to observed differences in the evolution of UDNs and river networks, including subnet heterogeneity and non-random branching.

Inferring transport characteristics in a fractured rock aquifer by combining single-hole GPR reflection monitoring and tracer test data

Abstract: Investigations of solute transport in fractured rock aquifers often rely on tracer test data acquired at a limited number of observation points. Such data do not, by themselves, allow detailed assessments of the spreading of the injected tracer plume. To better understand the transport behavior in a granitic aquifer, we combine tracer test data with single-hole ground-penetrating radar (GPR) reflection monitoring data. Five successful tracer tests were performed under various experimental conditions between two boreholes 6 m apart. For each experiment, saline tracer was injected into a previously identified packed-off transmissive fracture while repeatedly acquiring single-hole GPR reflection profiles together with electrical conductivity logs in the pumping borehole. By analyzing depth-migrated GPR difference images together with tracer breakthrough curves and associated simplified flow and transport modeling, we estimate (1) the number, the connectivity, and the geometry of fractures that contribute to tracer transport, (2) the velocity and the mass of tracer that was carried along each flow path, and (3) the effective transport parameters of the identified flow paths. We find a qualitative agreement when comparing the time evolution of GPR reflectivity strengths at strategic locations in the formation with those arising from simulated transport. The discrepancies are on the same order as those between observed and simulated breakthrough curves at the outflow locations. The rather subtle and repeatable GPR signals provide useful and complementary information to tracer test data acquired at the outflow locations and may help us to characterize transport phenomena in fractured rock aquifers.

Combining machine learning techniques, microanalyses and large geochemical datasets for tephrochronological studies in complex volcanic areas: new age constraints for the Pleistocene magmatism of Central Italy

Abstract: Characterization, correlation and provenance determination of tephra samples in sedimentary sections (tephrochronological studies) are powerful tools for establishing ages of depositional events, volcanic eruptions, and tephra dispersion. Despite the large literature and the advancements in this research field, the univocal attribution of tephra deposits to specific volcanic sources remains too often elusive. In this contribution, we test the application of a machine learning technique named Support Vector Machine to attempt shedding new light upon tephra deposits related to one of the most complex and debated volcanic regions on Earth: the Pliocene-Pleistocene magmatism in Italy. The machine learning algorithm was trained using one of the most comprehensive global petrological databases (GEOROC); 17 chemical elements including major and selected trace elements were chosen as input parameters. We first show the ability of support vector machines in discriminating among different Pliocene-Pleistocene volcanic provinces in Italy and then apply the same methodology to determine the volcanic source of tephra samples occurring in the Caio outcrop, an Early Pleistocene sedimentary section located in Central Italy.

Storage and dissipation of energy in Prabhakar viscoelasticity

Abstract: In this paper, after a brief review of the physical notion of quality factor in viscoelasticity, we present a complete discussion of the attenuation processes emerging in the Maxwell-Prabhakar model, recently developed by Giusti and Colombaro. Then, taking profit of some illuminating plots, we discuss some potential connections between the presented model and the modern mathematical modelling of seismic processes.

Forecasting the magnitude and onset of El Nino based on climate network

Abstract: El Nino is probably the most influential climate phenomenon on interannual time scales. It affects the global climate system and is associated with natural disasters and serious consequences in many aspects of human life. However, the forecasting of the onset and in particular the magnitude of El Nino are still not accurate, at least more than half a year in advance. Here, we introduce a new forecasting index based on network links representing the similarity of low frequency temporal temperature anomaly variations between different sites in the El Nino 3.4 region. We find that significant upward trends and peaks in this index forecast with high accuracy both the onset and magnitude of El Nino approximately 1 year ahead. The forecasting procedure we developed improves in particular the prediction of the magnitude of El Nino and is validated based on several, up to more than a century long, datasets.

Magmatic focusing to mid-ocean ridges: the role of grain size variability and non-Newtonian viscosity

Abstract: Melting beneath mid-ocean ridges occurs over a region that is much broader than the zone of magmatic emplacement to form the oceanic crust. Magma is focused into this zone by lateral transport. This focusing has typically been explained by dynamic pressure gradients associated with corner flow, or by a sub-lithospheric channel sloping upward toward the ridge axis. Here we discuss a novel mechanism for magmatic focusing: lateral transport driven by gradients in compaction pressure within the asthenosphere. These gradients arise from the co-variation of melting rate and compaction viscosity. The compaction viscosity, in previous models, was given as a function of melt fraction and temperature. In contrast, we show that the viscosity variations relevant to melt focusing arise from grain-size variability and non-Newtonian creep. The asthenospheric distribution of melt fraction predicted by our models provides an improved ex- planation of the electrical resistivity structure beneath one location on the East Pacific Rise. More generally, although grain size and non-Newtonian viscosity are properties of the solid phase, we find that in the context of mid-ocean ridges, their effect on melt transport is more profound than their effect on the mantle corner-flow.

Geochemical discrimination and characteristics of magmatic tectonic settings; a machine learning-based approach

Abstract: Geochemically discriminating between magmatism in different tectonic settings remains a fundamental part of understanding the processes of magma generation within the Earth's mantle. Here, we present an approach where machine-learning (ML) methods are used for quantitative tectonic discrimination and feature selection using global geochemical datasets containing data for volcanic rocks generated in eight different tectonic settings. This study uses support vector machine, random forest, and sparse multinomial regression (SMR) approaches. All these ML methods with data for 20 elements and 5 isotopic ratios allowed the successful geochemical discrimination between igneous rocks formed in eight different tectonic settings with a discriminant ratio better than 83% for all settings barring oceanic plateaus and back-arc basins. SMR is a particularly powerful and interpretable ML method because it quantitatively identifies geochemical signatures that characterize the tectonic settings of interest and the characteristics of each sample as a probability of the membership of the sample for each setting. We also present the most representative basalt composition for each tectonic setting. The new data provide reference points for future geochemical discussions. Our results indicate that at least 17 elements and isotopic ratios are required to characterize each tectonic setting, suggesting that geochemical tectonic discrimination cannot be achieved using only a small number of elemental compositions and/or isotopic ratios. The results show that volcanic rocks formed in different tectonic settings have unique geochemical signatures, indicating that both volcanic rock geochemistry and magma generation processes are closely connected to the tectonic setting.