pinecone/core-2020-05-10-deduplication huggingface dataset:
The mechanical coupling of fluid-filled granular material under shear

International audienceThe coupled mechanics of fluid-filled granular media controls the physics of many Earth systems, for example saturated soils, fault gouge, and landslide shear zones. It is well established that when the pore fluid pressure rises, the shear resistance of fluid-filled granular systems decreases, and, as a result, catastrophic events such as soil liquefaction, earthquakes, and accelerating landslides may be triggered. Alternatively, when the pore pressure drops, the shear resistance of these geosystems increases. Despite the great importance of the coupled mechanics of grain-fluid systems, the basic physics that controls this coupling is far from understood. Fundamental questions that must be addressed include: what are the processes that control pore fluid pressurization and depressurization in response to deformation of the granular skeleton? and how do variations of pore pressure affect the mechanical strength of the grains skeleton? To answer these questions, a formulation for the pore fluid pressure and flow has been developed from mass and momentum conservation, and is coupled with a granular dynamics algorithm that solves the grain dynamics, to form a fully coupled model. The pore fluid formulation reveals that the evolution of pore pressure obeys viscoelastic rheology in response to pore space variations. Under undrained conditions elastic-like behavior dominates and leads to a linear relationship between pore pressure and overall volumetric strain. Viscous-like behavior dominates under well-drained conditions and leads to a linear relationship between pore pressure and volumetric strain rate. Numerical simulations reveal the possibility of liquefaction under drained and initially over-compacted conditions, which were often believed to be resistant to liquefaction. Under such conditions liquefaction occurs during short compactive phases that punctuate the overall dilative trend. In addition, the previously recognized generation of elevated pore pressure under undrained compactive conditions is observed. Simulations also show that during liquefaction events stress chains are detached, the external load becomes completely supported by the pressurized pore fluid, and shear resistance vanishes
Effects of Surface Geology on Seismic Ground Motion Deduced from Ambient-Noise Measurements in the Town of Avellino, Irpinia Region (Italy)

The effects of surface geology on ground motion
provide an important tool in seismic hazard studies. It is well
known that the presence of soft sediments can cause amplification
of the ground motion at the surface, particularly when there is a
sharp impedance contrast at shallow depth. The town of Avellino is
located in an area characterised by high seismicity in Italy, about
30 km from the epicentre of the 23 November 1980, Irpinia
earthquake (M = 6.9). No earthquake recordings are available in
the area. The local geology is characterised by strong heterogeneity,
with impedance contrasts at depth. We present the results
from seismic noise measurements carried out in the urban area of
Avellino to evaluate the effects of local geology on the seismic
ground motion. We computed the horizontal-to-vertical (H/V)
noise spectral ratios at 16 selected sites in this urban area for which
drilling data are available within the first 40 m of depth. A Rayleigh
wave inversion technique using the peak frequencies of the
noise H/V spectral ratios is then presented for estimating Vs
models, assuming that the thicknesses of the shallow soil layers are
known. The results show a good correspondence between experimental
and theoretical peak frequencies, which are interpreted in
terms of sediment resonance. For one site, which is characterised
by a broad peak in the horizontal-to-vertical spectral-ratio curve,
simple one-dimensional modelling is not representative of the
resonance effects. Consistent variations in peak amplitudes are seen
among the sites. A site classification based on shear-wave velocity
characteristics, in terms of Vs30, cannot explain these data. The
differences observed are better correlated to the impedance contrast
between the sediments and basement. A more detailed investigation
of the physical parameters of the subsoil structure, together with
earthquake data, are desirable for future research, to confirm these
data in terms of site response
Tsunami Wave Analysis and Possibility of Splay Fault Rupture During the 2004 Indian Ocean Earthquake

