Microstructural Inhomogeneity and Mechanical Anisotropy Associated with Bedding in Rothbach Sandstone

This study present the result of conventional triaxial tests conducted on samples of Rothbach sandstone cored parallel, oblique (at 45 degrees) and perpendicular to the bedding at effective pressures ranging from 5 to 250 MPa. Mechanical and microstructural data were used to determine the role of the bedding on mechanical strength and failure mode. We find that samples cored at 45 degrees to the bedding yield at intermediate level of differential stress between the ones for parallel and perpendicular samples at all effective pressures. Strain localization at high confining pressure (i.e., in the compactive domain) is observed in samples perpendicular and oblique to the bedding but not in samples cored parallel to the bedding. However, porosity reduction is comparable whether compactive shear bands, compaction bands or homogeneous cataclastic flow develop. Microstructural data suggest that (1) mechanical anisotropy is controlled by a preferred intergranular contact alignment parallel to the bedding and that (2) localization of compaction is controlled by bedding laminations and grain scale heterogeneity, which both prevent the development of well localized compaction features.     

Buried paleo-channels on the New Jersey continental margin: channel porosity structures from electromagnetic surveying

We report on a marine electromagnetic (EM) survey across two portions of the New Jersey continental margin that have been previously shown to contain buried paleo-channels. The EM method used provides bulk porosity estimates to depths of around 20 m below the seafloor and is thus able to place porosity constraints on the nature of the channel infill and the contrast in physical properties across the channel boundaries. Our data show that a key condition for the channels to have an electrical signature is that they incise an underlying regional unconformity, R, thought to represent a subaerially eroded surface, exposed during the late Wisconsinan glaciation. Channels that cut R are seen through increases in apparent porosity. Another seismically imaged channel sequence, which lies within the outer-shelf sediment wedge sequence above R, does not have an electrical signature, indicating that the sediments above and below the channel boundaries have similar physical properties.     

Prehistoric soil and water detention structures (gabarbands) at Phang,Sindh Kohistan,Pakistan: An adaptation to environmental change?

Five gabarbands (dams), components of integrated soil and water conservation systems, were investigated in the arid region of Sindh Kohistan. The dams are associated with two prehistoric Kot Dijian settlements, Phang and Kohtrash (3200–2800 B. C.). The dams were constructed with relatively permeable materials. Three of the dams, South, East, and North, close off water gaps and a strike valley, respectively, to form a detention basin. The West Dam, located on Phang Nai, upstream of the detention basin removed the coarser sediment fraction from the flows. Spring Dam, located in a water gap to the north of the detention basin, acted as a low-head weir to distribute flow from a perennial spring onto the floodplain of Baran Nai, where double-cropping was probably practiced. Faulting eliminated the spring, which had probably supported the Kohtrash site, and led to construction of the detection basin for the purposes of subsurface storage of detained flows. Association of both modern and prehistoric sites with springs suggests that there has not been appreciable climate change in Sindh Kohistan in the last 5000 years. © 1993 John Wiley & Sons, Inc.     

Kinetic and equilibrium studies of copper-dissolved organic matter complexation in water column of the stratified Krka River estuary (Croatia)

An interaction of dissolved natural organic matter (DNOM) with copper ions in the water column of the stratified Krka River estuary (Croatia) was studied. The experimental methodology was based on the differential pulse anodic stripping voltammetric (DPASV) determination of labile copper species by titrating the sample using increments of copper additions uniformly distributed on the logarithmic scale. A classical at-equilibrium approach (determination of copper complexing capacity, CuCC) and a kinetic approach (tracing of equilibrium reconstitution) of copper complexation were considered and compared. A model of discrete distribution of organic ligands forming inert copper complexes was applied. For both approaches, a home-written fitting program was used for the determination of apparent stability constants (K_i^equ), total ligands concentration (L_iT) and association/dissociation rate constants (k_i¹,k_i^- 1).A non-conservative behaviour of dissolved organic matter (DOC) and total copper concentration in a water column was registered. An enhanced biological activity at the freshwater–seawater interface (FSI) triggered an increase of total copper concentration and total ligand concentration in this water layer. The copper complexation in fresh water of Krka River was characterised by one type of binding ligands, while in most of the estuarine and marine samples two classes of ligands were identified. The distribution of apparent stability constants (log K₁^equ: 11.2–13.0, log K₂^equ:8.8–10.0) showed increasing trend towards higher salinities, indicating stronger copper complexation by autochthonous seawater organic matter.Copper complexation parameters (ligand concentrations and apparent stability constants) obtained by at-equilibrium model are in very good accordance with those of kinetic model. Calculated association rate constants (k₁¹:6.1–20 × 10³ (M s)^− 1, k₂¹: 1.3–6.3 × 10³ (M s)^− 1) indicate that copper complexation by DNOM takes place relatively slowly. The time needed to achieve a new pseudo-equilibrium induced by an increase of copper concentration (which is common for Krka River estuary during summer period due to the nautical traffic), is estimated to be from 2 to 4 h.It is found that in such oligotrophic environment (dissolved organic carbon content under 83 µM_C, i.e. 1 mg_CL^− 1) an increase of the total copper concentration above 12 nM could enhance a free copper concentration exceeding the level considered as potentially toxic for microorganisms (10 pM).     

A study on mechanical properties and fracturing behavior of Carrara marble with the flat-jointed model

