Geophysical flows under location uncertainty,Part I Random transport and general models     

V. Resseguier  E. Mémin  B. Chapron 《地球物理与天体物理流体动力学》2017,111(3):149-176

A stochastic flow representation is considered with the Eulerian velocity decomposed between a smooth large scale component and a rough small-scale turbulent component. The latter is specified as a random field uncorrelated in time. Subsequently, the material derivative is modified and leads to a stochastic version of the material derivative to include a drift correction, an inhomogeneous and anisotropic diffusion, and a multiplicative noise. As derived, this stochastic transport exhibits a remarkable energy conservation property for any realizations. As demonstrated, this pivotal operator further provides elegant means to derive stochastic formulations of classical representations of geophysical flow dynamics.  相似文献   

Geophysical flows under location uncertainty,Part II Quasi-geostrophy and efficient ensemble spreading     

V. Resseguier  E. Mémin  B. Chapron 《地球物理与天体物理流体动力学》2017,111(3):177-208

Models under location uncertainty are derived assuming that a component of the velocity is uncorrelated in time. The material derivative is accordingly modified to include an advection correction, inhomogeneous and anisotropic diffusion terms and a multiplicative noise contribution. In this paper, simplified geophysical dynamics are derived from a Boussinesq model under location uncertainty. Invoking usual scaling approximations and a moderate influence of the subgrid terms, stochastic formulations are obtained for the stratified Quasi-Geostrophy and the Surface Quasi-Geostrophy models. Based on numerical simulations, benefits of the proposed stochastic formalism are demonstrated. A single realization of models under location uncertainty can restore small-scale structures. An ensemble of realizations further helps to assess model error prediction and outperforms perturbed deterministic models by one order of magnitude. Such a high uncertainty quantification skill is of primary interests for assimilation ensemble methods. MATLAB^® code examples are available online.  相似文献   

Three‐term inversion of prestack seismic data using a weighted l2, 1 mixed norm          下载免费PDF全文

Daniel O. Pérez  Danilo R. Velis  Mauricio D. Sacchi 《Geophysical Prospecting》2017,65(6):1477-1495

We present a new inversion method to estimate, from prestack seismic data, blocky P‐ and S‐wave velocity and density images and the associated sparse reflectivity levels. The method uses the three‐term Aki and Richards approximation to linearise the seismic inversion problem. To this end, we adopt a weighted mixed l_{2, 1}‐norm that promotes structured forms of sparsity, thus leading to blocky solutions in time. In addition, our algorithm incorporates a covariance or scale matrix to simultaneously constrain P‐ and S‐wave velocities and density. This a priori information is obtained by nearby well‐log data. We also include a term containing a low‐frequency background model. The l_{2, 1} mixed norm leads to a convex objective function that can be minimised using proximal algorithms. In particular, we use the fast iterative shrinkage‐thresholding algorithm. A key advantage of this algorithm is that it only requires matrix–vector multiplications and no direct matrix inversion. The latter makes our algorithm numerically stable, easy to apply, and economical in terms of computational cost. Tests on synthetic and field data show that the proposed method, contrarily to conventional l₂‐ or l₁‐norm regularised solutions, is able to provide consistent blocky and/or sparse estimators of P‐ and S‐wave velocities and density from a noisy and limited number of observations.  相似文献   

Estimating River Conductance from Prior Information to Improve Surface‐Subsurface Model Calibration          下载免费PDF全文

Yohann Cousquer  Alexandre Pryet  Nicolas Flipo  Célestine Delbart  Alain Dupuy 《Ground water》2017,55(3):408-418

Most groundwater models simulate stream‐aquifer interactions with a head‐dependent flux boundary condition based on a river conductance (CRIV). CRIV is usually calibrated with other parameters by history matching. However, the inverse problem of groundwater models is often ill‐posed and individual model parameters are likely to be poorly constrained. Ill‐posedness can be addressed by Tikhonov regularization with prior knowledge on parameter values. The difficulty with a lumped parameter like CRIV, which cannot be measured in the field, is to find suitable initial and regularization values. Several formulations have been proposed for the estimation of CRIV from physical parameters. However, these methods are either too simple to provide a reliable estimate of CRIV, or too complex to be easily implemented by groundwater modelers. This paper addresses the issue with a flexible and operational tool based on a 2D numerical model in a local vertical cross section, where the river conductance is computed from selected geometric and hydrodynamic parameters. Contrary to other approaches, the grid size of the regional model and the anisotropy of the aquifer hydraulic conductivity are also taken into account. A global sensitivity analysis indicates the strong sensitivity of CRIV to these parameters. This enhancement for the prior estimation of CRIV is a step forward for the calibration and uncertainty analysis of surface‐subsurface models. It is especially useful for modeling objectives that require CRIV to be well known such as conjunctive surface water‐groundwater use.  相似文献   

Simulating Stable Isotope Ratios in Plumes of Groundwater Pollutants with BIOSCREEN‐AT‐ISO          下载免费PDF全文

Patrick Höhener  Zhi M. Li  Maxime Julien  Pierrick Nun  Richard J. Robins  Gérald S. Remaud 《Ground water》2017,55(2):261-267

