Long-term soil moisture dynamics derived from GNSS interferometric reflectometry: a case study for Sutherland,South Africa

Soil moisture is a geophysical key observable for predicting floods and droughts, modeling weather and climate and optimizing agricultural management. Currently available in situ observations are limited to small sampling volumes and restricted number of sites, whereas measurements from satellites lack spatial resolution. Global navigation satellite system (GNSS) receivers can be used to estimate soil moisture time series at an intermediate scale of about 1000 m². In this study, GNSS signal-to-noise ratio (SNR) data at the station Sutherland, South Africa, are used to estimate soil moisture variations during 2008–2014. The results capture the wetting and drying cycles in response to rainfall. The GNSS Volumetric Water Content (VWC) is highly correlated (r ² = 0.8) with in situ observations by time-domain reflectometry sensors and is accurate to 0.05 m³/m³. The soil moisture estimates derived from the SNR of the L1 and L2P signals compared to the L2C show small differences with a RMSE of 0.03 m³/m³. A reduction in the SNR sampling rate from 1 to 30 s has very little impact on the accuracy of the soil moisture estimates (RMSE of the VWC difference 1–30 s is 0.01 m³/m³). The results show that the existing data of the global tracking network with continuous observations of the L1 and L2P signals with a 30-s sampling rate over the last two decades can provide valuable complementary soil moisture observations worldwide.     

Recursive algorithm for fast GNSS orbit fitting

Gaussian elimination is an efficient and numerically stable algorithm for estimating parameters and their precision. However, before estimating the parameters, it is often prudent to perform statistical tests to achieve the best fitting model. We use Gaussian elimination to select the best fitting model among candidate models. A succinct relationship between the weighted sum of squared residuals and the previous one is revealed by a volume formula. For quick parameter estimation and determination of weighted sum of squared residuals, a recursive elimination algorithm is proposed in the context of Gaussian elimination. In order to improve the model selection efficiency, the parameter estimation and the determination of the weighted sum of squared residuals are carried out in parallel using the proposed recursive elimination algorithm in which the improvement at each recursive stage is judged by the Bayesian information criterion. Ultimately, the computational complexity and numerical stability of the recursive elimination proposed are briefly discussed, and a GNSS orbit interpolation example is used to verify the results. It shows that the proposed recursive elimination algorithm inherits the numerical stability of the Gaussian elimination, and this algorithm can be used to examine the gain from the newly introduced parameter, dynamically assess the fitting model, and fix the optimal model efficiently. The optimal fitting model with the lowest information is very close to the real situation verified by checkpoints.     

A new method for specular and diffuse pseudorange multipath error extraction using wavelet analysis

Multipath remains one of the major challenges in Global Navigation Satellite System (GNSS) positioning because it is considered the dominant source of ranging errors, which can be classified into specular and diffuse types. We present a new method using wavelets to extract the pseudorange multipath in the time domain and breaking it down into the two components. The main idea is an analysis-reconstruction approach based on application of both continuous wavelet transform (CWT) and discrete wavelet transform (DWT). The proposed procedure involves the use of L1 code-minus-carrier (CMC) observable where higher-frequency terms are isolated as residuals. CMC residuals are analyzed by applying the CWT, and we propose the scalogram as a technique for discerning time–frequency variations of the multipath signal. Unlike Fourier transform, the potential of the CWT scalogram for examining the non-stationary and multifrequency nature of the multipath is confirmed as it simultaneously allows fine detection and time localization of the most representative frequencies of the signal. This interpretation of the CWT scalogram is relevant when choosing the levels of reconstruction with DWT, allowing accurate time domain extraction of both the specular and diffuse multipath. The performance and robustness of the method and its boundary applicability are assessed. The experiment was carried out using a receiver of Campania GNSS Network. The results are given in which specular multipath error is achieved using DWT level 7 approximation component and diffuse multipath error is achieved using DWT level 6 denoised detail component.     

Optimization and Parametrical Investigation to Assess the Reconstitution of Different Types of Indian Sand Using Portable Travelling Pluviator

By using the air pluviation technique, it is aimed to achieve the desired relative density with uniform void ratio throughout the specimen in order to maintain the homogeneity and to avoid the spatial variability. Further, in order to achieve the maximum deposition intensity, a systematic optimization study has been carried out rigorously in a test tank to determine the diameter of the orifice to be employed for the sieve plates of different porosity and the number of sieve plates to be installed in the diffuser sieve sets. The study has been conducted with four different patterns of sieves with different porosity to achieve a wide range of relative densities for four different uniformly graded Indian sands. The dynamic penetrometer which is considered to be one of the cost effective instruments has been efficiently used to determine the soil resistance at various locations of the test tank for every given height of fall in order to check the uniformity of placement density throughout the sand bed. The study reveals that the sand beds of different relative densities could be achieved using different patterns of diffuser sieves at optimum sand flow rate without compromising the uniformity. The effect of height of fall as well as porosity of diffuser sieves on the relative density of different sands has been studied in detail. The deposition intensity and the relative density obtained from the present study are compared with the values available in the literature.     

