Advancing process‐based watershed hydrological research using near‐surface geophysics: a vision for,and review of,electrical and magnetic geophysical methods

We want to develop a dialogue between geophysicists and hydrologists interested in synergistically advancing process based watershed research. We identify recent advances in geophysical instrumentation, and provide a vision for the use of electrical and magnetic geophysical instrumentation in watershed scale hydrology. The focus of the paper is to identify instrumentation that could significantly advance this vision for geophysics and hydrology during the next 3–5 years. We acknowledge that this is one of a number of possible ways forward and seek only to offer a relatively narrow and achievable vision. The vision focuses on the measurement of geological structure and identification of flow paths using electrical and magnetic methods. The paper identifies instruments, provides examples of their use, and describes how synergy between measurement and modelling could be achieved. Of specific interest are the airborne systems that can cover large areas and are appropriate for watershed studies. Although airborne geophysics has been around for some time, only in the last few years have systems designed exclusively for hydrological applications begun to emerge. These systems, such as airborne electromagnetic (EM) and transient electromagnetic (TEM), could revolutionize hydrogeological interpretations. Our vision centers on developing nested and cross scale electrical and magnetic measurements that can be used to construct a three‐dimensional (3D) electrical or magnetic model of the subsurface in watersheds. The methodological framework assumes a ‘top down’ approach using airborne methods to identify the large scale, dominant architecture of the subsurface. We recognize that the integration of geophysical measurement methods, and data, into watershed process characterization and modelling can only be achieved through dialogue. Especially, through the development of partnerships between geophysicists and hydrologists, partnerships that explore how the application of geophysics can answer critical hydrological science questions, and conversely provide an understanding of the limitations of geophysical measurements and interpretation. Copyright © 2008 John Wiley & Sons, Ltd.     

Observations and model estimates of diurnal water temperature dynamics in mosquito breeding sites in western Kenya

Water temperature is an important determinant of the growth and development of malaria mosquito immatures. To gain a better understanding of the daily temperature dynamics of malaria mosquito breeding sites and of the relationships between meteorological variables and water temperature, three clear water pools (diameter × depth: 0·16 × 0·04, 0·32 × 0·16 and 0·96 × 0·32 m) were created in Kenya. Continuous water temperature measurements at various depths were combined with weather data collections from a meteorological station. The water pools were homothermic, but the top water layer differed by up to about 2 °C in temperature, depending on weather conditions. Although the daily mean temperature of all water pools was similar (27·4–28·1 °C), the average recorded difference between the daily minimum and maximum temperature was 14·4 °C in the smallest versus 7·1 °C in the largest water pool. Average water temperature corresponded well with various meteorological variables. The temperature of each water pool was continuously higher than the air temperature. A model was developed that predicts the diurnal water temperature dynamics accurately, based on the estimated energy budget components of these water pools. The air–water interface appeared the most important boundary for energy exchange processes and on average 82–89% of the total energy was gained and lost at this boundary. Besides energy loss to longwave radiation, loss due to evaporation was high; the average estimated daily evaporation ranged from 4·2 mm in the smallest to 3·7 mm in the largest water pool. Copyright © 2008 John Wiley & Sons, Ltd.     

The use of inclined hemisphere projections for analyzing failure mechanisms in discontinuous rocks

This paper demonstrates the advantages of using inclined stereographic projections in kinematic analysis of rock blocks in discontinuous rock masses. Some examples of limiting cases are presented. The application of inclined projections is illustrated by its use in a mine slope in Brazil. It is clear from the discussion of these examples that inclined hemisphere projections provide better results than horizontal projections. It is also demonstrated that horizontal projections can lead to incorrect results in limiting cases.     

Stable isotopic compositions of Recent freshwater cyanobacterial carbonates from the British Isles: local and regional environmental controls

