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.     

Effect of macrofaunal bioturbation on bacterial distribution in marine sandy sediments, with special reference to sulphur-oxidising bacteria

We have studied the impact of the bioturbating macrofauna, in particular the lugworm Arenicola marina and the bivalve Cerastoderma edule, on abundances and distribution patterns of total bacteria and of bacteria of selected functional groups in sandy intertidal sediments. The selected groups comprised the colourless sulphur-oxidising bacteria and the anoxygenic phototrophic bacteria, which are expected to occupy small zones at the oxygen–sulphide interface in stable (non-bioturbated) sediments. The presence of a wooden wreck buried in the sediment at 10 cm depth within a large area of intertidal sand flat colonised by lugworms provided a unique opportunity to confront field observations with laboratory simulations. The site with the wooden wreck, which was used as control site, was devoid of both A. marina and C. edule, while the composition of the rest of the zoobenthic community was rather similar to that of the surrounding area. In the field, the density of total bacteria was approximately one order of magnitude higher in the control site than in the natural (bioturbated) site. This can be explained by the higher contents of silt and clay particles (higher surface-area/volume ratio) and higher total organic-carbon contents found at the control site. It appears that the presence of macrofauna affects sedimentation processes, which indirectly influence bacterial dynamics. Samples from the control site have been incubated in the laboratory with A. marina and C. edule added (bioturbated core), while an unamended core served as a control. The laboratory experiments contrasted with the field observations, because it was found that total bacteria were actually higher in the deeper layers of the bioturbated core. Moreover, the populations were more homogeneous (less stratified) and colourless sulphur bacteria were on average less numerous in the bioturbated core. In general, laboratory incubations resulted in a decrease of total bacteria with a concomitant increase of colourless and phototrophic sulphur-oxidising bacteria and thus in modifications of the bacterial community structure. Hence, our results demonstrate that care must be taken in extrapolating results from laboratory experiments (e.g. mesocosm research) to field situations.     

The Manila clam population in Arcachon Bay (SW France): Can it be kept sustainable?

The venerid clam Ruditapes philippinarum is the most prominent suspension-feeding bivalve inhabiting muddy intertidal seagrass beds in Arcachon Bay (SW France). It is exploited by fishermen, and Arcachon Bay ranks number one in France in terms of production and total biomass of this species. Previous studies revealed a decrease in the standing stock of R. philippinarum since 2003 and unbalanced length–frequency distributions with a lack of juveniles and of adults > 40 mm. Consequently, the population dynamics of this bivalve were studied at four intertidal sites and one oceanic site in Arcachon Bay. As clam size structure did not allow classical dynamics computations, field monitoring was coupled with field experiments (tagging–recapture) over two years. Monitoring of condition index and gonadal maturation stages highlighted a high variability in spawning number and intensity between sites. Recruitment events in the exploited area varied spatially but with uniformly low values. Von Bertalanffy Growth Function (VBGF) parameters (K, L_∞) were determined using Appeldoorn and ELEFAN methods. In the exploited sites in the inner lagoon, K was relatively high (mean = 0.72 yr^? 1) but L_∞ was low (mean = 41.1 mm) resulting in a moderate growth performance index (Φ′ = 2.99). Growth parameters were not correlated with immersion time and L_∞ was different between sites. Comparison of mortality coefficients (Z) between cage experiments and field monitoring suggested that fishing accounts for 65–75% of total adult mortality. Low recruitment, a low growth rate and a normal mortality rate led to low somatic production (4.1 and 8.7 g Shell-Free Dry Weight (SFDW) m^? 2 yr^? 1) and an annual P/B ratio from 0.44 to 0.92 yr^? 1. Under current conditions, the possibility of a sustainable population in Arcachon Bay will strongly depend on recruitment success and fishing management.     

Simulation of intrinsic random fields of order <Emphasis Type="Italic">k</Emphasis> with a continuous spectral algorithm

