A real‐space cellular automaton laboratory

Geomorphic investigations may benefit from computer modelling approaches that rely entirely on self‐organization principles. In the vast majority of numerical models, instead, points in space are characterized by a variety of physical variables (e.g. sediment transport rate, velocity, temperature) recalculated over time according to some predetermined set of laws. However, there is not always a satisfactory theoretical framework from which we can quantify the overall dynamics of the system. For these reasons, we prefer to concentrate on interaction patterns using a basic cellular automaton modelling framework. Here we present the Real‐Space Cellular Automaton Laboratory (ReSCAL), a powerful and versatile generator of 3D stochastic models. The objective of this software suite, released under a GNU licence, is to develop interdisciplinary research collaboration to investigate the dynamics of complex systems. The models in ReSCAL are essentially constructed from a small number of discrete states distributed on a cellular grid. An elementary cell is a real‐space representation of the physical environment and pairs of nearest‐neighbour cells are called doublets. Each individual physical process is associated with a set of doublet transitions and characteristic transition rates. Using a modular approach, we can simulate and combine a wide range of physical processes. We then describe different ingredients of ReSCAL leading to applications in geomorphology: dune morphodynamics and landscape evolution. We also discuss how ReSCAL can be applied and developed across many disciplines in natural and human sciences. Copyright © 2013 John Wiley & Sons, Ltd.     

ERORUN-STAFOR: A collaborative observatory for the multidisciplinary study of the critical zone processes in a tropical volcanic watershed including a Tropical Montane Cloud Forest

Tropical volcanic islands are biodiversity hotspots where the Critical Zone (CZ) still remains poorly studied. In such steep topographic environments associated with extreme climatic events (cyclones), deployment and maintenance of monitoring equipment is highly challenging. While a few Critical Zone Observatories (CZOS) are located in tropical volcanic regions, none of them includes a Tropical Montane Cloud Forest (TMCF) at the watershed scale. We present here the dataset of the first observatory from the French network of critical zone observatories (OZCAR) located in an insular tropical and volcanic context, integrating a ‘Tropical Montane Cloud Forest’: The ERORUN-STAFOR observatory. This collaborative observatory is located in the northern part of La Réunion island (Indian Ocean) within the 45.0 km² watershed of Rivière des Pluies (i.e., Rainfall river) which hosts the TMCF of Plaines des Fougères, one of the best preserved natural habitats in La Réunion Island. Since 2014, the ERORUN-STAFOR monitoring in collaboration with local partners collected a multidisciplinary dataset with a constant improvement of the instrumentation over time. At the watershed scale and in its vicinity, the ERORUN-STAFOR observatory includes 10 measurement stations covering the upstream, midstream and downstream part of the watershed. The stations record a total of 48 different variables through continuous (sensors) or periodic (sampling) monitoring. The dataset consists of continuous time series variables related to (i) meteorology, including precipitation, air temperature, relative humidity, wind speed and direction, net radiation, atmospheric pressure, cloud water flux, irradiance, leaf wetness and soil temperature, (ii) hydrology, including water level and temperature, discharge and electrical conductivity (EC) of stream, (iii) hydrogeology, including (ground)water level, water temperature and EC in two piezometers and one horizontally drilled groundwater gallery completed by soil moisture measurements under the canopy. The dataset is completed by periodic time series variables related to (iv) hydrogeochemistry, including field parameters and water analysis results. The periodic sampling survey provides chemical and isotopic compositions of rainfall, groundwater, and stream water at different locations of this watershed. The ERORUN-STAFOR monitoring dataset extends from 2014 to 2022 with an acquisition frequency from 10 min to hourly for the sensor variables and from weekly to monthly frequency for the sampling. Despite the frequent maintenance of the monitoring sites, several data gaps exist due to the remote location of some sites and instrument destruction by cyclones. Preliminary results show that the Rivière des Pluies watershed is characterized by high annual precipitation (>3000 mm y⁻¹) and a fast hydrologic response to precipitation (≈2 h basin lag time). The long-term evolution of the deep groundwater recharge is mainly driven by the occurrence of cyclone events with a seasonal groundwater response. The water chemical results support existing hydrogeological conceptual models suggesting a deep infiltration of the upstream infiltrated rainfall. The TMCF of Plaine des Fougères shows a high water storage capacity (>2000% for the Bryophytes) that makes this one a significant input of water to groundwater recharge which still needs to be quantified. This observatory is a unique research site in an insular volcanic tropical environment offering three windows of observation for the study of critical zone processes through upstream-midstream-downstream measurements sites. This high-resolution dataset is valuable to assess the response of volcanic tropical watersheds and aquifers at both event and long-term scales (i.e., global change). It will also provide insights in the hydrogeological conceptual model of volcanic islands, including the significant role of the TMCFs in the recharge processes as well as the watershed hydrosedimentary responses to extreme climatic events and their respective evolution under changing climatic conditions. All data sets are available at https://doi.org/10.5281/zenodo.7983138 .     

