Modelling soil moisture in a high‐latitude landscape using LiDAR and soil data

Soil moisture has a pronounced effect on earth surface processes. Global soil moisture is strongly driven by climate, whereas at finer scales, the role of non‐climatic drivers becomes more important. We provide insights into the significance of soil and land surface properties in landscape‐scale soil moisture variation by utilizing high‐resolution light detection and ranging (LiDAR) data and extensive field investigations. The data consist of 1200 study plots located in a high‐latitude landscape of mountain tundra in north‐western Finland. We measured the plots three times during growing season 2016 with a hand‐held time‐domain reflectometry sensor. To model soil moisture and its temporal variation, we used four statistical modelling methods: generalized linear models, generalized additive models, boosted regression trees, and random forests. The model fit of the soil moisture models were R² = 0.60 and root mean square error (RMSE) 8.04 VWC% on average, while the temporal variation models showed a lower fit of R² = 0.25 and RMSE 13.11 CV%. The predictive performances for the former were R² = 0.47 and RMSE 9.34 VWC%, and for the latter R² = 0.01 and RMSE 15.29 CV%. Results were similar across the modelling methods, demonstrating a consistent pattern. Soil moisture and its temporal variation showed strong heterogeneity over short distances; therefore, soil moisture modelling benefits from high‐resolution predictors, such as LiDAR based variables. In the soil moisture models, the strongest predictor was SAGA (System for Automated Geoscientific Analyses) wetness index (SWI), based on a 1 m² digital terrain model derived from LiDAR data, which outperformed soil predictors. Thus, our study supports the use of LiDAR based SWI in explaining fine‐scale soil moisture variation. In the temporal variation models, the strongest predictor was the field‐quantified organic layer depth variable. Our results show that spatial soil moisture predictions can be based on soil and land surface properties, yet the temporal models require further investigation. Copyright © 2017 John Wiley & Sons, Ltd.     

222Rn as Natural Tracer for LNAPL Recovery in a Crude Oil‐Contaminated Aquifer

The objective of this study was to investigate whether ²²²Rn in groundwater can be used as a tracer for light non‐aqueous phase liquid (LNAPL) quantification at a field site treated by dual‐phase LNAPL removal. After the break of a pipeline, 5 ha of soil in the nature reserve Coussouls de Crau in southern France was contaminated by 5100 m³ of crude oil. Part of this oil seeped into the underlying gravel aquifer and formed a floating oil body of about 3.9 ha. The remediation consists of plume management by hydraulic groundwater barriers and LNAPL extraction in the source zone. ²²²Rn measurements were performed in 21 wells in and outside the source zone during 15 months. In uncontaminated groundwater, the radon activity was relatively constant and remained always >11 Bq/L. The variability of radon activity measurements in wells affected by the pump‐and‐skim system was consistent with the measurements in wells that were not impacted by the system. The mean activities in wells in the source zone were, in general, significantly lower than in wells upgradient of the source zone, owing to partitioning of ²²²Rn into the oil phase. The lowest activities were found in zones with high non‐aqueous phase liquid (NAPL) recovery. LNAPL saturations around each recovery well were furthermore calculated during a period of high groundwater level, using a laboratory‐determined crude oil–water partitioning coefficient of 38.5 ± 2.9. This yielded an estimated volume of residual crude oil of 309 ± 93 m³ below the capillary fringe. We find that ²²²Rn is a useful and cheap groundwater tracer for finding zones of good LNAPL recovery in an aquifer treated by dual‐phase LNAPL removal, but that quantification of NAPL saturation using Rn is highly uncertain.     

Export and degassing of terrestrial carbon through watercourses draining a temperate podzolized catchment