Intrinsic random fields of order k, defined as random fields whose high-order increments (generalized increments of order k) are second-order stationary, are used in spatial statistics to model regionalized variables exhibiting spatial trends, a feature that is common in earth and environmental sciences applications. A continuous spectral algorithm is proposed to simulate such random fields in a d-dimensional Euclidean space, with given generalized covariance structure and with Gaussian generalized increments of order k. The only condition needed to run the algorithm is to know the spectral measure associated with the generalized covariance function (case of a scalar random field) or with the matrix of generalized direct and cross-covariances (case of a vector random field). The algorithm is applied to synthetic examples to simulate intrinsic random fields with power generalized direct and cross-covariances, as well as an intrinsic random field with power and spline generalized direct covariances and Matérn generalized cross-covariance.     

Admissible nested covariance models over spheres cross time

Nested covariance models, defined as linear combinations of basic covariance functions, are very popular in many branches of applied statistics, and in particular in geostatistics. A notorious limit of nested models is that the constants in the linear combination are bound to be nonnegative in order to preserve positive definiteness (admissibility). This paper studies nested models on d-dimensional spheres and spheres cross time. We show the exact interval of admissibility for the constants involved in the linear combinations. In particular, we show that at least one constant can be negative. One of the implications is that one can obtain a nested model attaining negative correlations. We provide characterization theorems for arbitrary linear combinations as well as for nonconvex combinations involving two covariance functions. We illustrate our findings through several examples involving nonconvex combinations of well-known parametric families of covariance functions.     

Western boundary upwelling dynamics off Oman

Despite its climatic and ecosystemic significance, the coastal upwelling that takes place off Oman is not well understood. A primitive-equation, regional model forced by climatological wind stress is used to investigate its dynamics and to compare it with the better-known Eastern Boundary Upwellings (EBUs). The solution compares favorably with existing observations, simulating well the seasonal cycles of thermal structure, surface circulation (mean and turbulent), and sea-surface temperature (SST). There is a 1.5-month lag between the maximum of the upwelling-favorable wind-stress-curl forcing and the oceanic response (minima in sea-surface height and SST), which we attribute to onshore-propagating Rossby waves. A southwestward-flowing undercurrent (opposite to the direction of the near-surface flow) is also simulated with a core depth of 1000 m, much deeper than found in EBUs (150–200 m). An EKE budget reveals that, in contrast to EBUs, the upwelling jet is more prone to barotropic than baroclinic instability and the contribution of locally-generated instabilities to EKE is higher by an order of magnitude. Advection and redistribution of EKE by standing mesoscale features also play a significant role in EKE budget.     

Assessing the accuracy of sequential gaussian simulation through statistical testing

Sequential Gaussian simulation is one of the most widespread algorithms for simulating regionalized variables in the earth sciences. Simplicity and flexibility of this algorithm are the most important reasons that make it popular, but its implementation is highly dependent on a screen effect approximation that allows users to use a moving neighborhood instead of a unique neighborhood. Because of this, the size of the moving neighborhood the number of conditioning data and the size of variogram range are important in the simulation process and should be chosen carefully. In this work, different synthetic and real case studies are presented to show the effect of the neighborhood size the number of conditioning data and the size of variogram range on the simulation result, with respect to the reproduction of the model first and second-order parameters. Results indicate that, in both conditional and non-conditional simulation cases, using a neighborhood with <50 conditioning data may lead to an inaccurate reproduction of the model statistics, and some cases require considering more than 200 conditioning data. It also can be understood from the result of example 3 that when the variogram range is beg compared to the simulation domain determination of inaccurate simulation program is harder.     

Late Quaternary biotic and abiotic controls on long-term sediment flux in a northern Australian tropical river system

