Multivariate statistical analysis of hydrochemical data to assess alluvial aquifer–stream connectivity during drought and flood : Cressbrook Creek, southeast Queensland, Australia

A catchment-scale multivariate statistical analysis of hydrochemistry enabled assessment of interactions between alluvial groundwater and Cressbrook Creek, an intermittent drainage system in southeast Queensland, Australia. Hierarchical cluster analyses and principal component analysis were applied to time-series data to evaluate the hydrochemical evolution of groundwater during periods of extreme drought and severe flooding. A simple three-dimensional geological model was developed to conceptualise the catchment morphology and the stratigraphic framework of the alluvium. The alluvium forms a two-layer system with a basal coarse-grained layer overlain by a clay-rich low-permeability unit. In the upper and middle catchment, alluvial groundwater is chemically similar to streamwater, particularly near the creek (reflected by high HCO3/Cl and K/Na ratios and low salinities), indicating a high degree of connectivity. In the lower catchment, groundwater is more saline with lower HCO3/Cl and K/Na ratios, notably during dry periods. Groundwater salinity substantially decreased following severe flooding in 2011, notably in the lower catchment, confirming that flooding is an important mechanism for both recharge and maintaining groundwater quality. The integrated approach used in this study enabled effective interpretation of hydrological processes and can be applied to a variety of hydrological settings to synthesise and evaluate large hydrochemical datasets.

Nitrous oxide emissions from fertilized, irrigated cotton (Gossypium hirsutum L.) in the Aral Sea Basin, Uzbekistan : influence of nitrogen applications and irrigation practices

Nitrous oxide emissions were monitored at three sites over a 2-year period in irrigated cotton fields in Khorezm, Uzbekistan, a region located in the arid deserts of the Aral Sea Basin. The fields were managed using different fertilizer management strategies and irrigation water regimes. N2O emissions varied widely between years, within 1 year throughout the vegetation season, and between the sites. The amount of irrigation water applied, the amount and type of N fertilizer used, and topsoil temperature had the greatest effect on these emissions. Very high N2O emissions of up to 3000 μg N2O-N m−2 h−1 were measured in periods following N-fertilizer application in combination with irrigation events. These “emission pulses” accounted for 80–95% of the total N2O emissions between April and September and varied from 0.9 to 6.5 kg N2O-N ha−1.. Emission factors (EF), uncorrected for background emission, ranged from 0.4% to 2.6% of total N applied, corresponding to an average EF of 1.48% of applied N fertilizer lost as N2O-N. This is in line with the default global average value of 1.25% of applied N used in calculations of N2O emissions by the Intergovernmental Panel on Climate Change. During the emission pulses, which were triggered by high soil moisture and high availability of mineral N, a clear diurnal pattern of N2O emissions was observed, driven by daily changes in topsoil temperature. For these periods, air sampling from 8:00 to 10:00 and from 18:00 to 20:00 was found to best represent the mean daily N2O flux rates. The wet topsoil conditions caused by irrigation favored the production of N2O from NO3− fertilizers, but not from NH4+ fertilizers, thus indicating that denitrification was the main process causing N2O emissions. It is therefore argued that there is scope for reducing N2O emission from irrigated cotton production; i.e. through the exclusive use of NH4+ fertilizers. Advanced application and irrigation techniques such as subsurface fertilizer application, drip irrigation and fertigation may also minimize N2O emission from this regionally dominant agro-ecosystem.

A case study in flood fatality : Beijing July 2012 flood

