Lacunarity of the Spatial Distributions of Soil Types in Europe.

Lacunarity as a means of quantifying textural properties of spatial distributions suggests a classification into three main classes of the most abundant soils that cover 92% of Europe. Soils with a well-defined self-similar structure of the linear class are related to widespread spatial patterns that are nondominant but ubiquitous at continental scale. Fractal techniques have been increasingly and successfully applied to identify and describe spatial patterns in natural sciences. However, objects with the same fractal dimension can show very different optical properties because of their spatial arrangement. This work focuses primary attention on the geometrical structure of the geographical patterns of soils in Europe. We made use of the European Soil Database to estimate lacunarity indexes of the most abundant soils that cover 92% of the surface of Europe and investigated textural properties of their spatial distribution. We observed three main classes corresponding to three different patterns that displayed the graphs of lacunarity functions, that is, linear, convex, and mixed. They correspond respectively to homogeneous or self-similar, heterogeneous or clustered and those in which behavior can change at different ranges of scales. Finally, we discuss the pedological implications of that classification.

Morphological Functions to Quantify Three-Dimensional Tomograms of Macropore Structure in a Vineyard Soil with Two Different Management Regimes

Soil tomography and morphological functions built over Minkowski functionals were used to describe the impact on pore structure of two soil management practices in a Mediterranean vineyard. Soil structure controls important physical and biological processes in soil–plant–microbial systems. Those processes are dominated by the geometry of soil pore structure, and a correct model of this geometry is critical for understanding them. Soil tomography has been shown to provide rich three-dimensional digital information on soil pore geometry. Recently, mathematical morphological techniques have been proposed as powerful tools to analyze and quantify the geometrical features of porous media. Minkowski functionals and morphological functions built over Minkowski functionals provide computationally efficient means to measure four fundamental geometrical features of three-dimensional geometrical objects, that is, volume, boundary surface, mean boundary surface curvature, and connectivity. We used the threshold and the dilation and erosion of three-dimensional images to generate morphological functions and explore the evolution of Minkowski functionals as the threshold and as the degree of dilation and erosion changes. We analyzed the three-dimensional geometry of soil pore space with X-ray computed tomography (CT) of intact soil columns from a Spanish Mediterranean vineyard by using two different management practices (conventional tillage versus permanent cover crop of resident vegetation). Our results suggested that morphological functions built over Minkowski functionals provide promising tools to characterize soil macropore structure and that the evolution of morphological features with dilation and erosion is more informative as an indicator of structure than moving threshold for both soil managements studied.

Lithologic Control on the Scaling Properties of the First-Order Streams of Drainage Networks: A Monofractal Analysis

The movement of water through the landscape can be investigated at different scales. This study dealt with the interrelation between bedrock lithology and the geometry of the overlying drainage systems. Parameters of fractal analysis, such as fractal dimension and lacunarity, were used to measure and quantify this relationship. The interrelation between bedrock lithology and the geometry of the drainage systems has been widely studied in the last decades. The quantification of this linkage has not yet been clearly established. Several studies have selected river basins or regularly shaped areas as study units, assuming them to be lithologically homogeneous. This study considered irregular distributions of rock types, establishing areas of the soil map (1:25,000) with the same lithologic information as study units. The tectonic stability and the low climatic variability of the study region allowed effective investigation of the lithologic controls on the drainage networks developed on the plutonic rocks, the metamorphic rocks, and the sedimentary materials existing in the study area. To exclude the effect of multiple in- and outflows in the lithologically homogeneous units, we focused this study on the first-order streams of the drainage networks. The geometry of the hydrologic features was quantified through traditional metrics of fluvial geomorphology and scaling parameters of fractal analysis, such as the fractal dimension, the reference density, and the lacunarity. The results demonstrate the scale invariance of both the drainage networks and the set of first-order streams at the study scale and a relationship between scaling in the lithology and the drainage network.

On the generative equations of fractal self-similarity in granular media and the related PSD models

The physical appearance of granular media suggests the existence of geometrical scale invariance. The paper discuss how this physico-empirical property can be mathematically encoded leading to different generative models: a smooth one encoded by a differential equation and another encoded by an equation coming from a measure theoretical property.

Scaling in Soil and Other Complex Porous Media