Overcoming two inherent limitations presented by the standard bonded-particle model (BPM), the flat-jointed model still lacks extensive application. This study hereby conducts a comprehensive investigation on the application of the flat-jointed model in rock mechanics study. First, after a careful examination on the mechanical behavior of the flat-joint contact, a systematic calibration has been conducted to build numerical models matching the mechanical properties of Carrara marble. Second, the validated model is used to simulate the fracturing behavior of marble specimens containing a single or en-echelon flaws. Finally, the appearance of particle flow code (PFC) models has been compared with the results obtained from optical observation in order to establish a possible correlation between PFC models and real marble regarding microcracking behavior. Discussions suggest that to remove the particle size effect, the model resolution (ratio of sample dimension to the average particle diameter) should be larger than 30 in the flat-jointed model and 150 in the BPM, respectively. Results from the above investigations suggest that the flat-jointed model is capable of matching both the mechanical properties and fracturing behavior of Carrara marble. The numerical model is more reliable to predict the fracturing path at the instant of initiation than at failure. Development of initial fractures is not only controlled by the magnitude of tension force but also influenced by the extent of distribution of tension force. The appearance and components of a fracture in flat-jointed models probably reveal not only its property (tensile or shear) but also the corresponding microcracking behavior of real rocks.     

Kernel Principal Component Analysis for Efficient,Differentiable Parameterization of Multipoint Geostatistics

This paper describes a novel approach for creating an efficient, general, and differentiable parameterization of large-scale non-Gaussian, non-stationary random fields (represented by multipoint geostatistics) that is capable of reproducing complex geological structures such as channels. Such parameterizations are appropriate for use with gradient-based algorithms applied to, for example, history-matching or uncertainty propagation. It is known that the standard Karhunen–Loeve (K–L) expansion, also called linear principal component analysis or PCA, can be used as a differentiable parameterization of input random fields defining the geological model. The standard K–L model is, however, limited in two respects. It requires an eigen-decomposition of the covariance matrix of the random field, which is prohibitively expensive for large models. In addition, it preserves only the two-point statistics of a random field, which is insufficient for reproducing complex structures. In this work, kernel PCA is applied to address the limitations associated with the standard K–L expansion. Although widely used in machine learning applications, it does not appear to have found any application for geological model parameterization. With kernel PCA, an eigen-decomposition of a small matrix called the kernel matrix is performed instead of the full covariance matrix. The method is much more efficient than the standard K–L procedure. Through use of higher order polynomial kernels, which implicitly define a high-dimensionality feature space, kernel PCA further enables the preservation of high-order statistics of the random field, instead of just two-point statistics as in the K–L method. The kernel PCA eigen-decomposition proceeds using a set of realizations created by geostatistical simulation (honoring two-point or multipoint statistics) rather than the analytical covariance function. We demonstrate that kernel PCA is capable of generating differentiable parameterizations that reproduce the essential features of complex geological structures represented by multipoint geostatistics. The kernel PCA representation is then applied to history match a water flooding problem. This example demonstrates that kernel PCA can be used with gradient-based history matching to provide models that match production history while maintaining multipoint geostatistics consistent with the underlying training image.     

gLucifer: next generation visualization framework for high-performance computational geodynamics

High-performance computing provides unprecedented capabilities to produce higher resolution 4-D models in a fraction of time. Thus, the need exists for a new generation of visualization systems able to maintain parity with the enormous volume of data generated. In attempting to write this much data to disk, each computational step introduces a significant performance bottleneck, yet most existing visualization software packages inherently rely on reading data in from a dump file. Available packages make this assumption of postprocessing at quite a fundamental level and are not very well suited for plotting very large numbers of specialized particles. This necessitates the creation of a new visualization system that meets the needs of large-scale geodynamic modeling. We have developed such a system, gLucifer, using a software framework approach that allows efficient reuse of our efforts in other areas of research. gLucifer is capable of producing movies of a 4-D data set “on the fly” (simultaneously with running the parallel scientific application) without creating a performance bottleneck. By eliminating most of the human efforts involved in visualizing results through postprocessing, gLucifer reconnects the scientist to the numerical experiment as it unfolds. Data sets that were previously very difficult to even manage may be efficiently explored and interrogated without writing to disk, and because this approach is based entirely on memory distributed across as many processors as are being utilized by the scientific application, the visualization solution is scalable into terabytes of data being rendered in real time.     

Submarine reworking of exhumed subcontinental mantle rocks: field evidence from the Lherz peridotites,French Pyrenees

In the Pyrenees, the lherzolites nowhere occur as continuous units. Rather, they always outcrop as restricted bodies, never more than 3 km wide, scattered across Mesozoic sedimentary units along the North Pyrenean Fault. We report the results of a detailed analysis of the geological setting of the Lherz massif (central Pyrenees), the type‐locality of lherzolites and one of the most studied occurrences of mantle rocks worldwide. The Lherz body is only 1.5 km long and belongs to a series of ultramafic bodies of restricted size (a few metres to some hundreds of metres), occurring within sedimentary formations composed mostly of carbonate breccias originating from the reworking of Mesozoic platform limestones and dolomites. The clastic formations also include numerous layers of polymictic breccias reworking lherzolitic clasts. These layers are found far from any lherzolitic body, implying that lherzolitic clasts cannot derive from the in situ fragmentation of an ultramafic body alone, but might also have been transported far away from their sources by sedimentary processes. A detailed analysis of the contacts between the Lherz ultramafic body and the surrounding limestones confirms that there is no fault contact and that sediments composed of ultramafic material have been emplaced into fissures within the brecciated carapace of the peridotites. These observations bear important constraints for the mode of emplacement of the lherzolite bodies. We infer that mantle exhumation may have occurred during Albian strike‐slip deformation linked to the rotation of Iberia along the proto‐North Pyrenean Fault.