Assessment and Prediction of Liquefaction Potential Using Different Artificial Neural Network Models: A Case Study

Soil liquefaction as a transformation of granular material from solid to liquid state is a type of ground failure commonly associated with moderate to large earthquakes and refers to the loss of strength in saturated, cohesionless soils due to the build-up of pore water pressures and reduction of the effective stress during dynamic loading. In this paper, assessment and prediction of liquefaction potential of soils subjected to earthquake using two different artificial neural network models based on mechanical and geotechnical related parameters (model A) and earthquake related parameters (model B) have been proposed. In model A the depth, unit weight, SPT-N value, shear wave velocity, soil type and fine contents and in model B the depth, stress reduction factor, cyclic stress ratio, cyclic resistance ratio, pore pressure, total and effective vertical stress were considered as network inputs. Among the numerous tested models, the 6-4-4-2-1 structure correspond to model A and 7-5-4-6-1 for model B due to minimum network root mean square errors were selected as optimized network architecture models in this study. The performance of the network models were controlled approved and evaluated using several statistical criteria, regression analysis as well as detailed comparison with known accepted procedures. The results represented that the model A satisfied almost all the employed criteria and showed better performance than model B. The sensitivity analysis in this study showed that depth, shear wave velocity and SPT-N value for model A and cyclic resistance ratio, cyclic stress ratio and effective vertical stress for model B are the three most effective parameters on liquefaction potential analysis. Moreover, the calculated absolute error for model A represented better performance than model B. The reasonable agreement of network output in comparison with the results from previously accepted methods indicate satisfactory network performance for prediction of liquefaction potential analysis.     

A multi-dimensional assessment of urban vulnerability to climate change in Sub-Saharan Africa

In this paper, we develop and apply a multi-dimensional vulnerability assessment framework for understanding the impacts of climate change-induced hazards in Sub-Saharan African cities. The research was carried out within the European/African FP7 project CLimate change and Urban Vulnerability in Africa, which investigated climate change-induced risks, assessed vulnerability and proposed policy initiatives in five African cities. Dar es Salaam (Tanzania) was used as a main case with a particular focus on urban flooding. The multi-dimensional assessment covered the physical, institutional, attitudinal and asset factors influencing urban vulnerability. Multiple methods were applied to cover the full range of vulnerabilities and to identify potential response strategies, including: model-based forecasts, spatial analyses, document studies, interviews and stakeholder workshops. We demonstrate the potential of the approach to assessing several dimensions of vulnerability and illustrate the complexity of urban vulnerability at different scales: households (e.g., lacking assets); communities (e.g., situated in low-lying areas, lacking urban services and green areas); and entire cities (e.g., facing encroachment on green and flood-prone land). Scenario modeling suggests that vulnerability will continue to increase strongly due to the expected loss of agricultural land at the urban fringes and loss of green space within the city. However, weak institutional commitment and capacity limit the potential for strategic coordination and action. To better adapt to urban flooding and thereby reduce vulnerability and build resilience, we suggest working across dimensions and scales, integrating climate change issues in city-level plans and strategies and enabling local actions to initiate a ‘learning-by-doing’ process of adaptation.     

Reconstruction of Holocene paleohydrological and sedimentological characteristics of Umm Ar-Rimam basin (northern Kuwait)

The paleohydrological and sedimentological characteristics of a playa lake in northern Kuwait (Arabian Gulf) are reconstructed using sedimentological, geochemical, and isotopic techniques. The sequence consists of up to 8 cycles of S-poor, alluvial sediments capped by a thin organic soil interbedded with gravity-fall calcrete sediments. The succession is locally derived from mainly Quaternary sediments and is regressive with upsection filling of the subsiding basin by cycles of sheetwash flow in response to climatic change. There is no natural, open-water lake water as indicated by low total organic carbon (TOC) data, but the presence of incised calcrete yardangs suggests that more extensive open-water conditions were operative in the past. Stable isotope (δ¹⁸O‰ and δ¹³C‰) values of the authigenic carbonates indicate the following three distinct processes: evaporation, meteoric fluid infiltration, and rapid per-descensum flow (rapid downward movement of water and playa sediment through pipes) through a porous, clastic sequence. Because evaporites are scarce, other factors besides evaporation action control chemical and isotopic compositions of the per-decensum lake fluids. Consequently, the isotopic composition cannot be interpreted exclusively as an indicator of salinity or evaporation ratio. The degassing of CO₂ during groundwater discharge may explain the enriched carbon isotope values for the authigenic carbonates precipitated in the sediments. Hydrologically closed lake water bodies tend to show low negative carbonate oxygen and carbon isotopic signatures. Isotopically negative δ¹³C values imply a strong input of soil-zone carbon to the groundwater of the top 60 cm of the sediment. Lakes that are hydrologically closed and evaporate or equilibrate with atmospheric CO₂ will tend to have low negative δ¹⁸O and δ¹³C values in the carbonates as reported by Talbot (Chem Geol: Isotope Geosci Sect 80(4):261–279, 1990). Biologically active lakes will tend toward lower δ¹³C of dissolved carbon due to the photosynthetic effects of ¹²C withdrawal as reported by Dunagan and Driese (J Sed Res 69:772–783, 1999). Increased biological activity during sedimentation may account for low carbon isotope values where plants were abundant, but in shrinkage-dominated systems (those of clay-rich soil subjecting to wet-dry conditions), carbon isotopes will be largely inherited from the calcretic limestones in the land extending landward of the coast and not influenced by coastal processes (known as hinterland), such as Umm Ar-Rimam depression. This basin does not fit the classic shallow playa-type basins of the Arabian Peninsula but rather the recharge playas of the southwestern USA.     

