Conditions for veining in the Barrandian Basin (Lower Palaeozoic), Czech Republic: evidence from fluid inclusion and apatite fission track analysis

The interplay between fracture propagation and fluid composition and circulation has been examined by deciphering vein sequences in Silurian and Devonian limestones and shales at Kosov quarry in the Barrandian Basin. Three successive vein generations were recognised that can be attributed to different stages of a basinal cycle. Almost all generations of fracture cements host abundant liquid hydrocarbon inclusions that indicate repeated episodes of petroleum migration through the strata during burial, tectonic compression and uplift.The earliest veins that propagated prior to folding were displacive fibrous “beef” calcite veins occurring parallel to the bedding of some shale beds. Hydrocarbon inclusions within calcite possess homogenisation temperatures between 58 and 68 °C and show that the “beef” calcites originated in the deeper burial environment, during early petroleum migration from overpressured shales.E–W-striking extension veins that postdate “beef” calcite formed in response to Variscan orogenic deformations. Based on apatite fission track analysis (AFTA) data and other geological evidence, the veins probably formed 380–315 Ma ago, roughly coinciding with peak burial heating of the strata, folding and the intrusion of Variscan synorogenic granites. The veins that crosscut diagenetic cements and low-amplitude stylolites in host limestones are oriented semi-vertically to the bedding plane and are filled with cloudy, twinned calcite, idiomorphic smoky quartz and residues of hardened bitumen. Calcite and quartz cements contain abundant blue and blue–green-fluorescing primary inclusions of liquid hydrocarbons that homogenise between 50 and 110 °C. Geochemical characteristics of the fluids as revealed by gas chromatography–mass spectrometry, particularly the presence of olefins and parent aromatic hydrocarbons (phenonthrene), suggest that the oil entrapped in the inclusions experienced intense but geologically fast heating that resulted in thermal pyrolysis of its hydrocarbons. This implies that the organic fluids in the fractures may have been partly influenced by heating associated with igneous intrusions that are hidden below the surface.Subvertical N–S-striking veins represent the most recent fracturing event(s). Some of these veins are only a few millimeters thick and sparsely mineralised with thin leaf-like quartz crystals that contain tiny blue and yellow–orange-fluorescing hydrocarbon inclusions. Most of the N–S veins, however, occur as thick calcite veins that generally crystallised at 70 °C or less from H₂O–NaCl solutions of variable salinity with admixture of petroleum. The origin of these fluids is interpreted in terms of deeply circulating meteoric waters that partially mixed with deep basinal fluids. Wider structural considerations combined with fission-track analysis of adjacent host sediments suggest that N–S veins formed during post-Mesozoic uplift of the area, probably in response to major Tertiary Alpine deformations transmitted far into the Bohemian Massif.     

The influence of porosity gradients on mixing coefficients in sediments

The surface layer of aquatic sediments is a zone characterized by both porosity gradients and intensive mixing. In the standard approach, porosity gradients are ignored when estimating mixing intensity. Here, model formulations with both constant and varying porosity are contrasted to estimate mixing coefficients D_b from tracer depth profiles. Complementing the well-known exponential solution of the constant-porosity model, we present a general solution to the variable-porosity model in terms of hypergeometric functions. When using these models in a forward way, the tracer activities predicted by the variable-porosity model are higher than those generated by the constant-porosity model. Similarly, when inverse modelling, D_b values estimated by the variable-porosity model are systematically higher than those derived from the constant-porosity model. Still, differences in D_b values remain relatively small. When applying both mixing models to excess ²¹⁰Pb data profiles from slope sediments, a maximal difference of 30% is obtained between D_b values, the average deviation being 16%. A systematic exploration of parameter space predicts a maximal underestimation of 60% when deriving D_b values from the constant-porosity mixing model. Given the uncertainty imposed by other model assumptions underlying the diffusive mixing model, the influence of porosity gradients on D_b values must be classified as rather modest. Hence, the current mixing coefficient database is not biased by the constant porosity approximation.     

Geomorphological research of large-scale slope instability at Machu Picchu, Peru

A multidisciplinary approach has been adopted to study the slope movements and landscape evolution at the archaeological site of Machu Picchu and its immediate surroundings. The basic event in the paleogeomorphological evolution of the area was the large-scale slope movement, which destroyed the originally higher ridge between Mt. Machupicchu and Mt. Huaynapicchu. Within remnants of that primary deformation, several younger generations of slope movements occurred. The laboratory analyses of granitoids revealed highly-strained zones on the slopes of Mt. Machupicchu, which strongly affect the largest slope deformation. The borders of the largest slope deformation are structurally predisposed by the existence of fault zones. The majority of various types of slope movements on the so-called Front Slope (E facing) and Back Slope (W facing) are influenced by the alignment between topography and joints. Along with slope movements, fluvial erosion and tectonic disturbance of the rocks have been affecting the evolution of the landscape. A monitoring network for dilatometric and extensometric measurements was used to detect the present-day activity of rock displacements within the archaeological site. In addition to standard mapping of surface hydrogeological phenomena, eleven express slug tests were conducted to verify the infiltration potential of precipitation. The results of these surveys indicate that recent large-scale slope movement as suggested by some previous studies is doubtful, and the detected movements can be explained by individual movements of rock blocks or several other mechanisms including sinking of archaeological structures, subsurface erosion and annual changes in the water content of the soils.     

