Understanding river dune splitting through flume experiments and analysis of a dune evolution model

Forecasts of water level during river floods require accurate predictions of the evolution of river dune dimensions, because the hydraulic roughness of the main channel is largely determined by the bed morphology. River dune dimensions are controlled by processes like merging and splitting of dunes. Particularly the process of dune splitting is still poorly understood and – as a result – not yet included in operational dune evolution models. In the current paper, the process of dune splitting is investigated by carrying out laboratory experiments and by means of a sensitivity analysis using a numerical dune evolution model. In the numerical model, we introduced superimposed TRIAS ripples (i.e. triangular asymmetric stoss side‐ripples) on the stoss sides of underlying dunes as soon as these stoss sides exceed a certain critical length. Simulations with the model including dune splitting showed that predictions of equilibrium dune characteristics were significantly improved compared to the model without dune splitting. As dune splitting is implemented in a parameterized way, the computational cost remains low which means that dune evolution can be calculated on the timescale of a flood wave. Subsequently, we used this model to study the mechanism of dune splitting. Literature showed that the initiation of a strong flow separation zone behind a superimposed bedform is one of the main mechanisms behind dune splitting. The flume experiments indicated that besides its height also the lee side slope of the superimposed bedform is an important factor to determine the strength of the flow separation zone and therefore is an important aspect in dune splitting. The sensitivity analysis of the dune evolution model showed that a minimum stoss side length was required to develop a strong flow separation zone. Copyright © 2014 John Wiley & Sons, Ltd.     

Anomalous radar backscatter from Titan’s surface?

Since Cassini arrived at Saturn in 2004, its moon Titan has been thoroughly mapped by the RADAR instrument at 2-cm wavelength, in both active and passive modes. Some regions on Titan, including Xanadu and various bright hummocky bright terrains, contain surfaces that are among the most radar-bright encountered in the Solar System. This high brightness has been generally attributed to volume scattering processes in the inhomogeneous, low-loss medium expected for a cold, icy satellite surface. We can test this assumption now that the emissivity has been obtained from the concurrent radiometric measurements for nearly all the surface, with unprecedented accuracy (Janssen et al., and the Cassini RADAR Team [2009]. Icarus 200, 222-239). Kirchhoff’s law of thermal radiation relates the radar and radiometric properties in a way that has never been fully exploited. In this paper we examine here how this law may be applied in this case to better understand the nature of Titan’s radar-bright regions. We develop a quantitative model that, when compared to the observational data, allows us to conclude that either the reflective characteristics of the putative volume scattering subsurface must be highly constrained, or, more likely, organized structure on or in the surface is present that enhances the backscatter.     

Spatial Sciences on Ice: 50 years of Australian activities on the Amery Ice Shelf,East Antarctica

The Antarctic climate system involves many complex interactions between the atmosphere, oceans and ice and is sensitive to variations in these components. Ice shelves represent the ice–ocean–atmosphere interface of the Antarctic continent and are therefore very important indicators of climate change in the region. The Amery Ice Shelf is the largest ice shelf in East Antarctica and has been the focus of many scientific research projects over the past 50 years. This paper presents a history of the use of spatial science techniques from basic survey equipment to satellite systems in Australian research projects based on the Amery Ice Shelf (and surrounding glaciers) since 1955. The application of these spatial data to projects based primarily in the fields of geodesy, glaciology, climatology, and oceanography has allowed the measurement and monitoring of the physical, dynamic and environmental characteristics of this large and remote region. This new information provides scientists with a better understanding of the ice shelf/ocean/atmosphere system allowing future monitoring to observe the effects of global climate change.     

Calibration of Parameter Fields Consisting of Multiple Statistical Populations

If a parameter field to be calibrated consists of more than one statistical population, usually not only the parameter values are uncertain, but the spatial distributions of the populations are uncertain as well. In this study, we demonstrate the potential of the multimodal calibration method we proposed recently for the calibration of such fields, as applied to real-world ground water models with several additional stochastic parameter fields. Our method enables the calibration of the spatial distribution of the statistical populations, as well as their spatially correlated parameterization, while honoring the complete prior geostatistical definition of the multimodal parameter field. We illustrate the implications of the method in terms of the reliability of the posterior model by comparing its performance to that of a "conventional" calibration approach in which the positions of the statistical populations are not allowed to change. Information from synthetic calibration runs is used to show how ignoring the uncertainty involved in the positions of the statistical populations not only denies the modeler the opportunity to use the measurement information to improve these positions but also unduly influences the posterior intrapopulation distributions, causes unjustified adjustments to the cocalibrated parameter fields, and results in poorer observation reproduction. The proposed multimodal calibration allows a more complete treatment of the relevant uncertainties, which prevents the abovementioned adverse effects and renders a more trustworthy posterior model.     

Percolation losses in paddy fields with a dynamic soil structure: model development and applications

