Methodology of hydrogeological characterization of deep carbonate aquifers as potential reservoirs of groundwater. Case of study: the Jurassic aquifer of <Emphasis Type="Italic">El Maestrazgo</Emphasis> (Castellón,Spain)

A methodology for the characterization of deep carbonate aquifers has been developed and applied to El Maestrazgo Jurassic aquifer in Castellón, Spain. Characterization of these aquifer formations, located at more than 300 m deep, consisted of a previous phase of compilation, analysis and synthesis of the existing information about the area, followed by a coordinated combination of different speciality studies: geology, stratigraphy, structural analysis, hydrogeology, hydrochemistry, geophysics and remote sensing. Geological studies included geological mapping, definition of stratigraphical units and facies and structural analysis. The aim of the hydrogeology study was to define aquifer formations, recharge area, aquifer points inventory and groundwater flow directions for the establishment of piezometric and water quality observation nets. Special techniques were applied, like thermal infrared aerial images and the evaluation of submarine groundwater discharge by means of natural radium isotopes. Hydrochemical techniques, including majority elements characterization and stable isotopes (¹⁸O, ²H and ³H) determination, allowed classifying hydrochemical facies and establishing a renewal pattern for water within the system. Geophysics was useful in determining the aquifer geometry, the features of the basement and the petrophysical characteristics of the geological formations. Preliminary results show an important tectonic complexity and the possibilities for groundwater uses in the area of study.     

Spatial and seasonal patterns of flood change across Brazil

Brazil has some of the largest rivers in the world and has the second greatest flood loss potential among the emergent countries. Despite that, flood studies in this area are still scarce. In this paper, we used flood seasonality and trend analysis at the annual and seasonal scales in order to describe flood regimes and changes across the whole of Brazil in the period 1976–2015. We identified a strong seasonality of floods and a well-defined spatio-temporal pattern for flood occurrence. There are positive trends in the frequency and magnitude of floods in the North, South and parts of Southeast Brazil; and negative trends in the North-east and the remainder of Southeast Brazil. Trends in the magnitude (frequency) were predominant in the winter (summer). Overall, floods are becoming more frequent and intense in Brazilian regions characterized by wet conditions, and less frequent and intense in drier regions.     

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.     

Effect of galvanic distortions on the series and parallel magnetotelluric impedances and comparison with other responses

The series and parallel impedances of the magnetotelluric tensor are appraised in relation to their relative immunity to galvanic electric distortions. The distorted responses are modeled using the Groom-Bailey decomposition of the tensor in terms of twist, shear, statics and strike direction. These four parameters and the undistorted responses are normally considered as unknowns, and are obtained from field data through the solution of an inverse problem. In the present work we use the decomposition as a forward model to simulate distorted sounding curves. Starting with undistorted 2-d te and tm responses, the tensor is distorted by assuming arbitrary values of twist, shear, static and strike direction. By default, both series and parallel responses are immune to the strike direction because they are invariants under rotation. In addition, series responses are immune to twist and shear and parallel responses only to twist. The dependence of the latter on shear is in the form of a real factor that shifts downwards the amplitude curves. On the other hand, the effect of statics on both series and parallel responses is more complicated than that on the impedance tensor because it cannot be accounted for by a simple shift of the curves. On the whole, there is a positive balance on the part of the series and parallel impedances over the te and tm responses because some of the distortions are filtered out by the invariants. It is shown that invariance is not sufficient to be immune to any of the distortions. The example chosen is Eggers’ eigenvalues, which are immune only to the by-the-fault strike direction. Invariance is not necessary either, as evidenced by the phase tensor, whose elements depend on strike but are immune to all distortions. The derivations are illustrated using soundings from the synthetic coprod2s1 and field-recorded coprod2 and bc87 data sets.     

Anthropogenic influence on trace element geochemistry of healing mud (peloid) from Makirina Cove (Croatia)

Due to their balneotherapeutic features, the organic-rich sediments in Makirina Cove are an important source of healing mud. An environmental geochemistry approach using normalization techniques was applied to evaluate the anthropogenic contribution of trace metals to sediments used as healing mud. Sediment geochemistry was found to be associated with land-use change and storm events, as well as with proximity of a road with heavy traffic in the summer months. Local valley topography preferentially channels lithogenic and pollutant transport to the cove. Concentrations and distribution of trace metals indicate lithogenic (Ni, Cr, Co) and anthropogenic (Pb, Cu, Zn and Se) contributions to the sediments. The calculation of enrichment factors indicates a moderate (EFs between 2–3.5) input of anthropogenic Cu and Pb in surficial sediments to a depth of 10 cm. Patients using the Makirina Cove sediments as healing mud could be to some extent exposed to enhanced uptake of metals from anthropogenic sources via dermal contact.     

The Barents Sea Ice Sheet — A model of its growth and decay during the last ice maximum

