Application of a mixing-ratios based formulation to model mixing-driven dissolution experiments

We address the question of how one can combine theoretical and numerical modeling approaches with limited measurements from laboratory flow cell experiments to realistically quantify salient features of complex mixing-driven multicomponent reactive transport problems in porous media. Flow cells are commonly used to examine processes affecting reactive transport through porous media, under controlled conditions. An advantage of flow cells is their suitability for relatively fast and reliable experiments, although measuring spatial distributions of a state variable within the cell is often difficult. In general, fluid is sampled only at the flow cell outlet, and concentration measurements are usually interpreted in terms of integrated reaction rates. In reactive transport problems, however, the spatial distribution of the reaction rates within the cell might be more important than the bulk integrated value. Recent advances in theoretical and numerical modeling of complex reactive transport problems [De Simoni M, Carrera J, Sanchez-Vila X, Guadagnini A. A procedure for the solution of multicomponent reactive transport problems. Water Resour Res 2005;41:W11410. doi: 10.1029/2005WR004056, De Simoni M, Sanchez-Vila X, Carrera J, Saaltink MW. A mixing ratios-based formulation for multicomponent reactive transport. Water Resour Res 2007;43:W07419. doi: 10.1029/2006WR005256] result in a methodology conducive to a simple exact expression for the space–time distribution of reaction rates in the presence of homogeneous or heterogeneous reactions in chemical equilibrium. The key points of the methodology are that a general reactive transport problem, involving a relatively high number of chemical species, can be formulated in terms of a set of decoupled partial differential equations, and the amount of reactants evolving into products depends on the rate at which solutions mix. The main objective of the current study is to show how this methodology can be used in conjunction with laboratory experiments to properly describe the key processes that occur in a complex, geochemically-active system under chemical equilibrium conditions. We model three CaCO₃ dissolution experiments reported in Singurindy et al. [Singurindy O, Berkowitz B, Lowell RP. Carbonate dissolution and precipitation in coastal environments: Laboratory analysis and theoretical consideration. Water Resour Res 2004;40:W04401. doi: 10.1029/2003WR002651, Singurindy O, Berkowitz B, Lowell RP. Correction to Carbonate dissolution and precipitation in coastal environments: laboratory analysis and theoretical consideration. Water Resour Res 2005;41:W11701. doi: 10.1029/2005WR004433], in which saltwater and freshwater were mixed in different proportions. The integrated reaction rate within the cell estimated from the experiments are modeled independently by means of (a) a state-of-the-art reactive transport code, and (b) the uncoupled methodology of [12, 13], both of which use dispersivity as a single, adjustable parameter. The good agreement between the results from both methodologies demonstrates the feasibility of using simple solutions to design and analyze laboratory experiments involving complex geochemical problems.     

Understanding fracture distribution and its relation to knickpoint evolution in the Rio Icacos watershed (Luquillo Critical Zone Observatory,Puerto Rico) using landscape-scale hydrogeophysics

The Rio Icacos watershed in the Luquillo Mountains (Puerto Rico) is unique due to its extremely rapid weathering rates. The watershed is incised into a quartz diorite that has developed a large knickzone defining the river profile. Regolith thickness within the watershed generally decreases from 20 to 30 m at the ridges to several meters in the quartz diorite-dominated valley to tens of centimeters near the major river knickpoint, as determined from previous studies. Above the knickzone, we observe spheroidal corestones, but below this weathering is much less apparent. Measured erosion rates from previous studies are also high in the knickzone compared with upper elevations within the river profile. A suite of near-surface geophysical methods (i.e. ground penetrating radar and terrain conductivity) capable of fast data acquisition in rugged landscapes, was deployed at kilometer scales to characterize critical zone structure. Concentrations of chaotic ground penetrating radar (GPR) reflections and diffraction hyperbolas with low electrical conductivity were observed in vertical zones that outcrop at the land surface as areas of intense fracturing and spheroidally weathered corestones. The width of these fractured and weathered zones showed an increase with proximity to the knickpoint, and was attributed to dilation of these sub-vertical fractures near the knickpoint, as postulated theoretically by a stress model calculated for the topographic variability across the knickzone in the Rio Icacos, and that shows a release of compressive stress near the knickpoint. We hypothesize that erosion rates increase in the knickzone because of this inferred dilation of fractures. Specifically, opened fractures could enhance access of water and in turn promote spalling, erosion, and spheroidal weathering. This study shows that ground-based hydrogeophysical methods used at the landscape-scale (traditionally applied at smaller scales) can be used to explore critical zone architecture at the scales needed to explain the extreme variability in erosion rates across river profiles. © 2018 John Wiley & Sons, Ltd.     

Nankai Trough and Zenisu Ridge: a deep-sea submersible survey

