Combination of multi-mission altimetry data along the Mekong River with spatio-temporal kriging

River water-level time series at fixed geographical locations, so-called virtual stations, have been computed from single altimeter crossings for many years. Their temporal resolution is limited by the repeat cycle of the individual altimetry missions. The combination of all altimetry measurements along a river enables computing a water-level time series with improved temporal and spatial resolutions. This study uses the geostatistical method of spatio-temporal ordinary kriging to link multi-mission altimetry data along the Mekong River. The required covariance models reflecting the water flow are estimated based on empirical covariance values between altimetry observations at various locations. In this study, two covariance models are developed and tested in the case of the Mekong River: a stationary and a non-stationary covariance model. The proposed approach predicts water-level time series at different locations along the Mekong River with a temporal resolution of 5 days. Validation is performed against in situ data from four gauging stations, yielding RMS differences between 0.82 and 1.29 m and squared correlation coefficients between 0.89 and 0.94. Both models produce comparable results when used for combining data from Envisat, Jason-1, and SARAL for the time period between 2002 and 2015. The quality of the predicted time series turns out to be robust against a possibly decreasing availability of altimetry mission data. This demonstrates that our method is able to close the data gap between the end of the Envisat and the launch of the SARAL mission with interpolated time series.     

A phase-field method for the direct simulation of two-phase flows in pore-scale media using a non-equilibrium wetting boundary condition

Advances in pore-scale imaging (e.g., μ-CT scanning), increasing availability of computational resources, and recent developments in numerical algorithms have started rendering direct pore-scale numerical simulations of multi-phase flow on pore structures feasible. Quasi-static methods, where the viscous and the capillary limit are iterated sequentially, fall short in rigorously capturing crucial flow phenomena at the pore scale. Direct simulation techniques are needed that account for the full coupling between capillary and viscous flow phenomena. Consequently, there is a strong demand for robust and effective numerical methods that can deliver high-accuracy, high-resolution solutions of pore-scale flow in a computationally efficient manner. Direct simulations of pore-scale flow on imaged volumes can yield important insights about physical phenomena taking place during multi-phase, multi-component displacements. Such simulations can be utilized for optimizing various enhanced oil recovery (EOR) schemes and permit the computation of effective properties for Darcy-scale multi-phase flows.We implement a phase-field model for the direct pore-scale simulation of incompressible flow of two immiscible fluids. The model naturally lends itself to the transport of fluids with large density and viscosity ratios. In the phase-field approach, the fluid-phase interfaces are expressed in terms of thin transition regions, the so-called diffuse interfaces, for increased computational efficiency. The conservation law of mass for binary mixtures leads to the advective Cahn–Hilliard equation and the condition that the velocity field is divergence free. Momentum balance, on the other hand, leads to the Navier–Stokes equations for Newtonian fluids modified for two-phase flow and coupled to the advective Cahn–Hilliard equation. Unlike the volume of fluid (VoF) and level-set methods, which rely on regularization techniques to describe the phase interfaces, the phase-field method facilitates a thermodynamic treatment of the phase interfaces, rendering it more physically consistent for the direct simulations of two-phase pore-scale flow. A novel geometric wetting (wall) boundary condition is implemented as part of the phase-field method for the simulation of two-fluid flows with moving contact lines. The geometric boundary condition accurately replicates the prescribed equilibrium contact angle and is extended to account for dynamic (non-equilibrium) effects. The coupled advective Cahn–Hilliard and modified Navier–Stokes (phase-field) system is solved by using a robust and accurate semi-implicit finite volume method. An extension of the momentum balance equations is also implemented for Herschel–Bulkley (non-Newtonian) fluids. Non-equilibrium-induced two-phase flow problems and dynamic two-phase flows in simple two-dimensional (2-D) and three-dimensional (3-D) geometries are investigated to validate the model and its numerical implementation. Quantitative comparisons are made for cases with analytical solutions. Two-phase flow in an idealized 2-D pore-scale conduit is simulated to demonstrate the viability of the proposed direct numerical simulation approach.     

Stratigraphy and Mesozoic–Cenozoic tectonic history of northern Sierra Los Ajos and adjacent areas,Sonora, Mexico

Geologic mapping in the northern Sierra Los Ajos reveals new stratigraphic and structural data relevant to deciphering the Mesozoic–Cenozoic tectonic evolution of the range. The northern Sierra Los Ajos is cored by Proterozoic, Cambrian, Devonian, Mississippian, and Pennsylvanian strata, equivalent respectively to the Pinal Schist, Bolsa Quartzite and Abrigo Limestone, Martin Formation, Escabrosa Limestone, and Horquilla Limestone. The Proterozoic–Paleozoic sequence is mantled by Upper Cretaceous rocks partly equivalent to the Fort Crittenden and Salero Formations in Arizona, and the Cabullona Group in Sonora, Mexico.Absence of the Upper Jurassic–Lower Cretaceous Bisbee Group below the Upper Cretaceous rocks and above the Proterozoic–Paleozoic rocks indicates that the Sierra Los Ajos was part of the Cananea high, a topographic highland during the Late Jurassic and Early Cretaceous. Deposition of Upper Cretaceous rocks directly on Paleozoic and Proterozoic rocks indicates that the Sierra Los Ajos area had subsided as part of the Laramide Cabullona basin during Late Cretaceous time. Basal beds of the Upper Cretaceous sequence are clast-supported conglomerate composed locally of basement (Paleozoic) clasts. The conglomerate represents erosion of Paleozoic basement in the Sierra Los Ajos area coincident with development of the Cabullona basin.The present-day Sierra Los Ajos reaches elevations of greater than 2600 m, and was uplifted during Tertiary basin-and-range extension. Upper Cretaceous rocks are exposed at higher elevations in the northern Sierra Los Ajos and represent an uplifted part of the inverted Cabullona basin. Tertiary uplift of the Sierra Los Ajos was largely accommodated by vertical movement along the north-to-northwest-striking Sierra Los Ajos fault zone flanking the west side of the range. This fault zone structurally controls the configuration of the headwaters of the San Pedro River basin, an important bi-national water resource in the US-Mexico border region.