Scaling is becoming an increasingly important topic in the earth and environmental sciences as researchers attempt to understand complex natural systems through the lens of an ever-increasing set of methods and scales. The guest editors introduce the papers in this issue’s special section and present an overview of some of the work being done. Scaling remains one of the most challenging topics in earth and environmental sciences, forming a basis for our understanding of process development across the multiple scales that make up the subsurface environment. Tremendous progress has been made in discovery, explanation, and applications of scaling. And yet much more needs to be done and is being done as part of the modern quest to quantify, analyze, and manage the complexity of natural systems. Understanding and succinct representation of scaling properties can unveil underlying relationships between system structure and response functions, improve parameterization of natural variability and heterogeneity, and help us address societal needs by effectively merging knowledge acquired at different scales.

Heated fiber optic distributed temperature sensing for measuring soil volumetric heat capacity and water content: A dual probe heat-pulse approach

The first feasibility study of using dual-probe heated fiber optics with distributed temperature sensing to measure soil volumetric heat capacity and soil water content is presented. Although results using different combinations of cables demonstrate feasibility, further work is needed to gain accuracy, including a model to account for the finite dimension and the thermal influence of the probes. Implementation of the dual-probe heat-pulse (DPHP) approach for measurement of volumetric heat capacity (C) and water content (θ) with distributed temperature sensing heated fiber optic (FO) systems presents an unprecedented opportunity for environmental monitoring (e.g., simultaneous measurement at thousands of points). We applied uniform heat pulses along a FO cable and monitored the thermal response at adjacent cables. We tested the DPHP method in the laboratory using multiple FO cables at a range of spacings. The amplitude and phase shift in the heat signal with distance was found to be a function of the soil volumetric heat capacity. Estimations of C at a range of moisture contents (θ = 0.09– 0.34 m3 m−3) suggest the feasibility of measurement via responsiveness to the changes in θ, although we observed error with decreasing soil water contents (up to 26% at θ = 0.09 m3 m−3). Optimization will require further models to account for the finite radius and thermal influence of the FO cables. Although the results indicate that the method shows great promise, further study is needed to quantify the effects of soil type, cable spacing, and jacket configurations on accuracy.

Modeling soil processes : Review, key challenges, and new perspectives

Peer reviewed

Comparison among monitoring strategies to assess water flow dynamic and soil hydraulic properties in agricultural soils

Irrigated agriculture is usually performed in semi-arid regions despite scarcity of water resources. Therefore, optimal irrigation management by monitoring the soil is essential, and assessing soil hydraulic properties and water flow dynamics is presented as a first measure. For this purpose, the control of volumetric water content, θ, and pressure head, h, is required. This study adopted two types of monitoring strategies in the same experimental plot to control θ and h in the vadose zone: i) non-automatic and more time-consuming; ii) automatic connected to a datalogger. Water flux was modelled with Hydrus-1D using the data collected from both acquisition strategies independently (3820 daily values for the automatic; less than 1000 for the non-automatic). Goodness-of-fit results reported a better adjustment in case of automatic sensors. Both model outputs adequately predicted the general trend of θ and h, but with slight differences in computed annual drainage (711 mm and 774 mm). Soil hydraulic properties were inversely estimated from both data acquisition systems. Major differences were obtained in the saturated volumetric water content, θs, and the n and α van Genuchten model shape parameters. Saturated hydraulic conductivity, Ks, shown lower variability with a coefficient of variation range from 0.13 to 0.24 for the soil layers defined. Soil hydraulic properties were better assessed through automatic data acquisition as data variability was lower and accuracy was higher.

Quantifying effects of soil heterogeneity on groundwater pollution in four sites in the USA

Four sites located in the north-eastern region of the United States of America have been chosen to investigate the impacts of soil heterogeneity in the transport of solutes (bromide and chloride) through the vadose zone (the zone in the soil that lies below the root zone and above the permanent saturated groundwater). A recently proposed mathematical model based on the cumulative beta distribution has been deployed to compare and contrast the regions' heterogeneity from multiple sample percolation experiments. Significant differences in patterns of solute leaching were observed even over a small spatial scale, indicating that traditional sampling methods for solute transport, for example the gravity pan or suction lysimeters, or more recent inventions such as the multiple sample percolation systems may not be effective in estimating solute fluxes in soils when a significant degree of soil heterogeneity is present. Consequently, ignoring soil heterogeneity in solute transport studies will likely result in under- or overprediction of leached fluxes and potentially lead to serious pollution of soils and/or groundwater. The cumulative beta distribution technique is found to be a versatile and simple technique of gaining valuable information regarding soil heterogeneity effects on solute transport. It is also an excellent tool for guiding future decisions of experimental designs particularly in regard to the number of samples within one site and the number of sampling locations between sites required to obtain a representative estimate of field solute or drainage flux.