Sand spit dynamics in a large tidal-range environment: Insight from multiple LiDAR,UAV and hydrodynamic measurements on multiple spit hook development,breaching, reconstruction,and shoreline changes

Sand spits with distal hooks have been well documented from coasts with low to moderate tidal ranges, unlike high tidal-range environments. Datasets from 15 LiDAR and 3 UAV surveys between 2009 and 2019 on the Agon spit in Normandy (France), a setting with one of the largest tidal ranges in the world (mean spring tidal range: 11 m), combined with in-situ hydrodynamic records between 2013 and 2017, highlight a three-stage pattern of spit hook evolution. Stage 1 (2009–2013) commenced with the onshore migration and attachment of a swash bar, followed by persistent spit accretion updrift of the bar and erosion downdrift because of the slow speed of bar migration in this large tidal-range environment. In stage 2 (2013–2016), three overwash events and a 220 m-wide breach culminating in the total destruction of the spit during winter 2015–2016 involved the landward mobilization of thousands of cubic metres of sand. These events occurred during short durations (a few hours) when spring high tides coincided with relatively energetic waves, underscoring the importance of storms in rapid spit morphological change. Strong spring tidal currents maintained the breach. Stage 3 (2016–2019) has involved new hook construction through welding of a swash bar and spit longshore extension, highlighting the resilience of the spit over the 10-year period, and involving a positive sediment balance of 244 000 m³. The three stages bring out, by virtue of the temporal density of LiDAR and UAV data used, a high detail of spit evolution relative to earlier studies in this macrotidal setting. The large tidal range strongly modulates the role of waves and wave-generated longshore currents, the main process drivers of spit evolution, by favouring long periods of inertia in the course of the spring–neap tidal cycle, but also brief episodes of significant morphological change when storm waves coincide with spring high tides. © 2020 John Wiley & Sons, Ltd.     

The morphology of craters on Mercury: Results from MESSENGER flybys

Topographic data measured from the Mercury Laser Altimeter (MLA) and the Mercury Dual Imaging System (MDIS) aboard the MESSENGER spacecraft were used for investigations of the relationship between depth and diameter for impact craters on Mercury. Results using data from the MESSENGER flybys of the innermost planet indicate that most of the craters measured with MLA are shallower than those previously measured by using Mariner 10 images. MDIS images of these same MLA-measured craters show that they have been modified. The use of shadow measurement techniques, which were found to be accurate relative to the MLA results, indicate that both small bowl-shaped and large complex craters that are fresh possess depth-to-diameter ratios that are in good agreement with those measured from Mariner 10 images. The preliminary data also show that the depths of modified craters are shallower relative to fresh ones, and might provide quantitative estimates of crater in-filling by subsequent volcanic or impact processes. The diameter that defines the transition from simple to complex craters on Mercury based on MESSENGER data is consistent with that reported from Mariner 10 data.     

Path planning in uncertain flow fields using ensemble method

An ensemble-based approach is developed to conduct optimal path planning in unsteady ocean currents under uncertainty. We focus our attention on two-dimensional steady and unsteady uncertain flows, and adopt a sampling methodology that is well suited to operational forecasts, where an ensemble of deterministic predictions is used to model and quantify uncertainty. In an operational setting, much about dynamics, topography, and forcing of the ocean environment is uncertain. To address this uncertainty, the flow field is parametrized using a finite number of independent canonical random variables with known densities, and the ensemble is generated by sampling these variables. For each of the resulting realizations of the uncertain current field, we predict the path that minimizes the travel time by solving a boundary value problem (BVP), based on the Pontryagin maximum principle. A family of backward-in-time trajectories starting at the end position is used to generate suitable initial values for the BVP solver. This allows us to examine and analyze the performance of the sampling strategy and to develop insight into extensions dealing with general circulation ocean models. In particular, the ensemble method enables us to perform a statistical analysis of travel times and consequently develop a path planning approach that accounts for these statistics. The proposed methodology is tested for a number of scenarios. We first validate our algorithms by reproducing simple canonical solutions, and then demonstrate our approach in more complex flow fields, including idealized, steady and unsteady double-gyre flows.     

Observations of coastal cliff base waves,sand levels,and cliff top shaking

