The application of GPR and ERI in combination with exposure logging and retrodeformation analysis to characterize sinkholes and reconstruct their impact on fluvial sedimentation

This work illustrates the practicality of investigating sinkholes integrating data gathered by ground penetrating radar (GPR), electrical resistivity imaging (ERI) and trenching or direct logging of the subsidence‐affected sediments in combination with retrodeformation analysis. This mutidisciplinary approach has been tested in a large paleosinkhole developed during the deposition of a Quaternary terrace on salt‐bearing evaporites. The subsidence structure, exposed in an artificial excavation, is located next to Puilatos, a village that was abandoned in the 1970s due to severe subsidence damage. Detailed logging of the exposure revealed that the subsidence structure corresponds to an asymmetric sagging and collapse paleosinkhole with no clear evidence of recent activity. The sedimentological and structural relationships together with the retrodeformation analysis indicate that synsedimentary subsidence controlled channel location, the development of a palustrine environment and local changes in the channel pattern. GPR profiles were acquired using an array of systems with different antenna frequencies, including some recently developed shielded antennas with improved vertical resolution and penetration depth. Although radargrams imaged the faulted sagging structure and provided valuable data on fault throw, they did not satisfactorily image the complex architecture of the fluvial deposit. ERI showed lower resolution but higher penetration depth when compared to GPR, roughly capturing the subsidence structure and yielding information on the thickness of the high‐resistivity alluvium and the nature of the underlying low‐resistivity karstic residue developed on top of the halite‐bearing evaporitic bedrock. Data comparison allows the assessment of the advantages and limitations of these complementary techniques, highly useful for site‐specific sinkhole risk management. Copyright © 2016 John Wiley & Sons, Ltd.     

Evaluation of trenching,ground penetrating radar (GPR) and electrical resistivity tomography (ERT) for sinkhole characterization

This paper explores the suitability and advantages of combining the trenching technique with geophysical surveys [ground penetrating radar (GPR) and electrical resistivity tomography (ERT)] for sinkhole characterization in a mantled karst area. The approach is applied to two active sinkholes concealed by anthropogenic deposits and formed by contrasting subsidence mechanisms; collapse and sagging. The ERT section acquired across the collapse sinkhole images the clayey fill of the depression as an obvious low resistivity area, showing the approximate location of the sinkhole edges. Spatially dense GPR surveys provide information on the position of the boundaries of the concealed subsidence structures and their three‐dimensional (3D) internal geometry, revealing the dominant subsidence mechanism. We illustrate the impact of several factors on the quality of the GPR data such as sinkhole size, nominal frequency of the antennas, antenna shielding, and the presence of backfilled excavations and above‐surface objects. Trenches provided detailed information on the subsurface structure of the sinkhole, subsidence magnitude, partitioning of the strain, and the position of the sinkhole edges, especially when they are deep enough and excavated across the central sector and perpendicular to the boundaries. The stratigraphic and structural relationships observed in the trench were then used to infer the spatial evolution of the sinkholes (e.g. enlargement), their kinematic behavior (episodic versus progressive), and to differentiate discrete subsidence events and their associated magnitude. Numerical dates were used to estimate average subsidence rates and the recurrence of subsidence events. Such integrated data sets may be used as an objective basis to forecast the future behavior of potentially damaging sinkholes and to assess the associated hazard and risk. Copyright © 2013 John Wiley & Sons, Ltd.     

Detection, modelling and estimation of non-linear trends by using a non-homogeneous Vasicek stochastic diffusion. Application to CO2 emissions in Morocco

This paper proposes a new stochastic model, based on a Vasicek non-homogeneous diffusion process, in which the non-linear trend coefficient (drift) depends on deterministic functions that describe the dynamic evolution of certain exogenous variables. After studying its probabilistic characteristics, and in particular the transition probability density and trend function, the associated stochastic inference based on discrete sampling in time is established using maximum likelihood methodology. This model is applied to detect, estimate and model the non-linear trend present in data corresponding to CO₂ emissions in Morocco. Energy and financial variables that affect the behaviour of this trend are also detected, and substantial improvement provided by this non-homogeneous model with respect to its corresponding homogeneous version, is confirmed.     

