SiSeRHMap v1.0: A simulator for mapped seismic response using a hybrid model

SiSeRHMap is a computerized methodology capable of drawing up prediction maps of seismic response. It was realized on the basis of a hybrid model which combines different approaches and models in a new and non-conventional way. These approaches 5 and models are organized in a code-architecture composed of five interdependent modules. A GIS (Geographic Information System) Cubic Model (GCM), which is a layered computational structure based on the concept of lithodynamic units and zones, aims at reproducing a parameterized layered subsoil model. A metamodeling process confers a hybrid nature to the methodology. In this process, the one-dimensional linear 10 equivalent analysis produces acceleration response spectra of shear wave velocitythickness profiles, defined as trainers, which are randomly selected in each zone. Subsequently, a numerical adaptive simulation model (Spectra) is optimized on the above trainer acceleration response spectra by means of a dedicated Evolutionary Algorithm (EA) and the Levenberg–Marquardt Algorithm (LMA) as the final optimizer. In the fi15 nal step, the GCM Maps Executor module produces a serial map-set of a stratigraphic seismic response at different periods, grid-solving the calibrated Spectra model. In addition, the spectra topographic amplification is also computed by means of a numerical prediction model. This latter is built to match the results of the numerical simulations related to isolate reliefs using GIS topographic attributes. In this way, different sets 20 of seismic response maps are developed, on which, also maps of seismic design response spectra are defined by means of an enveloping technique

Adjoint models of mantle convection with seismic, plate motion, and stratigraphic constraints: North America since the Late Cretaceous

We apply adjoint models of mantle convection to North America since the Late Cretaceous. The present-day mantle structure is constrained by seismic tomography and the time-dependent evolution by plate motions and stratigraphic data (paleoshorelines, borehole tectonic subsidence, and sediment isopachs). We infer values of average upper and lower mantle viscosities, provide a synthesis of North American vertical motions (relative sea level) from the Late Cretaceous to the present, and reconstruct the geometry of the Farallon slab back to the Late Cretaceous. In order to fit Late Cretaceous marine inundation and borehole subsidence, the adjoint model requires a viscosity ratio across 660 km discontinuity of 15:1 (reference viscosity of 10^(21) Pa s), which is consistent with values previously inferred by postglacial rebound studies. The dynamic topography associated with subduction of the Farallon slab is localized in western North America over Late Cretaceous, representing the primary factor controlling the widespread flooding. The east coast of the United States is not stable; rather, it has been experiencing continuous dynamic subsidence over the Cenozoic, coincident with an overall eustatic fall, explaining a discrepancy between sea level derived from the New Jersey coastal plain and global curves. The east coast subsidence further constrains the mantle viscosity structure and requires an uppermost mantle viscosity of 10^(20) Pa s. Imposed constraints require that the Farallon slab was flat lying during Late Cretaceous, with an extensive zone of shallow dipping Farallon subduction extending beyond the flat-lying slab farther east and north by up to 1000 km than previously suggested

Simultaneous inversion of mantle properties and initial conditions using an adjoint of mantle convection

Through the assimilation of present-day mantle seismic structure, adjoint methods can be used to constrain the structure of the mantle at earlier times, i.e., mantle initial conditions. However, the application to geophysical problems is restricted through both the high computational expense from repeated iteration between forward and adjoint models and the need to know mantle properties (such as viscosity and the absolute magnitude of temperature or density) a priori. We propose that an optimal first guess to the initial condition can be obtained through a simple backward integration (SBI) of the governing equations, thus lessening the computational expense. Given a model with known mantle properties, we show that a solution based on an SBI-generated first guess has smaller residuals than arbitrary guesses. Mantle viscosity and the effective Rayleigh number are crucial for mantle convection models, neither of which is exactly known. We place additional constraints on these basic mantle properties when the convection-induced dynamic topography on Earth's surface is considered within an adjoint inverse method. Besides assimilating present-day seismic structure as a constraint, we use dynamic topography and its rate of change in an inverse method that allows simultaneous inversion of the absolute upper and lower mantle viscosities, scaling between seismic velocity and thermal anomalies, and initial condition. The theory is derived from the governing equations of mantle convection and validated by synthetic experiments for both one-layer viscosity and two-layer viscosity regionally bounded spherical shells. For the one-layer model, at any instant of time, the magnitude of dynamic topography is controlled by the temperature scaling while the rate of change of topography is controlled by the absolute value of viscosity. For the two-layer case, the rate of change of topography constrains upper mantle viscosity while the magnitude of dynamic topography determines the temperature scaling (lower mantle viscosity) when upper-mantle (lower-mantle) density anomaly dominates the flow field; this two-stage scheme minimizes the tradeoff between temperature and lower mantle viscosity. For both cases, we show that the theory can constrain mantle properties with errors arising through the adjoint recovery of the initial condition; for the two-layer model, this error is manifest as a tradeoff between the temperature scaling and lower mantle viscosity

