Estimating water balance components in irrigated agriculture using a combined approach of soil moisture and energy balance monitoring,and numerical modelling

Information on water balance components such as evapotranspiration and groundwater recharge are crucial for water management. Due to differences in physical conditions, but also due to limited budgets, there is not one universal best practice, but a wide range of different methods with specific advantages and disadvantages. In this study, we propose an approach to quantify actual evapotranspiration, groundwater recharge and water inflow, i.e. precipitation and irrigation, that considers the specific conditions of irrigated agriculture in warm, arid environments. This approach does not require direct measurements of precipitation or irrigation quantities and is therefore suitable for sites with an uncertain data basis. For this purpose, we combine soil moisture and energy balance monitoring, remote sensing data analysis and numerical modelling using Hydrus. Energy balance data and routine weather data serve to estimate ET₀. Surface reflectance data from satellite images (Sentinel-2) are used to derive leaf area indices, which help to partition ET₀ into energy limited evaporation and transpiration. Subsequently, first approximations of water inflow are derived based on observed soil moisture changes. These inflow estimates are used in a series of forward simulations that produce initial estimates of drainage and ET_act, which in turn help improve the estimate of water inflow. Finally, the improved inflow estimates are incorporated into the model and then a parameter optimization is performed using the observed soil moisture as the reference figure. Forward simulations with calibrated soil parameters result in final estimates for ET_act and groundwater recharge. The presented method is applied to an agricultural test site with a crop rotation of cotton and wheat in Punjab, Pakistan. The final model results, with an RMSE of 2.2% in volumetric water content, suggest a cumulative ET_act and groundwater recharge of 769 and 297 mm over a period of 281 days, respectively. The total estimated water inflow accounts for 946 mm, of which 77% originates from irrigation.     

Aspects of Mechanical Behavior and Modeling of a Tropical Unsaturated Soil

The research studies the applicability of two elastoplastic models for the collapse prediction of the lateritic soil profile from Southeastern Brazil. These tropical soils have peculiar geotechnical behavior, due to their mineralogical composition and porous structure coming from intense process of formation. Two elastoplastic models were analyzed: the Barcelona Basic Model (BBM) and another one based on BBM, however developed for tropical soils. Oedometric tests with suction control were performed at three distinct depths of the soil profile. The BBM was not suitable for the upper layer of the soil profile, because BBM considers the compressible behavior of the soil in function of the reduction of the elastoplastic compressibility index with the increase of the matric suction. The model developed for tropical soils showed better suited to the compressible behavior of the soil profile, resulting in good prediction of the collapse potential, mainly by accepting increasing values of the elastoplastic compressibility index of the soil profile with the matric suction rise.     

Angular momentum transport in quasi-Keplerian accretion disks

We reexamine arguments advanced by Hayashi & Matsuda (2001), who claim that several simple, physically motivated derivations based on mean free path theory for calculating the viscous torque in a quasi-Keplerian accretion disk yield results that are inconsistent with the generally accepted model. If correct, the ideas proposed by Hayashi & Matsuda would radically alter our understanding of the nature of the angular momentum transport in the disk, which is a central feature of accretion disk theory. However, in this paper we point out several fallacies in their arguments and show that there indeed exists a simple derivation based on mean free path theory that yields an expression for the viscous torque that is proportional to the radial derivative of the angular velocity in the accretion disk, as expected. The derivation is based on the analysis of the epicyclic motion of gas parcels in adjacent eddies in the disk.     

An 8-year record of gas geochemistry and isotopic composition of methane during baseline sampling at a groundwater observation well in Alberta (Canada)

Variability in baseline groundwater methane concentrations and isotopic compositions was assessed while comparing free and dissolved gas sampling approaches for a groundwater monitoring well in Alberta (Canada) over an 8-year period. Methane concentrations in dissolved gas samples (n?=?12) were on average 4,380?±?2,452 μg/L, yielding a coefficient of variation (CV) >50 %. Methane concentrations in free gas samples (n?=?12) were on average 228,756?±?62,498 ppm by volume, yielding a CV of 27 %. Quantification of combined sampling, sample handling and analytical uncertainties was assessed via triplicate sampling (CV of 19 % and 12 % for free gas and dissolved gas methane concentrations, respectively). Free and dissolved gas samples yielded comparable methane concentration patterns and there was evidence that sampling operations and pumping rates had a marked influence on the obtained methane concentrations in free gas. δ¹³C_CH4 and δ²H_CH4 values of methane were essentially constant (?78.6?±?1.3 and ?300?±?3?‰, respectively) throughout the observation period, suggesting that methane was derived from the same biogenic source irrespective of methane concentration variations. The isotopic composition of methane constitutes a robust and highly valuable baseline parameter and increasing δ¹³C_CH4 and δ²H_CH4 values during repeat sampling may indicate influx of thermogenic methane. Careful sampling and analytical procedures with identical and repeatable approaches are required in baseline-monitoring programs to generate methane concentration and isotope data for groundwater that can be reliably compared to repeat measurements once potential impact from oil and gas development, for example, may occur.     

