Earth Model Research Articles

Mechanical properties of hcp Fe at high pressures and temperatures from large-scale molecular dynamics simulations

We study the elastic behavior of hexagonal close-packed (hcp) Fe at the high temperature and pressure conditions of the Earth Core, using an embedded-atom method interatomic potential adjusted to those conditions. We calculate diffusivity, elastic constants, density, bulk modulus, shear modulus, and sound velocities vs temperature. We obtain reasonable agreement with ab initio simulations and with other empirical potential simulations. Our densities and shear modulus are slightly higher than those in the preliminary reference earth model for the core. Phase stability is discussed in terms of the Born criteria and free energies, finding that hcp is mechanically stable and that the free energy difference between hcp and body-centered cubic (bcc) is very small compared to the thermal energy. We compare our simulated shear modulus G to several analytical models, obtaining excellent agreement with the Atom in Jelium model by Swift and co-workers. Assuming that the yield strength Y is equal to the shear modulus G, Y=G/30, we find reasonable agreement with a recent parametrization of the Steinberg–Guinan model. These results can lead to future large-scale, multi-million simulations of Fe under core conditions for samples with microstructure like grain boundaries and twins, which might be present under those conditions.

Time-lapse Full Waveform Inversion with a new gradient-weighting strategy applied to marine seismic data

Time-lapse full waveform inversion (FWI) has the potential to reveal high-resolution reservoir production-related velocity changes. However, time-lapse FWI can introduce inversion artifacts, which might mask the true time-lapse signature within the reservoir zone, especially for the inversion of field data examples. We implement sequential bootstrap time-lapse FWI on marine seismic data from North Australia and observe noticeable inversion artifacts. To better understand these time-lapse FWI artifacts, we implement time-lapse FWI using the seismic data with different frequency bands and observe the distribution of the potential true time-lapse signature and inversion artifacts.On the inverted δVp (P wave velocity change) models using data with different frequency bands, the distribution of inversion artifacts exhibits inconsistencies, while the inverted δVp at the reservoir remains consistent. In order to comprehend these observations, we analyze the origin of these time-lapse artifacts from a theoretical perspective and find that a portion of these artifacts may arise from inversion null space and differences in data residuals between baseline and monitoring inversions. Building on these insights, we develop a novel time-lapse FWI method to suppress these inversion artifacts. We use the energy of the inverted δVp utilizing data with a lower dominant frequency band as a gradient-weighting term for time-lapse FWI using data with a higher dominant frequency band. The test of the novel time-lapse FWI method on marine data demonstrates its capability to suppress inversion artifacts effectively. The inverted δVp primarily reflects the 4D response, which is interpreted to be the softening of the reservoir caused by gas coming out of the solution within the oil column due to a pressure drop within the reservoir. To validate this data-driven research solution, we conduct tests on the SEG Advanced Modeling 4D (SEAM4D) model and synthetic 4D data where we know the true underlying earth models and time-lapse reservoir changes.

Consequences of Ignoring the Curvature of the Earth in Nineteenth-Century Large Surveys: A Case Study of Geometrical Geodesy and the Survey of the Texas and Pacific Railway Company Eighty-Mile Reserve

Consequences of Ignoring the Curvature of the Earth in Nineteenth-Century Large Surveys: A Case Study of Geometrical Geodesy and the Survey of the Texas and Pacific Railway Company Eighty-Mile Reserve

Erratum for “Consequences of Ignoring the Curvature of the Earth in Nineteenth-Century Large Surveys: A Case Study of Geometrical Geodesy and the Survey of the Texas and Pacific Railway Company Eighty-Mile Reserve”

Erratum for “Consequences of Ignoring the Curvature of the Earth in Nineteenth-Century Large Surveys: A Case Study of Geometrical Geodesy and the Survey of the Texas and Pacific Railway Company Eighty-Mile Reserve”

Detection of Earth’s free oscillation and analysis of the non-synchronous oscillation phenomenon of normal modes