Recent (<50 years old) freshwater cyanobacterial carbonates from diverse environments (streams, lakes, waterfalls) throughout Britain and Ireland were analysed for their stable carbon and oxygen isotope compositions. The mean δ¹⁸O value of ?5–9‰ PDB for river and stream data represents calcite precipitation in equilibrium with the mean oxygen isotopic composition of precipitation in central Britain (?7–5‰SMOW) assuming a mean water temperature of 9°C. The mean δ¹⁸O of lake data, ?4–5‰ PDB, is statistically different, reflecting the effects of residence time and/or variations in the oxygen isotopic composition of rainfall. Carbon isotopes have wide variations in both fluviatile and lake data sets (+ 3 to ?12‰ PDB). These variations are principally controlled in the fluviatile samples by contribution of isotopically light ‘soil zone’ carbon relative to isotopically heavier carbon from limestone aquifer rock dissolution. Lake samples have the heaviest carbon isotope values, reflecting a trend toward isotopic equilibrium between atmospheric CO₂ and aqueous HCO^?₃. We infer that isotopic compositions of ancient cyanobacterial carbonates should also record environmental information, although the effects of stabilization and diagenesis on primary δ¹⁸O values will need careful consideration. Primary carbon isotope compositions should be well preserved, although in marine samples values will be buffered by the isotopic composition of aqueous marine bicarbonate.     

Temperature evolution between 50 K and 320 K of the thermal expansion tensor of gypsum derived from neutron powder diffraction data

The unit cell parameters, extracted from Rietveld analysis of neutron powder diffraction data collected between 4.2 K and 320 K, have been used to calculate the temperature evolution of the thermal expansion tensor for gypsum for 50 ≤ T ≤ 320 K. At 300 K the magnitudes of the principal axes are α ₁₁  = 1.2(6) × 10⁻⁶ K⁻¹, α ₂₂  = 36.82(1) × 10⁻⁶ K⁻¹ and α ₃₃  = 25.1(5) × 10⁻⁶ K⁻¹. The maximum axis, α ₂₂ , is parallel to b, and using Institution of Radio Engineers (IRE) convention for the tensor orthonormal basis, the axes α ₁₁ and α ₃₃ have directions equal to (−0.979, 0, 0.201) and (0.201, 0, 0.979) respectively. The orientation and temperature dependent behaviour of the thermal expansion tensor is related to the crystal structure in the I2/a setting. Received 12 February 1998 / Revised, accepted 19 October 1998     

Soil moisture potential and water content in the unsaturated zone within the arid Ejina Oasis in Northwest China

Three soil profiles were selected in the Ejina Oasis, northwest China, to determine water content profiles and evolution of soil moisture potentials in the unsaturated zone within the arid area. The total soil moisture potentials have been monitored for about 3 months in 2001 at different depths in the soil profiles. The occurrence and movement of water in the unsaturated zone was analyzed using the zero flux plane (ZFP) method. It is shown that convergent ZFPs and divergent ZFPs may occur at depths between 0.5 and 3.0 m, and that the depth of the ZFPs was controlled by the root zone of plants growing on the land surface. Profiles of the total soil moisture potentials were observed to be coincident with those of the water contents at the three experimental sites. The total soil moisture potential showed a slight increasing trend and the ZFPs tend to vanish from summer to winter as the water extraction by roots decreased. Evapotranspiration through vegetation has an important bearing on the water content and the total potential in the unsaturated zone.     

Chemical assembly of Archean anorthosites from amphibolite- and granulite-facies terranes, SW Greenland

We report whole-rock, major- and trace-element compositions (obtained by XRF and INA methods) for the amphibolite-facies Buksefjorden and granulite-facies Nordland anorthosites, SW Greenland. In a previous petrologic study on the same sample suite, we documented differences in texture, mineralogy, and mineral compositions between these two anorthosite bodies. Chemical analyses confirm differences in composition between the two bodies, but these differences cannot be explained by variations in metamorphic conditions, and point towards differences in the nature of their protoliths. Analyzed Nordland samples are anorthosites and leucogabbros with 88-98% normative plagioclase, whereas those from Buksefjorden include anorthosite, leucogabbro, and gabbro with ~55-95% normative plagioclase. Two or more compositional groupings can be recognized at each site, which correspond to differences in color and mineralogy of the hand samples. Samples from Buksefjorden are mainly quartz-normative, whereas those from Nordland are olivine (- nepheline) normative. Other differences include higher Ni/Co ratio and REE contents in the granulite-facies anorthosites from Nordland. REE pattern shapes are similar, however, being moderately fractionated at ~0.5-102 chondrites with positive Eu-anomalies. Calculated equilibrium melt patterns are similar for both anorthosites, being relatively flat at ~50-1502 chondrites, suggesting unfractionated (but evolved) parental magmas. Olivine must have been present in the protoliths of the Nordland rocks compared with Buksefjorden. Otherwise, the protoliths contained plagioclase with variable An-content (~An₆₂-An₉₂) and a mafic component with variable Fe/Mg (mg ~0.3-0.8). This mafic component was either hornblende or a combination of ortho- and clinopyroxene in fixed proportions, plus a small amount of magnetite. Mixing calculations demonstrate that some Buksefjorden anorthosites contain two varieties of plagioclase: a calcic type that may correspond to cumulus crystals, and a sodic-type that may correspond to a trapped-melt component. On plots of normative whole-rock An versus mg, compositions of the Buksefjorden and Nordland anorthosites form crude negative arrays that differ from the generally positive trends of mafic layered intrusions (Kiglapait, Skaergaard) and from the generally flats trends of plagioclase-rich cumulate rocks (St. Urbain and Stillwater anorthosites). This difference provides further evidence for the distinctive nature of Archean calcic-anorthosite complexes compared with other types of mafic intrusions. Moreover, this distribution of data points is consistent with the assembly of the protolith of the SW Greenland anorthosites mainly as mixtures of plagioclase and hornblende. Finally, the field for the Buksefjorden and Nordland anorthosites overlaps only slightly with that for the Fiskenaesset Complex, thus extending the known range of compositions for Archean anorthosites in West Greenland.     