Eight submersible dives between 3000 and 4200 m water depth were made off southern Japan in the eastern Nankai subduction zone. Benthic communities associated with chemosynthetic processes were discovered along the 800 m wide active tectonic zone, at the toe of the accretionary prism. A benthic community was also discovered along a zone of active compression, at the foot of Zenisu Ridge, 30 km south of Nankai Trough. Temperature measurements within the sediments below the benthic communities confirm that upward motion of interstitial water occurs there. Studies of water samples indicate advection of methane and light hydrocarbons. Specimens of the benthic community have been shown to have included in their shells carbonate resulting from methane consumption. Thus the benthic communities are related to overpressure-driven fluid advection along tectonic zones with active surface deformation. A 300 m high active scarp at the toe of the accretionary prism is related to relative motion in a 280° direction which is close to the 305° average direction of subduction in this area. The dives establish further that compressive deformation is presently occurring at the foot of Zenisu Ridge. The previous interpretation of the Zenisu Ridge as a zone of recent north-south intraplate shortening, 40 km south of the Nankai Trench, is confirmed. We conclude that tectonic evolution might well lead to future detachment of the Zenisu Ridge and overthrusting of this large piece of oceanic crust over the continental margin. Such a process might be an efficient one to emplace ophiolites over continents.     

Dynamics of intrathermocline vortices in a gyre flow over a seamount chain

The interaction of meddies with a complex distribution of seamounts is studied in a three-layer quasi-geostrophic model on the f-plane. This study aims at understanding if and how this seamount chain can represent a barrier to the propagation of these eddies and how it can be involved in their decay. The eddies are idealized as vortex patches in the middle layer, interacting with a regional cyclonic current and with ten idealized seamounts. The numerical code is based on the contour surgery technique. The initial position, radius, shape, number and polarity of the eddies are varied. The main results are the following: (1) Though they do not describe the unsteady flow, the streamlines of the regional and topographic flow provide a useful estimate of the vortex trajectories, in particular towards the major seamounts, where stronger velocity shears are expected. (2) The tallest and widest seamounts which have the largest vorticity reservoir are able to considerably erode the vortices, but also to draw anticyclones towards the seamount top. The ability of narrower seamounts to erode vortices is related to their multiplicity. (3) Only 1/3 of the anticyclones with about 30-km radius reach the southern boundary of the seamount chain, and their erosion is larger than 50 %. The other anticyclones are either completely eroded or trapped over a wide seamount top. Cyclones are less affected by seamounts because they oppose the topographic draft towards the seamount top and they drift along the side of the seamount. (4) Large vortices resist topographic erosion more efficiently. The rate of erosion grows from a few percent to about 35–50 % as the vortex radius decreases from about 60 to 30 km. Small cyclones are not eroded, contrary to small anticyclones (which completely decay), in relation with the different trajectories of these eddies in the vicinity of the seamounts. (5) The detailed vortex shape does not appear critical for their evolution, if they are close enough to the seamount chain initially. The interaction between a group of vortices initially north of the seamount chain can modify their trajectory to such an extent that they finally avoid collision with seamounts. (6) Finally, meddy trajectories across the Horseshoe Seamounts (data from the AMUSE experiment) show qualitative similarity with the vortex paths in the model. Several events of vortex decay also occur at comparable locations (in particular over the wide and tall seamounts) in the model and observations.     

Seagrass colonization: Knock-on effects on zoobenthic community,populations and individual health

This study provided evidence that Zostera noltii presence affects macrofauna community structure independently from median sediment grain-size and that the notion of ecosystem health is rather subjective: in the present case, we recorded “good health” in terms of seagrass development, “no impact” in terms of macrobenthic biotic indices and “negative effect” for a given key-population. The occurrence and development of a Z. noltii seagrass bed was surveyed at Banc d’Arguin, Arcachon Bay (France), to estimate the modification of the macrozoobenthic community and of the dynamics of a key-population for the local ecosystem, – the cockle Cerastoderma edule. Even though median grain-size of the sediment decreased only at the very end of the survey, i.e. when seagrass totally invaded the area, most of the macrofauna community characteristics (such as abundance and biomass) increased as soon as Z. noltii patches appeared. The structure of the macrofauna community also immediately diverged between sand and seagrass habitats, without however modifying the tested biotic indices (BENTIX, BOPA, AMBI). The health of the cockle population (growth, abundance, recruitment) was impacted by seagrass development. Related parasite communities slowly diverged between habitats, with more parasites in the cockles from seagrass areas. However, the number of parasites per cockle was always insufficient to alter cockle fitness.     

Information infrastructure for Australia’s Integrated Marine Observing System