Earth’s free oscillation can provide essential constraints for refining Earth models, inverting seismic source mechanisms, and studying the deep internal structure of the Earth. Large earthquakes can simultaneously excite numerous normal modes. Due to the Earth’s ellipticity, rotation, and internal heterogeneities, these normal modes undergo splitting, with the frequencies of singlets of normal modes becoming very close (only a few µHz apart). This imposes greater demands on the detection of normal modes. This paper introduces a novel method for normal mode detection based on the normal time–frequency transform (NTFT). Compared to classical FT spectrum methods and recent optimal sequence estimation (OSE), the proposed method not only detects more weak normal modes but also reveals the spatial distribution of the phase of each normal mode. Taking the detection of 0S2 as an example, the phase measurements of each singlet are spatially inconsistent. This phenomenon can provide prior information for other methods, such as product spectrum analysis (PSA), spherical harmonic stacking (SHS), multistation experiments (MSE), and OSE. Additionally, understanding the phase distribution patterns contributes to further study of geological structures, offering crucial foundational data and observational support.

Chapter 5. Gravity surveying and the ‘Figure of the Earth’ from Newton to CHAMP

Chapter 5. Gravity surveying and the ‘Figure of the Earth’ from Newton to CHAMP

Georecon: a coarse-to-fine visual 3D reconstruction approach for high-resolution images with neural matching priors

ABSTRACT Visual 3D reconstruction enables rebuilding 3D scenes from captured images, serving as a fundamental data source for digital earth modeling and intelligent cities. In the foundational step, recent methods leverage learning-based descriptors for image registration and achieve tremendous advances in precision and robustness. However, these methods inevitably execute down sampling towards high-resolution images to fit the needs of neural networks, which leads to precision degradation of feature localization and matching. Thus, we propose GeoRecon: a novel coarse-to-fine visual 3D reconstruction method that optimally utilizes high-resolution images for high-quality visual 3D reconstruction. Firstly, the coarse stage conducts coarse reconstruction from downsampled images by performing neural matching with geometric priors. Secondly, we define the fine-grained stage, proposing a GPU-based algorithm for generating image-patch correspondences based on the neural matching priors to perform fine-grained image registration. Finally, based on the optimized camera poses under this coarse-to-fine paradigm, progressive dense reconstruction leveraging efficient neural radiance fields is proposed to accomplish the high-quality MVS reconstruction. Comparative experiments across various scenarios demonstrate the proposed method’s superior precision, robustness, and reconstruction quality.

Structure and Non‐Ideal Mixing of Fe‐Ni‐S Liquid at High Temperature and Pressure and Its Implication for the Earth's Outer Core Composition

AbstractThe effect of light elements (LEs) such as sulfur on the physical properties of liquid iron‐nickel alloy under the earth's outer core conditions is critical for understanding the core composition and dynamics. First‐principles molecular dynamics simulations were employed to model Fe‐Ni‐S liquid with S concentrations in the range of (0–25) atomic percent (at%) at 4050 K and (0–33.33) at% at 5530 K and pressures relevant to the core‐mantle boundary (CMB) and inner core boundary (ICB), respectively. The thermodynamic mixing properties of Fe‐Ni‐S liquid were calculated, showing that the excess volume for Fe‐Ni‐S alloys deviates negatively from ideal mixing by −0.33% at 12.5 at% S at the CMB and −0.35% at 17 at% S at the ICB. Similarly, the excess enthalpy negatively deviated from the ideal mixing by −3.4 kJ/mole and −13 kJ/mole at the similar S concentrations at CMB and ICB, respectively, indicating non‐ideal mixing throughout the outer core. Similar behaviors are observed for isothermal bulk modulus (KT) and seismic velocity. The short‐ and intermediate‐range structures were analyzed and used to explain the non‐ideal mixing behaviors. The results suggest that extrapolations using ideal mixing underestimates the sound velocity by ∼0.14 km/s near CMB and ∼0.10 km/s near ICB, which is significant for constraining the core composition. If S is the only LE, the density at 10–12 wt% S matches the preliminary reference earth model (PREM). The seismic velocity at 12–15 wt% S matches PREM. These results suggest the presence of other LEs in the outer core.

Simulation-Based Optimization Workflow of CO2-EOR for Hydraulic Fractured Wells in Wolfcamp A Formation