The 2004 Indian Ocean tsunami was observed by two satellites, close in space and time, that traversed the Indian ocean two hours after the Sumatra-Andaman earthquake, but which observed different tsunami lead wave morphologies. The earlier satellite, Jason-1, recorded a lead wave with two peaks of similar amplitude and wavelength, while the later satellite, TOPEX/Poseidon, recorded a lead wave with only one longer wavelength uplift. To resolve this disparity, we examine the travel paths of long wavelength waves over the seafloor bathymetry. Waves traveling from the margin will traverse significantly different paths to arrive at the two satellite transects. The result is that the satellites are sensitive to different parts of the margin; Jason-1 is highly sensitive to the margin in the area of the epicenter, while TOPEX is sensitive to a more northerly section. By developing solutions of the ocean gravity wave equations, accounting for dispersion, we show that the double peak of the Jason-1 satellite observations are consistent with coseismic rupture of a splay fault of limited along-strike extent, located north of Simeulue Island. The doubly peaked morphology can be reproduced with co-activation of the subduction zone interface and the splay fault, which creates a seafloor uplift pattern with two distinct areas of uplift. The Jason-1 satellite is sensitive to a splay fault in this portion of the margin, whereas the TOPEX satellite would not be significantly affected by this uplift pattern. By back-projecting satellite observation points to the margin, we constrain the location of the proposed splay fault and find that it correlates with a bathymetric high. The aftershock locations, uplift of corals on Simeulue Island and a fault scarp on Pulau Salaut Besar are also consistent with the activation of a splay fault in the area delimited by the back-projection. Our work also shows that it is critical to fully capture gravity wave dispersion in order to represent features of the lead wave profile that may not be as well characterized by the shallow water (long-wavelength) model. It is also necessary to account for dispersion so as to precisely assess wavefront travel times; this leads us to conclude that the rupture must have reached very near to the trench and propagated with an updip rupture velocity of order 2.0 km/s or more.Earth and Planetary SciencesEngineering and Applied Science
Modeling the Phase Transition Associated with Melting Snow in a 1D Kinematic Framework: Sensitivity to the Microphysics

A simple 1D kinematic cloud model coupled to a two-moment bulk microphysics scheme is used to perform quasi-idealized simulations of snow, with a prescribed upper boundary snow field based on observed radar reflectivity and temperature, falling into a low-level melting layer. The model realistically simulates the formation of a nearly isothermal layer below the melting level, the surface precipitation rate, and the phase transition from liquid to solid, consistent with observations for this case. A series of test runs is performed to examine the sensitivity of modeling the timing and duration of the phase transition period to details of specific parameterization aspects related to snow in the microphysics scheme. The sensitivity tests include varying the number of prognostic moments, the mass–diameter relation, the fall velocity–diameter relation, the treatment of aggregation, and the lower limit for the slope of the size distribution. It is shown that the simulated transition period, for such a case with the initial melting level being close to the surface, can be quite sensitive to model parameters specified within realistic ranges and/or ranges within our physical understanding
Coupled Simulation of Seismic Wave Propagation and Failure Phenomena by Use of an MPS Method

The failure of brittle materials, for example glasses and rock masses, is commonly observed to be discontinuous. It is, however, difficult to simulate these phenomena by use of conventional numerical simulation methods, for example the finite difference method or the finite element method, because of the presence of computational grids or elements artificially introduced before the simulation. It is, therefore, important for research on such discontinuous failures in science and engineering to analyze the phenomena seamlessly. This study deals with the coupled simulation of elastic wave propagation and failure phenomena by use of a moving particle semi-implicit (MPS) method. It is simple to model the objects of analysis because no grid or lattice structure is necessary. In addition, lack of a grid or lattice structure makes it simple to simulate large deformations and failure phenomena at the same time. We first compare analytical and MPS solutions by use of Lamb’s problem with different offset distances, material properties, and source frequencies. Our results show that analytical and numerical seismograms are in good agreement with each other for 20 particles in a minimum wavelength. Finally, we focus our attention on the Hopkinson effect as an example of failure induced by elastic wave propagation. In the application of the MPS, the algorithm is basically the same as in the previous calculation except for the introduction of a failure criterion. The failure criterion applied in this study is that particle connectivity must be disconnected when the distance between the particles exceeds a failure threshold. We applied the developed algorithm to a suspended specimen that was modeled as a long bar consisting of thousands of particles. A compressional wave in the bar is generated by an abrupt pressure change on one edge. The compressional wave propagates along the interior of the specimen and is visualized clearly. At the other end of the bar, the spalling of the bar is reproduced numerically, and a broken piece of the bar is formed and falls away from the main body of the bar. Consequently, these results show that the MPS method effectively reproduces wave propagation and failure phenomena at the same time
Dynamic Development of Hydrofracture