Aeration for plant root respiration in a tidal marsh

This paper investigates the tidal effects on aeration conditions for plant root respiration in a tidal marsh. We extend the work of Ursino et al. ( 2004) by using a two-phase model for air and water flows in the marsh. Simulations have been conducted to examine directly the link between the airflow dynamics and the aeration condition in the marsh soil. The results show that the effects of entrapped air on water movement in the vadose zone are significant in certain circumstances. Single-phase models based on Richards' equation, which neglect such effects, may not be adequate for quantifying the aeration condition in tidal marsh. The optimal aeration condition, represented by the maximum of the integral magnitude of tidally advected air mass ( TAAM) flux, is found to occur near the tidal creek for the four soil textures simulated. This may explain the observation that some salt marsh plant species grow better near tidal creeks than in the inner marsh areas. Our analyses, based on the two-phase model and predicted TAAM flux magnitude, provide further insight into the positive feedback'' mechanism proposed by Ursino et al. ( 2004). That is, pioneer plants may grow successfully near the creek where the root aeration condition is optimal. The roots of the pioneer plants can soften and loosen the rhizosphere soil, which increases the evapotranspiration rate, the soil porosity, and absolute permeability and weakens the capillary effects. These, in turn, improve further the root aeration conditions and may lead to colonization by plants less resistant to anaerobic conditions.

Artificial pumping errors in the Kool-Parker scaling model of soil moisture hysteresis

Hysteresis models that eliminate the artificial pumping errors associated with the Kool-Parker (KP) soil moisture hysteresis model, such as the Parker-Lenhard (PL) method, can be computationally demanding in unsaturated transport models since they need to retain the wetting-drying history of the system. The pumping errors in these models need to be eliminated for correct simulation of cyclical systems (e.g. transport above a tidally forced watertable, infiltration and redistribution under periodic irrigation) if the soils exhibit significant hysteresis. A modification is made here to the PL method that allows it to be more readily applied to numerical models by eliminating the need to store a large number of soil moisture reversal points. The modified-PL method largely eliminates any artificial pumping error and so essentially retains the accuracy of the original PL approach. The modified-PL method is implemented in HYDRUS-1D (version 2.0), which is then used to simulate cyclic capillary fringe dynamics to show the influence of removing artificial pumping errors and to demonstrate the ease of implementation. Artificial pumping errors are shown to be significant for the soils and system characteristics used here in numerical experiments of transport above a fluctuating watertable. (c) 2005 Elsevier B.V. All rights reserved.

Sea level rise and South Florida coastal forests

Coastal ecosystems lie at the forefront of sea level rise. We posit that before the onset of actual inundation, sea level rise will influence the species composition of coastal hardwood hammocks and buttonwood (Conocarpus erectus L.) forests of the Everglades National Park based on tolerance to drought and salinity. Precipitation is the major water source in coastal hammocks and is stored in the soil vadose zone, but vadose water will diminish with the rising water table as a consequence of sea level rise, thereby subjecting plants to salt water stress. A model is used to demonstrate that the constraining effect of salinity on transpiration limits the distribution of freshwater-dependent communities. Field data collected in hardwood hammocks and coastal buttonwood forests over 11 years show that halophytes have replaced glycophytes. We establish that sea level rise threatens 21 rare coastal species in Everglades National Park and estimate the relative risk to each species using basic life history and population traits. We review salinity conditions in the estuarine region over 1999–2009 and associate wide variability in the extent of the annual seawater intrusion to variation in freshwater inflows and precipitation. We also examine species composition in coastal and inland hammocks in connection with distance from the coast, depth to water table, and groundwater salinity. Though this study focuses on coastal forests and rare species of South Florida, it has implications for coastal forests threatened by saltwater intrusion across the globe.

Water source utilization and foliar nutrient status differs between upland and flooded plant communities in wetland tree islands