Hydraulic fracturing has enabled production from unconventional reservoirs in the U.S., but production rates often decline sharply, limiting recovery factors to under 10%. This study proposes an optimization workflow for the CO2 huff-n-puff process for multistage-fractured horizontal wells in the Wolfcamp A formation in the Delaware Basin. The potential for enhanced oil recovery and CO2 sequestration simultaneously was addressed using a coupled geomechanics–reservoir simulation. Geomechanical properties were derived from a 1D mechanical earth model and integrated into reservoir simulation to replicate hydraulic fracture geometries. The fracture model was validated using a robust production history matching. A fluid phase behavior analysis refined the equation of state, and 1D slim tube simulations determined a minimum miscibility pressure of 4300 psi for CO2 injection. After the primary production phase, various CO2 injection rates were tested from 1 to 25 MMSCFD/well, resulting in incremental oil recovery ranging from 6.3% to 69.3%. Different injection, soaking and production cycles were analyzed to determine the ideal operating condition. The optimal scenario improved cumulative oil recovery by 68.8% while keeping the highest CO2 storage efficiency. The simulation approach proposed by this study provides a comprehensive and systematic workflow for evaluating and optimizing CO2 huff-n-puff in hydraulically fractured wells, enhancing the recovery factor of unconventional reservoirs.

Application of first- and second-order adjoint methods to glacial isostatic adjustment incorporating rotational feedbacks

SUMMARY This paper revisits and extends the adjoint theory for glacial isostatic adjustment (GIA) of Crawford et al. (2018). Rotational feedbacks are now incorporated, and the application of the second-order adjoint method is described for the first time. The first-order adjoint method provides an efficient means for computing sensitivity kernels for a chosen objective functional, while the second-order adjoint method provides second-derivative information in the form of Hessian kernels. These latter kernels are required by efficient Newton-type optimization schemes and within methods for quantifying uncertainty for non-linear inverse problems. Most importantly, the entire theory has been reformulated so as to simplify its implementation by others within the GIA community. In particular, the rate-formulation for the GIA forward problem introduced by Crawford et al. (2018) has been replaced with the conventional equations for modelling GIA in laterally heterogeneous earth models. The implementation of the first- and second-order adjoint problems should be relatively easy within both existing and new GIA codes, with only the inclusions of more general force terms being required.

Analytical solution for the skip trajectory of hypersonic vehicles with time-varying aerodynamic coefficients

Abstract An analytical method is proposed for calculating the altitude and flight path angle of hypersonic vehicles’ skip trajectories under the influence of arbitrary time-varying aerodynamic coefficients. First, a longitudinal plane dynamics equation is established based on the inertial Earth model and the specific kinetic energy is introduced as independent variable to simplify the derivation. Subsequently, the solution is obtained through the perturbation method, where the zero-order analytical solution is based on existing literature, and the first-order linear time-varying system is decoupled through matrix linear transformation. Finally, by comparing numerical integration and analytical solutions with constant aerodynamic coefficients, the precision of this method in dealing with time-varying aerodynamic coefficients is verified, demonstrating its advantages over analytical solutions with constant aerodynamic coefficients.

Prediction of shear wave velocity in three sedimentary rocks in East Baghdad oilfield using multiple regression analysis

Shear wave is a crucial parameter for assessing the wellbore stability, the stress response, and rock deformation. It is essential for constructing the mechanical earth model (MEM) for many applications related to reservoir geomechanics including wellbore stability, sand production, hydraulic fracturing, and fault reactivation. However, shear sonic data is often omitted during the well-logging measurements for cost and saving purposes. To overcome this challenge, recent research has been focused on determining shear wave velocity through the use of core plugs, empirical correlations, artificial intelligence techniques, and multiple regression to quantify and evaluate the mechanical properties of subsurface formations without performing direct measurements at the wellbore. The greatest difference between this study and the literature is to predict the shear wave velocities for three sedimentary rocks based on conventional well logs. This study has been conducted on datasets of two wells drilled in the East Baghdad oilfield, for which there is a lack of shear wave data. Two formations (Tanuma and Zubair formations) within the production section of this field were conducted to develop new models for determining the shear wave velocity using multiple regression analysis. These two formations primarily consist of three lithologies: limestone, sandstone, and shale. Before the model development, data analysis on the selected data was applied to figure the most influential parameter(s) in determining the shear wave velocity. The results of the developed models are then compared with the previous models in the literature. The results showed that the multiple regression analysis technique is a powerful technique in determining shear wave velocity with high-performance capacity. The correlation coefficient ( ) and the root mean square error (RMSE) were 0.84 and 0.092 for limestone, 0.84 and 0.0972 for sandstone, and 0.86 and 0.0796 for shale respectively. Furthermore, the performance of the developed models is well matched to the actual shear wave data rather than the Castagna correlations. The findings of this study are effective in determining shear wave velocity for future applications related to reservoir geomechanics without needing costly well-log or core measurements.

Neutrino tomography of the earth: the earth total mass, moment of inertia and hydrostatic equilibrium constraints