Background The aim of this study is to examine the flood fatality with a view to identifying risks which may inform public policy responses to future flood. On July 21st, 2012, Beijing suffered the heaviest rain since 1963. The average rainfall was 215 mm over a 24 hour period in the central city (301mm in Fangshan District). The rain resulted in a flood that caused severe health, social and financial impact. Results This flood caused 79 deaths. Of the 71 deaths for which a specific cause could be identified, 5 were rescue team members, 42 were killed by drowning (11 in the car), and the others by electricity shock, fallen house, falling items and lightning. The total financial cost was estimated to be US$ 1.7 billion. The causations of the deaths inform the risks associated with the flood. Discussion This flood had a catastrophic impact on Beijing, mainly due to the intensity of the rain (the rain was the heaviest in the modern Beijing history; possibly due to global warming and urban heat island effect), the vulnerability of the infrastructure (poor standards of drainage, disorganized water management systems and decreased permeability of the earth as a result of the city’s rapid development), and the capacity of the response system (mainly dependent on the awareness of the citizens, warning systems and the capacity of the emergency rescue). Implication Many risk management measures have been implemented as a result of this flood, including water level warning marks, flood safety education and warnings sent to mobile phones, a project to move about 74,500 farmers away from the flood-prone areas within 5 years. However, further measures targeted at the fundamental issues identified by this analysis are necessary, especially those targeting at health issues. These may include better planning, improved drainage systems and ecological development to increase permeability etc..

Spatial and temporal variability in a residual soil profile

It is well understood that that there is variation inherent in all testing techniques, and that all soil and rock materials also contain some degree of natural variability. Less consideration is normally given to variation associated with natural material heterogeneity within a site, or the relative condition of the material at the time of testing. This paper assesses the impact of spatial and temporal variability upon repeated insitu testing of a residual soil and rock profile present within a single residential site over a full calendar year, and thus range of seasonal conditions. From this repeated testing, the magnitude of spatial and temporal variation due to seasonal conditions has demonstrated that, depending on the selected location and moisture content of the subsurface at the time of testing, up to a 35% variation within the test results can be expected. The results have also demonstrated that the completed insitu test technique has a similarly large measurement and inherent variability error and, for the investigated site, up to a 60% variation in normalised results was observed. From these results, it is recommended that the frequency and timing of insitu tests should be considered when deriving geotechnical design parameters from a limited data set.

Effects of magnetite on high-frequency ground-penetrating radar

Large concentrations of magnetite in sedimentary deposits and soils with igneous parent material have been reported to affect geophysical sensor performance. We have undertaken the first systematic experimental effort to understand the effects of magnetite for ground-penetrating radar (GPR) characterization of the shallow subsurface. Laboratory experiments were conducted to study how homogeneous magnetite-sand mixtures and magnetite concentrated in layers affect the propagation behavior (velocity, attenuation) of high-frequency GPR waves and the reflection characteristics of a buried target. Important observations were that magnetite had a strong effect on signal velocity and reflection, at magnitudes comparable to what has been observed in small-scale laboratory experiments that measured electromagnetic properties of magnetite-silica mixtures. Magnetite also altered signal attenuation and affected the reflection characteristics of buried targets. Our results indicated important implications for several fields, including land mine detection, Martian exploration, engineering, and moisture mapping using satellite remote sensing and radiometers.

Landform characterization using geophysics—Recent advances, applications, and emerging tools

This paper presents an overview of the strengths and limitations of existing and emerging geophysical tools for landform studies. The objectives are to discuss recent technical developments and to provide a review of relevant recent literature, with a focus on propagating field methods with terrestrial applications. For various methods in this category, including ground-penetrating radar (GPR), electrical resistivity (ER), seismics, and electromagnetic (EM) induction, the technical backgrounds are introduced, followed by section on novel developments relevant to landform characterization. For several decades, GPR has been popular for characterization of the shallow subsurface and in particular sedimentary systems. Novel developments in GPR include the use of multi-offset systems to improve signal-to-noise ratios and data collection efficiency, amongst others, and the increased use of 3D data. Multi-electrode ER systems have become popular in recent years as they allow for relatively fast and detailed mapping. Novel developments include time-lapse monitoring of dynamic processes as well as the use of capacitively-coupled systems for fast, non-invasive surveys. EM induction methods are especially popular for fast mapping of spatial variation, but can also be used to obtain information on the vertical variation in subsurface electrical conductivity. In recent years several examples of the use of plane wave EM for characterization of landforms have been published. Seismic methods for landform characterization include seismic reflection and refraction techniques and the use of surface waves. A recent development is the use of passive sensing approaches. The use of multiple geophysical methods, which can benefit from the sensitivity to different subsurface parameters, is becoming more common. Strategies for coupled and joint inversion of complementary datasets will, once more widely available, benefit the geophysical study of landforms.Three cases studies are presented on the use of electrical and GPR methods for characterization of landforms in the range of meters to 100. s of meters in dimension. In a study of polygonal patterned ground in the Saginaw Lowlands, Michigan, USA, electrical resistivity tomography was used to characterize differences in subsurface texture and water content associated with polygon-swale topography. Also, a sand-filled thermokarst feature was identified using electrical resistivity data. The second example is on the use of constant spread traversing (CST) for characterization of large-scale glaciotectonic deformation in the Ludington Ridge, Michigan. Multiple CST surveys parallel to an ~. 60. m high cliff, where broad (~. 100. m) synclines and narrow clay-rich anticlines are visible, illustrated that at least one of the narrow structures extended inland. A third case study discusses internal structures of an eolian dune on a coastal spit in New Zealand. Both 35 and 200. MHz GPR data, which clearly identified a paleosol and internal sedimentary structures of the dune, were used to improve understanding of the development of the dune, which may shed light on paleo-wind directions.

