Flexural deformation controls on Late Quaternary sediment dispersal in the Garo-Rajmahal Gap,NW Bengal Basin

Subsurface deformation is a driver for river path selection when deformation rates become comparable to the autogenic mobility rate of rivers. Here we combine geomorphology, soil and sediment facies analyses, and geophysical data of the Late Quaternary sediments of the central Garo-Rajmahal Gap in Northwest Bengal to link subsurface deformation with surface processes. We show variable sedimentation characteristics, from slow rates (<0.8 mm/year) in the Tista megafan at the foot of the Himalaya to nondeposition at the exposed surface of the Barind Tract to the south, enabling the development of mature soils. Combined subsidence in the Tista fan and uplift of the Barind Tract are consistent with a N-S flexural response of the Indian plate to loading of the Himalaya Mountains given a low value of elastic thickness (15–25 km). Provenance analysis based on bulk strontium concentration suggests a dispersal of sediment consistent with this flexural deformation—in particular the abandonment of the Barind Tract by a Pleistocene Brahmaputra River and the current extents of the Tista megafan lobes. Overall, these results highlight the control by deeply rooted deformation patterns on the routing of sediment by large rivers in foreland settings.     

The distribution of chlorophyll-a and dominant planktonic components in the coastal transition zone off Concepción, central Chile, during different oceanographic conditions

The oceanographic setting and the planktonic distribution in the coastal transition zone off Concepción (∼35-38°S, ∼73-77°W), an area characterized by its high biological production, were assessed during two different seasons: austral spring with equatorward upwelling favorable winds and austral winter with predominately northerly winds. Oceanographic and biological data (total chlorophyll-a, particulate organic carbon, microplankton, large mesozooplankton >500 μm as potential consumers of microplankton) were obtained during two cruises (October 1998, July 1999) together with satellite imagery for wind stress, geostrophic flow, surface temperature, and chlorophyll-a data. The physical environment during the spring sampling was typical of the upwelling period in this region, with a well-defined density front in the shelf-break area and high concentrations of surface chlorophyll-a (>5 mg m⁻³) on the shelf over the Itata terrace. During the winter sampling, highly variable though weakly upwelling-favorable winds were observed along with lower surface chlorophyll-a values (<2 mg m⁻³) on the shelf. In the oceanic area (>100 km from the coast), cyclonic and anti-cyclonic eddies were evident in the flow field during both periods, the former coinciding with higher chlorophyll-a contents (∼1 mg m⁻³) than in the surrounding waters. Also, a cold, chlorophyll-a rich filament was well defined during the spring sampling, extending from the shelf out to 350-400 km offshore. Along a cross-shelf transect, the micro- and meso-planktonic assemblages displayed higher coastal abundances during the spring cruise but secondary peaks appeared in the oceanic area during the winter cruise, coinciding with the distribution of the eddies. These results suggest that the mesoscale features in this region, in combination with upwelling, play a role in potentially increasing the biological productivity of the coastal transition zone off Concepción.     

Groundwater flow systems in turbidites of the Northern Apennines (Italy): natural discharge and high speed railway tunnel drainage

Turbidites crop out extensively in the Northern Apennine mountains (Italy). The huge amounts of groundwater drained by tunnels, built for the high speed railway connection between Bologna and Florence, demonstrate the aquifer-like behaviour of these units, up to now considered as aquitards. A conceptual model of groundwater flow systems (GFS) in fractured aquifers of turbidites is proposed, taking into account both system natural state and the perturbation induced by tunnel drainage. Analysis of hydrological data (springs, streams and tunnel discharge), collected over 10 years, was integrated with analysis of hydrochemical and isotopic data and a stream-tunnel tracer test. Hydrologic recession analysis of undisturbed conditions is a key tool in studying turbiditic aquifer hydrogeology, permitting the discrimination of GFS, the estimation of recharge relative to the upstream reach portion and the identification of springs most vulnerable to tunnel drainage impacts. The groundwater budgeting analysis provides evidence that the natural aquifer discharge was stream-focused through GFS, developed downslope or connected to main extensional tectonic lineaments intersecting stream beds; now tunnels drain mainly active recharge groundwater and so cause a relevant stream baseflow deplenishment (approximately two-thirds of the natural value), possibly resulting in adverse effects on local ecosystems.     

World-wide synthetic tide parameters for gravity and vertical and horizontal displacements