We investigate the implications of the constraints following from the precise knowledge of the total Earth mass, M⊕\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$M_\\oplus $$\\end{document}, and moment of inertia, I⊕\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$I_\\oplus $$\\end{document}, and from the requirement that Earth be in hydrostatic equilibrium (EHE), in the neutrino tomography studies of the Earth density structure. In order to estimate the sensitivity of a given neutrino detector to possible deviations of the inner core (IC), outer core (OC), core (IC + OC) and mantle Earth densities from those obtained using geophysical and seismological data and described by the preliminary reference Earth model (PREM), in the statistical analyses performed within the neutrino tomography studies one typically varies the density of each of these structures. These variations, however, must respect the M⊕\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$M_\\oplus $$\\end{document}, I⊕\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$I_\\oplus $$\\end{document} and EHE constraints. Working with PREM average densities we derive the M⊕\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$M_\\oplus $$\\end{document}, I⊕\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$I_\\oplus $$\\end{document} and EHE constraints on the possible density variations when one approximates the Earth density structure with (i) three layers – mantle, outer core and inner core, and (ii) four layers – upper mantle, lower mantle, outer core and inner core. We get drastically different results in the two cases.

Sand Production Prediction and Management Using a One-Dimensional Geomechanical Model: A Case Study in NahrUmr Formation, Subba Oil Field

Sand production in sandstone reservoirs is a complex problem for oil and gas companies. This problem can damage downhole equipment and surface production facilities. This paper presents a sand production case and quantifies sand risks for the NahrUmr formation located in Subbaoilfied.The risk of sanding or rock failure is commonly observed during production phases as well as drilling operations. A (1D) mechanical earth model was created using wireline log data and core data for estimating the formations mechanical properties.The extrapolation method is utilized for vertical stress prediction, whereas pore pressure and static Young’s modulus were calculated using Slowness Eaton method and John Fuller correlation respectively. The poro-elastic model was used to determine horizontal principal stresses and hydrocarbon intervals susceptible to sand production issues.The purpose of this study was achieved using Tech-Log software. Sand onset intervals are firstly predicted as an open hole using ten methods: Loading factor (LF), Critical drawdown pressure (CDDP), B-index, schlumberger (Sc), Elastic combination index (Ec-index), Interval Transit-Time (DTc), Total porosity method (PHIT), unconfined compressive strength (UCS), Shear modulus to bulk compressibility ratio (G/Cb), and shear failure method. Approximately ten methods refer to the same intervals that can produce sand. Secondly, sand depths were predicted by representing the well according to what actually existed as a cased hole by calculating the critical drawdown pressure (CDDP), which is calculated at various depletion rates (0%, 15%, 25%, and 35%) in order to measure the effect of reservoir pressure depletion on sand production. The results show that as the rock strength is increased, the amount of sand-free drawdown and depletion becomes larger. Also, a single-depth analysis was conducted for problem intervals, indicatg a sandwillproducefromdepth of 2527.7 m under certain conditions. Therefore, it is necessary to install sand control equipment in this well.

Three-Dimensional Coupled Temporal Geomechanical Model for Fault-Reactivation and Surface-Deformation Evaluation during Reservoir Depletion and CO2 Sequestration, Securing Long-Term Reservoir Sustainability

Assessing reservoir subsidence due to depletion involves understanding the geological and geophysical processes that lead to ground subsidence as a result of reservoir fluid extraction. Subsidence is a gradual sinking or settling of the Earth’s surface, and it can occur when hydrocarbons are extracted from underground reservoirs. In this study, a time-integrated 3D coupled geomechanical modeling incorporating the fourth dimension—time—into traditional 3D geomechanical models has been constructed utilizing seismic inversion volumes and a one-dimensional mechanical Earth model (1D MEM). The 3D geomechanical model was calibrated to the 1D MEM results. Geomechanical rock properties were derived from the density and sonic log data that was distributed with conditioning to the seismic inversion volumes obtained from running pre-stack inversion. The standard elastic parameter equations were used to generate estimates of the elastic moduli. These properties are dynamic but have been converted to static values using additional equations used in the 1D MEM study. This included estimating the Unconfined Compressive Strength. In situ stresses were matched using different minimum horizontal principal stress gradients and horizontal principal stress ratios. The match is good except where the weak carbonate faults are close to the wells, where the Shmin magnitudes tend to decrease. The SHmax orientations were assessed from image log data and indicated to be 110° in the reservoir section. A time-integrated 3D coupled simulation was created using the finite-element method (FEM). The effective stresses increase while there is depletion in all directions, especially in the Z direction. The predicted compaction in the reservoir and overburden was 350 mm. Most of the compaction occurs at the reservoir level and dissipates towards the surface (seabed). Furthermore, the case displayed no shear failure that might cause or fault reactivation in the reservoir interval (Kangan–Dalan Formations) located in the simulated area. In this study, we applied an integrated and comprehensive geomechanical approach to evaluate subsidence, fault reactivation and stress alteration, while reservoir depletion was assessed using seismic inversion, well logs, and experiment data. The deformation monitoring of geological reservoirs, whether for gas storage or hazardous gas disposal, is essential due to the economic value of the stored assets and the hazardous nature of the disposed materials. This monitoring is vital for ensuring the sustainability of the reservoir by maintaining operational success and detecting integrity issues.