International audienceMany natural examples of complex joint and vein networks in layered sedimentary rocks are hydro-fractures that form by a combination of pore fluid overpressure and tectonic stresses. In this paper, a two-dimensional hybrid hydro-mechanical formulation is proposed to model the dynamic development of natural hydrofractures. The numerical scheme combines a Discrete Element Model (DEM) framework that represents a porous solid medium with a supplementary Darcy based pore-pressure diffusion as continuum description for the fluid. This combination yields a porosity controlled coupling between an evolving fracture network and the associated hydraulic field. The model is tested on some basic cases of hydro-driven fracturing commonly found in nature i.e., fracturing due to local fluid overpressure in rocks subjected to hydrostatic and nonhydrostatic tectonic loadings. In our models we find that seepage forces created by hydraulic pressure gradients together with poroelastic feedback upon discrete fracturing play a significant role in subsurface rock deformation. These forces manipulate the growth and geometry of hydrofractures in addition to tectonic stresses and the mechanical properties of the porous rocks. Our results show characteristic failure patterns that reflect different tectonic and lithological conditions and are qualitatively consistent with existing analogue and numerical studies as well as field observations. The applied scheme is numerically efficient, can be applied at various scales and is computational cost effective with the least involvement of sophisticated mathematical computation of hydrodynamic flow between the solid grains
Comparative simulations of the 2011 Tohoku tsunami with MOST and Cliffs

Wet-dry interface which appreciably increases accuracy of tsunami simulations
with the MOST model (Method of Splitting Tsunamis, adapted by the NOAA for
tsunami forecasting operations) is highlighted here with a few comparative
simulations. The new solution, termed Cliffs, exceeds all MOST versions in
accuracy of computing later waves. This is demonstrated with simulation of the
Tohoku-2011 tsunami to Monterey Bay, CA, and into fiords, bays, and inlets of
southeastern Alaska, followed by comparison with tide gage records.Comment: 6 pages, 4 figures; highlights one aspect of a paper in the proces
The seismotectonics of the Po Plain (northern Italy): tectonic diversity in a blind faulting domain

We present a systematic and updated overview of a seismotectonic model for the Po
Plain (northern Italy). This flat and apparently quiet tectonic domain is in fact rather
active as it comprises the shortened foreland and foredeep of both the Southern Alps
and the Northern Apennines. Assessing its seismic hazard is crucial due to the
concentration of population, industrial activities and critical infrastructures, but it is also
complicated because a) the region is geologically very diverse, and b) nearly all
potential seismogenic faults are buried beneath a thick blanket of Pliocene-Pleistocene
sediments, and hence can be investigated only indirectly.
Identifying and parameterizing the potential seismogenic faults of the Po Plain requires
proper consideration of their depth, geometry, kinematics, earthquake potential and
location with respect to the two confronting orogens. To this end we subdivided them
into four main homogeneous groups. Over the past 15 years we developed new
strategies for coping with this diversity, resorting to different data and modeling
approaches as required by each individual fault group.
The most significant faults occur beneath the thrust fronts of the Ferrara-Romagna and
Emilia arcs, which correspond to the most advanced and buried portions of the
Northern Apennines and were the locus of the destructive May 2012 earthquake
sequence. The largest known Po Plain earthquake, however, occurred on an elusive
reactivated fault cutting the Alpine foreland south of Verona. Significant earthquakes
are expected to be generated also by a set of transverse structures segmenting the
thrust system, and by the deeper ramps of the Apennines thrusts.
The new dataset is intended to be included in the next version of the Database of
Seismogenic Sources (DISS; http://diss.rm.ingv.it/diss/, version 3.2.0, developed and
maintained by INGV) to improve completeness of potential sources for seismic hazard
assessment
Complexity Phenomena and ROMA of the Earth’s Magnetospheric Cusp, Hydrodynamic Turbulence, and the Cosmic Web

