Improving the representation of soil moisture by using a semi‐analytical infiltration model

Soil moisture is widely recognized as a fundamental variable governing the mass and energy fluxes between the land surface and the atmosphere. In this study, the soil moisture modelling at sub‐daily timescale is addressed by using an accurate representation of the infiltration component. For that, the semi‐analytical infiltration model proposed by Corradini et al. (1997) has been incorporated into a soil water balance model to simulate the evolution in time of surface and profile soil moisture. The performances of this new soil moisture model [soil water balance module‐semi‐analytical (SWBM‐SA)] are compared with those of a precedent version [SWBM‐Green–Ampt (GA)] where the GA approach was employed. Their capability to reproduce in situ soil moisture observations at three sites in Italy, Spain and France is analysed. Hourly observations of quality‐checked rainfall, temperature and soil moisture data for a 2‐year period are used for testing the modelling approaches. Specifically, different configurations for the calibration and validation of the models are adopted by varying a single parameter, that is, the saturated hydraulic conductivity. Results indicate that both SWBMs are able to reproduce satisfactorily the hourly soil moisture temporal pattern for the three sites with root mean square errors lower than 0.024 m³/m³ both in the calibration and validation periods. For all sites, the SWBM‐SA model outperforms the SWBM‐GA with an average reduction of the root mean square error of ~20%. Specifically, the higher improvement is observed for the French site for which in situ observations are measured at 30 cm depth, and this is attributed to the capability of the SA infiltration model to simulate the time evolution of the whole soil moisture profile. The reasonable models performance coupled with the need to calibrate only a single parameter makes them useful tools for soil moisture simulation in different regions worldwide, also in scarcely gauged areas. Copyright © 2013 John Wiley & Sons, Ltd.     

Isotopic fingerprints of Milankovitch cycles in Pennsylvanian carbonate platform‐top deposits: the Valdorria record,Northern Spain

Sedimentary cyclic sequences deposited during the Late Palaeozoic Ice Age are widespread. Glacio‐eustatic control of the cyclic patterns is commonly accepted, and the durations of the cyclothems generally match the short‐ and long‐eccentricity Milankovitch orbital parameters. Nevertheless, geochemical fingerprints of orbital parameters are poorly documented in deep‐time sedimentary records. Here, we report on well‐exposed Bashkirian cyclothems of c. 123 ka and c. 400 ka durations from the Valdorria platform. The shorter‐term cyclothems can be grouped into longer‐term composite sequences that are consistent with generally accepted durations of c. 125 ka and c. 400 ka for Milankovitch eccentricity cycles. The stratigraphic pattern is mirrored by the isotope geochemical signals, which show distinct recurring trends. These trends are confirmed by statistical tests. Whereas intrinsic factors and/or subaerial exposure related to sea‐level lowstands might have truncated cycle patterns in tectonically stable basins, rapid subsidence of the Valdorria platform's foreland basin appears to have contributed to a faithful recording of cyclothems of different orders. The patterns and biostratigraphic constraint revealed in this study demonstrate the power of orbital forcing in imprinting sedimentary and geochemical signals in the rock record.     

Importance of porosity and transfer of matter in the rock weathering processes: two examples in central Spain

Some physical properties (bulk and free porosity, pore size distribution), and the chemical composition and mass balance of two deeply weathered profiles one developed on Hercynian granodiorite and the other on pre-Cambrian slates were studied. Hydric and mercury porosimetry, nitrogen adsorption techniques, chemical analyses and XRD techniques were used. On granodiorite, weathering has created increased porosity with a pore diameter <5 μm, whereas on slates the weathering has produced of ca. 1 μm in diameter. These pore sizes have played an important role in the weathering processes. Assuming that weathering preserves volumes, except in the uppermost part of the profiles, it brought about a loss of matter of more than 12% (~300 kg/m³) on granodiorite and ca. 30% (~800 kg/m³) on slates. These changes are related to shifts in the mineralogical evolution, with the appearance of new 2:1 and 1:1 phyllosilicates and Fe oxy-hydroxides as the main authigenic minerals. The release of matter, at least since the upper Neogene until the present, has led to the lowering of relief in a more or less homogeneous way, giving rise to gentle hillsides and flat surfaces below which the current river networks are incised. Porosity studies have the potential to explain several specific landforms as well as affecting landscape development in general.     

Radar interferometry techniques for the study of ground subsidence phenomena: a review of practical issues through cases in Spain

Subsidence related to multiple natural and human-induced processes affects an increasing number of areas worldwide. Although this phenomenon may involve surface deformation with 3D displacement components, negative vertical movement, either progressive or episodic, tends to dominate. Over the last decades, differential SAR interferometry (DInSAR) has become a very useful remote sensing tool for accurately measuring the spatial and temporal evolution of surface displacements over broad areas. This work discusses the main advantages and limitations of addressing active subsidence phenomena by means of DInSAR techniques from an end-user point of view. Special attention is paid to the spatial and temporal resolution, the precision of the measurements, and the usefulness of the data. The presented analysis is focused on DInSAR results exploitation of various ground subsidence phenomena (groundwater withdrawal, soil compaction, mining subsidence, evaporite dissolution subsidence, and volcanic deformation) with different displacement patterns in a selection of subsidence areas in Spain. Finally, a cost comparative study is performed for the different techniques applied.     