The origin and implications of primordial helium depletion in the Afar mantle plume

Mantle plumes are responsible for the Earth’s largest volcanic provinces. In the prevailing paradigm, the deep mantle is less degassed than convecting shallow mantle, implying that plume-derived lavas have higher concentrations of primordial volatiles such as helium (He). Demonstrating this has led to explanations that question the established Earth model. Here, we show that the 3He/4He of basalts from the Red Sea display coherent relationships with trace elements, allowing the helium concentration of the Afar plume to be calculated. Contrary to the prevailing model it appears the helium concentration of the Afar plume is 10-25% of the upper mantle. This contradiction is resolved if the plume material itself is a mixture of helium-rich high-3He/4He deep mantle with helium-depleted low-3He/4He recently subducted oceanic crust. This implies that helium-depleted domains may exist in convecting mantle and that moderately high 3He/4He plumes likely do not contain a notable contribution of the deep mantle.

Feedback mechanisms controlling Antarctic glacial-cycle dynamics simulated with a coupled ice sheet–solid Earth model

Abstract. The dynamics of the ice sheets on glacial timescales are highly controlled by interactions with the solid Earth, i.e., the glacial isostatic adjustment (GIA). Particularly at marine ice sheets, competing feedback mechanisms govern the migration of the ice sheet's grounding line (GL) and hence the ice sheet stability. For this study, we developed a coupling scheme and performed a suite of coupled ice sheet–solid Earth simulations over the last two glacial cycles. To represent ice sheet dynamics we apply the Parallel Ice Sheet Model (PISM), and to represent the solid Earth response we apply the 3D VIscoelastic Lithosphere and MAntle model (VILMA), which, in addition to load deformation and rotation changes, considers the gravitationally consistent redistribution of water (the sea-level equation). We decided on an offline coupling between the two model components. By convergence of trajectories of the Antarctic Ice Sheet deglaciation we determine optimal coupling time step and spatial resolution of the GIA model and compare patterns of inferred relative sea-level change since the Last Glacial Maximum with the results from previous studies. With our coupling setup we evaluate the relevance of feedback mechanisms for the glaciation and deglaciation phases in Antarctica considering different 3D Earth structures resulting in a range of load-response timescales. For rather long timescales, in a glacial climate associated with the far-field sea-level low stand, we find GL advance up to the edge of the continental shelf mainly in West Antarctica, dominated by a self-amplifying GIA feedback, which we call the “forebulge feedback”. For the much shorter timescale of deglaciation, dominated by the marine ice sheet instability, our simulations suggest that the stabilizing sea-level feedback can significantly slow down GL retreat in the Ross sector, which is dominated by a very weak Earth structure (i.e., low mantle viscosity and thin lithosphere). This delaying effect prevents a Holocene GL retreat beyond its present-day position, which is discussed in the scientific community and supported by observational evidence at the Siple Coast and by previous model simulations. The applied coupled framework, PISM–VILMA, allows for defining restart states to run multiple sensitivity simulations from. It can be easily implemented in Earth system models (ESMs) and provides the tools to gain a better understanding of ice sheet stability on glacial timescales as well as in a warmer future climate.

Kinematic vorticity of shear zones that accommodate vertical crustal advection: Implications for metamorphic core complexes and pluton emplacement