The hydraulic characteristics of the plough pan of paddy fields provide continuous ponding conditions during the growing season and control the water use efficiency in wet rice production. Its saturated hydraulic conductivity K_s, however, exhibits a large spatiotemporal variability as a consequence of a highly dynamic soil structure involving temporary shrinkage cracks. Water flow through the earthen bunds surrounding the fields further contributes to the uncertainty in water flux calculations. The objective of this study was to develop a simple deterministic model with stochastic elements (‘PADDY‐FLUX’) for depiction of deep percolation, and to assess the effect of different water management scenarios on percolation in two channel command areas. Darcy's law is used as the fundamental equation for water flow calculations with the ponding water depth h as a time‐dependent variable. Flux uncertainty is estimated by a Monte‐Carlo‐type implementation. K_s is treated as a random variable of a bimodal probability density function (PDF), which is the weighted sum of two Gaussian PDFs (accounting for a matrix and a preferential flow domain). The weighing factor α is a function of h, reflecting an increasing risk for preferential flow situations after desiccation and the development of shrinkage cracks. Under‐bund percolation is calculated using transfer functions. The results demonstrate that percolation losses increase in the following order: continuous soil saturation < continuous flooding (CF) < mid‐season drainage and intermittent irrigation (MD + II) < mid‐season drainage and continuous flooding. The bunds contribute up to 54 and 17% to total fluxes under CF and MD + II, respectively. Preferential water fluxes are responsible for the major part of water losses as soon as desiccation causes the formation of shrinkage cracks. As a conclusion, continuous soil saturation should be promoted as the least water‐intensive irrigation regime, while intermittent irrigation is recommended only in case that irreversible shrinkage cracks have already developed. Copyright © 2010 John Wiley & Sons, Ltd.     

Occurrence and distribution of microplastics in marine sediments along the Belgian coast

Plastic debris is known to undergo fragmentation at sea, which leads to the formation of microscopic particles of plastic; the so called ‘microplastics’. Due to their buoyant and persistent properties, these microplastics have the potential to become widely dispersed in the marine environment through hydrodynamic processes and ocean currents. In this study, the occurrence and distribution of microplastics was investigated in Belgian marine sediments from different locations (coastal harbours, beaches and sublittoral areas).Particles were found in large numbers in all samples, showing the wide distribution of microplastics in Belgian coastal waters. The highest concentrations were found in the harbours where total microplastic concentrations of up to 390 particles kg⁻¹ dry sediment were observed, which is 15-50 times higher than reported maximum concentrations of other, similar study areas.The depth profile of sediment cores suggested that microplastic concentrations on the beaches reflect the global plastic production increase.     

Impact of the water salinity on the hydraulic conductivity of fen peat

Coastal peatlands represent an interface between marine and terrestrial ecosystems; their hydrology is affected by salt and fresh water inflow alike. Previous studies on bog peat have shown that pore water salinity can have an impact on the saturated hydraulic conductivity (K_s) of peat because of chemical pore dilation effects. In this study, we aimed at quantifying the impact of higher salinities (up to 3.5% NaCl) on K_s of fen peat. Two experiments employing a constant‐head upward‐flow permeameter and differing in measurement and salinity change duration were conducted. Additionally, a third experiment to determine the impact of water salinity on the release of dissolved organic carbon (DOC) of the studied peat type was carried out. The results show a decrease of K_s with time, which does not depend on the water salinity but is differently shaped for different peat types. We assume pore clogging due to a conglomerate of physical, chemical, and biological processes, which rather depend on water movement rate and time than on water salinity. However, an increased water salinity did increase the DOC release. We conclude that salinity‐dependent behaviour of K_s is a function of peat chemistry and that for some peat types, salinity may only affect the DOC release without having a pronounced impact on water flow.     

Combining satellite and seismic images to analyse the shallow structure of the Dead Sea Transform near the DESERT transect

The left-lateral Dead Sea Transform (DST) in the Middle East is one of the largest continental strike-slip faults of the world. The southern segment of the DST in the Arava/Araba Valley between the Dead Sea and the Red Sea, called Arava/Araba Fault (AF), has been studied in detail in the multidisciplinary DESERT (DEad SEa Rift Transect) project. Based on these results, here, the interpretations of multi-spectral (ASTER) satellite images and seismic reflection studies have been combined to analyse geologic structures. Whereas satellite images reveal neotectonic activity in shallow young sediments, reflection seismic image deep faults that are possibly inactive at present. The combination of the two methods allows putting some age constraint on the activity of individual fault strands. Although the AF is clearly the main active fault segment of the southern DST, we propose that it has accommodated only a limited (up to 60 km) part of the overall 105 km of sinistral plate motion since Miocene times. There is evidence for sinistral displacement along other faults, based on geological studies, including satellite image interpretation. Furthermore, a subsurface fault is revealed ≈4 km west of the AF on two ≈E–W running seismic reflection profiles. Whereas these seismic data show a flower structure typical for strike-slip faults, on the satellite image this fault is not expressed in the post-Miocene sediments, implying that it has been inactive for the last few million years. About 1 km to the east of the AF another, now buried fault, was detected in seismic, magnetotelluric and gravity studies of DESERT. Taking together various evidences, we suggest that at the beginning of transform motion deformation occurred in a rather wide belt, possibly with the reactivation of older ≈N–S striking structures. Later, deformation became concentrated in the region of today’s Arava Valley. Till ≈5 Ma ago there might have been other, now inactive fault traces in the vicinity of the present day AF that took up lateral motion. Together with a rearrangement of plates ≈5 Ma ago, the main fault trace shifted then to the position of today’s AF.     

Spatial variability and interpolation of daily precipitation amount

Daily precipitation amounts show spatial variation over sub-continential regions. Point measurements, representative for regions of land, have to be interpolated towards unobserved locations. In this study four days in 1984 were selected to investigate the spatial variability of daily precipitation amount in North-western Europe in relation to the meteorological conditions. Data were interpolated using Kriging. Crossvalidation was used to compare interpolated values with measured values. Large differences in the spatial structure of daily precipitation amount are obsered as a result of different meterological conditions. Stratification of the study area into a coastal, a mountainous and an interior stratum proved to be successful, reducing the Mean Squared Error of Prediction with up to 55%.