“Dynamic complexity” is a phenomenon observed for a nonlinearly interacting system within which multitudes of different sizes of large scale coherent structures emerge, resulting in a globally nonlinear stochastic behavior vastly different from that which could be surmised from the underlying equations of interaction. A characteristic of such nonlinear, complex phenomena is the appearance of intermittent fluctuating events with the mixing and distribution of correlated structures on all scales. We briefly review here a relatively recent method, ROMA (rank-ordered multifractal analysis), explicitly developed for analysis of the intricate details of the distribution and scaling of such types of intermittent structure. This method is then used for analysis of selected examples related to the dynamic plasmas of the cusp region of the Earth’s magnetosphere, velocity fluctuations of classical hydrodynamic turbulence, and the distribution of the structures of the cosmic gas obtained by use of large-scale, moving mesh simulations. Differences and similarities of the analyzed results among these complex systems will be contrasted and highlighted. The first two examples have direct relevance to the Earth’s environment (i.e., geoscience) and are summaries of previously reported findings. The third example, although involving phenomena with much larger spatiotemporal scales, with its highly compressible turbulent behavior and the unique simulation technique employed in generating the data, provides direct motivation for applying such analysis to studies of similar multifractal processes in extreme environments of near-Earth surroundings. These new results are both exciting and intriguing.National Science Foundation (U.S.)Seventh Framework Programme (European Commission) (FP7/2007-2013 Grant Agreement 313038/STORM
Earth’s Rotation: A Challenging Problem in Mathematics and Physics

A suitable knowledge of the orientation and motion of the Earth in space is a common need in various fields. That knowledge has been ever necessary to carry out astronomical observations, but with the advent of the space age, it became essential for making observations of satellites and predicting and determining their orbits, and for observing the Earth from space as well. Given the relevant role it plays in Space Geodesy, Earth rotation is considered as one of the three pillars of Geodesy, the other two being geometry and gravity. Besides, research on Earth rotation has fostered advances in many fields, such as Mathematics, Astronomy and Geophysics, for centuries. One remarkable feature of the problem is in the extreme requirements of accuracy that must be fulfilled in the near future, about a millimetre on the tangent plane to the planet surface, roughly speaking. That challenges all of the theories that have been devised and used to-date; the paper makes a short review of some of the most relevant methods, which can be envisaged as milestones in Earth rotation research, emphasizing the Hamiltonian approach developed by the authors. Some contemporary problems are presented, as well as the main lines of future research prospected by the International Astronomical Union/International Association of Geodesy Joint Working Group on Theory of Earth Rotation, created in 2013.This work has been partially supported by the Spanish government under Grants AYA2010-22039-C02-01 and AYA2010-22039-C02-02 from Ministerio de Economía y Competitividad (MINECO), the University of Alicante under Grant GRE11-08 and the Generalitat Valenciana, Grant GV/2014/072
Disperse two-phase flows, with applications to geophysical problems

In this paper we study the motion of a fluid with several dispersed particles
whose concentration is very small (smaller than $10^{-3}$), with possible
applications to problems coming from geophysics, meteorology, and oceanography.
We consider a very dilute suspension of heavy particles in a
quasi-incompressible fluid (low Mach number). In our case the Stokes number is
small and --as pointed out in the theory of multiphase turbulence-- we can use
an Eulerian model instead of a Lagrangian one. The assumption of low
concentration allows us to disregard particle--particle interactions, but we
take into account the effect of particles on the fluid (two-way coupling). In
this way we can study the physical effect of particle inertia (and not only
passive tracers), with a model similar to the Boussinesq equations. The
resulting model is used in both direct numerical simulations and large eddy
simulations of a dam-break (lock-exchange) problem, which is a well-known
academic test case. Keywords: Dilute suspensions, Eulerian models, direct and
large eddy simulations, slightly compressible flows, dam-break (lock-exchange)
problem.Comment: 28 pages, 7 figure
Geosystemics: a systemic view of the Earth’s magnetic field and possibilities for an imminent geomagnetic transition

Geosytemics is a way to see and study the Earth in its wholeness, together with the eventual couplings among the subsystems composing our planet. This paper will provide this view for the
Earth’s magnetic field, reviewing most of the results obtained in our recent works. The main tools used by geosystemics are some nonlinear quantities, such as some kinds of entropy. Through them, it is possible to: a) establish the chaoticity and ergodicity of the recent geomagnetic field in a direct
and simple way; b) indentify the most extreme events in its history, as the most rapid and the slowest ones, i.e. jerks and polarity changes (reversals or excursions). In particular, regarding the latter phenomena, with the help of these entropic concepts and together with the use of the theory of
critical transitions, some clues can be given for a possible imminent change of the geomagnetic field dynamical regime
An Overview of Geodetic Volcano Research in the Canary Islands

