Fitting the Linear Model of Coregionalization by Generalized Least Squares   总被引：2，自引：0，他引：2  

Bernard Pelletier  Pierre Dutilleul  Guillaume Larocque  James W. Fyles 《Mathematical Geology》2004,36(3):323-343

In geostatistical studies, the fitting of the linear model of coregionalization (LMC) to direct and cross experimental semivariograms is usually performed with a weighted least-squares (WLS) procedure based on the number of pairs of observations at each lag. So far, no study has investigated the efficiency of other least-squares procedures, such as ordinary least squares (OLS), generalized least squares (GLS), and WLS with other weighing functions, in the context of the LMC. In this article, we compare the statistical properties of the sill estimators obtained with eight least-squares procedures for fitting the LMC: OLS, four WLS, and three GLS. The WLS procedures are based on approximations of the variance of semivariogram estimates at each distance lag. The GLS procedures use a variance–covariance matrix of semivariogram estimates that is (i) estimated using the fourth-order moments with sill estimates (GLS₁), (ii) calculated using the fourth-order moments with the theoretical sills (GLS₂), and (iii) based on an approximation using the correlation between semivariogram estimates in the case of spatial independence of the observations (GLS₃). The current algorithm for fitting the LMC by WLS while ensuring the positive semidefiniteness of sill matrix estimates is modified to include any least-squares procedure. A Monte Carlo study is performed for 16 scenarios corresponding to different combinations of the number of variables, number of spatial structures, values of ranges, and scale dependence of the correlations among variables. Simulation results show that the mean square error is accounted for mostly by the variance of the sill estimators instead of their squared bias. Overall, the estimated GLS₁ and theoretical GLS₂ are the most efficient, followed by the WLS procedure that is based on the number of pairs of observations and the average distance at each lag. On that basis, GLS₁ can be recommended for future studies using the LMC.  相似文献   

Fabric studies within the Cascade Lake shear zone, Sierra Nevada, California   总被引：1，自引：0，他引：1  

Basil Tikoff  Matthew R. Davis  Christian Teyssier  Michel de St. Blanquat  Guillaume Habert  Sven Morgan   《Tectonophysics》2005,400(1-4):209-226

The Cascade Lake shear zone occurs on the eastern margin of the Tuolumne Intrusive Suite, Sierra Nevada Batholith, California. Foliation in the zone is NNW trending and subvertical, and lineation is moderately south plunging. Deformation is syn-tectonic with emplacement of the Cathedral Peak granodiorite. A deformation gradient exists toward the NE margin of this pluton, with higher strains and lower temperatures of deformation found near the contact. We compare fabric data collected very densely in this shear zone using several techniques: field fabrics, 3D orientation of K-feldspar megacrysts, and AMS (anisotropy of magnetic susceptibility) analysis. In general, the results from the three different methods are in agreement. Deformation in this shear zone is part of a larger pattern of deformation within the Cathedral Peak granodiorite, as recorded by AMS analysis, and dextral shearing associated within the last stage of plutonism within the Sierra Nevada magmatic arc.  相似文献   

Backward modeling of the rifting kinematics in the Upper Rhine Graben: insights from an elastic-perfect contact law on the restoration of a multi-bloc domain     

Tristan G. M. Cornu  Guillaume Bertrand 《International Journal of Earth Sciences》2005,94(4):751-757

This paper addresses the problem of volume restoration for 3-D sedimentary basin kinematic deformation. The primary purpose is methodological and concerns the use of contact mechanics with the finite element method, in order to deform a geological multi-bloc domain. This approach is applied to backward model the later stage of rifting of a segment of the southern Upper Rhine Graben (France–Germany border). Preliminary results from our modeling demonstrate the ability of the method not only to handle complex geometries, but also to successfully perform retro-deformation of a complex geological domain. In addition, they provide or confirm crucial information on the rifting evolution and tectonic features of this segment of the Upper Rhine Graben, such as the distribution of deformation, the asymmetry of the graben and a significant left-lateral strike-slip component of displacement.  相似文献   

Native and artificial radiation-induced defects in montmorillonite. An EPR study     

S. Sorieul  T. Allard  G. Morin  B. Boizot  G. Calas 《Physics and Chemistry of Minerals》2005,32(1):1-7