The response of the Earth’s crust to the direct effect of lunisolar gravitational forcing is known as the body tide. The body tide is superimposed by surface-loading forces due to the pressure of the periodically varying ocean tide acting on the Earth, called ocean tide loading (OTL). Both body tide and OTL can be decomposed into components of the same frequency known as tidal parameters. However, OTL is more complicated than body tides because of the dynamic effects of the ocean. Estimating OTL requires a model of the ocean tides and knowledge of the elastic properties of the solid Earth. Thus, synthetic tide parameters (amplitude factors and phase leads) have been developed here on a world-wide grid for gravity and positional displacements. The body tide contributions were added to the oceanic contribution to provide the Earth tide response. The accuracy and reliability of the synthetic tidal parameters have been estimated by comparing observed gravity and vertical-displacement tide parameters with those interpolated from our synthetic model, which shows good agreement. Tests also indicate that the synthetic tide parameters provide realistic gravimetric and displacements for practical use in tidal prediction.     

Uplift and subsidence from oblique slip: the Ganos–Marmara bend of the North Anatolian Transform, western Turkey

Bends that locally violate plate-motion-parallel geometry are common structural elements of continental transform faults. We relate the vertical component of crustal motion in the western Marmara Sea region to the NNW-pointing 18° bend on the northern branch of the North Anatolian Fault (NAF-N) between the Ganos segment, which ruptured in 1912, and the central Marmara segment, a seismic gap. Crustal shortening and uplift on the transpressive west side of the bend results in the Ganos Mountain; crustal extension and subsidence on the transtensional east side produce the Tekirdağ Basin. We propose that this vertical component of deformation is controlled by oblique slip on the non-vertical north-dipping Ganos and Tekirdağ segments of the North Anatolian Fault. We compare Holocene with Quaternary structure across the bend using new and recently published data and conclude the following. First, bend-related vertical motion is occurring primarily north of the NAF-N. This suggests that this bend is fixed to the Anatolian side of the fault. Second, current deformation is consistent with an antisymmetric pattern centered at the bend, up on the west and down on the east. Accumulated deformation is shifted to the east along the right-lateral NAF-N, however, leading to locally opposite vertical components of long- and short-term motion. Uplift has started as far west as the landward extension of the Saros trough. Current subsidence is most intense close to the bend and to the Ganos Mountain, while the basin deepens gradually from the bend eastward for 28 km along the fault. The pattern of deformation is time-transgressive if referenced to the material, but is stable if referenced to the bend. The lag between motion and structure implies a 1.1–1.4 Ma age for the basin at current dextral slip rate (2.0–2.5 cm/year). Third, the Tekirdağ is an asymmetric basin progressively tilted down toward the NAF-N, which serves as the border fault. Progressive tilt suggests that the steep northward dip of the fault decreases with depth in a listric geometry at the scale of the upper crust and is consistent with reactivation of Paleogene suture-related thrust faults. Fourth, similar thrust-fault geometry west of the bend can account for the Ganos Mountain anticline/monocline as hanging-wall-block folding and back tilting. Oblique slip on a non-vertical master fault may accommodate transtension and transpression associated with other bends along the NAF and other continental transforms.     

Salinity and periodic inundation controls on the soil‐plant‐atmosphere continuum of gray mangroves

Salinity and periodic inundation are both known to have a major role in shaping the ecohydrology of mangroves through their controls on water uptake, photosynthesis, stomatal conductance, gas exchanges, and nutrient availability. Salinity, in particular, can be considered one of the main abiotic regulating factors for halophytes and salt‐tolerant species, due to its influence on water use patterns and growth rate. Ecohydrological literature has rarely focused on the effects of salinity on plant transpiration, based on the fact that the terrestrial plants mostly thrive in low‐saline, unsaturated soils where the role of osmotic potential can be considered negligible. However, the effect of salinity cannot be neglected in the case of tidal species like mangroves, which have to cope with hyperosmotic conditions and waterlogging. We introduce here a first‐order ecohydrological model of the soil/plant‐atmosphere continuum of Avicennia marina—also known as gray mangrove—a highly salt‐tolerant pioneer species able to adapt to hyperarid intertidal zones and characterized by unique morphological and ecophysiological traits. The A. marina's soil‐plant‐atmosphere continuum takes explicitly into account the role of water head, osmotic water potential, and water salinity in governing plant water fluxes. A. marina's transpiration is thus modeled as a function of salinity based on a simple parameterization of salt exclusion mechanisms at the root level and a modified Jarvis' expression accounting for the effects of salinity on stomatal conductance. Consistently with previous studies investigating the physiology of mangroves in response to different environmental drivers, our results highlight the major influence of salinity on mangrove transpiration when contrasted with other potential stressors such as waterlogging and water stress.     

Segmentation and kinematics of the North America‐Caribbean plate boundary offshore Hispaniola