BIOSCREEN is a well‐known simple tool for evaluating the transport of dissolved contaminants in groundwater, ideal for rapid screening and teaching. This work extends the BIOSCREEN model for the calculation of stable isotope ratios in contaminants. A three‐dimensional exact solution of the reactive transport from a patch source, accounting for fractionation by first‐order decay and/or sorption, is used. The results match those from a previously published isotope model but are much simpler to obtain. Two different isotopes may be computed, and dual isotope plots can be viewed. The dual isotope assessment is a rapidly emerging new approach for identifying process mechanisms in aquifers. Furthermore, deviations of isotope ratios at specific reactive positions with respect to “bulk” ratios in the whole compound can be simulated. This model is named BIOSCREEN‐AT‐ISO and will be downloadable from the journal homepage.  相似文献   

Simulating Head Recovery and Intermittent Discharge in Free Flowing Wells          下载免费PDF全文

Ferenc Székely 《Ground water》2017,55(6):879-884

Modeling of the head recovery (buildup) in and intermittent operation of flowing wells lacks rigorous analytical tools. Presently different methods are applied in discharging and recovery phases. The presented reliable approximate solution considers a continuous free flow process and time variant screen loss parameters. The latter are temporarily set to large value to prevent inflow into the well bore during the single recovery or in repeated shut down periods. This method is called as screen loss control or SLC technique and is verified by means of comparative analysis with the constant rate pumping‐recovery test simulated in leaky aquifer. The evaluation of the free flow and recovery test in the free flowing well 28 (Artesia Heights) concluded with calibrated parameters close to those obtained in the earlier analyses based on split well flow process. The simulation of intermittent flowing well operation also underwent successful validation resulting in reasonable values of the Q_w and s response functions.  相似文献   

Occurrence and distribution of bacteria indicators,chemical tracers and pathogenic vibrios in Singapore coastal waters     

Shin Giek Goh  Stéphane Bayen  David Burger  Barry C. Kelly  Ping Han  Vladan Babovic  Karina Yew-Hoong Gin 《Marine pollution bulletin》2017,114(1):627-634

Water quality in Singapore's coastal area was evaluated with microbial indicators, pathogenic vibrios, chemical tracers and physico-chemical parameters. Sampling sites were grouped into two clusters (coastal sites at (i) northern and (ii) southern part of Singapore). The coastal sites located at northern part of Singapore along the Johor Straits exhibited greater pollution. Principal component analysis revealed that sampling sites at Johor Straits have greater loading on carbamazepine, while turbidity poses greater influence on sampling sites at Singapore Straits. Detection of pathogenic vibrios was also more prominent at Johor Straits than the Singapore Straits. This study examined the spatial variations in Singapore's coastal water quality and provided the baseline information for health risk assessment and future pollution management.  相似文献   

Paleomagnetism and magnetic fabrics of Mio-Pliocene hypabyssal rocks of the Combia event,Colombia: tectonic implications     

Victor A. Piedrahita  Roberto S. Molina-Garza  Gloria M. Sierra  José F. Duque-Trujillo 《Studia Geophysica et Geodaetica》2017,61(4):772-800

Mio-Pliocene hypabyssal rocks of the Combia event in the Amagá basin (NW Andes-Colombia), contain a deformational record of the activity of the Cauca-Romeral fault system, and the interaction of terranes within the Choco and northern Andean blocks. Previous paleomagnetic studies interpreted coherent counterclockwise rotations and noncoherent modes of rotation about horizontal axes for the Combia intrusives. However, rotations were determined from in-situ paleomagnetic directions and the existing data set is small. In order to better understand the deformational features of these rocks, we collected new paleomagnetic, structural, petrographic and magnetic fabric data from well exposed hypabyssal rocks of the Combia event. The magnetizations of these rocks are controlled by a low-coercivity ferromagnetic phase. Samples respond well to alternatingfield demagnetization isolating a magnetization component of moderate coercivity. These rocks do not have ductile deformation features. Anisotropy of magnetic susceptibility and morphotectonic analysis indicate that rotation about horizontal axes is consistently to the south-east, suggesting the need to apply a structural correction to the paleomagnetic data. The relationships between magnetic foliations and host-rock bedding planes indicate tectonic activity initiated before ~10 Ma. We present a mean paleomagnetic direction (declination D = 342.8°, inclination I = 12.1°, 95% confidence interval α₉₅ = 12.5°, precision parameter k = 8.6, number of specimens n = 18) that incorporates structural corrections. The dispersion S = 27° of site means cannot be explained by secular variation alone, but it indicates a counterclockwise rotation of 14.8° ± 12.7° relative to stable South America. Paleomagnetic data within a block bounded by the Sabanalarga and Cascajosa faults forms a more coherent data set (D = 336.5°, I = 17.4°, α₉₅ = 11.7°, k = 12.5, n = 14), which differs from sites west of the Sabanalarga fault and shows a rotation about a vertical axis of 20.2° ± 10.7°. Deformation in the Amagá basin may be tentatively explained by the obduction of the Cañas Gordas terrane over the northwestern margin of the northern Andean block. However, it can also be related to the local effects of the Cauca-Romeral fault system.  相似文献   

Breaking walls towards fully open source hydrological modeling     

Kazi Rahman  Nicolas Ray  Grégory Giuliani  Chetan Maringanti  Chris George  Anthony Lehmann 《Water Resources》2017,44(1):23-30

  相似文献