Jonkershoek: Africa's oldest catchment experiment - 80 years and counting

The Jonkershoek Forestry Research Station was established in the Mediterranean climate region of South Africa in 1935 to implement a multiple catchment experiment to determine the effects of afforestation on water yield. The experiment consists of six neighbouring catchments previously supporting indigenous fynbos shrublands, five of which were sequentially afforested with Pinus radiata plantations every 8 years from 1940 to 1980 and one kept as the control. They conclusively demonstrated a significant impact of afforestation on streamflow. The treatment catchments have seen subsequent plantation rotations since 1980 and rainfall and streamflow observations have been continued to date. Here we describe the site, experimental design, rainfall and streamflow records, instrumentation, and how to access the data. We also provide details of recently installed instrumentation, including full weather stations, fog gauges, and an eddy covariance flux tower. The Jonkershoek catchment experiment was the core of a globally significant interdisciplinary research programme (1935–1992) that shaped most environmental policies and practices relating to biodiversity, forestry, fire and water in South Africa and beyond. The South African Environmental Observation Network (SAEON) inherited the experiment in 2010 and is maintaining it as a long-term eco-hydrological research platform and global change observatory.     

Mapping and analysing virtual outcrops

Laser scanning is a very efficient way to generate realistic, high-resolution digital models of 3-D geological outcrops. This paper discusses the methodologies involved in the creation and analysis of virtual outcrops, based on laser scanner data. The visualisation of the laser scanner data as a photorealistic 3-D object is described. Geological features picked out on the virtual outcrop (e.g. fractures, faults or bedding planes) can be extrapolated outward, into space, and inward, into the subsurface, using tension surfaces.Electronic Supplementary Material Supplementary material is available in the online version of this article at Reviewed by: J.D. Clemens, D. Yuen     

A modified light transmission visualization method for DNAPL saturation measurements in 2-D models

In this research, a light transmission visualization (LTV) method was used to quantify dense non-aqueous phase liquids (DNAPL) saturation in two-dimensional (2-D), two fluid phase systems. The method is an expansion of earlier LTV methods and takes into account both absorption and refraction light theories. Based on this method, DNAPL and water saturations can rapidly be obtained point wise across sand-packed 2-D flow chambers without the need to develop a calibration curve. A single point calibration step is, however, needed when dyed DNAPL is used to account for the change in the transmission factor at the dyed DNAPL–water interface. The method was applied to measure, for the first time, undyed DNAPL saturation in small 2-D chambers. Known amounts of DNAPL, modeled by tetrachloroethylene (PCE), were added to the chamber and these amounts were compared to results obtained by this LTV method. Strong correlation existed between results obtained based on this method and the known PCE amounts with an R² value of 0.993. Similar experiments conducted using dyed PCE showed a stronger correlation between results obtained by this LTV method and the known amounts of dyed PCE added to the chamber with an R² value of 0.999. The method was also used to measure dyed PCE saturation in a large 2-D model following sparging experiments. Results obtained from image analyses following each sparging event were compared to results obtained by two independent techniques, namely gas chromatography–mass spectrometry (GC/MS) analyses and carbon column extraction. There was a good agreement between the results obtained by this LTV method and those obtained by the two independent techniques when experiments were carried out under stable light source conditions and errors in mass balance were minor. The method presented here can be expanded to measure fluid contents in three fluid phase systems and provide a non-destructive, non-intrusive tool to investigate changes in DNAPL architecture and flow characteristics in laboratory experiments.     

The influence of autocorrelation on the ability to detect trend in hydrological series

This study investigated using Monte Carlo simulation the interaction between a linear trend and a lag‐one autoregressive (AR(1)) process when both exist in a time series. Simulation experiments demonstrated that the existence of serial correlation alters the variance of the estimate of the Mann–Kendall (MK) statistic; and the presence of a trend alters the estimate of the magnitude of serial correlation. Furthermore, it was shown that removal of a positive serial correlation component from time series by pre‐whitening resulted in a reduction in the magnitude of the existing trend; and the removal of a trend component from a time series as a first step prior to pre‐whitening eliminates the influence of the trend on the serial correlation and does not seriously affect the estimate of the true AR(1). These results indicate that the commonly used pre‐whitening procedure for eliminating the effect of serial correlation on the MK test leads to potentially inaccurate assessments of the significance of a trend; and certain procedures will be more appropriate for eliminating the impact of serial correlation on the MK test. In essence, it was advocated that a trend first be removed in a series prior to ascertaining the magnitude of serial correlation. This alternative approach and the previously existing approaches were employed to assess the significance of a trend in serially correlated annual mean and annual minimum streamflow data of some pristine river basins in Ontario, Canada. Results indicate that, with the previously existing procedures, researchers and practitioners may have incorrectly identified the possibility of significant trends. Copyright © Environment Canada. Published by John Wiley & Sons, Ltd.     

On the effect of line current width and relative position on the multi-spacecraft curlometer technique

The response of the multi-spacecraft curlometer technique to variations in the size and relative position of infinitely long line currents with radially varying current density is systematically investigated for spacecraft in a regular tetrahedral formation. It is shown that, for line currents with a width less than the spacecraft separation, there is significant variation in the returned current with position of that current within the tetrahedron. For infinitely thin line currents, the curlometer tends to detect approximately 20% of the input current. For increasingly wide line currents there is less variation of the curlometer results with position of the current and the percentage of current magnitude detected increases. When the width of the current system is half the spacecraft separation, the curlometer tends to detect approximately 80% of the input current. These results are discussed in the context of multi-scale, multi-spacecraft missions.     

Spherical harmonic synthesis and least squares computations in satellite gravity gradiometry

The computational requirements in the simulations of geopotential estimation from satellite gravity gradiometry are discussed. Fast algorithms for spherical harmonic synthesis and least squares accumulation on a vectorizing supercomputers are presented. Using these methods, in a test case estimation of 2595 coefficients of a degree and order 50 gravity field, sustained program execution speeds of 275 Mflops (87 % peak machine speed) on a single processor of a CRAY Y-MP were achieved, with spherical harmonics computation accounting for less than 1 % of total cost. From the results, it appears that brute-force estimation of a degree and order 180 field would require 537 Million Words of memory and 85 hours of CPU time, assuming mission duration of 1 month, and execution speed of 1 Gflops. Both memory size and execution speed requirements are within the capabilities of modern multi-processor supercomputers.