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

  相似文献   

Geodynamic Evolution of Northeastern Tunisia During the Maastrichtian–Paleocene Time: Insights from Integrated Seismic Stratigraphic Analysis     

Oussama Abidi  Mohamed Hédi Inoubli  Kawthar Sebei  Adnen Amiri  Haifa Boussiga  Imen Hamdi Nasr  Abdelhamid Ben Salem  Mahmoud Elabed 《Surveys in Geophysics》2017,38(3):617-649

The Maastrichtian–Paleocene El Haria formation was studied and defined in Tunisia on the basis of outcrops and borehole data; few studies were interested in its three-dimensional extent. In this paper, the El Haria formation is reviewed in the context of a tectono-stratigraphic interval using an integrated seismic stratigraphic analysis based on borehole lithology logs, electrical well logging, well shots, vertical seismic profiles and post-stack surface data. Seismic analysis benefits from appropriate calibration with borehole data, conventional interpretation, velocity mapping, seismic attributes and post-stack model-based inversion. The applied methodology proved to be powerful for charactering the marly Maastrichtian–Paleocene interval of the El Haria formation. Migrated seismic sections together with borehole measurements are used to detail the three-dimensional changes in thickness, facies and depositional environment in the Cap Bon and Gulf of Hammamet regions during the Maastrichtian–Paleocene time. Furthermore, dating based on their microfossil content divulges local and multiple internal hiatuses within the El Haria formation which are related to the geodynamic evolution of the depositional floor since the Campanian stage. Interpreted seismic sections display concordance, unconformities, pinchouts, sedimentary gaps, incised valleys and syn-sedimentary normal faulting. Based on the seismic reflection geometry and terminations, seven sequences are delineated. These sequences are related to base-level changes as the combination of depositional floor paleo-topography, tectonic forces, subsidence and the developed accommodation space. These factors controlled the occurrence of the various parts of the Maastrichtian–Paleocene interval. Detailed examinations of these deposits together with the analysis of the structural deformation at different time periods allowed us to obtain a better understanding of the sediment architecture in depth and the delineation of the geodynamic evolution of the region.  相似文献   

A wind tunnel investigation of particle segregation,ripple formation and armouring within sand beds of systematically varied texture          下载免费PDF全文

Cheryl McKenna Neuman  ottO Bédard 《地球表面变化过程与地形》2017,42(5):749-762

This paper reports on a wind tunnel investigation of particle segregation, ripple formation and surface armouring within sand beds of systematically varied particle size distribution, from coarsely skewed to bimodal. By design, the system was closed with no external inputs of mass from an external particle feed. Particles too coarse to travel in saltation for the given range in wind speed were dyed red in order to distinguish them in optical images from finer sand particles, which could be entrained into the unidirectional airflow. A 3D laser scanner measured the changing bed topography at regular time intervals during 18 experiments involving varied combinations of wind speed and bed texture. Image classification techniques were used to investigate the coincident self‐organization of the two populations of particles, as distinguished by their colour. As soon as saltation commenced, some of the red particles segregated into thin discontinuous patches. Particle trapping and sheltering on these rough patches was strongly favoured, causing them to grow preferentially. During the earliest stages of formation, bedform growth coincided with: (i) rapid coarsening of the surface texture; and (ii) the merging of proto‐ripple ‘crests’ to generate larger rhythmic bedforms of lower frequency. Consistent with previous work, ripple size was observed to increase under stronger winds when not exceeding the threshold for entrainment of the coarse‐mode or red particles from the crest. With declining rates of mass transport and particle segregation as the bed surface armoured, and the consequent deceleration of ripple propagation through to the end of each experiment, all surfaces eventually attained a steady‐state morphometry. At saturation, the largest ripples developed on beds having the lowest initial concentration of red particles. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献