Geomechanical model for the post-Variscan evolution of the Permocarboniferous–Mesozoic basins in Northeast Germany

The Permocarboniferous basins in Northeast Germany formed on the heterogeneous and eroded parts of the Variscan orogene and its deformed northern foreland. Transtensional tectonic movements and thermal re-equilibration lead to medium-scale crustal fragmentation, fast subsidence rates and regional emplacement of large amounts of mostly acidic volcanics. The later basin formation and differentiation was triggered by reversals of the large-scale stress field and reactivation of prominent zones of weakness like the Elbe Fault System and the Rhenohercynian/Saxothuringian boundary that separate different Variscan basement domains in the area. The geomechanical behaviour of the latter plays an important role for the geodynamic evolution of the medium to large-scale structural units, which we can observe today in three dimensions on structural maps, geophysical recordings and digital models. This study concentrates on an area that comprises the southern Northeast German Basin, the Saale Basin, the Flechtingen High, the Harz Mountains High and the Subhercynian Basin. The presented data include re-evaluations of special geological and structural maps, the most recent interpretation of the DEKORP BASIN 9601 seismic profile and observations of exposed rock sections in Northeast Germany. On the basis of different structural inventories and different basement properties, we distinguish two structural units to the south and one structural unit to the north of the Elbe Fault System. For each unit, we propose a geomechanical model of basin formation and basin inversion, and show that the Rhenohercynian Fold and Thrust Belt domain is deformed in a thin-skinned manner, while the Mid-German Crystalline Rise Domain, which is the western part of the Saxothuringian Zone, rather shows a thick-skinned deformation pattern. The geomechanical model for the unit north to the Elbe Fault System takes account to the fact that the base of the Zechstein beneath the present Northeast German basin shows hardly any evidence for brittle deformation, which indicates a relative stable basement. Our geomechanical model suggests that the Permocarboniferous deposits may have contributed to the structural stiffness by covering small to medium scale structures of the upper parts of the brittle basement. It is further suggested that the pre-Zechstein successions underneath the present Northeast German basin were possibly strengthening during the Cretaceous basin inversion, which resulted in stress transfer to the long-lived master faults, as indicated for example by the shape of the salt domes in the vicinity of the latter faults. Contrary to this, post-Zechstein successions deformed in a different and rather complex way that was strongly biased by intensive salt tectonic movements.     

Influence of compound bedforms on hydraulic roughness in a tidal environment

The effect exerted by the seabed morphology on the flow is commonly expressed by the hydraulic roughness, a fundamental parameter in the understanding and simulation of hydro- and sediment dynamics in coastal areas. This study quantifies the hydraulic roughness of large compound bedforms throughout a tidal cycle and investigates its relationship to averaged bedform dimensions. Consecutive measurements with an acoustic Doppler current profiler and a multibeam echosounder were carried out in the Jade tidal channel (North Sea, Germany) along large compound bedforms comprising ebb-oriented primary bedforms with superimposed smaller secondary bedforms. Spatially averaged velocity profiles produced log-linear relationships which were used to calculate roughness lengths. During the flood phase, the velocity profiles were best described by a single log-linear fit related to the roughness created by the secondary bedforms. During the ebb phase, the velocity profiles were segmented, showing the existence of at least two boundary layers: a lower one scaling with the superimposed secondary bedforms and an upper one scaling with the ebb-oriented primary bedforms. The drag induced by the primary bedform during the ebb phase is suggested to be related to flow expansion, separation, and recirculation on the downstream side of the bedform. Three existing formulas were tested to predict roughness lengths from the local bedform dimensions. All three predicted the right order of magnitude for the average roughness length but failed to predict its variation over the tidal cycle.     

The fractal geometry of flow paths in natural fractures in rock and the approach to percolation

