Equilibrium modeling of the beach profile on a macrotidal embayed low tide terrace beach

Eleven-year long time series of monthly beach profile surveys and hourly incident wave conditions are analyzed for a macrotidal Low Tide Terrace beach. The lower intertidal zone of the beach has a pluriannual cycle, whereas the upper beach profile has a predominantly seasonal cycle. An equilibrium model is applied to study the variation of the contour elevation positions in the intertidal zone as a function of the wave energy, wave power, and water level. When forcing the model with wave energy, the predictive ability of the equilibrium model is around 60% in the upper intertidal zone but decreases to 40% in the lower intertidal zone. Using wave power increases the predictive ability up to 70% in both the upper and lower intertidal zones. However, changes around the inflection point are not well predicted. The equilibrium model is then extended to take into account the effects of the tide level. The initial results do not show an increase in the predictive capacity of the model, but do allow the model free parameters to represent more accurately the values expected in a macrotidal environment. This allows comparing the empirical model calibration in different tidal environment. The interpretation of the model free parameter variation across the intertidal zone highlights the behavior of the different zones along the intertidal beach profile. This contributes to a global interpretation of the four model parameters for beaches with different tidal ranges, and therefore to a global model applicable at a wide variety sites.     

Formation events of shoreline sand waves on a gravel beach

Kilometric-scale shoreline sand waves (KSSW) have been observed in the north-east flank of the Dungeness Cuspate Foreland (southeastern coast of the UK). They consist of two bumps separated by embayments with a 350–450-m spacing. We have analysed 36 shoreline surveys of 2-km length using the Discrete Fourier Transformation (DFT), from 2005 to 2016, and seven topographic surveys encompassing the intertidal zone, from 2010 to 2016. The data set shows two clear formation events. In order to test the role of high-angle waves on the KSSW formation, the 10-year wave series is propagated from the wave buoy located at 43 m depth up to a location in front of the undulations at 4 m depth using the SWAN wave model. The dominating SW waves arrive with a very high incidence angle (～ 80°) while the NE waves arrive almost shore normal. The ratio R, which measures the degree of dominance of high-angle waves with respect to low-angle waves, correlates well with the shoreline DFT magnitude values of the observed wavelength undulations. In particular, the highest R values coincide with the formation events. Finally, a linear stability model based on the one-line approximation is applied to the Dungeness profile and the 10-year propagated wave series. It predicts accurately the formation moments, with positive growth rates in the correct order of magnitude for wavelengths similar to the observed ones. All these results confirm that the shoreline undulations in Dungeness are self-organized and that the underlying formation mechanism is the high-angle wave instability. The two detected formation events provide a unique opportunity to validate the existing morphodynamic models that include such instability.     

The Wadden Sea in transition - consequences of sea level rise

The impact of sea level rise (SLR) on the future morphological development of the Wadden Sea (North Sea) is investigated by means of extensive process-resolving numerical simulations. A new sediment and morphodynamic module was implemented in the well-established 3D circulation model GETM. A number of different validations are presented, ranging from an idealized 1D channel over a semi-idealized 2D Wadden Sea basin to a fully coupled realistic 40-year hindcast without morphological amplification of the Sylt-Rømøbight, a semi-enclosed subsystem of the Wadden Sea. Based on the results of the hindcast, four distinct future scenarios covering the period 2010–2100 are simulated. While these scenarios differ in the strength of SLR and wind forcing, they also account for an expected increase of tidal range over the coming century. The results of the future projections indicate a transition from a tidal-flat-dominated system toward a lagoon-like system, in which large fractions of the Sylt-Rømøbight will remain permanently covered by water. This has potentially dramatic implications for the unique ecosystem of the Wadden Sea. Although the simulations also predict an increased accumulation of sediment in the back-barrier basin, this accumulation is far too weak to compensate for the rise in mean sea level.     

Dense collection areas and terrestrial alteration of meteorites in the Atacama Desert

In the last 15 years, more than 2700 meteorites have been recovered and officially classified from the Atacama Desert. Although the number of meteorites collected in the Atacama has risen, the physical and climatic properties of the dense collection areas (DCAs) have not been fully characterized. In this article, we compiled the published data of all classified meteorites found in the Atacama Desert to (i) describe the distribution by meteorite groups, (ii) compare the weathering degree of chondrites among different Atacama DCAs and other hot and cold deserts, and (iii) determine the preservation conditions of chondrites in the main Atacama DCAs in relation with the local climatic conditions. The 35 DCAs so far identified in the Atacama Desert are located in three main morphotectonic units: The Coastal Range (CR), Central Depression (CD), and Pre-Andean Range/Basement. A comparison with reported weathering data from other cold and hot deserts indicates that the mean terrestrial weathering of Atacama chondrites (W1–2), displays less alteration than other hot deserts (W2–3) and resembles the weathering distribution of the Antarctic meteorites (W1–2). The highest abundance of Atacama chondrites with low weathering (≤W2) is localized in the CD (78.8%, N = 1435), which is protected from the coastal fog influence and seasonal rainfalls and displays the oldest surfaces in the Atacama Desert. The morphogenetic classification based on present-day temperatures and precipitations of the main Atacama DCAs reveals similar regional/subregional climatic conditions in the most productive areas and a truly productive surface for meteorite recovery between 5% and 58% of the quadrangles formally defined for each Atacama DCA. Our morphogenetic classification lacks consideration of some meteorological parameters such as the coastal fog, so it cannot fully explain the differences in weathering patterns among CR chondrites. Future studies of chondrite preservation in the Atacama DCAs should consider other meteorological variables such as relative humidity, specific humidity, or dew point, in combination with exposure ages of meteorites and its surfaces.     