A natural montmorillonite containing radiation-induced defects was studied with Electron Paramagnetic Resonance (EPR) spectroscopy (X- and Q-band). A first dominant native defect, namely native defect 1, is identified. It gives rise to an orthorhombic spectrum with g_x= 2.004 ± 0.005 g_y= 2.010 ± 0.003, g_z= 2.065 ± 0.002. Simulation of the EPR spectrum at X- and Q-band reveals a second native defect with isotropic spectrum at g= 2.019 ± 0.005. Both are electron holes trapped on oxygen atoms of the structure. The native defect 1 is located on an oxygen-silicon bond or a non-bonding orbital parallel to the c* axis. These defects are annealed at 500°C and the half-life determined for native defect 1 is circa 3,000 years. Irradiations with beta rays produced two additional hole centers of lower stability and distinct EPR parameters. Artificial irradiations show that montmorillonite can be used as a dosimeter in a large dose range.  相似文献   

Chronology and stratigraphy of the Magdalen Islands archipelago from the last glaciation to the early Holocene: new insights into the glacial and sea‐level history of eastern Canada          下载免费PDF全文

Audrey M. Rémillard  Guillaume St‐Onge  Pascal Bernatchez  Bernard Hétu  Jan‐Pieter Buylaert  Andrew S. Murray  Benoit Vigneault 《Boreas: An International Journal of Quaternary Research》2016,45(4):604-628

  相似文献   

Lateral extrusion along the Altyn Tagh Fault,Qilian Shan (NE Tibet): insight from a 3D crustal budget          下载免费PDF全文

Feng Cheng  Marc Jolivet  Guillaume Dupont‐Nivet  Lin Wang  Xiangjiang Yu  Zhaojie Guo 《地学学报》2015,27(6):416-425

The lithospheric strike‐slip Altyn Tagh Fault has accommodated hundreds of kilometres of displacement between the Qaidam and Tarim blocks since its Eocene reactivation. However, the way the deformation is accommodated in the Qilian Shan and further east remains uncertain. Based on 360 km of north‐eastward migration of the relatively rigid Qaidam block along the Altyn Tagh Fault and 3D isovolumetric balancing of the crustal deformation within the Altyn Tagh Fault–Qilian Shan system, we demonstrate that 250 ± 28 km (43.8–49.4%) of N20E directed crustal shortening and an additional ~250–370 km of eastward motion of the Qilian Shan crust must be accounted for by strike‐slip faulting in the Qilian Shan and crustal thickening in the Qinling area, as well as by extension in the adjoining North China block graben systems.  相似文献   

Modeling Channel Forms and Related Sedimentary Objects Using a Boundary Representation Based on Non-uniform Rational B-Splines     

Jeremy Ruiu  Guillaume Caumon  Sophie Viseur 《Mathematical Geosciences》2016,48(3):259-284

  相似文献   

Impacts of a large Sahelian city on groundwater hydrodynamics and quality: example of Niamey (Niger)     

Aïssata B. Hassane  Christian Leduc  Guillaume Favreau  Barbara A. Bekins  Thomas Margueron 《Hydrogeology Journal》2016,24(2):407-423

The management of groundwater resources is very important in the semiarid Sahel region, which is experiencing rapid urban development. Impacts of urbanization on groundwater resources were investigated in the unconfined aquifer of the Continental Terminal beneath the city of Niamey, Niger, using water level and chemical data. Hydrodynamic and chemical changes are best described by a combination of factors including the historical development of the city, current land use, water-table depth and topography. Seasonal groundwater recharge occurs with high spatial variability, as indicated by water-level monitoring in all wells, but there was no interannual trend over the 5-year study period. Groundwater salinity shows high spatial variability and a minor rising trend. The highest salinity is in the old city centre, with Na–NO₃ dominant, and it increases seasonally with recharge. Salinity is much lower and more variable in the suburbs (Ca–HCO₃, Ca–NO₃, and Na–NO₃ dominant). Nitrate is the main ionic contaminant and is seasonally or permanently above the international guidelines for drinking water quality in 36 % of sampled wells, with a peak value of 112 mg L^?1 NO₃–N (8 meq L^?1). Comparison of urban and rural sites indicates a long-term increase in groundwater recharge and nitrate enrichment in the urban area with serious implications for groundwater management in the region.  相似文献   

Indirect unstructured hex-dominant mesh generation using tetrahedra recombination     

Arnaud Botella  Bruno Lévy  Guillaume Caumon 《Computational Geosciences》2016,20(3):437-451