The distributions of contact areas in single, natural fractures in quartz monzonite (Stripa granite) are found to have fractal dimensions which decrease fromD=2.00 to values nearD=1.96 as stress normal to the fractures is increased from 3 MPa up to 85 MPa. The effect of stress on fluid flow is studied in the same samples. Fluid transport through a fracture depends on two properties of the fracture void space geometry. the void aperture; and the tortuosity of the flow paths, determined through the distribution of contact area. Each of these quantities change under stress and contribute to changes observed in the flow rate. A general flow law is presented which separates these different effects. The effects of tortuosity on flow are largely governed by the proximity of the flow path distribution to a percolation threshold. A fractal model of correlated continuum percolation is presented which quantitatively reproduces the flow path geometries. The fractal dimension in this model is fit to the measured fractal dimensions of the flow systems to determine how far the flow systems are above the percolation threshold.     

Restructuring of polygalacturonate on alumina upon hydration—Effect on phosphate sorption kinetics

Hydration of organic coatings in soils is expected to affect the sorption of oxyanions onto hydrous Fe and Al oxides. We hypothesized that the hydration of polygalacturonate (PGA) coatings on alumina (Al₂O₃) increases their permeability for phosphate. Pure and PGA-coated alumina were equilibrated in deionized water for 2 and 170 h at pH 5 and 20 °C before studying (i) their porosity with N₂ gas adsorption and ¹H NMR relaxometry, (ii) structural changes of PGA-coatings with differential scanning calorimetry (DSC), and (iii) the kinetics of phosphate sorption and PGA desorption in batch experiments. Scanning electron micrographs revealed that PGA molecules formed three-dimensional networks with pores ranging in size from <10 to several hundred nanometers. Our NMR results showed that the water content of intraparticle alumina pores decreased upon PGA sorption, indicating a displacement of pore water by PGA. The amount of water in interparticle alumina pores increased strongly after PGA addition, however, and was attributed to water in pores of PGA and/or in pores at the PGA-alumina interface. The flexibility of PGA molecules and the fraction of a PGA gel phase increased within one week of hydration, implying restructuring of PGA. Hydration of PGA coatings increased the amount of phosphate defined as instantaneously sorbed by 84%, showing that restructuring of PGA enhanced the accessibility of phosphate to external alumina surfaces. Despite the fact that the efficacy of phosphate to displace PGA was higher after 170 h than after 2 h, a higher phosphate surface loading was required after 170 h to set off PGA desorption. Our findings imply that the number of PGA chain segments directly attached to the alumina surface decreased with time. We conclude that hydration/dehydration of polymeric surface coatings affects the sorption kinetics of oxyanions, and may thus control the sorption and transport of solutes in soils.     

Linking tracers,water age and conceptual models to identify dominant runoff processes in a sparsely monitored humid tropical catchment

Assessing catchment runoff response remains a key research frontier because of limitations in current observational techniques to fully characterize water source areas and transit times in diverse geographical environments. Here, we report a study that combines empirical data with modelling to identify dominant runoff processes in a sparsely monitored humid tropical catchment. The analysis integrated isotope tracers into conceptual rainfall–runoff models of varying complexity (from 5 to 11 calibrated parameters) that are able to simulate discharge and tracer concentrations and track the evolving age of stream water exiting the catchment. The model structures can be seen as competing hypotheses of catchment functioning and were simultaneously calibrated against uncertain streamflow gaugings and a 2‐year daily isotope rainfall–runoff record. Comparison of the models was facilitated using global parameter sensitivity analysis and the resulting effect on calibration. We show that a variety of tested model structures reproduced water and tracer dynamics in stream, but the simpler models failed to adequately reproduce both. The resulting water age distributions of the tested models varied significantly with little similarity between the stream water age and stored water age distributions. The sensitivity analysis revealed that only some of the more complex models (from eight parameters) could be better constrained to infer more plausible water age distributions and catchment storage estimates. These models indicated that the age of water stored in the catchment is generally older compared with the age of water fluxes, with evapotranspiration age being younger compared with streamflow. However, the water age distributions followed a similar temporal behaviour dominated by climatic seasonality. Stream water ages increased during the dry season (greater than 1 year) and decreased with increased streamflow (a few weeks old) during the wet season. We further show that the ratios of the streamwater age to stored water age distribution and the water age distribution of actual evapotranspiration to the stored water age distribution from constrained models could potentially serve as useful hydrological indicators of catchment functioning. Copyright © 2016 John Wiley & Sons, Ltd.     