Hydrological consequences of land-cover change: Quantifying the influence of plants on soil moisture with time-lapse electrical resistivity

Electrical resistivity of soils and sediments is strongly influenced by the presence of interstitial water. Taking advantage of this dependency, electrical-resistivity imaging (ERI) can be effectively utilized to estimate subsurface soil-moisture distributions. The ability to obtain spatially extensive data combined with time-lapse measurements provides further opportunities to understand links between land use and climate processes. In natural settings, spatial and temporal changes in temperature and porewater salinity influence the relationship between soil moisture and electrical resistivity. Apart from environmental factors, technical, theoretical, and methodological ambiguities may also interfere with accurate estimation of soil moisture from ERI data. We have examined several of these complicating factors using data from a two-year study at a forest-grassland ecotone, a boundary between neighboring but different plant communities.At this site, temperature variability accounts for approximately 20-45 of resistivity changes from cold winter to warm summer months. Temporal changes in groundwater conductivity (mean=650 S/cm =57.7) and a roughly 100-S/cm spatial difference between the forest and grassland had only a minor influence on the moisture estimates. Significant seasonal fluctuations in temperature and precipitation had negligible influence on the basic measurement errors in data sets. Extracting accurate temporal changes from ERI can be hindered by nonuniqueness of the inversion process and uncertainties related to time-lapse inversion schemes. The accuracy of soil moisture obtained from ERI depends on all of these factors, in addition to empirical parameters that define the petrophysical soil-moisture/resistivity relationship. Many of the complicating factors and modifying variables to accurately quantify soil moisture changes with ERI can be accounted for using field and theoretical principles.

Characterization and Mapping of Patterned Ground in the Saginaw Lowlands, Michigan: Possible Evidence for Late-Wisconsin Permafrost

We identified, mapped, and characterized a widespread area (gt;1,020 km2) of patterned ground in the Saginaw Lowlands of Michigan, a wet, flat plain composed of waterlain tills, lacustrine deposits, or both. The polygonal patterned ground is interpreted as a possible relict permafrost feature, formed in the Late Wisconsin when this area was proximal to the Laurentide ice sheet. Cold-air drainage off the ice sheet might have pooled in the Saginaw Lowlands, which sloped toward the ice margin, possibly creating widespread but short-lived permafrost on this glacial lake plain. The majority of the polygons occur between the Glacial Lake Warren strandline (~14.8 cal. ka) and the shoreline of Glacial Lake Elkton (~14.3 cal. ka), providing a relative age bracket for the patterned ground. Most of the polygons formed in dense, wet, silt loam soils on flat-lying sites and take the form of reticulate nets with polygon long axes of 150 to 160 m and short axes of 60 to 90 m. Interpolygon swales, often shown as dark curvilinears on aerial photographs, are typically slightly lower than are the polygon centers they bound. Some portions of these interpolygon swales are infilled with gravel-free, sandy loam sediments. The subtle morphology and sedimentological characteristics of the patterned ground in the Saginaw Lowlands suggest that thermokarst erosion, rather than ice-wedge replacement, was the dominant geomorphic process associated with the degradation of the Late-Wisconsin permafrost in the study area and, therefore, was primarily responsible for the soil patterns seen there today.

