Fluid entrapment and reequilibration during subduction and exhumation: A case study from the high-grade Nestos shear zone, Central Rhodope, Greece

A fluid inclusion study on metamorphic minerals of successive growth stages was performed on highly deformed paragneisses from the Nestos Shear Zone at Xanthi (Central Rhodope), in which microdiamonds provide unequivocal evidence for ultrahigh-pressure (UHP) metamorphism. The correlation of fluid inclusion density isochores and fluid inclusion reequilibration textures with geothermobarometric data and the relative chronology of micro- and macro-scale deformation stages allow a better understanding of both the fluid and metamorphic evolution along the P–T–d path. Textural evidence for subduction towards the NE is recorded by the orientation of intragranular NE-oriented fluid inclusion planes and the presence of single, annular fluid inclusion decrepitation textures. These textures occur within quartz “foam” structures enclosed in an earlier generation of garnets with prolate geometries and rarely within recrystallized matrix quartz, and reequilibrated both in composition and density during later stages of exhumation. No fluid inclusions pertaining to the postulated ultrahigh-pressure stage for microdiamond-bearing garnet–kyanite–gneisses have yet been found. The prolate shape of garnets developed during the earliest stages of exhumation that is recorded structurally by (L  S) tectonites, which subsequently accommodated progressive ductile SW shearing and folding up to shallow crustal levels. The majority of matrix kyanite and a later generation of garnet were formed during SW-directed shear under plane-strain conditions. Fluid inclusions entrapped in quartz during this stage of deformation underwent density loss and transformed to almost pure CO₂ inclusions by preferential loss of H₂O. Those inclusions armoured within garnet retained their primary 3-phase H₂O–CO₂ compositions. Reequilibration of fluid inclusions in quartz aggregates is most likely the result of recrystallization along with stress-induced, preferential H₂O leakage along dislocations and planar lattice defects which results in the predominance of CO₂ inclusions with supercritical densities. Carbonic fluid inclusions from adjacent kyanite–corundum-bearing pegmatoids and, the presence of shear-plane-parallel fluid inclusion planes within late quartz boudin structures consisting of pure CO₂-fluid inclusions with negative crystal shapes, bear witness of the latest stage of deformation by NE-directed extensional shear.This study shows that the textures of early fluid inclusions that formed already during the prograde metamorphic path can be preserved and used to derive information about the kinematics of subduction that is difficult to obtain from other sources. The textures of early inclusions, together with later generations of unaltered primary and secondary inclusions in metamorphic index minerals that can be linked to specific deformation stages and even P–T conditions, are a welcome supplement for the reconstruction of a rather detailed P–T–d path.     

Strontium isotope evolution of Late Permian and Triassic seawater

The ⁸⁷Sr/⁸⁶Sr values based on brachiopods and conodonts define a nearly continuous record for the Late Permian and Triassic intervals. Minor gaps in measurements exist only for the uppermost Brahmanian, lower part of the Upper Olenekian, and Middle Norian, and only sparse data are available for the Late Permian. These 219 measurements include 67 brachiopods and 114 conodont samples from the Tethyan realm as well as 37 brachiopods and one conodont sample from the mid-European Middle Triassic Muschelkalk Sea. The Late Permian/Lower Triassic interval is characterized by a steep 1.3 × 10⁻³ rise, from 0.7070 at the base of the Dzhulfian to 0.7082 in the late Olenekian, a rate of change comparable to that in the Cenozoic. In the mid-Triassic (Anisian and Ladinian), the isotope values fall to 0.7075, followed again by a rise to 0.7081 in the Middle/Late Norian. The ⁸⁷Sr/⁸⁶Sr values decline again in the Late Norian (Sevatian) and Rhaetian to 0.7076.The sharp rise in the ⁸⁷Sr/⁸⁶Sr values during the Late Permian/Early Triassic was coincident with widespread clastic sedimentation. Because of the paucity of tectonic uplifts, the enhanced erosion may have been due to intermittent humid phases, during mainly an arid interval, coupled with the absence of a dense protective land plant cover following the mass extinction during the latest Permian. The apex of the ⁸⁷Sr/⁸⁶Sr curve at the Olenekian/Anisian boundary coincides with cessation of the large-scale clastic sedimentation and also marks the final recovery of land vegetation, as indicated by the renewed onset of coal formation in the Middle Triassic. The rising ⁸⁷Sr/⁸⁶Sr values from the Middle Carnian to the Late Norian coincide with the uplift and erosion of the Cimmeride-Indosinian orogens marking the closure of the Palaeotethys. The subsequent Rhaetian decline that continues into Jurassic (Pliensbachian/Toarcian boundary), on the other hand, coincides with the opening of the Vardar Ocean and its eastern continuation in the Izmir-Ankara Ophiolitic Belt.Samples from the Upper Muschelkalk are more radiogenic than the global trend. This may reflect separation of the basin from the open ocean. Due to strong meteoric influx from a large land mass in the north, the Germanic Basin became increasing brackish up section in the north and east, but because of the high evaporation rates, the salt content was not much reduced in the southern and central basin where a rich, but increasingly endemic, marine fauna survived.     