We measured spatial and temporal variations in carbon concentrations, isotopic compositions and exports during a complete hydrological cycle in nine watercourses draining a lowland forested podzolized catchment, flowing into the Arcachon lagoon (France). In addition, integrated fluxes of CO₂ across the water-atmosphere interface were estimated to assess the relative importance of CO₂ evasion versus lateral carbon transport at the catchment scale. Watercourse similarities and specificities linked to the local catchment characteristics are discussed and compared with other riverine systems. Low concentrations of suspended particulate matter and particulate organic carbon (POC) were generally measured in all the watercourses (8.4 ± 3.4 and 1.6 ± 0.6 mg L^?1, respectively), reflecting limited mechanical soil erosion. The generally high POC content in the suspended matter (20 %), low Chl a concentrations (1.3 ± 1.4 μg L^?1) and the relatively constant δ¹³C-POC value (near ?28 ‰) throughout the year reveal this POC originates from terrestrial C₃ plant and soil detritus. The presence of podzols leads to high levels of dissolved organic carbon (DOC; 6.6 ± 2.2 mg L^?1). Similarly, high dissolved inorganic carbon (DIC) concentrations were measured in the Arcachon lagoon catchment (5.9 ± 2.2 mg L^?1). The δ¹³C-DIC value around ?20 ‰ throughout the year in many small watercourses reveals the predominance of terrestrial carbon mineralisation and silicate rock weathering in soils as the major DIC source. With pCO₂ between 1,000 and 10,000 ppmv, all watercourses were a source of CO₂ to the atmosphere, particularly during the low river stage. Organic carbon parameters remained relatively stable throughout the year, whereas DIC parameters showed strong seasonal contrasts closely linked to the hydrological regime and hyporheic flows. In total, the carbon export from the Arcachon watershed was estimated at 15,870 t C year^?1 or 6 t C km^?2 year^?1, mostly exported to the lagoon as DOC (35 %), DIC (24 %) and lost as CO₂ degassing to the atmosphere (34 %).     

Dynamique et modélisation des crues dans le bassin du Mono à Nangbéto (Togo/Bénin)

Résumé

Les inondations des dernières années en Afrique de l’Ouest rendent nécessaire la prévision hydrologique pour réduire la vulnérabilité des populations. Nous avons étudié la dynamique des crues dans le bassin du Mono en amont du barrage de Nangbéto (Togo et Bénin) par des méthodes statistiques et déterministes. L’augmentation des précipitations au cours de la décennie 1990 et la mise en eau du barrage de Nangbéto en 1988 ont modifié le régime du Mono, dont les débits d’étiages ont augmenté. Les débits entrants et sortants du barrage sont équivalents, même lors des crues. Les maximums annuels des débits entrants peuvent être modélisés par une distribution des valeurs extrêmes (GEV) de type Fréchet, tandis que ceux des débits sortants, suivent une loi GEV de type Weibull (bornée). Le mécanisme dominant de génération des crues dans le bassin est la montée progressive des eaux associée à l’augmentation de la saturation des sols. L’analyse des débits de crue entre 1988 et 2010 n’a pas montré d’augmentation de la magnitude ou de la fréquence des crues. Le modèle GR4J reproduit bien les débits de crues en période calage, mais ces critères de performance se dégradent lorsque l’on s’en éloigne. Dans ces conditions, l’utilisation de ce modèle pour évaluer les risques futurs de crues serait hasardeuse et il est nécessaire d’envisager d’autres types d’approches.     

Effect of hysteretic smoothness on inelastic response spectra with constant‐ductility

This paper presents an efficient methodology for computing constant‐ductility inelastic response spectra. The computation of constant‐ductility spectra involves numerical root‐finding algorithms to find the strongest structure providing a desired ductility response. Smooth inelastic structural behavior is modeled using a first‐order nonlinear differential equation and the transient structural response is solved using an implicit algorithm requiring Newton iterations at each time step. For structural models with smooth hysteretic behavior (not piece‐wise linear), a simple root‐finding method involving a combination of hyperbolic fits, linear interpolation, and Newton's method converges upon the highest strength (conservative) solution with a small number of iterations. The effect of the hysteretic smoothness on the occurrence of multiple roots is examined for two near‐fault and two far‐fault earthquake records, and for two measures of ductility and for normalized hysteretic energy. The results indicate how the smoothness of the hysteretic behavior affects ductility demand and constant‐ductility response spectra. Copyright © 2010 John Wiley & Sons, Ltd.     

Alkaline fluid circulation in ultramafic rocks and formation of nucleotide constituents: a hypothesis

Seawater is constantly circulating through oceanic basement as a low-temperature hydrothermal fluid (<150°C). In cases when ultramafic rocks are exposed to the fluids, for instance during the initial phase of subduction, ferromagnesian minerals are altered in contact with the water, leading to high pH and formation of secondary magnesium hydroxide, among other – brucite, that may scavenge borate and phosphate from seawater. The high pH may promote abiotic formation of pentoses, particularly ribose. Pentoses are stabilized by borate, since cyclic pentoses form a less reactive complex with borate. Analyses have shown that borate occupies the 2' and 3' positions of ribose, thus leaving the 5' position available for reactions like phosphorylation. The purine coding elements (adenine, in particular) of RNA may be formed in the same general hydrothermal environments of the seafloor.     