Relationship between footwall composition,crustal contamination,and fluid–rock interaction in the Platreef,Bushveld Complex,South Africa

In the northern limb of the 2.06-Ga Bushveld Complex, the Platreef is a platinum group elements (PGE)-, Cu-, and Ni-mineralized zone of pyroxenite that developed at the intrusion margin. From north to south, the footwall rocks of the Platreef change from Archaean granite to dolomite, hornfels, and quartzite. Where the footwall is granite, the Sr-isotope system is more strongly perturbed than where the footwall is Sr-poor dolomite, in which samples show an approximate isochron relationship. The Nd-isotope system for samples of pyroxenite and hanging wall norite shows an approximate isochron relationship with an implied age of 2.17 ± 0.2 Ga and initial Nd-isotope ratio of 0.5095. Assuming an age of 2.06 Ga, the ɛNd values range from −6.2 to −9.6 (ave. −7.8, n = 17) and on average are slightly more negative than the Main Zone of the Bushveld. These data are consistent with local contamination of an already contaminated magma of Main Zone composition. The similarity in isotope composition between the Platreef pyroxenites and the hanging wall norites suggests a common origin. Where the country rock is dolomite, the Platreef has generally higher plagioclase and pyroxene δ ¹⁸O values, and this indicates assimilation of the immediate footwall. Throughout the Platreef, there is considerable petrographic evidence for sub-solidus interaction with fluids, and the Δ _{plagioclase–pyroxene} values range from −2 to +6, which indicates interaction at both high and low temperatures. Whole-rock and mineral δD values suggest that the Platreef interacted with both magmatic and meteoric water, and the lack of disturbance to the Sr-isotope system suggests that fluid–rock interaction took place soon after emplacement. Where the footwall is granite, less negative δD values suggest a greater involvement of meteoric water. Consistently higher values of Δ _{plagioclase–pyroxene} in the Platreef pyroxenites and hanging wall norites in contact with dolomite suggest prolonged interaction with CO₂-rich fluid derived from decarbonation of the footwall rocks. The overprint of post crystallization fluid–rock interaction is the probable cause of the previously documented lack of correlation between PGE and sulfide content on the small scale. The Platreef in contact with dolomite is the focus of the highest PGE grades, and this suggests that dolomite contamination played a role in PGE concentration and deposition, but the exact link remains obscure. It is a possibility that the CO₂ produced by decarbonation of assimilated dolomite enhanced the process of PGE scavenging by sulfide precipitation.     

Efficient and accurate high-degree spherical harmonic synthesis of gravity field functionals at the Earth’s surface using the gradient approach

Spherical harmonic synthesis (SHS) of gravity field functionals at the Earth’s surface requires the use of heights. The present study investigates the gradient approach as an efficient yet accurate strategy to incorporate height information in SHS at densely spaced multiple points. Taylor series expansions of commonly used functionals quasigeoid heights, gravity disturbances and vertical deflections are formulated, and expressions of their radial derivatives are presented to arbitrary order. Numerical tests show that first-order gradients, as introduced by Rapp (J Geod 71(5):282–289, 1997) for degree 360 models, produce cm- to dm-level RMS approximation errors over rugged terrain when applied with EGM2008 to degree 2190. Instead, higher-order Taylor expansions are recommended that are capable of reducing approximation errors to insignificance for practical applications. Because the height information is separated from the actual synthesis, the gradient approach can be applied along with existing highly efficient SHS routines to compute surface functionals at arbitrarily dense grid points. This confers considerable computational savings (above or well above one order of magnitude) over conventional point-by-point SHS. As an application example, an ultra-high resolution model of surface gravity functionals (EurAlpGM2011) is constructed over the entire European Alps that incorporates height information in the SHS at 12,000,000 surface points. Based on EGM2008 and residual topography data, quasigeoid heights, gravity disturbances and vertical deflections are estimated at ~200m resolution. As a conclusion, the gradient approach is efficient and accurate for high-degree SHS at multiple points at the Earth’s surface.     