The 1988–2003 Greenland ice sheet melt extent using passive microwave satellite data and a regional climate model

Measurements from ETH-Camp and JAR1 AWS (West Greenland) as well as coupled atmosphere-snow regional climate simulations have highlighted flaws in the cross-polarized gradient ratio (XPGR) technique used to identify melt from passive microwave satellite data. It was found that dense clouds (causing notably rainfall) on the ice sheet severely perturb the XPGR melt signal. Therefore, the original XPGR melt detection algorithm has been adapted to better incorporate atmospheric variability over the ice sheet and an updated melt trend for the 1988–2003 period has been calculated. Compared to the original algorithm, the melt zone area increase is eight times higher (from 0.2 to 1.7% year⁻¹). The increase is higher with the improved XPGR technique because rainfall also increased during this period. It is correlated to higher atmospheric temperatures. Finally, the model shows that the total ice sheet runoff is directly proportional to the melt extent surface detected by satellites. These results are important for the understanding of the effect of Greenland melting on the stability of the thermohaline circulation.     

The influence of extreme floods from the River Danube in 2006 on phytoplankton communities in a floodplain lake: Shift to a clear state

The influence of extreme floods from the River Danube in 2006 on the species composition and vertical distributions of phytoplankton was studied in a shallow floodplain lake, Lake Sakadaš (Kopa?ki Rit Nature Park, Croatia) which in the last few decades was in a turbid state characterised by high phytoplankton concentrations. As a consequence of extremely high floods, the whole floodplain area (approximately 16 km²) became one lentic habitat with well developed macrophyte vegetation. Seasonal dynamics of chlorophyll a (Chl a) concentration in the lake had a characteristic pattern for the shallow lakes with dense macrophyte vegetation. Extremely low mean phytoplankton abundance and biomass were found in the conditions of very high nutrient concentrations. Dominant phytoplankton species were diatoms and chlorococcal green algae from the functional groups characteristic for a mixed environment. The canonical correspondence analysis (CCA) demonstrated that nutrients and temperature were significant environmental variables for their development. The sequence of phytoplankton seasonality, vertical distribution of phytoplankton, as well as the domination of rapidly acclimating phytoplankton forms (R-strategists) indicated clear, well-mixed conditions and a highly disturbed environment. Our results suggest that the occurrence of extreme flooding can be a stressor high enough for the transition from a turbid to a clear state of the floodplain lake. Possibly, cyclic shifts between alternative stable states in floodplain ecosystems can be expected as a consequence of the impact of extreme hydrological events induced by a climate change.     

Modeling reactive transport in sediments subject to bioturbation and compaction

Current bioturbation models are marked by confusion in their treatment of porosity. Different equations appear to be needed for different biodiffusion mechanisms, i.e., interphase mixing, where biological activity causes bulk mixing of sediment affecting both tracer and porosity profiles, versus intraphase mixing, where the solid components are intermixed, but the porosity is left unchanged. Another issue is whether the model depends upon the particle type with which tracers are associated, e.g., ¹³⁷Cs on small clay particles versus ²¹⁰Pb on larger grains. This uncertainty has lead to conflicting conservation equations for radiotracers, and in particular, to the question whether the porosity should be placed inside or outside of the differential term that governs the biodiffusive flux. We have reexamined this situation in the context of multiphase, multicomponent continuum theory. Most importantly, we prove that under the assumption of steady-state porosity, there exists only one correct form of the steady-state conservation equation for a radiotracer, regardless of biodiffusion mechanism and particle type, i.e.,

     

Greenland surface mass balance simulated by a regional climate model and comparison with satellite-derived data in 1990–1991

The 1990 and 1991 ablation seasons over Greenland are simulated with a coupled atmosphere-snow regional climate model with a 25-km horizontal resolution. The simulated snow water content allows a direct comparison with the satellite-derived melt signal. The model is forced with 6-hourly ERA-40 reanalysis at its boundaries. An evaluation of the simulated precipitation and a comparison of the modelled melt zone and the surface albedo with remote sensing observations are presented. Both the distribution and quantity of the simulated precipitation agree with observations from coastal weather stations, estimates from other models and the ERA-40 reanalysis. There are overestimations along the steep eastern coast, which are most likely due to the “topographic barrier effect”. The simulated extent and time evolution of the wet snow zone compare generally well with satellite-derived data, except during rainfall events on the ice sheet and because of a bias in the passive microwave retrieved melt signal. Although satellite-based surface albedo retrieval is only valid in the case of clear sky, the interpolation and the correction of these data enable us to validate the simulated albedo on the scale of the whole Greenland. These two comparisons highlight a large sensitivity of the remote sensing observations to weather conditions. Our high-resolution climate model was used to improve the retrieval algorithms by taking more fully into account the atmosphere variability. Finally, the good agreement of the simulated melting surface with the improved satellite signal allows a detailed estimation of the melting volume from the simulation.