Delineating depth to bedrock beneath shallow unconfined aquifers: a gravity transect across the Palmer River Basin

In view of the increasing demand on ground water supplies in the northeastern United States, it is imperative to develop appropriate methods to geophysically characterize the most widely used sources of ground water in the region: shallow unconfined aquifers consisting of well-sorted, stratified glacial deposits laid down in bedrock valleys and channels. The gravity method, despite its proven value in delineating buried bedrock valleys elsewhere, is seldom used by geophysical contractors in this region. To demonstrate the method's effectiveness for evaluating such aquifers, a pilot study was undertaken in the Palmer River Basin in southeastern Massachusetts. Because bedrock is so shallow beneath this aquifer (maximum depth is 30 m), the depth-integrated mass deficiency of the overlying unconsolidated material was small, so that the observed gravity anomaly was on the order of 1 milligal (mGal) or less. Thus data uncertainties were significant. Moreover, unlike previous gravity studies elsewhere, we had no a priori information on the density of the sediment. Under such circumstances, it is essential to include model constraints and weighted least-squares in the inversion procedure. Among the model constraints were water table configuration, bedrock outcrops, and depth to bedrock from five water wells. Our procedure allowed us to delineate depth to bedrock along a 3.5 km profile with a confidence interval of 1.8 m at a nominal depth of 17 m. Moreover, we obtained a porosity estimate in the range of 39% to 44%. Thus the gravity method, with appropriate refinements, is an effective tool for the reconnaissance of shallow unconfined aquifers.     

Preliminary report on radionuclide concentrations in the bottom waters at the entrance of Wakasa Bay with special reference to the Japan Sea Proper Water

Marked fluctuation of concentrations of⁹⁰Sr and¹³⁷Cs was observed in the bottom waters at the entrance of Wakasa Bay during 1987–1992, and the cause was investigated. The concentrations of⁹⁰Sr and¹³⁷Cs in the bottom waters were significantly low when the upper level of the Japan Sea Proper Water (JSPW) was high and covered the sampling depth, but high when the upper level of the JSPW was low. The cause of the fluctuation observed in the bottom waters is, therefore, suggested to be the vertical fluctuation of the upper level of the JSPW on the shelf slope, which has been little described before.     

Simulation of Tsaoling landslide, Taiwan, based on Saint Venant equations over general topography

The 1999 Chi-Chi earthquake triggered the catastrophic Tsaoling landslide in Taiwan. The geomorphological change measured from the data of the 1989 and 2000 aerial photos reveals that the scar and deposit volumes are about 0.126 km³ and 0.15 km³ respectively. The debris material ran over a distance of 1.6 km with 500 m descent in elevation. In this paper, we use the continuum model of hydraulic flow, SHALTOP^2D, based on the equations of Bouchut and Westdickenberg to simulate numerically the landslide dynamics. When the mass is moving, the flow is assumed hydrostatic with a basal Coulomb friction. The best fit is obtained using the basal friction angle equal to 6°, the only parameter of the simulation, uniformly applied in the calculation domain. The landslide front reaches the Chinshui river valley, the foothill of the slope, within only about 25 s after initiation and the motion settles in about 113 s. The maximum speed is estimated 75 m/s. The spreading of the deposit agrees well with the field measurement.