Tectono-sedimentary evolution of the Miocene Oued Dayr Formation (Ghomaride Complex,Internal Domain of the Maghrebian Rif Belt,Morocco)

This work deals with sedimentological, petrographic, and structural analyses of a middle Miocene late-orogenic sedimentary cycle, denoted Oued Dayr Formation, recognized in the Rifian sector of the Maghrebian Chain (Morocco). The analyzed Formation (75 m thick) starts with 15–20 m of light colored polymict conglomerates, with minor sandstone beds, lying on the Paleozoic basement and Mesozoic cover of the Ghomaride Nappe. Facies analysis indicates a fining-upward deposition in a marine environment characterized by increasing deepening, reflecting a subsidence rate that exceeds sedimentary supply. Petrographic analysis points out that sandstones are represented by litharenites originated by erosion of recycled orogen. The conglomerates pebbles and cobbles consist of Alpine low- to high-grade metamorphic rocks as metarenites, phyllites, mylonitic quartzites, micaschists, augen gneisses deriving from the exhumed deep metamorphic basement, the overlying metasedimentary of the Sebtide Nappes and of sedimentary rocks as sandstones, jaspes, limestones, and shales deriving from the Ghomaride Nappes and their sedimentary cover. Data reveal mixed provenance indicating that the Oued Dayr Formation was fed by the Internal Nappes stack of the Maghrebian Chain. Structural analysis shows that the Oued Dayr Formation accumulated in a Thrust-Top basin, during an early extension (D₀ phase), recorded by synsedimentary normal faults within middle Langhian deposits on the rear of the Internal Nappes stack. Subsequent ductile and brittle compressional (D₁, D₂, D₃) and extensional (D₄) deformation phases occurred during and/or after the stacking, exhumation, and early unroofing of Sebtide Complex coeval with the opening of the western Mediterranean back-arc basins since middle Miocene time.     

Sensitivity analysis for the GALDIT method based on the assessment of vulnerability to pollution in the northern Sfax coastal aquifer,Tunisia

Coastal aquifer of northern Sfax (Tunisia) suffers from the high risk to seawater intrusion and the water quality degradation due to the overexploitation. Hence, assessing the study area vulnerability to pollution is highly crucial so as to protect the groundwater resources. The assessment has been performed by applying the GALDIT method using Geographic Information System (GIS) software and multi-criteria evaluation techniques, and the sensitivity analysis approach to evaluate the effect of each GALDIT parameter on the vulnerability assessment. The GALDIT vulnerability map classifies the study area into three vulnerability classes: low vulnerability (30–50), moderate vulnerability (50–70), and high vulnerability (70–90), which represent 5, 30, and 65 % of the study area, respectively. The map illustrates that the coastal zones of the aquifer are the most threatened areas. The sensitivity analysis results show that the aquifer hydraulic conductivity (A) and the thickness of the aquifer (T) represent the determining factors in the modified vulnerability model. The real weight was used to elaborate the modified GALDIT model which was correlated with resistivity values for validation. This study could serve as a scientific basis for sustainable land planning and groundwater management in the study area.     

<Emphasis>H</Emphasis>ydrological variability of the <Emphasis>S</Emphasis>oummam watershed <Emphasis>(N</Emphasis>ortheastern <Emphasis>A</Emphasis>lgeria<Emphasis>)</Emphasis> and the possible links to climate fluctuations

The long-term variability of rainfall in the Soummam watershed (NE Algeria) has been analysed over the past 108 years using continuous wavelet method in order to identify the interannual modes controlling the rainfall variability. Statistical analyses of rainfall timeseries have shown its distribution following five periods of time, limited by a series of discontinuities around 1935, 1950, 1970 and 1990. The continuous wavelet transform have demonstrated different low frequency modes: 2–4, 4–8, 8–16 and 16–32 years.The annual band is expanded during the full study period with some pics around 1905, 1920–1935 and 1960; it shows a negative long-term trend, in particular since the period 1970–1990 when a major change has been identified. Then, the relationships between climate patterns of North Atlantic Oscillation (NAO) and Southern Oscillation Index (SOI) and the hydrological variability in the frequency domain have been investigated; they have shown a mean explained variance of 40 and 24 %, respectively. Such variances are less obvious for the annual mode and increase for the interannual frequencies. The coherence suffer from high perturbations since the period 1970–1990 when the NAO (SOI) shifts from negative (positive) phases to positive (negative) ones. Such anomalies are responsible for significant changes of rainfall variability, emphasising the global warming effects.