High‐pressure metamorphic evolution of eclogite and associated metapelite from the Chuacús complex (Guatemala Suture Zone): Constraints from phase equilibria modelling coupled with Lu‐Hf and U‐Pb geochronology

As is common in suture zones, widespread high‐pressure rocks in the Caribbean region reached eclogite facies conditions close to ultrahigh‐pressure metamorphism. Besides eclogite lenses, abundant metapelitic rocks in the Chuacús complex (Guatemala Suture Zone) also preserve evidence for high‐pressure metamorphism. A comprehensive petrological and geochronological study was undertaken to constrain the tectonometamorphic evolution of eclogite and associated metapelite from this area in central Guatemala. The integration of field and petrological data allows the reconstruction of a previously unknown segment of the prograde P–T path and shows that these contrasting rock types share a common high‐pressure evolution. An early stage of high‐pressure/low‐temperature metamorphism at 18–20 kbar and 530–580°C is indicated by garnet core compositions as well as the nature and composition of mineral inclusions in garnet, including kyanite–jadeite–paragonite in an eclogite, and chloritoid–paragonite–rutile in a pelitic schist. Peak high‐pressure conditions are constrained at 23–25 kbar and 620–690°C by combining mineral assemblages, isopleth thermobarometry and Zr‐in‐rutile thermometry. A garnet/whole‐rock Lu‐Hf date of 101.8 ± 3.1 Ma in the kyanite‐bearing eclogite indicates the timing of final garnet growth at eclogite facies conditions, while a Lu‐Hf date of 95.5 ± 2.1 Ma in the pelitic schist reflects the average age of garnet growth spanning from an early eclogite facies evolution to a final amphibolite facies stage. Concordant U‐Pb LA‐ICP‐MS zircon data from the pelitic schist, in contrast, yield a mean age of 74.0 ± 0.5 Ma, which is equivalent to a U‐Pb monazite lower‐intercept age of 73.6 ± 2.0 Ma in the same sample, and comparable within errors with a less precise U‐Pb lower‐intercept age of 80 ± 13 Ma obtained in post‐eclogitic titanite from the kyanite‐bearing eclogite. These U‐Pb metamorphic ages are interpreted as dating an amphibolite facies overprint. Protolith U‐Pb zircon ages of 167.1 ± 4.2 Ma and 424.6 ± 5.0 Ma from two eclogite samples reveal that mafic precursors in the Chuacús complex originated in multiple tectonotemporal settings from the Silurian to Jurassic. The integration of petrological and geochronological data suggests that subduction of the continental margin of the North American plate (Chuacús complex) beneath the Greater Antilles arc occurred during an Albian‐Cenomanian pre‐collisional stage, and that a subsequent Campanian collisional stage is probably responsible of the amphibolite facies overprint and late syncollisional exhumation.     

On the ecohydrology of the Yucatan Peninsula: Evapotranspiration and carbon intake dynamics across an eco‐climatic gradient

The hydrology and productivity of the ecosystems of the Yucatan Peninsula (YP) are highly constrained by two factors: (a) the lack of surface drainage networks due to the existence of a highly permeable and connected karstic aquifer roughly the size of the peninsula and (b) a climatic gradient that leads to a transition from seasonally dry deciduous and sub‐deciduous tropical forests, in the north‐western and central parts of the Peninsula, to evergreen forests, in the southern and eastern parts. As a result, surface water fluxes of the YP are restricted to evapotranspiration (ET) that are tightly coupled to ecosystems health and gross primary productivity (GPP). The magnitude and seasonal variation of these fluxes are sensitive to climatic variability and perturbations caused by extreme events such as droughts and tropical storms that are frequent in the YP. In this study, we assess the spatio‐temporal dynamics of ET and GPP above average dry and wet conditions through time series analyses of 15 years of remotely sensed data from both Moderate Resolution Imaging Spectroradiometer and Tropical Rainfall Measuring Mission satellite products. Our results show that ET and GPP follow a regional moisture and temperature gradient that highly controls the distribution of ecosystems within the peninsula. We observe that ET and GPP are in phase with the rainy season in the deciduous forests, but for the evergreen forests, only the GPP is in phase. Additionally, and with the exception of droughts on deciduous ecosystems of the northern part of the YP, the productivity of these ecosystems shows a legacy effect, responding more to a defined trajectory (wetting or drying on the previous years), rather than to punctual extreme climatic events. This has implications on the resilience of these ecosystems to natural perturbations of climate. Comparisons between deciduous and evergreen forest indicate that both types of ecosystems have different plant water use strategies in response to hydrologic variability.     