Highly siderophile and chalcogen element constraints on the origin of components of the Allende and Murchison meteorites

¹⁸⁷Re‐¹⁸⁷Os systematics, abundances of highly siderophile elements (HSE: Re, PGE, and Au), chalcogen elements (Te, Se, and S), and some major and minor elements were determined in physically separated components of the Allende (CV3) and Murchison (CM2) carbonaceous chondrites. Substantial differences exist in the absolute and relative abundances of elements in the components, but the similarity of calculated and literature bulk rock abundances of HSE and chalcogens indicate that chemical complementarity exists among the components, with CI chondrite‐like ratios for many elements. Despite subsequent alteration and oxidation, the overall cosmochemical behavior of most moderately to highly siderophile elements during high‐temperature processing has been preserved in components of Allende at the sampling scale of the present study. The ¹⁸⁷Re‐¹⁸⁷Os systematics and element variations of Allende are less disturbed compared with Murchison, which reflects different degrees of oxidation and alteration of these meteorites. The HSE systematics (with the exception of Au) is controlled by two types of materials: Pd‐depleted condensates and CI chondrite‐like material. Enrichment and heterogeneous distribution of Au among the components is likely the result of hydrothermal alteration. Chalcogen elements are depleted compared with HSE in all components, presumably due to their higher volatility. Small systematic variations of S, Se, and Te in components bear the signature of fractional condensation/partial evaporation and metal–sulfide–silicate partitioning.     

A Rapid Decision Support Tool for Estimating Impacts of a Vadose Zone Volatile Organic Compound Source on Groundwater and Soil Gas

Diminishing rates of subsurface volatile contaminant removal by soil vapor extraction (SVE) oftentimes warrants an in-depth performance assessment to guide remedy decision-making processes. Such a performance assessment must include quantitative approaches to better understand the impact of remaining vadose zone contamination on soil gas and groundwater concentrations. The spreadsheet-based Soil Vapor Extraction Endstate Tool (SVEET) software functionality has recently been expanded to facilitate quantitative performance assessments. The updated version, referred to as SVEET2, includes expansion of the input parameter ranges for describing a site (site geometry, source characteristics, etc.), an expanded list of contaminants, and incorporation of elements of the Vapor Intrusion Estimation Tool for Unsaturated-zone Sources software to provide soil gas concentration estimates for use in vapor intrusion evaluation. As part of the update, SVEET2 was used to estimate the impact of a tetrachloroethene (PCE) vadose zone source on groundwater concentrations, comparing SVEET2 results to field-observed values at an undisclosed site where SVE was recently terminated. PCE concentrations from three separate monitoring wells were estimated by SVEET2 to be within the range of 6.0–6.7 μg/L, as compared to actual field concentrations that ranged from 3 to 11 μg/L PCE. These data demonstrate that SVEET2 can rapidly provide representative quantitative estimates of impacts from a vadose zone contaminant source at field sites. In the context of the SVE performance assessment, such quantitative estimates provide a basis to support remedial and/or regulatory decisions regarding the continued need for vadose zone volatile organic compound remediation or technical justification for SVE termination, which can significantly reduce the cost to complete for a site.     

New Compositional and Structural Constraints on the Smithsonian Microanalytical Reference Materials: Amphiboles from Kakanui and Arenal

The amphiboles from Kakanui and Arenal are two natural minerals that have been used worldwide as microanalytical reference materials, but their compositions and crystal structures are still poorly constrained. In this paper, we report new data on H₂O and trace element mass fractions and single-crystal structural refinement of these two amphiboles. H₂O mass fractions of the Kakanui and Arenal amphiboles determined via Karl-Fischer titration are 0.92 ± 0.18 (2s) and 1.56 ± 0.22% m/m (2s), respectively; these values estimated based on crystal-structure refinement are 0.86 and 1.46% m/m, respectively. Trace element mass fractions measured via LA-ICP-MS in two laboratories are in good agreement, and spots from five fragments for both Kakanui and Arenal amphiboles are generally consistent within reproducibility precision (2s). Our measurements indicate a better homogeneity for the amphiboles from Kakanui than that from Arenal. According to the latest scheme for amphibole classification and nomenclature (Hawthorne et al. 2012), the sample from Arenal is a (partially dehydrogenated) pargasite, and that from Kakanui is a kaersutite. The significant amount of oxo-component and that ^CTi⁴⁺ content is strongly ordered at the M(1) site for both amphiboles indicate crystallisation under high fO₂ conditions.     