Australia’s Integrated Marine Observing System (IMOS, imos.org.au) is research infrastructure to establish an enduring observing program for Australian oceanic waters and shelf seas. The observations cover physical, biological, and chemical variables to address themes of multi-decadal ocean change, climate variability and weather extremes, boundary currents and inter-basin flows, continental shelf processes and ecosystem responses.IMOS observations are collected by national facilities based on various platform types and operated by partner institutions around the country. In this paper we describe the infrastructure and workflows developed to manage and distribute the data to the public. We highlight the existing standards and open-source software we have adopted, and the contributions we have made. To demonstrate the value of this infrastructure we provide some illustrations of use and uptake.All IMOS data are freely and openly available to the public via the Ocean Portal (https://imos.aodn.org.au). All IMOS-developed software is open-source and accessible at https://github.com/aodn.     

Ongoing research on seismic safety assessment

Research on seismic safety assessment has been the centre of great interest among the scientific community in recent years. Although the devastating impact of earthquakes on current society should be incentive enough to increase research, the development of more realistic mechanical behaviour models and the continuous enhancement of computation capabilities are paramount factors contributing a great deal to the increase of such interest. In this context, three research areas can be identified as currently leading to important developments: code related research, especially in Europe where new design codes are in the implementation process; risk analysis, namely concerning the definition of methodologies for safety assessment that involve the evaluation of the failure probability and could be included in future generations of codes; and experimental characterization of constitutive laws which provides support for the development and calibration of accurate and realistic numerical models for seismic analysis and for the adequate characterization of limit state capacities. The paper will present some of the current scientific research trends on these three seismic safety assessment related topics. Studies addressing the seismic safety assessment of structures, either probabilistically or based on code provisions, and that consider realistic nonlinear mechanical behaviour models will be focussed. Reference will also be made to experimental research on the seismic behaviour of structural elements, emphasizing its crucial role to support the development of numerical models simulating the effects of different seismic strengthening techniques. Finally, given the development of studies leading to new trends and perspectives for performance based earthquake engineering, a possible scenario for seismic design in the future is presented, emphasizing the key issues for its implementation.     

Organic matter in sediments of canyons and open slopes of the Portuguese,Catalan, Southern Adriatic and Cretan Sea margins

We describe the quantitative and compositional (phytopigment, protein, carbohydrate and lipid) patterns of sedimentary organic matter along bathymetric gradients in seven submarine canyons and adjacent open slopes located at four European regions: one along the NE Atlantic and three along the Mediterranean continental margins. The investigated areas are distributed along a putative longitudinal gradient of decreasing primary production from the Portuguese (northeastern Atlantic Ocean), to the Catalan (western Mediterranean Sea), Southern Adriatic (central Mediterranean Sea) and Southern Cretan (eastern Mediterranean Sea) margins. Sediment concentrations of organic matter differed significantly between the Portuguese margin and the Mediterranean regions and also from one study area to the other within the Mediterranean Sea. Differences in quantity and composition of sediment organic matter between canyons and open slopes were limited and significant only in the eutrophic Portuguese margin, where the differences were as large as those observed between regions (i.e. at the mesoscale). These results suggest that the overall trophic status of deep margin sediments is controlled mostly by the primary productivity of the overlying waters rather than by the local topography. Moreover, we also report that the quantity and nutritional quality of sediment organic matter in canyons and adjacent open slopes do not show any consistent depth-related pattern. Only the Nazaré and Cascais canyons in the Portuguese margin, at depths deeper than 500 m, displayed a significant accumulation of labile organic matter. The results of our study underline the need of further investigations of deep margins through sampling strategies accounting for adequate temporal and spatial scales of variability.     

Microstructure measurements in natural waters: Methodology and applications

Modern approaches to microstructure data processing, including wavelet denoising, are discussed. The wavelet procedure is applied to small-scale shear signals before estimating the dissipation rate ε and to the temperature/density profiles used to calculate Thorpe scales. Microstructure data obtained on the Mediterranean shelf of Catalonia are used to illustrate various approaches to the Thorpe displacement calculations. It is suggested that the Weibull probability function is an appropriate model for the Thorpe scale distribution. Microstructure measurements from the upper layer of the Boadella reservoir (Catalonia, Spain) support this finding.A new analytical approximation for the 1D Panchev–Kesich spectrum is deduced and the results of ε computation are compared with spectral fitting by the widely used Nasmyth spectrum. Applying the Kraichnan spectral model to compute ε from temperature spectra in the convective-viscous sub-range is examined as an alternative to the Batchelor spectrum. Microstructure measurements taken in Lake Banyoles (Catalonia, Spain) and in the North Atlantic were used for spectral calculations.Statistical analysis of eddy K_b and thermal K_θ diffusivities measured on a shallow shelf of the Black Sea shows the importance of process-orientated domain averaging of the diffusivities in obtaining good correspondence between K_b and K_θ in active turbulent regions. In weakly turbulent, stratified interior layers, the averaged K_b and K_θ differ significantly, which may point to the inapplicability of isotropic formulae used for ε and temperature dissipation χ_θ estimates, as well as to a dependence of the mixing efficiency γ on the Richardson number or in some cases on regions of fossil turbulence.