Concurrent observations of waves at the base of a southern California coastal cliff and seismic cliff motion were used to explore wave–cliff interaction and test proxies for wave forcing on coastal cliffs. Time series of waves and sand levels at the cliff base were extracted from pressure sensor observations programmatically and used to compute various wave impact metrics (e.g. significant cliff base wave height). Wave–cliff interaction was controlled by tide, incident waves, and beach sand levels, and varied from low tides with no wave–cliff impacts, to high tides with continuous wave–cliff interaction. Observed cliff base wave heights differed from standard Normal and Rayleigh distributions. Cliff base wave spectra levels were elevated at sea swell and infragravity frequencies. Coastal cliff top response to wave impacts was characterized using microseismic shaking in a frequency band (20–45 Hz) sensitive to wave breaking and cliff impacts. Response in the 20–45 Hz band was well correlated with wave–cliff impact metrics including cliff base significant wave height and hourly maximum water depth at the cliff base (r² = 0.75). With site‐specific calibration relating wave impacts and shaking, and acceptable anthropogenic (traffic) noise levels, cliff top seismic observations are a viable proxy for cliff base wave conditions. The methods presented here are applicable to other coastal settings and can provide coastal managers with real time coastal conditions. Copyright © 2016 John Wiley & Sons, Ltd.     

Input of particulate heavy metals from rivers and associated sedimentary deposits on the Gulf of Lion continental shelf

Fluxes of the heavy metals chromium (Cr), cobalt (Co), nickel (Ni), copper (Cu), lead (Pb), cadmium (Cd) and zinc (Zn) delivered by rivers to the Gulf of Lion (NW Mediterranean Sea) were estimated over a three year study of the River Rhone and its smaller tributaries. Most of the particulate metal fluxes (80–90%) delivered by these rivers occurred within a very short period of time (less than 12%), a typical trend for the Mediterranean environment, where highly contrasting hydrological regimes were observed over the year. Temporal and spatial variations in the fluxes of these particulate metals were driven by the fluxes in both water discharge and suspended particulate matter load. On the shelf, these particulate metal fluxes, largely arising from the Rhone watershed, were two to ten times more important than those resulting from atmospheric deposition. Co, Cr and Ni in the rivers and on the shelf surface sediments were mainly natural and associated with the finest particles. Cd and Phosphorus appeared to be associated with the silt fraction and to be enriched in the prodelta areas. Pb, Zn and Cu were more closely associated with the organic matter content and also showed enrichment in the organic rich prodeltaic sediments. Anthropogenic influences diminished offshore, except for Pb and Zn which could be supplied from the atmosphere by man-made aerosols. Although most of the metals tended to be enriched in the prodelta areas, these did not constitute a permanent sink due to resuspension processes affecting these shallow depths. A resuspension experiment conducted on sediment cores from the Rhone prodelta demonstrated that metal deposited on the surface layer, especially those associated with the organic matter, may be resuspended; this should be taken into account for a complete understanding of the biogeochemical cycle of these metals.     

TIMS measurements of Ra and Ra in the Gulf of Lion, an attempt to quantify submarine groundwater discharge

Submarine groundwater discharge (SGD) is now recognized as an important pathway for water and chemical species fluxes to the coastal ocean. In order to determinate SGD to the Gulf of Lion (France), we measured the activities of ²²⁶Ra and ²²⁸Ra by thermal ionization mass spectrometry (TIMS) in coastal waters and in the deep aquifer waters of the Rhone deltaic plain after pre-concentration of radium by MnO₂. Compared to conventional counting techniques, TIMS requires lower quantities of water for the analyses, and leads to higher analytical precision. Radium isotopes were thus measured on 0.25–2 L water samples containing as little as 20 fg of ²²⁶Ra and 0.2–0.4 fg of ²²⁸Ra with precision equal to 2%. We demonstrate that coastal surface waters samples are enriched in ²²⁶Ra and ²²⁸Ra compared to the samples further offshore. The high precision radium measurements display a small but significant ²²⁶Ra and ²²⁸Ra enrichment within a strip of circa 30 km from the coast. Radium activities decrease beyond this region, entrained in the northern current along the shelf break or controlled by eddy diffusion. The radium excess in the first 30 km cannot be accounted for by the river nor by the early diagenesis. The primary source of the radium enrichment must therefore be ascribed to the discharge of submarine groundwater. Using a mass-balance model, we estimated the advective fluxes of ²²⁶Ra and ²²⁸Ra through SGD to be 5.2 × 10¹⁰ and 21 × 10¹⁰ dpm/d respectively. The ²²⁶Ra activities measured in the groundwater from the Rhone deltaic plain aquifer are comparable to those from other coastal groundwater studies throughout the world. By contrast, ²²⁸Ra activities are higher by up to one order of magnitude. Taking those groundwater radium activities as typical of the submarine groundwater end-member, a minimum volume of 0.24–4.5 × 10¹⁰ l/d is required to support the excess radium isotopes on the inner shelf. This has to be compared with the average rivers water runoff of 15.4 × 10¹⁰ l/d during the study period (1.6 to 29% of the river flow).