On the basis of geomorphological and sedimentological data, we believe that the entire Barents Sea was covered by grounded ice during the last glacial maximum. ¹⁴C dates on shells embedded in tills suggest marine conditions in the Barents Sea as late as 22 ka BP; and models of the deglaciation history based on uplift data from the northern Norwegian coast suggest that significant parts of the Barents Sea Ice Sheet calved off as early as 15 ka BP. The growth of the ice sheet is related to glacioeustatic fall and the exposure of shallow banks in the central Barents Sea, where ice caps may develop and expand to finally coalesce with the expanding ice masses from Svalbard and Fennoscandia.The outlined model for growth and decay of the Barents Sea Ice Sheet suggests a system which developed and existed under periods of maximum climatic deterioration, and where its growth and decay were strongly related to the fall and rise of sea level.     

Cenozoic erosion and the preglacial uplift of the Svalbard–Barents Sea region

The creation of the huge fans observed in the western Barents Sea margin can only be explained by assuming extremely high glacial erosion rates in the Barents Sea area. Glacial processes capable of producing such high erosion rates have been proposed, but require the largest part of the preglacial Barents Sea to be subaerial. To investigate the validity of these proposals we have attempted to reconstruct the western preglacial Barents Sea. Our approach was to combine erosion maps based on prepublished data into a single mean valued erosion map covering the whole western Barents Sea and consequently use it together with a simple Airy isostatic model to obtain a first rough estimate of the preglacial topography and bathymetry of the western Barents Sea margin. The mean valued erosion map presented herein is in good volumetric agreement with the sediments deposited in the western Barents Sea margin areas, and as a direct consequence of the averaging procedures employed in its construction we can safely assume that it is the most reliable erosion map based on the available information. By comparing the preglacial sequences with the glacial sequences in the fans we have concluded that 1/2 to 2/3 of the total Cenozoic erosion was glacial in origin and therefore a rough reconstruction of the preglacial relief of the western Barents Sea could be obtained. The results show a subaerial preglacial Barents Sea. Thus, during interglacials and interstadials the area may have been partly glaciated and intensively eroded up to 1 mm/y, while during relatively brief periods of peak glaciation with grounded ice extending to the shelf edge, sediments have been evacuated and deposited at the margins at high rates. The interplay between erosion and uplift represents a typical chicken and egg problem; initial uplift is followed by intensive glacial erosion, compensated by isostatic uplift, which in turn leads to the maintenance of an elevated, and glaciated, terrain. The information we have on the initial tectonic uplift suggests that the most likely mechanism to cause an uplift of the dimensions and magnitude of the one observed in the Barents Sea is a thermal mechanism.     

Sources and distribution of organic matter in a river-dominated estuary (Winyah Bay, SC, USA)

The sources and distribution of organic matter (OM) in surface waters and sediments from Winyah Bay (South Carolina, USA) were investigated using a variety of analytical techniques, including elemental, stable isotope and organic biomarker analyses. Several locations along the estuary salinity gradient were sampled during four different periods of contrasting river discharge and tidal range. The dissolved organic carbon (DOC) concentrations of surface waters ranged from 7 mg l⁻¹ in the lower bay stations closest to the ocean to 20 mg l⁻¹ in the river and upper bay samples. There was a general linear relationship between DOC concentrations and salinity in three of the four sampling periods. In contrast, particulate organic carbon (POC) concentrations were significantly lower (0.1–3 mg l⁻¹) and showed no relationship with salinity. The high molecular weight dissolved OM (HMW DOM) isolated from selected water samples collected along the bay displayed atomic carbon:nitrogen ratios ([C/N]a) and stable carbon isotopic compositions of organic carbon (δ¹³C_OC) that ranged from 10 to 30 and from −28 to −25‰, respectively. Combined, such compositions indicate that in most HMW DOM samples, the majority of the OM originates from terrigenous sources, with smaller contributions from riverine and estuarine phytoplankton. In contrast, the [C/N]a ratios of particulate OM (POM) samples varied significantly among the collection periods, ranging from low values of 5 to high values of >20. Overall, the trends in [C/N]a ratios indicated that algal sources of POM were most important during the early and late summer, whereas terrigenous sources dominated in the winter and early spring.In Winyah Bay bottom sediments, the concentrations of the mineral-associated OM were positively correlated with sediment surface area. The [C/N]a ratios and δ¹³C_OC compositions of the bulk sedimentary OM ranged from 5 to 45 and from −28 to −23‰, respectively. These compositions were consistent with predominant contributions of terrigenous sources and lesser (but significant) inputs of freshwater, estuarine and marine phytoplankton. The highest terrigenous contents were found in sediments from the river and upper bay sites, with smaller contributions to the lower parts of the estuary. The yields of lignin-derived CuO oxidation products from Winyah Bay sediments indicated that the terrigenous OM in these samples was composed of variable mixtures of relatively fresh vascular plant detritus and moderately altered soil OM. Based on the lignin phenol compositions, most of this material appeared to be derived from angiosperm and gymnosperm vascular plant sources similar to those found in the upland coastal forests in this region. A few samples displayed lignin compositions that suggested a more significant contribution from marsh C₃ grasses. However, there was no evidence of inputs of Spartina alterniflora (a C₄ grass) remains from the salt marshes that surround the lower sections of Winyah Bay.