We explored the submarine portions of the Enriquillo–Plantain Garden Fault zone (EPGFZ) and the Septentrional–Oriente Fault zone (SOFZ) along the Northern Caribbean plate boundary using high‐resolution multibeam echo‐sounding and shallow seismic reflection. The bathymetric data shed light on poorly documented or previously unknown submarine fault zones running over 200 km between Haiti and Jamaica (EPGFZ) and 300 km between the Dominican Republic and Cuba (SOFZ). The primary plate‐boundary structures are a series of strike‐slip fault segments associated with pressure ridges, restraining bends, step overs and dogleg offsets indicating very active tectonics. Several distinct segments 50–100 km long cut across pre‐existing structures inherited from former tectonic regimes or bypass recent morphologies formed under the current strike‐slip regime. Along the most recent trace of the SOFZ, we measured a strike‐slip offset of 16.5 km, which indicates steady activity for the past ~1.8 Ma if its current GPS‐derived motion of 9.8 ± 2 mm a^?1 has remained stable during the entire Quaternary.     

Hierarchical simulation of aquifer heterogeneity: implications of different simulation settings on solute-transport modeling

The fine-scale heterogeneity of porous media affects the large-scale transport of solutes and contaminants in groundwater and it can be reproduced by means of several geostatistical simulation tools. However, including the available geological information in these tools is often cumbersome. A hierarchical simulation procedure based on a binary tree is proposed and tested on two real-world blocks of alluvial sediments, of a few cubic meters volume, that represent small-scale aquifer analogs. The procedure is implemented using the sequential indicator simulation, but it is so general that it can be adapted to various geostatistical simulation tools, improving their capability to incorporate geological information, i.e., the sedimentological and architectural characterization of heterogeneity. When compared with a standard sequential indicator approach on bi-dimensional simulations, in terms of proportions and connectivity indicators, the proposed procedure yields reliable results, closer to the reference observations. Different ensembles of three-dimensional simulations based on different hierarchical sequences are used to perform numerical experiments of conservative solute transport and to obtain ensembles of equivalent pore velocity and dispersion coefficient at the scale length of the blocks (meter). Their statistics are used to estimate the impact of the variability of the transport properties of the simulated blocks on contaminant transport modeled on bigger domains (hectometer). This is investigated with a one-dimensional transport modeling based on the Kolmogorov-Dmitriev theory of branching stochastic processes. Applying the proposed approach with diverse binary trees and different simulation settings provides a great flexibility, which is revealed by the differences in the breakthrough curves.     

High-dimensional geostatistical history matching

Hydrocarbon reservoir modelling and characterisation is a challenging subject within the oil and gas industry due to the lack of well data and the natural heterogeneities of the Earth’s subsurface. Integrating historical production data into the geo-modelling workflow, commonly designated by history matching, allows better reservoir characterisation and the possibility of predicting the reservoir behaviour. We present herein a geostatistical-based multi-objective history matching methodology. It starts with the generation of an initial ensemble of the subsurface petrophysical property of interest through stochastic sequential simulation. Each model is then ranked according the match between its dynamic response, after fluid flow simulation, and the observed available historical production data. This enables building regionalised Pareto fronts and the definition of a large ensemble of optimal subsurface Earth models that fit all the observed production data without compromising the exploration of the uncertainty space. The proposed geostatistical multi-objective history matching technique is successfully implemented in a benchmark synthetic reservoir dataset, the PUNQ-S3, where 12 objectives are targeted.     

An orthotropic interface damage model for simulating drying processes in soils

The study of drying process in soils has received an increased attention in the last few years. This is very complex phenomenon that generally leads to the formation and propagation of desiccation cracks in the soil mass. In recent engineering applications, high aspect ratio elements have proved to be well suited to tackle this type of problem using finite elements. However, the modeling of interfaces between materials with orthotropic properties that generally exist in this type of problem using standard (isotropic) constitutive model is very complex and challenging in terms of the mesh generation, leading to very fine meshes that are intensive CPU demanding. A novel orthotropic interface mechanical model based on damage mechanics and capable of dealing with interfaces between materials in which the strength depends on the direction of analysis is proposed in this paper. The complete mathematical formulation is presented together with the algorithm suggested for its numerical implementation. Some simple yet challenging synthetic benchmarks are analyzed to explore the model capabilities. Laboratory tests using different textures at the contact surface between materials were conducted to evaluate the strengths of the interface in different directions. These experiments were then used to validate the proposed model. Finally, the approach is applied to simulate an actual desiccation test involving an orthotropic contact surface. In all the application cases the performance of the model was very satisfactory.