Joint multifractal description of the relationship between wind patterns and land surface air temperature

Land surface air temperature (SAT), registered at 1.25–2 m above the ground, is influenced by wind patterns. As a consequence, some phenomena such as urban heat islands and the formation of ground-level air pollutants are affected. Detailed understanding of the effects of wind circulations on SAT is convenient to improve the knowledge of these phenomena. Thus, the joint multifractal analysis has been applied to describe time series recorded at Cordoba (southern Spain) from 2001 to 2008 revealing the presence of seasonal patterns related with warm and cold winds blowing from the SW (later spring and summer) and NE (later autumn and winter), respectively, that provoke different heterogeneity in SAT values. In addition, the extreme high SAT values seem to be related with summer SW winds. However, the presence of rare low SAT values produced by NE winds is less relevant.     

Active Deformation along the Southern End of the Tosco-Abreojos Fault System: New Insights from Multibeam Swath Bathymetry

The relative motion of the Pacific plate with respect to the North America plate is partitioned between transcurrent faults located along the western margin of Baja California and transform faults and spreading ridges in the Gulf of California. However, the amount of right lateral offset along the Baja California western margin is still debated. We revisited multibeam swath bathymetry data along the southern end of the Tosco-Abreojos fault system. In this area the depths are less than 1,000?m and allow a finer gridding at 60?m cell spacing. This improved resolution unveils several transcurrent right lateral faults offsetting the seafloor and canyons, which can be used as markers to quantify local offsets. The seafloor of the southern end of the Tosco-Abreojos fault system (south of 24°N) displays NW?CSE elongated bathymetric highs and lows, suggesting a transtensional tectonic regime associated with the formation of pull-apart basins. In such an active tectonic context, submarine canyon networks are unstable. Using the deformation rate inferred from kinematic predictions and pull-apart geometry, we suggest a minimum age for the reorganization of the canyon network.     

The versatility of petrological modeling: Thermobarometry of high-pressure metabasites from the Renge and Sanbagawa belts and phase evolution during warm subduction at Nankai

Petrological modeling is a powerful technique to address different types of geological problems via phase-equilibria predictions at different pressure–temperature-composition conditions. Here, we show the versatility of this technique by (1) performing thermobarometrical calculations using phase equilibrium diagrams to explore the petrological evolution of high-pressure (HP) metabasites from the Renge and Sanbagawa belts, Japan and (2) forward-modeling the mineral–melt evolution of the subducted fresh and altered oceanic crust along the Nankai subduction zone geotherm at the Kii peninsula, Japan. In the first case, we selected three representative samples from these metamorphic belts: a glaucophane eclogite and a garnet glaucophane schist from the Renge belt (Omi area) and a quartz eclogite from the Sanbagawa belt (Besshi area). We calculated the peak metamorphic conditions at ~2.0–2.3 GPa and ~550–630 °C for the HP metabasites from the Renge belt, whereas for the quartz eclogite, the peak equilibrium conditions were calculated at 2.5–2.8 GPa and ~640–750 °C. According to our models, the quartz eclogite experienced partial melting after peak metamorphism. In terms of the petrological evolution of the subducted uppermost portion of the oceanic crust along the warm Nankai geotherm, our models show that fluid release occurs at ~20–60 km, likely promoting high pore-fluid pressure, and thus, seismicity at these depths; dehydration is controlled by chlorite breakdown. Our petrological models predict partial melting at >60 km, mainly driven by phengite and amphibole breakdown. According to our models, the melt proportion is relatively small, suggesting that slab anatexis is not an efficient mechanism for generating voluminous magmatism at these conditions. Modeled melt compositions correspond to high-SiO₂ adakites; these are similar to compositions found in the Daisen and Sambe volcanoes, in southwest Japan, suggesting that the modeled melts may serve as an analog to explain adakite petrogenesis.