Modeling Variably Saturated Subsurface Solute Transport with MODFLOW‐UZF and MT3DMS

The MT3DMS groundwater solute transport model was modified to simulate solute transport in the unsaturated zone by incorporating the unsaturated‐zone flow (UZF1) package developed for MODFLOW. The modified MT3DMS code uses a volume‐averaged approach in which Lagrangian‐based UZF1 fluid fluxes and storage changes are mapped onto a fixed grid. Referred to as UZF‐MT3DMS, the linked model was tested against published benchmarks solved analytically as well as against other published codes, most frequently the U.S. Geological Survey's Variably‐Saturated Two‐Dimensional Flow and Transport Model. Results from a suite of test cases demonstrate that the modified code accurately simulates solute advection, dispersion, and reaction in the unsaturated zone. Two‐ and three‐dimensional simulations also were investigated to ensure unsaturated‐saturated zone interaction was simulated correctly. Because the UZF1 solution is analytical, large‐scale flow and transport investigations can be performed free from the computational and data burdens required by numerical solutions to Richards' equation. Results demonstrate that significant simulation runtime savings can be achieved with UZF‐MT3DMS, an important development when hundreds or thousands of model runs are required during parameter estimation and uncertainty analysis. Three‐dimensional variably saturated flow and transport simulations revealed UZF‐MT3DMS to have runtimes that are less than one tenth of the time required by models that rely on Richards' equation. Given its accuracy and efficiency, and the wide‐spread use of both MODFLOW and MT3DMS, the added capability of unsaturated‐zone transport in this familiar modeling framework stands to benefit a broad user‐ship.     

Tracing the exhumation of the Eclogite Zone (Tauern Window,Eastern Alps) by <Superscript>40</Superscript>Ar/<Superscript>39</Superscript>Ar dating of white mica in eclogites

New radiometric ages from the Subpenninic nappes (Eclogite Zone and Rote Wand – Modereck Nappe, Tauern Window) show that phengites formed under eclogite-facies metamorphic conditions retain their initial isotopic signature, even when associated lithologies were overprinted by greenschist- to amphibolite-facies metamorphism. Different stages of the eclogite-facies evolution can be dated provided ⁴⁰Ar/³⁹Ar dating is combined with micro-structural analyses. An age of 39 Ma from the Rote Wand – Modereck Nappe is interpreted to be close to the burial age of this unit. Eclogite deformation within the Eclogite Zone started at the pressure peak along distinct shear zones, and prevailed along the exhumation path. An age of ca. 38 Ma is only observed for eclogites not affected by subsequent deformation and is interpreted as maximum age due to the possible influence of homogenously distributed excess argon. During exhumation deformation was localised along distinct mylonitic shear zones. This stage is mainly characterised by the formation of dynamically recrystallized omphacite2 and phengite. Deformation resulted in the resetting of the Ar isotopic system within the recrystallized white mica. Flat argon release spectra showing ages of 32 Ma within mylonites record the timing of cooling along the exhumation path, and the emplacement onto the Venediger Nappe. Ar-release patterns and ³⁶Ar/⁴⁰Ar vs.³⁹Ar/⁴⁰Ar isotope correlation analyses indicate no significant ⁴⁰Ar-loss after initial closure, and only a negligible incorporation of excess argon. From the pressure peak onwards, eclogitic conditions prevailed for almost 8–10 Ma.