Stretching and contraction of extensional basins with pre-rift salt: a numerical modelling approach

We present a series of 2D thermo-mechanical numerical experiments of thick-skinned crustal extension including a pre-rift salt horizon and subsequent thin-, thick-skinned, or mixed styles of convergence accompanied by surface processes. Extension localization along steep basement faults produces half-graben structures and leads to variations in the original distribution of pre-rift salt. Thick-skinned extension rate and salt rheology control hanging wall accommodation space as well as the locus and timing of minibasin grounding. Upon shortening, extension-related basement steps hinder forward propagation of evolving shallow thrust systems; conversely, if full basin inversion takes place along every individual fault, the regional salt layer is placed back to its pre-extensional configuration, constituting a regionally continuous décollement. Continued shortening and basement involvement deform the shallow fold-thrust structures and locally breaches the shallow décollement. We aim at obtaining a series of structural, stratigraphic and kinematic templates of fold-and-thrust belts involving rift basins with an intervening pre-rift salt horizon. Numerical results are compared to natural cases of salt-related inversion tectonics to better understand their structural evolution

Constrained Parameterization of the Multichannel Analysis of Surface Waves Approach with Application at Yuma Proving Ground, Arizona

Field data from Yuma Proving Ground, Arizona was used to test the feasibility of merging common multichannel analysis of surface waves (MASW) processing routines with mode- consistent shear-wave refraction traveltime tomography and synthetic modeling to optimize and constrain inversion results. Shear-wave first-arrival refraction tomography was used to enhance layer-model resolution and refine the MASW layer model with independent body-wave information. Shear-wave tomograms suggested a high-velocity layer, not found in initial `smooth' MASW velocity sections that were used as initial models for tomographic inversion. Increasing the stratification of the MASW layer model, to generally match tomogram structure, resulted in a higher-resolution MASW model constrained through joint analysis. This mutual analysis of shear-wave velocity (Vs) provided multiplicity to the structural interpretation of the site. Constrained-parameterization MASW results, compressional-wave tomography (Vp:Vs ratio), and density well logs populated a 2D model for numerical modeling, which was manually updated over several iterations to converge upon the site's first-arrival and dispersion characteristics. Further evaluation of the synthetic seismograms gave insight into the relationship between acquisition geometry (offset selection) and the associated dispersion-image character. Furthermore, modeling gave a secondary measurement on depth to half-space, velocity structure, and relative Vp:Vs ratios, which formulated a final MASW profile. The gradual change of the earth model, given an evolving hierarchy of constraint, is seen as the main finding of this thesis. The calculated movement towards a higher-resolution inversion based on joint geophysical measurements, analysis, and interpretation, engenders a constrained-parameterization solution with highest confidence

Computational modelling of iron-ore mineralisation with stratigraphic permeability anisotropy

This study develops a computational framework to model fluid transport in sedimentary basins, targeting iron ore deposit formation. It offers a simplified flow model, accounting for geological features and permeability anisotropy as driving factors. A new finite element method lessens computational effort, facilitating robust predictions and cost-effective exploration. This methodology, applicable to other mineral commodities, enhances understanding of genetic models, supporting the search for new mineral deposits amid the global energy transition