What’s best for our lymphoedema patients? Have we lost the patient as an individual in the quest for good science?

In a recent issue of the Journal of Lymphoedema, Nickolaidis and Karlsson (2013) indicated that most of the standard treatments for lymphoedema patients were explored and developed early last century, and suggested that holistic assessment of the individual is critical for good outcomes, but that perhaps “less emphasis should be placed on manual lymphatic drainage (MLD) and more on compression, exercise and weight reduction.”

The vertebral venous system in healthy and scoliotic adolescent spines - a 3D MRI investigation

Introduction. The venous drainage system within vertebral bodies (VBs) has been well documented previously in cadaveric specimens. Advances in 3D imaging and image processing now allow for in vivo quantification of larger venous vessels, such as the basivertebral vein. Differences between healthy and scoliotic VB veins can therefore be investigated. Methods. 20 healthy adolescent controls and 21 AIS patients were recruited (with ethics approval) to undergo 3D MRI, using a 3 Tesla, T1-weighted 3D gradient echo sequence, resulting in 512 slices across the thoraco-lumbar spine, with a voxel size of 0.5x0.5x0.5mm. Using Amira Filament Editor, five transverse slices through the VB were examined simultaneously and the resulting observable vascular network traced. Each VB was assessed, and a vascular network recorded when observable. A local coordinate system was created in the centre of each VB and the vascular networks aligned to this. The length of the vascular network on the left and right sides (with a small central region) of the VB was calculated, and the spatial patterning of the networks assessed level-by-level within each subject. Results. An average of 6 (range 4-10) vascular networks, consistent with descriptions of the basivertebral vein, were identifiable within each subject, most commonly between T10-L1. Differences were seen in the left/right distribution of vessels in the control and AIS subjects. Healthy controls saw a percentage distribution of 29:18:53 across the left:centre:right regions respectively, whereas the AIS subjects had a slightly shifted distribution of 33:25:42. The control group showed consistent spatial patterning of the vascular networks across most levels, but this was not seen in the AIS group. Conclusion. Observation and quantification of the basivertebral vein in vivo is possible using 3D MRI. The AIS group lacked the spatial pattern repetition seen in the control group and minor differences were seen in the left/right distribution of vessels.

Variations in vertebral venous vasculature

INTRODUCTION Cadaveric studies have previously documented typical patterns of venous drainage within vertebral bodies (VBs) [1,2,3], comprised primarily of the basivertebral vein, a planar tree like structure at the mid-height of the VB. These studies, however, are limited in the number of samples available, and so have not examined any potential differences in this anatomy in conditions such as scoliosis. MRI is able to create 3D images of soft tissue structures in the spine, including the basivertebral vein without the use of contrast. As a non-invasive imaging technique this opens up the possibility of examining the venous network in multiple VBs within the same subject, in healthy controls as well as in subjects with abnormal anatomy such as adolescent idiopathic scoliosis (AIS). CONCLUSIONS High resolution MRI scans allow in vivo quantification of the vertebral venous system at multiple levels on healthy and scoliotic populations for the first time. The length of the basivertebral vein was seen to have a significant bias to the right hand side of the VB in both healthy and AIS adolescents. The spatial pattern of this vein showed large variations in branching both within and across individuals.

Temporal trend of pesticide concentrations in the Chikugo River (Japan) with changes in environmental regulation and field infrastructure

Contamination of pesticides, which are applied to rice paddy fields, in river water has been a major problem in Japan for decades. A prolonged water holding period after pesticide application in paddy fields is expected to reduce the concentration of rice pesticides in river water. Therefore, a long monitoring campaign was conducted from 2004 to 2010 to measure the concentrations of pesticides in water samples collected from several points along the Chikugo River (Japan) including tributaries and the main stream to see if there was any reduction in the level of pesticide contamination after the extension of the water holding period (from 3–4 days to 7 days) was introduced in 2007 by the new water management regulation. No significant difference (p > 0.05) was found in pesticide concentrations between the periods before and after 2007 in all monitoring points, except in one tributary where the pesticide concentrations after 2007 were even higher than that of the previous period. A detailed study in one of the tributaries also revealed that the renovated infrastructure did not reduce the pesticide concentrations in the drainage canals. Neither the introduction of the new regulation nor the improved infrastructure had any significant effect on reducing the contamination of pesticides in water of the Chikugo River. It is probably because most farmers did not properly implement the new requirement of holding paddy water within the field for 7 days after the application of pesticides. Only tightening the regulation would not be sufficient and more actions should be taken to enforce/provide extension support for the new water management regulation in order to reduce the level of residual pesticides in river water in Japan.