Influence of inclination and permeability of solitary woody riparian plants on local hydraulic and sedimentary processes

The role of solitary woody riparian plants with respect to local erosion and deposition of sediments is investigated. A focus is laid on the characteristics ‘inclination’ and ‘permeability’ of the plant's projected frontal area. Therefore, two experimental studies using cylindrical obstacles were carried out in a laboratory flume, one aiming at inclination, the other at permeability. The first series revealed that the total amount of mobilized sediment around the cylinder on average decreased by 8–10% per 5° increasing inclination in streamwise direction. Locations of maximum scour depth simultaneously shifted downstream. A horseshoe vortex system, causing the frontal and lateral scouring, ceased to exist below inclinations of 25–30°. The second series revealed that with increasing permeability, frontal scour incision is delayed, and the eroded sediment volume is significantly reduced. With permeable obstacles, two system states were observed: first, frontal scouring with leeside deposition at higher flow velocities and, second, moderate leeside scouring at lower flow velocities. For up‐scaling and comparison, a field study focussing on fluvial obstacle marks at poplars and willows in secondary channels of the River Loire was additionally conducted. A modified analytical model enabled us to quantify the amount of deposited sediments leeside of the plants. Leeside sediment ridges are significantly stabilized and have a higher preservation potential when covered by pioneer vegetation. Under such conditions, they may indeed induce the development of stable islands. Eventually, ‘sediment ridge width’ turned out to be a suitable indicator for leeside deposited sediment volume, irrespective of spatial scale. Copyright © 2012 John Wiley & Sons, Ltd.     

Diopatra neapolitana and Diopatra marocensis from the Portuguese coast: Morphological and genetic comparison

This paper reports the presence of Diopatra marocensis in European waters, for which Diopatra neapolitana was the only species recognized until recently. Both species coexist in transitional waters, where D. marocensis may be mistaken for young specimens of D. neapolitana. The population of D. marocensis studied in the coastal shelf can be traced back to 1997 and is increasing in density, apparently benefiting from a local anthropogenic organic enrichment source.     

In situ experimental study of reed leaf decomposition along a full salinity gradient

An experimental study on Phragmites australis leaf litter decomposition was conducted in the estuarine environment, Ria de Aveiro, Western Portugal, using the leaf-bag technique, with fine- (1 mm) and coarse-mesh (5 mm) bags. The leaf bags were placed in the field sites at day 0, covering a complete salinity gradient, and replicates were collected over time, at days 3 (leaching), 7, 15, 30 and 60. The biomass loss through the leaching phase, about 20% of the initial leaf mass, was independent of both the salinity and the bag mesh size. The biomass decay pattern along the salinity gradient varied through time and presented strong similarities between the two mesh sizes, with the remaining biomass always lower in the 5 mm mesh-size bags. At days 7 and 15, the lowest remaining biomass was observed at the head of the estuary, the preferential distribution area of P. australis. At day 30, the remaining biomass was higher in the marine area and diminished under a direct relationship with salinity, reaching the lowest value in the freshwater environment, with values ranging from 66% to 44% of the initial weight in 5 mm bags, and from 79% to 51% in 1 mm bags. The largest heterogeneity in the remaining biomass among the study areas positioned along the salinity gradient was found close to days 30 (5 mm) and 40 (1 mm). The overall results indicate that the relationship between leaf decay rate and salinity depends on the decay time considered (k₁₅, k₃₀ or k₆₀) and, for the later stages (k₆₀), also on the leaf-bag mesh size. This implies that the use of leaf litter decay rates as a functional indicator in transitional waters will need to take into consideration the factor location in the salinity gradient and leaf litter stage at which the decay rate is determined. The differences between the decay rates with the mesh size acted mainly at the level of the absolute k value and not at the level of the pattern along the salinity gradient. Even so, the data obtained at the mouth of the estuary, in the area closest to a fully marine environment, indicated that after the initial biomass loss through leaching, P. australis decayed either very slowly, in the 5 mm, or not at all, in the 1 mm mesh bags.