Corner-point gridding is widely used in reservoir and basin modeling but generally yields approximations in the representation of geological interfaces. This paper introduces an indirect method to generate a hex-dominant mesh conformal to 3D geological surfaces and well paths suitable for finite-element and control-volume finite-element simulations. By indirect, we mean that the method first generates an unstructured tetrahedral mesh whose tetrahedra are then merged into primitives (hexahedra, prisms, and pyramids). More specifically, we focus on determining the optimal set of primitives that can be recombined from a given tetrahedral mesh. First, we detect in the tetrahedral mesh all the feasible volumetric primitives using a pattern-matching algorithm (Meshkat and Talmor Int. J. Numer. Meth. Eng. 49(1-2), 17–30 2000) that we re-visit and extend with configurations that account for degenerated tetrahedra (slivers). Then, we observe that selecting the optimal set of primitives among the feasible ones can be formalized as a maximum weighted independent set problem (Bomze et al. 1999), known to be \(\mathcal {N}\mathcal {P}\)-Complete. We propose several heuristic optimizations to find a reasonable set of primitives in a practical time. All the tetrahedra of each selected primitive are then merged to build the final unstructured hex-dominant mesh. This method is demonstrated on 3D geological models including a faulted and folded model and a discrete fracture network.  相似文献   

Paleogene evolution and demise of the proto‐Paratethys Sea in Central Asia (Tarim and Tajik basins): Role of intensified tectonic activity at ca. 41 Ma     

Mustafa Y. Kaya  Guillaume Dupont‐Nivet  Jean‐Noël Proust  Pierrick Roperch  Laurie Bougeois  Niels Meijer  Joost Frieling  Chiara Fioroni  Sevin zkan Alt&#x;ner  Ezgi Vardar  Natasha Barbolini  Marius Stoica  Jovid Aminov  Mehmut Mamtimin  Guo Zhaojie 《Basin Research》2019,31(3):461-486

The proto‐Paratethys Sea covered a vast area extending from the Mediterranean Tethys to the Tarim Basin in western China during Cretaceous and early Paleogene. Climate modelling and proxy studies suggest that Asian aridification has been governed by westerly moisture modulated by fluctuations of the proto‐Paratethys Sea. Transgressive and regressive episodes of the proto‐Paratethys Sea have been previously recognized but their timing, extent and depositional environments remain poorly constrained. This hampers understanding of their driving mechanisms (tectonic and/or eustatic) and their contribution to Asian aridification. Here, we present a new chronostratigraphic framework based on biostratigraphy and magnetostratigraphy as well as a detailed palaeoenvironmental analysis for the Paleogene proto‐Paratethys Sea incursions in the Tajik and Tarim basins. This enables us to identify the major drivers of marine fluctuations and their potential consequences on Asian aridification. A major regional restriction event, marked by the exceptionally thick (≤ 400 m) shelf evaporites is assigned a Danian‐Selandian age (ca. 63–59 Ma) in the Aertashi Formation. This is followed by the largest recorded proto‐Paratethys Sea incursion with a transgression estimated as early Thanetian (ca. 59–57 Ma) and a regression within the Ypresian (ca. 53–52 Ma), both within the Qimugen Formation. The transgression of the next incursion in the Kalatar and Wulagen formations is now constrained as early Lutetian (ca. 47–46 Ma), whereas its regression in the Bashibulake Formation is constrained as late Lutetian (ca. 41 Ma) and is associated with a drastic increase in both tectonic subsidence and basin infilling. The age of the final and least pronounced sea incursion restricted to the westernmost margin of the Tarim Basin is assigned as Bartonian–Priabonian (ca. 39.7–36.7 Ma). We interpret the long‐term westward retreat of the proto‐Paratethys Sea starting at ca. 41 Ma to be associated with far‐field tectonic effects of the Indo‐Asia collision and Pamir/Tibetan plateau uplift. Short‐term eustatic sea level transgressions are superimposed on this long‐term regression and seem coeval with the transgression events in the other northern Peri‐Tethyan sedimentary provinces for the 1st and 2nd sea incursions. However, the 3rd sea incursion is interpreted as related to tectonism. The transgressive and regressive intervals of the proto‐Paratethys Sea correlate well with the reported humid and arid phases, respectively in the Qaidam and Xining basins, thus demonstrating the role of the proto‐Paratethys Sea as an important moisture source for the Asian interior and its regression as a contributor to Asian aridification.  相似文献