Late Miocene to Quaternary slip history across the Qiulitag anticline in the southern Tianshan piedmont

In this study, we reconstruct the Miocene to Quaternary shortening history across the Qiulitag anticline, a complex fault‐bend fold located in southern Tianshan. We studied the Yaha and Kuche sections, where we combined surface structural measurements and seismic imaging to model the stratigraphic horizons. The history of folding was reconstructed based on magnetostratigraphic analyses and eight cosmogenic burial ages in Kuche. Pleistocene deformation rates were also quantified in Yaha based on a deformed fluvial terrace that we dated to ~67 ka using a cosmogenic depth profile. Our results suggest that the fold grew at a mean slip rate of 0.9–1.3 mm/a in both sections but accelerated to ~2.5 mm/a during the Pleistocene in Kuche. These results support a migration of the deformation towards the basin during the Pleistocene and suggest that most of the present deformation in the Tianshan is accommodated across the external structures of the range.     

Zodiacal light gathered along the line of sight: Retrieval of the local scattering coefficient from photometric surveys of the ecliptic plane

A clue towards a retrieval of the zodiacal brightness gathering along a line of sight in the ecliptic plane consists in introducing the other intersection of that line with the terrestrial orbit (Fig. 1). The distribution of the elemental contribution to the brightness, or of the local quantity $\text{D}$ [directional scattering coefficient, i.e. cross-section of the unit-volume, which gives very simple expressions (1), (2) for the brightness integral] can then be approached with reduced uncertainty. The assumptions-steady state of the zodiacal cloud; smooth distribution of $\text{D}$—are strongly suggested by the observations, and are much less controversial than the classical assumption of uniform composition and size everywhere.The scattering coefficient may vary along the line of sight as seen in Fig. 3 : an uncertainty bar highly dependent of the abscissa, and considerably reduced in the vicinity of two “nodes”. Both in abscissae and in ordinates, these nodes are conspicuously insensitive to the arbitrary choice of a mathematical model (Table 1).The node exterior to the Earth's orbit (“martian node”) remains at r ? 1.5 a.u. from the Sun (Fig. 4). It gives access to a range of the scattering phase function near Mars' orbit, deconvolved from any radial dependence of that function (Fig. 5). The backscattering effect obtained is a new confirmation of the non-terrestrial origin of the gegenschein.The node interior to the Earth's orbit remains located not far from the middle of each chord (“quasi-radial node”. Fig. 4). It allows to retrieve the radial dependence of $\text{D}$, partly deconvolved from its angular dependence, between 0.5 and 1 a.u. (Fig. 6 and Table 4).The uncertainty bars on $\text{D}$ at the two observing locations yield two uncertainty bars of the phase function σ(θ) at 1 a.u. (Fig. 7). At θ = 30°, the forward scattering efficiency (normalized to θ = 90°) cannot exceed 6 and more likely 4. This disagrees with higher values obtained assuming spherical particles, and even obtained in part of the more realistic studies (assuming irregularly shaped particles, or mainly observational) reviewed in Table 5.All of these results are derived, with fair agreement, from three independent observational sources.     

The provenance of Australian uranium ore concentrates by elemental and isotopic analysis

Elemental and isotopic ratio analyses of U ore concentrate samples, from the 3 operating U mining facilities in Australia, were carried out to determine if significant variations exist between their products, thereby allowing the U ore concentrate’s origin to be identified. Elemental analyses were conducted using inductively coupled plasma mass spectrometry (ICP-MS) and X-ray fluorescence spectrometry (XRF). Lead isotope ratios were measured using ICP-MS and U isotope analyses were conducted using thermal ionisation mass spectrometry (TIMS). Minute quantities of sample, such as that obtained from a swipe, were also examined for elemental concentrations using secondary ion mass spectrometry (SIMS). The results of multivariate statistical analysis show clear patterns in the trace elemental composition of the processed U ores, indicating that it is possible to use this feature as a unique identifier of an Australian U ore concentrate’s source. Secondary ion mass spectrometry analyses also allow individual particles to be differentiated using this ‘fingerprinting’ technique. Isotope ratios determined using TIMS reveal that there is a significant difference in the n(²³⁴U)/n(²³⁸U) isotope ratio between the U ore concentrate from each mine.