The modern distribution of monsoonal rainforest in the Australian tropics is patchy and is mainly associated with river corridors and groundwater springs, which indicates a strong dependence on hydrologic and geomorphic conditions. While their present distribution is well known, very little data exists on past spatial and temporal dynamics of these ecosystems, or their medium- to longer-term controls. Factors such as (i) fire frequency and type, and/or (ii) hydroclimatic conditions (e.g. droughts) have been proposed to control riverine corridor rainforest extent. Recent observations, however, also suggest an additional (iii) geomorphic control induced by alluvial knickpoint migration. Sediment sequences provide valuable archives for the reconstruction of longer-term (a) floodplain sedimentary dynamics, (b) local vegetation history, and (c) catchment-wide fire histories. This study investigates such a sediment sequence at Wangi Creek, and shows that a phase of aggradation, lasting ~4000 years, was recently disrupted by channel incision and floodplain erosion. The aggradational phase is characterized by sand deposition with average vertical floodplain accretion rates of 0.8 cm/yr and includes phases of soil development. The recent incisional phase has changed hydro-geomorphic conditions and caused widespread degradation of vegetation, erosion and lowering of the macro-channel surface. While there is no evidence in our data for an erosional event of similar magnitude since the onset of late Holocene floodplain aggradation, Wangi Creek experienced significant erosion and incision immediately before ~4000 years, providing the first evidence for a tropical cut-and-fill river system. We hence argue that phases of aggradation mainly controlled by biotic processes alternate and depend on feedbacks with incision phases controlled mainly by abiotic processes. The results show that eco-hydro-geomorphic feedbacks may play a crucial role in the medium- to longer-term history of tropical fluvial systems and need to be considered when interpreting fluvial archives with regards to climate, fire or human induced change. © 2019 John Wiley & Sons, Ltd.     

Low densities of sea urchins influence the structure of algal assemblages in the western Mediterranean

Numerous studies of interactions between urchins and algae in temperate areas have shown an important structuring effect of sea urchin populations. These studies focused almost wholly on the effect of high urchin densities on laminarian forests. In contrast, algal communities below 5–6 m depth in the northwestern Mediterranean are characterised by low sea urchin densities (<5 ind m⁻²) and the absence of laminarian forests. No previous research has addressed sea urchin/algal interactions in this type of community. To determine the effect of the most abundant echinoid species, Paracentrotus lividus, on well-established algal communities in this area, we performed a removal–reintroduction experiment in rocky patches located between 13 and 16 m depth in the northwestern Mediterranean, where sea urchin densities ranged between 0.9 and 3.4 ind m⁻². After 6 months, the cover of non-crustose algae was significantly higher in the plots from which sea urchins had been removed than in control plots (84 vs 67% cover). These removal plots reverted to their original state upon reintroduction of sea urchins. The non-crustose algae consisted of turfing and frondose forms, with the former representing some 70% of the non-crustose algal cover. Change in the cover of turfing algae was responsible for the significant increase in algal development in the sea urchin removal plots. The response of frondose algae to the treatment varied between algal species. It is concluded that grazing by P. lividus exerts a significant effect on habitat structure, even in communities with low sea urchin densities, such as those found in vast areas of the Mediterranean sublittoral.     

Modelling of drainage dynamics influence on sediment routing system in a fold‐and‐thrust belt

Drainage networks link erosional landscapes and sedimentary basins in a source‐to‐sink system, controlling the spatial and temporal distribution of sediment flux at the outlets. Variations of accumulation rates in a sedimentary basin have been classically interpreted as changes in erosion rates driven by tectonics and/or climate. We studied the interactions between deformation, rainfall rate and the intrinsic dynamics of drainage basins in an experimental fold‐and‐thrust belt subjected to erosion and sedimentation under constant rainfall and shortening rates. The emergence of thrust sheets at the front of a prism may divert antecedent transverse channels (perpendicular to the structural grain) leading to the formation of longitudinal reaches, later uplifted and incorporated in the prism by the ongoing deformation. In the experiments, transverse incisions appear in the external slopes of the emerging thrust sheets. Headward erosion in these transverse channels results in divide migration and capture of the uplifted longitudinal channels located in the inner parts of the prism, leading to drainage network reorganization and modification of the sediment routing system. We show that the rate of drainage reorganization increases with the rainfall rate. It also increases in a nonlinear way with the rate of uplift. We explain this behaviour by an exponent > 1 on the slope variable in the framework of the stream power erosion model. Our results confirm the view that early longitudinal‐dominated networks are progressively replaced by transverse‐dominated rivers during mountain building. We show that drainage network dynamics modulate the distribution of sedimentary fluxes at the outlets of experimental wedges. We propose that under constant shortening and rainfall rates the drainage network reorganization can also modulate the composition and the spatial distribution of clastic fluxes in foreland basins.