Mapping the hydraulic connection between a coalbed and adjacent aquifer: example of the coal-seam gas resource area,north Galilee Basin,Australia

Coal-seam gas production requires groundwater extraction from coal-bearing formations to reduce the hydraulic pressure and improve gas recovery. In layered sedimentary basins, the coalbeds are often separated from freshwater aquifers by low-permeability aquitards. However, hydraulic connection between the coalbed and aquifers is possible due to the heterogeneity in the aquitard such as the existence of conductive faults or sandy channel deposits. For coal-seam gas extraction operations, it is desirable to identify areas in a basin where the probability of hydraulic connection between the coalbed and aquifers is low in order to avoid unnecessary loss of groundwater from aquifers and gas production problems. A connection indicator, the groundwater age indictor (GAI), is proposed, to quantify the degree of hydraulic connection. The spatial distribution of GAI can indicate the optimum positions for gas/water extraction in the coalbed. Depressurizing the coalbed at locations with a low GAI would result in little or no interaction with the aquifer when compared to the other positions. The concept of GAI is validated on synthetic cases and is then applied to the north Galilee Basin, Australia, to assess the degree of hydraulic connection between the Aramac Coal Measure and the water-bearing formations in the Great Artesian Basin, which are separated by an aquitard, the Betts Creek Beds. It is found that the GAI is higher in the western part of the basin, indicating a higher risk to depressurization of the coalbed in this region due to the strong hydraulic connection between the coalbed and the overlying aquifer.     

Factors affecting the Nd3+ (REE3+) luminescence of minerals

In this paper, possibilities and limits of the application of REE³⁺ luminescence (especially the Nd^{3+ 4}F_3/2 → ⁴I_9/2 emission) as structural probe are evaluated. Important factors controlling the Nd³⁺ luminescence signal are discussed, including effects of the crystal-field, crystal orientation, structural state, and temperature. Particular attention was paid to the study of the accessory minerals zircon (ZrSiO₄), xenotime–(Y) (YPO₄), monazite–(Ce) (CePO₄) and their synthetic analogues. Based on these examples we review in short that (1) REE³⁺ luminescence can be used as non-destructive phase identification method, (2) the intensities of certain luminescence bands are strongly influenced by crystal orientation effects, and (3) increased widths of REE³⁺-related emission bands are a strong indicator for structural disorder. We discuss the potential of luminescence spectroscopy, complementary to Raman spectroscopy, for the quantitative estimation of chemical (and potentially also radiation-induced) disorder. For the latter, emissions of Nd³⁺-related centres are found to be promising candidates.     

Alkaline leaching characteristics of biogenic opal in Eocene sediments from the central Arctic Ocean: a case study in the ACEX cores

Accurate knowledge of the extent of biogenic opal preservation in marine sediment cores is important for paleoceanographic reconstructions. The alkaline leaching method is widely employed for %biogenic opal analysis due to its ease and speed. In this study, a revised method for measuring %biogenic opal in sediment from arctic coring expedition samples was suggested. The studied middle Eocene sediments from the central Arctic Ocean presented a problem in insufficiently leaching biogenic opal with a Na₂CO₃ solution. Based on XRD analysis, it was suggested that such an alkaline resistance results from slight diagenesis of biogenic opal. In order to solve this problem, an alkaline leaching method utilizing a 2 M NaOH solution was suggested for the accurate measurement of %biogenic opal in the Eocene sediments from the central Arctic. Furthermore, dissolution rates from lithogenic matter by NaOH solution were measured in order to correct the %biogenic opal values.     

Glacier inventory of the Gran Campo Nevado Ice Cap in the Southern Andes and glacier changes observed during recent decades

The Gran Campo Nevado (GCN) forms an isolated ice cap on the Península Muñoz Gamero (PMG) located 200 km to the south of the Southern Patagonia Icefield (SPI). We present a glacier inventory of the GCN made up by 27 drainage basins (in total 199.5 km²) and other small cirque and valley glaciers of the southern part of PMG (in total 53 km²). The glacier inventory is based on a digital elevation model (DEM) and ortho-photos. Contour lines from maps, relief information derived from Landsat TM satellite imagery from 1986 and 2002 and stereoscopic data from aerial photos were combined in a knowledge-based scheme to obtain a DEM of the area. A digital ortho-photo map based on aerial photos from 1998 and several ortho-photos based on aerial photos from 1942 and 1984 could be produced from the initial DEM. A geographical information system (GIS) served to outline the extent of the present glaciation. All major glaciers of the GCN show a significant glacier retreat during the last 60 yr. Some of the outlet glaciers lost more than 20% of their total area during this period. Overall glacier retreat amounts to 2.8% of glacier length per decade and the glacier area loss is 2.4% per decade in the period from 1942 to 2002. We hypothesise that GCN glaciers may have reacted faster and more synchronously with the observed warming trend during recent decades when compared with the SPI.     