Effect of water management practice on pesticide behavior in paddy water

The fate and transport of three herbicides commonly used in rice production in Japan were compared using two water management practices. The herbicides were simetryn, thiobencarb and mefenacet. The first management practice was an intermittent irrigation scheme using an automatic irrigation system (AI) with a high drainage gate and the second one was a continuous irrigation and overflow drainage scheme (CI) in experimental paddy fields. Dissipation of the herbicides appeared to follow first order kinetics with the half-lives (DT₅₀) of 1.6-3.4 days and the DT₉₀ (90% dissipation) of 7.4-9.8 days. The AI scheme had little drainage even during large rainfall events thus resulting in losses of less than 4% of each applied herbicide through runoff. Meanwhile the CI scheme resulted in losses of about 37%, 12% and 35% of the applied masses of simetryn, thiobencarb and mefenacet, respectively. The intermittent irrigation scheme using an automatic irrigation system with a high drainage gate saved irrigation water and prevented herbicide runoff whereas the continuous irrigation and overflow scheme resulted in significant losses of water as well as the herbicides. Maintaining the excess water storage is important for preventing paddy water runoff during significant rainfall events. The organic carbon partition coefficient K_oc seems to be a strong indicator of the aquatic fate of the herbicide as compared to the water solubility (S_W). However, further investigations are required to understand the relation between K_oc and the agricultural practices upon the pesticide fate and transport. An extension of the water holding period up to 10 days after herbicide application based on the DT₉₀ from the currently specified period of 3-4 days in Japan is recommended to be a good agricultural practice for controlling the herbicide runoff from paddy fields. Also, the best water management practice, which can be recommended for use during the water holding period, is the intermittent irrigation scheme using an automatic irrigation system with a high drainage gate. © 2006 Elsevier B.V. All rights reserved.

Applicability of ELISA in pesticide monitoring to control runoff of bensulfuron-methyl and simetryn from paddy fields

The applicability of ELISA kits was evaluated as an alternative to monitor bensulfuron-methyl and simetryn behavior in paddy water under intermittent (Plot 1) and continuous (Plot 2) irrigation schemes. Simetryn concentrations in both plots decreased exponentially from the peak of the first day. However, the simetryn kit systematically underestimated by a factor of 0.79 as compared to the GC method. Bensulfuron-methyl concentrations exhibited similar dissipation kinetics in paddy water and the drainage water. The bensulfuron-methyl kit was capable of distinguishing spatial variations of concentrations in the paddy field. The ELISA kits clearly indicated differences in the loss of both herbicides between the two plots and therefore may be useful for evaluating the water management practice of pesticide runoff control in paddy fields.

Simulation of pesticide behavior in a paddy block by a pesticide fate and transport model

A simulation model (PCPF-B) was developed based on the PCPF-1 model to predict the runoff of pesticides from paddy plots to a drainage canal in a paddy block. The block-scale model now comprises three modules: (1) a module for pesticide application, (2) a module for pesticide behavior in paddy fields, and (3) a module for pesticide concentration in the drainage canal. The PCPF-B model was first evaluated by published data in a single plot and then was applied to predict the concentration of bensulfuron-methyl in one paddy block in the Sakura river basin, Ibaraki, Japan, where a detailed field survey was conducted. The PCPF-B model simulated well the behavior of bensulfuron-methyl in individual paddy plots. It also reflected the runoff pattern of bensulfuron-methyl at the block outlet, although overestimation of bensulfuronmethyl concentrations occurred due to uncertainty in water balance estimation. Application of water management practice such as water-holding period and seepage control also affected the performance of the model. A probabilistic approach may be necessary for a comprehensive risk assessment in large-scale paddy areas.