Uniform and lateral preferential flows under flood irrigation at field scale

Flood irrigation is globally one of the most used irrigation methods. Typically, not all water that is applied during flood irrigation is consumed by plants or lost to evaporation. Return flow, the portion of applied water from flood irrigation that returns back to streams either via surface or subsurface flow, can constitute a large part of the water balance. Few studies have addressed the connection between vertical and lateral subsurface flows and its potential role in determining return flow pathways due to the difficulty in observing and quantifying these processes at plot or field scale. We employed a novel approach, combining induced polarization, time‐lapse electrical resistivity tomography, and time‐lapse borehole nuclear magnetic resonance, to identify flow paths and quantify changes in soil hydrological conditions under nonuniform application of flood irrigation water. We developed and tested a new method to track the wetting front in the subsurface using the full range of inverted resistivity values. Antecedent soil moisture conditions did not play an important role in preferential flow path activation. More importantly, boundaries between lithological zones in the soil profile were observed to control preferential flow pathways with subsurface run‐off occurring at these boundaries when saturation occurred. Using the new method to analyse time‐lapse resistivity measurements, we were able to track the wetting front and identify subsurface flow paths. Both uniform infiltration and preferential lateral flows were observed. Combining three geophysical methods, we documented the influence of lithology on subsurface flow processes. This study highlights the importance of characterizing the subsurface when the objective is to identify and quantify subsurface return flow pathways under flood irrigation.     

Conditioning geological surfaces to horizontal wells

Computational Geosciences - Kriging is a standard method for conditioning surfaces to observations. Kriging works for vertical wells, but may produce surfaces that cross horizontal wells between...     

Airborne laser altimetry survey of Glaciar Tyndall, Patagonia

The first airborne laser altimetry measurements of a glacier in South America are presented. Data were collected in November of 2001 over Glaciar Tyndall, Torres del Paine National Park, Chilean Patagonia, onboard a Twin Otter airplane of the Chilean Air Force. A laser scanner with a rotating polygon-mirror system together with an Inertial Navigation System (INS) were fixed to the floor of the aircraft, and used in combination with two dual-frequency GPS receivers. Together, the laser–INS–GPS system had a nominal accuracy of 30 cm after data processing. On November 23rd, a total of 235 km were flown over the ablation area of Glaciar Tyndall, with 5 longitudinal tracks with a mean swath width of 300 m, which results in a point spacing of approximately 2 m both along and across track. A digital elevation model (DEM) generated using the laser altimetry data was compared with a DEM produced from a 1975 map (1:50,000 scale — Instituto Geográfico Militar (IGM), Chile). A mean thinning of − 3.1 ± 1.0 m a^− 1 was calculated for the ablation area of Glaciar Tyndall, with a maximum value of − 7.7 ± 1.0 m a^− 1 at the calving front at 50 m a.s.l. and minimum values of between − 1.0 and − 2.0 ± 1.0 m a^− 1 at altitudes close to the equilibrium line altitude (900 m a.s.l.). The thinning rates derived from the airborne survey were similar to the results obtained by means of ground survey carried out at  600 m of altitude on Glaciar Tyndall between 1975 and 2002, yielding a mean thinning of − 3.2 m a^− 1 [Raymond, C., Neumann, T.A., Rignot, E., Echelmeyer, K.A., Rivera, A., Casassa, G., 2005. Retreat of Tyndall Glacier, Patagonia, over the last half century. Journal of Glaciology 173 (51), 239–247.]. A good agreement was also found between ice elevation changes measured with laser data and previous results obtained with Shuttle Radar Topography Mission (SRTM) data. We conclude that airborne laser altimetry is an effective means for accurately detecting glacier elevation changes in Patagonia, where an ice thinning acceleration trend has been observed during recent years, presumably in response to warming and possibly also drier conditions.     

Leaky savannas: the significance of lateral carbon fluxes in the seasonal tropics