Approaching water management at watershed scale using distributed water balances in the Yanuncay River basins (Ecuador)

In this study, the excess water generated during a 10-year period (1998–2008) in the upper and middle Yanuncay River basins is estimated. The distributed water balance method, which analyses the interaction between all the different parameters that form part of the water cycle in nature, was developed. To create the model, basic parameters such as rainfall, temperature, soil type and surface cover are required. Moreover, by using computer software such as Microsoft Office Excel and geographic information systems, it is possible to obtain monthly data showing the water excess and generate thematic maps which allow for an effective monitoring of the behaviour of the middle and upper Yanuncay River basins. This study takes the spatial variability of the various factors that influence the behaviour of a basin into account. This allows determining zones with greater water excess, which are areas that would need to be protected. The model can easily incorporate new data such as land use and surface cover (based on time frames) which would facilitate the comparison of different scenarios.     

Simultaneous Solving of Three-Dimensional Gravity Anomalies Caused by Pumping Tests in Unconfined Aquifers

The contribution of geophysical techniques to groundwater characterization has largely evolved during the last two decades. As gravity is a geophysical technique sensitive to underground mass variations and due to the improved resolution of modern gravity meters, it can provide information on changes in the phreatic level caused by water being pumped from unconfined aquifers. Previous studies simulated the hydraulic head and the gravimetric anomaly using independent codes. The mass change associated to the pumping well-simulated drawdown had to be externally transferred to the gravity code used to simulate the gravimetrical anomaly, which has severe drawbacks. This article describes how to solve the forward coupled hydro-gravity problem using a unique finite-element code. To illustrate it, it is shown the case of a two-dimensional hydraulic model coupled to its three-dimensional gravity anomaly and also a more complex case where both related domains are three-dimensional. Both are compared against analytical solutions and discussed. The methodology is very flexible, general and amenable to extensions like including heterogeneous domains or coupling with the inverse problem in the same loop.     

Stalactite drip rate variations controlled by air pressure changes: an example of non‐linear infiltration processes in the ‘Cueva del Agua’ (Spain)

A study of drip water from a stalactite in the Cueva del Agua (Granada, southern Spain) over four hydrological years has enabled a detailed characterization of infiltration through the non‐saturated zone of this cave. The most significant aspects are: (1) The drip water regime is not seasonal, but is linked instead to slow infiltration. Sudden changes in the drip water regime are detected, due to infiltration along the preferential flow paths and the draining of water of supersaturated water reserves from the microfissure and pore system. (2) The accumulation of excess rain in the unsaturated zone dislodges the reserve of supersaturated water from the matrix of microfissures and pores, giving rise to an average increase in drip intensity of 2·2 ± 0·5 mm h^?1 and 37 ± 13 µS cm^?1 in the electrical conductivity of the water. (3) Time‐series analysis of the drip water demonstrates a lack of linearity over time of the drip rate, indicating a chaotic regime. (4) When the dripping is constant, barometric oscillation of the air is the principal factor that causes the lack of linearity in the drip flow. In this way, over an average of 2–3 days, a mean variation of the air pressure inside the cave of 10 ± 3·7 mbar causes a mean oscillation in the drip rate of 0·5 ± 0·2 mm h^?1. This increase in air pressure is translated as an increase in the relative thickness of the gaseous phase of the drip water at the cost of the aqueous phase, so reducing the drip intensity of the stalactite. Copyright © 2006 John Wiley & Sons, Ltd.     

Martian Dust Storms: A Review

A review of the dust storms observed on Mars is made. This includes the seasonal and interannual variability of planet- encircling and regional dust storms. Although there is a significant interannual variability, planet-encircling dust storms have been observed to form during the southern spring and summer seasons, while regional dust storms tend to occur more frequently. Some aspects of possible mechanisms associated with the origin, maintenance and decay of the dust storms are also discussed. This revised version was published online in July 2006 with corrections to the Cover Date.     

Infrared Observations Of Dust Emission From Comet Hale-Bopp

We present infrared imaging and photometry of the bright, giant comet C/1995 O1 (Hale-Bopp). The comet was observed in an extended infrared and optical observing campaign in 1996–1997. The infrared morphology of the comet was observed to change from the 6 to 8 jet “porcupine” structure in 1996 to the “pinwheel” structure seen in 1997; this has implications for the position of the rotational angular momentum vector. Long term light curves taken at 11.3 μm indicate a dust production rate that varies with heliocentric distance as ∶ r^−1.4. Short term light curves taken at perihelion indicate a rotational periodicity of 11.3 hours and a projected dust outflow speed of ∶ 0.4 km s⁻¹. The spectral energy distribution of the dust on October 31, 1996 is well modeled by a mixture of 70% silicaceous and 30% carbonaceous non-porous grains, with a small particle dominated size distribution like that seen for comet P/Halley (McDonnell et al., 1991), an overall dust production rate of 2 × 10⁵ kg s⁻¹, a dust-to-gas ratio of ∶5, and an albedo of 39%. This revised version was published online in July 2006 with corrections to the Cover Date.