Tree islands in the Everglades wetlands are centers of biodiversity and targets of restoration, yet little is known about the pattern of water source utilization by the constituent woody plant communities: upland hammocks and flooded swamp forests. Two potential water sources exist: (1) entrapped rainwater in the vadose zone of the organic soil (referred to as upland soil water), that becomes enriched in phosphorus, and (2) phosphorus-poor groundwater/surface water (referred to as regional water). Using natural stable isotope abundance as a tracer, we observed that hammock plants used upland soil water in the wet season and shifted to regional water uptake in the dry season, while swamp forest plants used regional water throughout the year. Consistent with the previously observed phosphorus concentrations of the two water sources, hammock plants had a greater annual mean foliar phosphorus concentration over swamp forest plants, thereby supporting the idea that tree island hammocks are islands of high phosphorus concentrations in the oligotrophic Everglades. Foliar nitrogen levels in swamp forest plants were higher than those of hammock plants. Linking water sources with foliar nutrient concentrations can indicate nutrient sources and periods of nutrient uptake, thereby linking hydrology with the nutrient regimes of different plant communities in wetland ecosystems. Our results are consistent with the hypotheses that (1) over long periods, upland tree island communities incrementally increase their nutrient concentration by incorporating marsh nutrients through transpiration seasonally, and (2) small differences in micro-topography in a wetland ecosystem can lead to large differences in water and nutrient cycles.

The Influence of Vegetation on the Hydrodynamics and Geomorphology of a Tree Island in Everglades National Park (Florida, United States)

Transpiration-driven nutrient accumulation has been identified as a potential mechanism governing the creation and maintenance of wetland vegetation patterning. This process may contribute to the formation of nutrient-rich tree islands within the expansive oligotrophic marshes of the Everglades (Florida, United States). This study presents hydrogeochemical data indicating that tree root water uptake is a primary driver of groundwater ion accumulation across one of these islands. Sap flow, soil moisture, water level, water chemistry, and rainfall were measured to identify the relationships between climate, transpiration, and groundwater uptake by phreatophytes and to examine the effect this uptake has on groundwater chemistry and mineral formation in three woody plant communities of differing elevations. During the dry season, trees relied more on groundwater for transpiration, which led to a depressed water table and the advective movement of groundwater and dissolved ions, including phosphorus, from the surrounding marsh towards the centre of the island. Ion exclusion during root water uptake led to elevated concentrations of all major dissolved ions in the tree island groundwater compared with the adjacent marsh. Groundwater was predominately supersaturated with respect to aragonite and calcite in the lower-elevation woody communities, indicating the potential for soil formation. Elevated groundwater phosphorous concentrations detected in the highest-elevation woody community were associated with the leaching of inorganic sediments (i.e. hydroxyapatite) in the vadose zone. Understanding the complex feedback mechanisms regulating plant/groundwater/surface water interactions, nutrient dynamics, and potential soil formation is necessary to manage and restore patterned wetlands such as the Everglades.

Espessura da zona de transição do manto sob a Província Borborema, NE Brasil, através de conversões Ps telessísmicas

The mantle transition zone is defined by two seismic discontinuities, nominally at 410 and 660 km depth, which result from transformations in the mineral olivine. The topography of these discontinuities provides information about lateral temperature changes in the transition zone. In this work, P-to-S conversions from teleseismic events recorded at 32 broadband stations in the Borborema Province were used to determine the transition zone thickness beneath this region and to investigate whether there are lateral temperature changes within this depth range. For this analysis, stacking and migration of receiver functions was performed. In the Borborema Province, geophysical studies have revealed a geoid anomaly which could reflect the presence of a thermal anomaly related to the origin of intraplate volcanism and uplift that marked the evolution of the Province in the Cenozoic. Several models have been proposed to explain these phenomena, which include those invoking the presence of a deep-seated mantle plume and those invoking shallower sources, such as small-scale convection cells. The results of this work show that no thermal anomalies are present at transition zone depths, as significant variations in the transition zone thickness were not observed. However, regions of depressed topography for both discontinuities (410 and 660 km) that approximately overlap in space were identified, suggesting that lower-thanaverage, lateral variations in seismic velocity above 410 km depth may exist below the the Borborema Province. This is consistent with the presence of a thermally-induced, low-density body independently inferred from analysis of geoid anomalies. Therefore, the magma source responsible for the Cenozoic intraplate volcanism and related uplift in the Province, is likely to be confined above the upper mantle transition zone.