Globally, dissolved inorganic carbon (DIC) accounts for more than half the annual flux of carbon exported from terrestrial ecosystems via rivers. Here, we assess the relative influences of biogeochemical and hydrological processes on DIC fluxes exported from a tropical river catchment characterized by distinct land cover, climate and geology transition from the wet tropical mountains to the low‐lying savanna plains. Processes controlling changes in river DIC were investigated using dissolved organic carbon, particulate organic carbon and DIC concentrations and stable isotope ratios of DIC (δ¹³C_DIC) at two time scales: seasonal and diel. The recently developed Isotopic Continuous Dissolved Inorganic Carbon Analyser was used to measure diel DIC concentration and δ¹³C_DIC changes at a 15‐min temporal resolution. Results highlight the predominance of biologically mediated processes (photosynthesis and respiration) controlling diel changes in DIC. These resulted in DIC concentrations varying between 3.55 and 3.82 mg/l and δ¹³C_DIC values ranging from ?19.7 ± 0.31‰ to ?17.1 ± 0.08‰. In contrast, at the seasonal scale, we observed wet season DIC variations predominantly from mixing processes and dry season DIC variations due to both mixing processes and biological processes. The observed wet season increases in DIC concentrations (by 6.81 mg/l) and δ¹³C_DIC values of river water (by 5.4‰) largely result from proportional increases in subsurface inflows from the savanna plains (C₄ vegetation) region relative to inflows from the rainforest (C₃ vegetation) highlands. The high DIC river load during the wet season resulted in the transfer of 97% of the annual river carbon load. Therefore, in this gaining river, there are significant seasonal variations in both the hydrological and carbon cycles, and there is evidence of substantial coupling between the carbon cycles of the terrestrial and the fluvial environments. Recent identification of a substantial carbon sink in the savannas of northern Australia during wetter years in the recent past does not take into account the possibility of a substantial, rapid, lateral flux of carbon to rivers and back to the atmosphere. Copyright © 2015 John Wiley & Sons, Ltd.     

Kriging with Inequality Constraints

A Gaussian random field with an unknown linear trend for the mean is considered. Methods for obtaining the distribution of the trend coefficients given exact data and inequality constraints are established. Moreover, the conditional distribution for the random field at any location is calculated so that predictions using e.g. the expectation, the mode, or the median can be evaluated and prediction error estimates using quantiles or variance can be obtained. Conditional simulation techniques are also provided.     

Indication of a geomagnetic low-inclination excursion in supposed Middle Weichselian interstadial marine clay at Rubjerg,Denmark

Palaeomagnetic measurements have been made on specimens from Late Pleistocene sediments of a piston boring at Rubjerg in Vendsyssel, northern Jutland. The stratigraphy of the deposits is based on content of foraminifera. A total of 70 relatively oriented specimens were investigated palaeomagnetically. Normal steep inclinations close to that of the axial dipole field were found in the Upper Saxicava Sand and in the Younger Yoldia Clay (radiocarbon dated at 14,650 ± 190?12,650 ± 180 B.P.), and a secular variation with an amplitude of 10–12° in the inclincation and a “period” roughly estimated at about 350–400 years was found in the Younger Yoldia Clay.Seventeen relatively oriented specimens from undisturbed older marine deposits revealed a stable low inclination of 11° with α₉₅ = 3°. The age of this apparent geomagnetic excursion falls somewhere between 23,000 and 40,000 B.P. (Older Yoldia Clay). Among other known geomagnetic excursions and events within this interval are Laschamp in France, Mono Lake in California and Lake Mungo in Australia. Until more definite ages have been obtained, the excursion is provisionally named the “Rubjerg Excursion”.     

Petro-and paleomagnetic investigations of Tuzla section sediments (krasnodarsk territory)

Petro-and paleomagnetic methods are applied to the study of the upper part of the Late Pleistocene Tuzla section (Azov coast of the Taman Peninsula) composed of continental sediments and dated at 50–10 ka. The detailed curves of the angular components of the geomagnetic field obtained in this study display an anomalous direction coinciding in time (～25–35 ka) with an anomalous horizon discovered in rocks of the Roxolany section (Ukraine). According to the world time scale of geomagnetic excursions, the anomalous direction correlates with the Mono Lake excursion. A significant correlation between the time series NRM_0.015/SIRM (Tuzla section) and NRM₂₅₀/KB (Roxolany section) in the interval 50–10 ka and the world composite curves VADM-21 and Sint-800 implies that, in this time interval, the curve NRM_0.015/SIRM reflects the variation in the relative paleointensity of the geomagnetic field.     

Transient heat flow in a stratified medium

The equation of heat conduction is solved for a horizontally stratified medium initially at constant temperature and subject to a step-function surface temperature change. The problem is solved by Laplace transformation and by applying a concept of further stratifications of the medium into unitary layers of constant ratio of thickness to the square root of thermal diffusivity. The solution takes the form of an infinite sum of complementary error functions, with coefficients given by recursion relations, is suitable for numerical applications and offers an attractive alternative to harmonic and quasi transient approaches in calculating the penetration of transient surface temperature variations into a layered medium. Use of the theorem of superposition yields a general expression for an arbitrary surface temperature function. The utility of the theory is illustrated by modelling examples of palaeoclimatically induced subsurface temperature and heat flow perturbations.