Modelling regional effects of artificial groundwater recharge in a multilayer aquifer characterized by perched water tables

This paper presents an approach to estimate the effects of a managed recharge experiment in a multilayer aquifer characterized by the presence of perched water tables in the Medina del Campo groundwater body, Douro basin, central Spain. A numerical model was developed to evaluate the effect of artificial recharge on the shallow sector of a regional-scale aquifer and on formerly active wetlands. The model was developed in the Visual MODFLOW Pro v.2011.1 environment in order to represent and analyse the regional impact of this artificial recharge event. Results suggest that the assumption of a single perched system may prove useful in regional contexts where data is limited. From a study site perspective, managed recharge is observed to increase shallow storage along the riverbanks, which is considered valuable for environmental purposes. However, downstream wetlands are unlikely to experience a significant recovery. Furthermore, only a small percentage of artificial recharge is expected to reach the deep regional aquifer. This method can be exported to settings characterized by the presence of perched aquifers and associated groundwater dependent ecosystems.     

Role of time variant and invariant terms on soil moisture spatio-temporal variability in a subhumid coastal wetland

The knowledge of soil moisture spatio-temporal variability is highly relevant for water resources management. This paper reports an analysis of the spatial–temporal variability of soil moisture data for a small to medium-scale soil-sensors network in a coastal wetland of southwestern Spain. Measurements were taken from five sites located in the Doñana National Park over the time-period of one hydrological year from September 2017 to September 2018. The total area of the soil-sensors network shows an extension about 25 × 3 km. Soil moisture data was separated into time invariant (the temporal mean of the whole period at each site) and time-variant terms (the deviations of soil moisture from the mean, or anomalies). The time-invariant component was generally the main contributor to the total spatial variance of soil moisture and it was mostly controlled by the groundwater levels in the area. Nevertheless, the time variant terms have a huge effect on soil moisture variability in very dry states. Characteristic convex time-dependent patterns for this field site were found between spatially averaged soil moisture and its variability. This information could be used for the up and downscaling of soil moisture from satellite data. Those patterns of relation between spatial mean and variability of soil moisture were only affected by heavy rainfalls giving rise to hysteretic behaviour. This study shows that even though groundwater level is a time-variant variable, it significantly affects soil moisture's time-variant but also time-invariant terms due to the different average groundwater level depths at the different sites.     

A first look at the effects of ionospheric signal bending on a globally processed GPS network

This study provides a first attempt at quantifying potential signal bending effects on the GPS reference frame, coordinates and zenith tropospheric delays (ZTDs). To do this, we homogeneously reanalysed data from a global network of GPS sites spanning 14 years (1995.0–2009.0). Satellite, Earth orientation, tropospheric and ground station coordinate parameters were all estimated. We tested the effect of geometric bending and dTEC bending corrections, which were modelled at the observation level based, in part, on parameters from the International Reference Ionosphere 2007 model. Combined, the two bending corrections appear to have a minimal effect on site coordinates and ZTDs except for low latitude sites. Considering five days (DOY 301–305, 28 October–1 November 2001) near ionospheric maximum in detail, they affect mean ZTDs by up to ~1.7 mm at low latitudes, reducing to negligible levels at high latitudes. Examining the effect on coordinates in terms of power-spectra revealed the difference to be almost entirely white noise, with noise amplitude ranging from 0.3 mm (high latitudes) to 2.4 mm (low latitudes). The limited effect on station coordinates is probably due to the similarity in the elevation dependence of the bending term with that of tropospheric mapping functions. The smoothed z-translation from the GPS reference frame to ITRF2005 changes by less than 2 mm, though the effect combines positively with that from the second order ionospheric refractive index term. We conclude that, at the present time, and for most practical purposes, the geometric and dTEC bending corrections are probably negligible at current GPS/reference frame precisions.     

Analysis of soil drying incorporating a constitutive model for curling

This paper focuses on the shrinkage behavior of soil specimens involving sand, kaolinite, and kaolinite/sand mixtures subjected to desiccation under controlled conditions. Both, free and restrained shrinkage conditions are studied. The experiments show that pure soils do not curl upon unrestrained shrinkage; however, (under the same conditions) kaolinite/sand mixtures exhibited a marked curling. Furthermore, the mixture with the higher sand content broke through the middle of the sample after displaying a significant curling. Soils subjected to restricted shrinkage developed cracks with slight curling. To simulate the observed behavior, a mechanical model able to reproduce the detachment of the soil sample from the mold is proposed in this work and implemented in a fully coupled hydro-mechanical finite-element code. It is concluded that suction and differential shrinkage are key factors influencing the curling behavior of soils. The proposed framework was able to satisfactorily explain and reproduce the different stages and features of soil behavior observed in the experiments.

     

The Radial Basis Functions Method for Improved Numerical Approximations of Geological Processes in Heterogeneous Systems

Mathematical Geosciences - A robust, high order modeling approach is introduced, based on the finite difference-based radial basis functions method, for solving the groundwater flow equation in the...