The Canary Islands are mostly characterized by diffuse and scattered volcanism affecting a large area, with only one active stratovolcano, the Teide?Pico Viejo complex (Tenerife). More than 2 million people live and work in the 7,447 km2 of the archipelago, resulting in an average population density three times greater than the rest of Spain. This fact, together with the growth of exposure during the past 40 years, increases volcanic risk with respect previous eruptions, as witnessed during the recent 2011?2012 El Hierro submarine eruption. Therefore, in addition to purely scientific reasons there are economic and population-security reasons for developing and maintaining an efficient volcano monitoring system. In this scenario geodetic monitoring represents an important part of the monitoring system. We describe volcano geodetic monitoring research carried out in the Canary Islands and the results obtained. We consider for each epoch the two main existing constraints: the level of volcanic activity in the archipelago, and the limitations of the techniques available at the time. Theoretical and observational aspects are considered, as well as the implications for operational volcano surveillance. Current challenges of and future perspectives in geodetic volcano monitoring in the Canaries are also presented
Non-dipole and regional effects on the geomagnetic dipole moment estimation

The study of the temporal evolution of the dipole moment variations is a forefront research topic in Earth sciences. It constrains geodynamo simulations and is used to correct cosmogenic isotope production, which is evidence of past solar activity, and it is used to study possible correlations between the geomagnetic field and the climate. In this work, we have analysed the main error sources in the geomagnetic dipole moment computation from palaeomagnetic data: the influence of the non-dipole terms in the average approach, the inhomogeneous distribution of the current palaeomagnetic database, and the averaging procedure used to obtain the evolution of the dipole moment. To evaluate and quantify these effects, we have used synthetic data from a global model based on instrumental and satellite data, the International Geomagnetic Reference Field: 11th generation. Results indicate that the non-dipole terms contribute on a global scale of < 6 % in the averaged dipole moment, whereas the regional non-dipole contribution can show deviations of up to 35 % in some regions such as Oceania, and different temporal trends with respect to the global dipole moment evolution in other ones, such as Europe and Asia. A regional weighting scheme seems the best option to mitigate these effects in the dipole moment average approach. But when directional and intensity palaeomagnetic information is available on a global scale, and in spite of the inhomogeneity of the database, global modelling presents more reliable values of the geomagnetic dipole moment