Characterizing hydrological processes within kettle holes using stable water isotopes in the Uckermark of northern Brandenburg,Germany

Understanding the hydrologic connectivity between kettle holes and shallow groundwater, particularly in reaction to the highly variable local meteorological conditions, is of paramount importance for tracing water in a hydro(geo)logically complex landscape and thus for integrated water resource management. This article is aimed at identifying the dominant hydrological processes affecting the kettle holes' water balance and their interactions with the shallow groundwater domain in the Uckermark region, located in the north-east of Germany. For this reason, based on the stable isotopes of oxygen (δ¹⁸O ) and hydrogen (δ²H ), an isotopic mass balance model was employed to compute the evaporative loss of water from the kettle holes from February to August 2017. Results demonstrated that shallow groundwater inflow may play the pivotal role in the processes taking part in the hydrology of the kettle holes in the Uckermark region. Based on the calculated evaporation/inflow (E/I) ratios, most of the kettle holes (86.7%) were ascertained to have a partially open, flow-through-dominated system. Moreover, we identified an inverse correlation between E/I ratios and the altitudes of the kettle holes. The same holds for electrical conductivity (EC) and the altitudes of the kettle holes. In accordance with the findings obtained from this study, a conceptual model explaining the interaction between the shallow groundwater and the kettle holes of Uckermark was developed. The model exhibited that across the highest altitudes, the recharge kettle holes are dominant, where a lower ratio of E/I and a lower EC was detected. By contrast, the lowest topographical depressions represent the discharge kettle holes, where a higher ratio of E/I and EC could be identified. The kettle holes existing in between were categorized as flow-through kettle holes through which the recharge takes place from one side and discharge from the other side.     

Conceptual examination of climate–soil controls upon rainfall partitioning in an open-fractured soil II: Response to a population of storms

A simple two-domain bucket model of fractured soil was coupled with a stochastic model of rainfall variability, in order to investigate the climate and soil controls upon the stochastic properties of the triggering of fracture flow and surface runoff, and the partitioning of rainfall between the matrix and fracture domains and surface runoff. Conventionally, soils are regarded as time domain filters between rainfall and hydrological response. This investigation highlights an additional type of threshold filtering especially important in understanding the infiltration behaviour of fractured soils, for which an event-based characterisation of rainfall in modelling is crucial. A priori-definable indices were derived which are capable of describing elements of this threshold filtering, by allowing the statistical properties of fracture flow- and surface runoff-triggering storms (i.e., mean and variance of storm duration, intensity and effective inter-storm period, as well as cumulative partitioning of rainfall), to be inferred directly from average storm and soil properties. Using these indices, the long-term response of fractured soils, including the long-term hydrological importance of fractures, can be estimated without simulation.     

Lunar noise correlation,imaging and monitoring

Passive seismic techniques have revolutionarised seismology,leading for example to increased resolution in surface wave tomography,to the possibility to monitor changes in the propagation medium,and to many new processing strategies in seismic exploration.Here we review applications of the new techniques to a very particular dataset,namely data from the Apollo 17 lunar network.The special conditions of the lunar noise environment are investigated,illustrating the interplay between the properties of the noise and the ability to reconstruct Green's functions.With a dispersion analysis of reconstructed Rayleigh waves new information about the shallow shear velocity structure of the Moon are obtained.Passive image interferometry is used to study the effect of temperature changes in the subsurface on the seismic velocities providing direct observation of a dynamic process in the lunar environment.These applications highlight the potential of passive techniques for terrestrial and planetary seismology.     

Flux contribution of coherent structures and its implications for the exchange of energy and matter in a tall spruce canopy

The flux contribution of coherent structures to the total exchange of energy and matter is investigated in a spruce canopy of moderate density in heterogeneous, complex terrain. The study deploys two methods of analysis to estimate the coherent exchange: conditional averages in combination with wavelet analysis, and quadrant analysis. The data were obtained by high-frequency single-point measurements using sonic anemometers and gas analysers at five observation heights above and within the canopy and subcanopy, and represent a period of up to 2.5 months. The study mainly addresses the momentum transfer and exchange of sensible heat throughout the roughness sublayer, while results are provided for the exchange of carbon dioxide and water vapour above the canopy. The magnitude of the flux contribution of coherent structures largely depends on the method of analysis, and it is demonstrated that these differences are attributed to differences in the sampling strategy between the two methods. Despite the differences, relational properties such as sweep and ejection ratios and the variation of the flux contribution with height were in agreement for both methods. The sweep phase of coherent structures is the dominant process close to and within the canopy, whereas the ejections gain importance with increasing distance to the canopy. The efficiency of the coherent exchange in transporting scalars exceeds that for momentum by a factor of two. The occurrence of coherent structures results in a flux error less than 4% for the eddy-covariance method. Based on the physical processes identified from the analysis of the ejection and sweep phases along the vertical profile in the roughness sublayer, a classification scheme for the identification of exchange regimes is developed. This scheme allows one to estimate the region of the canopy participating in the exchange of energy and matter with the above-canopy air under varying environmental conditions.