Kinematic analysis of ductile shear zones is an important method to interpret the dynamic evolution of many tectonic and magmatic processes on Earth, such as orogeny, rifting, and plutonism. Despite decades of study, kinematic vorticity analysis is an underutilized tool to interpret the dynamic drivers of shear zones. Determination of the kinematic vorticity number (Wk) quantifies the relative contributions of pure- versus simple-shear strain in shear zones. Here we systematically investigated Wk from mylonitic shear zones associated with metamorphic core complexes (MCCs) developed across the central and southern North American Cordillera using electron backscatter diffraction (EBSD) data that allows consistent and reproducible results from a large number of recrystallized grains. We investigated samples from four distinct MCC systems and the structural aureole of a Cordilleran pluton to investigate their kinematics. We find that most MCC samples display pure-to-general-shear strain (average 70 % pure shear), consistent with significant bulk shear-zone shortening (>80 % shortening) observed in many of the MCC systems. Such strain patterns are remarkably similar to deformation observed around plutons that were forcefully emplaced as diapirs. To further validate and investigate these results, we constructed the Wk field of these geologic processes using numerical simulations to highlight that pure shear kinematics are more common with diapirism and coupled wallrock deformation compared with discrete detachment-involved normal fault systems. These observations support that buoyant doming can be a viable end-member process to form the investigated MCCs. Our results also suggest that pure shear deformation may be a diagnostic strain characteristic for diapir-like crustal processes, and therefore Wk analyses could be used to test similar processes like dome-and-keel models for Early Earth.

Rejuvenating exploration opportunities in Egypt's Mediterranean region using basin-scale 3D seismic reprocessing and subsurface evaluations

The Eastern Mediterranean region has proven to be a key province for natural gas resources over the last two decades. Gas discoveries in the Nile Delta Basin, Levantine Basin, and over the Eratosthenes platform have proven that an efficient working petroleum system extends to the deepwater areas and into the distal parts of the Nile Delta and Levantine basins. While the suprasalt biogenic post-Messinian play is adequately imaged in various areas in the Eastern Mediterranean region, seismic imaging of the deeper subsalt pre-Messinian plays is challenged by legacy seismic acquisition and processing techniques. Moreover, the complexity of the Messinian salt structures and the heterogeneity of the lithology of the Messinian layer across the Mediterranean basins represent a key subsurface challenge. In phase 1 of the study presented in this paper, we reprocessed 12 legacy 3D seismic data sets over 12,000 km2 and constructed a consistent basin-scale earth model from the shelf areas offshore Sinai to the deepwater areas around the Zohr carbonate platform. The data processing and imaging sequences incorporate full-waveform inversion and tomography in the earth model building process followed by a basin-to-lead evaluation to delineate exploration opportunities in phase 1. The evaluations resulted in delineating Cretaceous and Oligo-Miocene untested structures. Furthermore, potential deepwater Lower Miocene sand fairways were outlined, supported by improved regional seismic imaging and consistent seismic attributes over multiple legacy surveys. The authors of this paper expect to update the results of their work after completing the earth model building and the prestack depth migration workflows of an additional 13 legacy surveys over 35,000 km2 in phase 2 of the study.

Estimating the Ebro river discharge at 1 km/daily resolution using indirect satellite observations

Estimating river discharge Q at global scale from satellite observations is not yet fully satisfactory in part because of limited space/time resolution. Furthermore, on highly anthropized basins, it is essential to anchor the analysis to reliable Q measurements. Gauge networks are however very sparse and limited in time, and SWOT (Surface Water Ocean Topography) river discharge estimates at global scale are not yet available. The method proposed here is able to obtain continuous daily Q estimates at 1 km/daily resolution, using indirect satellite data and ground-based estimates. We focus here on the Ebro. Over such an anthropized basin (e.g. change of land use, irrigation), the exploitation of 205 available gauges at their nominal resolution (i.e., daily point measurements) is a necessity. The hydrological Continuum model is used to help interpolate spatially and temporally the observations into our optimal interpolation scheme. The proposed Q-mapping is similar to an assimilation scheme were Earth observations (precipitation, evapotranspiration and total water storage change) and model simulations are constrained by in situ gauge measurements. The Q estimates are evaluated using a rigorous leave-one-out experiment, showing a good agreement with the in situ data: a correlation of 0.72 (median), and a 75th percentile of Nash-Sutcliffe Efficiency up to 0.62. Our spatio-temporal continuous Q estimates at high spatial/temporal resolution can describe complex continental water dynamics, including extreme events. SWOT estimates will soon be available, at the global scale but with irregular space/time sampling: our method should help exploit them to obtain a regular space-temporal description of the water cycle at high resolution.