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.     

Jurassic‐to‐present thermal history of the central High Atlas (Morocco) assessed by low‐temperature thermochronology

Apatite fission track (AFT) and (U–Th)/He data from the High Atlas have been obtained for the first time to constrain the tectono‐thermal evolution of the central part of the chain. Results from Palaeozoic basement massifs indicate long residence at low temperatures, consistently with their original location out of the deepest Mesozoic rift troughs and indicating minor exhumation. The best rocks for extracting the Alpine history of the Atlas Mountains are Jurassic intrusives, which yield AFT ages centred on c. 80 Ma; thermal models based on AFT data and constrained by (U–Th)/He suggest that these ages are included in a slow cooling trend from intrusion age to c. 50 Ma ago that we attribute to post‐rift thermal relaxation. This is followed by a stability period of c. 30 Ma and then by a final exhumational cooling until present exposure. Eocene intrusives yield AFT ages similar to those of Rb–Sr and K–Ar suggesting rapid emplacement in the uppermost crust.     

Uniform Permian 40Ar/39Ar detrital mica ages in the eastern Qaidam Basin (NW China): where is the source?

New single‐grain ⁴⁰Ar/³⁹Ar detrital white‐mica ages from the Lulehe section at the eastern Qaidam Basin yield uniform Permian ages between 250 ± 3 and 279 ± 3 Ma throughout the whole Cenozoic sequence. This is inconsistent with the present hinterland, which is composed of early Palaeozoic metamorphic units with subordinate early Palaeozoic and few Permian granites. The new data indicate that Permian tectonic units are likely more widespread at the north‐eastern margin of the Tibetan plateau as known at present, particularly within the Qilian Mountains. The preferred explanation is that the Qaidam block represents a rigid indenter, which indented during late Tertiary times into early Palaeozoic orogenic units. This is consistent with recent findings of a NW‐trending sinistral Permian ductile shear zone and a dextral, NW‐trending Tertiary fault system close to the north‐eastern margin of the Qaidam Basin.     

Scattering of elastic wave by a partially sealed cylindrical shell embedded in poroelastic medium

Scattering of an elastic wave by a cylindrical shell embedded in poroelastic medium is investigated theoretically with the assumption that the shell material is also a porous elastic medium. The porous medium is modellized via Biot's theory and the scattering by cylindrical shell is expressed by the definition of scattering matrix. The normal mode expansion technique is employed for analysing the scattering field, and the asymptotic solutions of displacements, stresses and pore pressure are derived. Two limiting cases‐scattering by a poroelastic cylinder in Biot medium and a elastic cylindrical shell in elastic medium are obtained from the general solutions. The dispersion curves of displacement amplitude at the interface of shell and medium is compared with the case of elastic shell. The scattering amplitude associated with the fast, slow and transverse waves are identified by numerical simulation. Furthermore, the influence of the poroelastic property of shell material on scattering amplitude is analysed. Copyright © 2005 John Wiley & Sons, Ltd.     

Analysis of geomorphologic and hydrological characteristics in watershed saturated areas using topographic‐index threshold and geomorphology‐based runoff model

The objective of this paper is to investigate the variation of geomorphology and runoff characteristics in saturated areas under different partial contributing area (PCA) conditions. Geomorphologic information and hydrologic records from two mid‐size watersheds in northern Taiwan were selected for analysis. The PCA ratio in the watershed during a storm was assumed equal to the ratio of the surface‐flow volume to the direct runoff volume from measured hydrologic data. The extents of PCA regions were then determined by using a topographic‐index threshold. Consequently, the geomorphologic factors in saturated and unsaturated areas could be calculated using a digital elevation model, and these factors could then be linked to a geomorphology‐based IUH model for runoff simulation, which can consider both the surface‐ and subsurface‐flow processes in saturated and unsaturated areas, respectively. The results show that geomorphologic characteristics in the saturated areas vary significantly with different PCA ratios especially for higher order streams. A large PCA ratio results in a sharp hydrograph because the quick surface flow dominates the runoff process, whereas the hydrologic response in a low PCA case is dominated by the delayed subsurface flow. Copyright © 2007 John Wiley & Sons, Ltd.     

Water temperature and conductivity variability as indicators of groundwater behaviour in complex aquifer systems in the south‐east of Spain

Groundwater temperature and electrical conductivity logs were applied in this study to gain an understanding of the hydrogeological behaviour of a complex system. This system, compartmentalized by a large number of fractures, which may even cause geothermal anomalies, is that of the Campo de Dalías in the south‐east of the Iberian Peninsula. In addition to this complexity one must consider the different levels with aquifer properties and their varying lithologies, which in some cases are in direct contact with the sea. Vertical sections based upon the various logs allow lateral correlation to be made between the anomalies detected, and provide information on the existence of both ascending and descending vertical flows, and the inflow of sea water into the system. Analysis of environmental isotopes (¹⁸O and D) indicate the areas of recharge with the greatest flow velocity. No seasonal variation in the isotopes was observed, which suggests that the residence time of the groundwater in the aquifer is several years. Radiocarbon dating gave estimates of generally long residence times, although in some sectors it is possible to detect the current recharge based on its tritium content, from rainwater infiltration or seawater intrusion. Copyright © 2002 John Wiley & Sons, Ltd.     

利用黄河泥沙围出“渤海湖”的初步设想

本文是一个关于将渤海变成淡水湖的水利工程的构想。让黄河的大量泥沙通过人工导流的方式沿莱州湾南岸与渤海海峡沉积延伸,产生的陆地将渤海封闭并逐渐形成淡水湖。这一工程将陆地连接山东半岛与辽东半岛,并为华北提供巨大淡水库,具有显著的经济、生态效益。具体实施上,可通过两道渔网式围栏促使黄河泥沙定点沉积并阻挡浪潮侵蚀,通过加强汛期水量增加黄河口泥沙。     

Biodegradation and transport of benzene,toluene, and xylenes in a simulated aquifer: comparison of modelled and experimental results

Both laboratory experiments and numerical modelling were conducted to study the biodegradation and transport of benzene–toluene–xylenes (BTX) in a simulated semi‐confined aquifer. The factors incorporated into the numerical model include advection, hydrodynamic dispersion, adsorption, and biodegradation. The various physico‐chemical parameters required by the numerical model were measured experimentally. In the experimental portion of the study, BTX compounds were introduced into the aquifer sand. After the contaminants had been transported through the system, BTX concentrations were measured at 12 equally spaced wells. Subsequently, microorganisms obtained from the activated sludge of a sewage treatment plant and cultured in BTX mixtures were introduced into the aquifer through the 12 sampling wells. The distribution data for BTX adsorption by the aquifer sand form a nonlinear isotherm. The degree of adsorption by the sand varies, depending on the composition of the solute. The degradation time, measured from the time since the bacteria were added to the aquifer until a specific contaminant was no longer detectable, was 35–42 h for BTX. The dissolved oxygen, after degradation by BTX compounds and bacteria, was consumed by about 40–60% in the entire simulated aquifer; thus the aerobic conditions were maintained. This study provides insights for the biodegradation and transport of BTX in aquifers by numerical modelling and laboratory experiments. Experimental and numerical comparisons indicate that the results by Monod degradation kinetics are more accurate than those by the first‐order degradation kinetics. Copyright © 2002 John Wiley & Sons, Ltd.     

Simultaneous identification of parameter,initial condition,and boundary condition in groundwater modelling

In this paper, we perform an inverse method to simultaneously estimate aquifer parameters, initial condition, and boundary conditions in groundwater modelling. The parameter estimation is extended to a complete inverse problem that makes the calibrated groundwater flow model more realistic. The adjoint state method, the gradient search method, and the least square error algorithm are combined to build the optimization procedure. Horizontal two‐dimensional groundwater flow in a confined aquifer is exemplified to demonstrate the correlation between unknowns, the contribution of observation, as well as the suitability of applying the inverse method. The correlation analysis shows the connection between storage coefficient and initial condition. Besides, transmissivity and boundary conditions are also highly correlated. More observations at different location and time are necessary to provide sufficient information. A time series of unsteady head is requested for estimation of storage coefficient and initial condition. Observation near boundary is very effective for boundary condition estimation. The observation at pumping well mostly contributes to the estimation of transmissivity. According to all observations, it is possible to identify parameters, initial condition, and boundary condition simultaneously. Furthermore, the results not only illustrate the traditional assumption of known boundary condition but also initial condition, which may cause an incorrect estimation. Copyright © 2009 John Wiley & Sons, Ltd.     

Towards a climate event stratigraphy for New Zealand over the past 30 000 years (NZ‐INTIMATE project)

It is widely recognised that the acquisition of high‐resolution palaeoclimate records from southern mid‐latitude sites is essential for establishing a coherent picture of inter‐hemispheric climate change and for better understanding of the role of Antarctic climate dynamics in the global climate system. New Zealand is considered to be a sensitive monitor of climate change because it is one of a few sizeable landmasses in the Southern Hemisphere westerly circulation zone, a critical transition zone between subtropical and Antarctic influences. New Zealand has mountainous axial ranges that amplify the climate signals and, consequently, the environmental gradients are highly sensitive to subtle changes in atmospheric and oceanic conditions. Since 1995, INTIMATE has, through a series of international workshops, sought ways to improve procedures for establishing the precise ages of climate events, and to correlate them with high precision, for the last 30 000 calendar years. The NZ‐INTIMATE project commenced in late 2003, and has involved virtually the entire New Zealand palaeoclimate community. Its aim is to develop an event stratigraphy for the New Zealand region over the past 30 000 years, and to reconcile these events against the established climatostratigraphy of the last glacial cycle which has largely been developed from Northern Hemisphere records (e.g. Last Glacial Maximum (LGM), Termination I, Younger Dryas). An initial outcome of NZ‐INTIMATE has been the identification of a series of well‐dated, high‐resolution onshore and offshore proxy records from a variety of latitudes and elevations on a common calendar timescale from 30 000 cal. yr BP to the present day. High‐resolution records for the last glacial coldest period (LGCP) (including the LGM sensu stricto) and last glacial–interglacial transition (LGIT) from Auckland maars, Kaipo and Otamangakau wetlands on eastern and central North Island, marine core MD97‐2121 east of southern North Island, speleothems on northwest South Island, Okarito wetland on southwestern South Island, are presented. Discontinuous (fragmentary) records comprising compilations of glacial sequences, fluvial sequences, loess accumulation, and aeolian quartz accumulation in an andesitic terrain are described. Comparisons with ice‐core records from Antarctica (EPICA Dome C) and Greenland (GISP2) are discussed. A major advantage immediately evident from these records apart from the speleothem record, is that they are linked precisely by one or more tephra layers. Based on these New Zealand terrestrial and marine records, a reasonably coherent, regionally applicable, sequence of climatically linked stratigraphic events over the past 30 000 cal. yr is emerging. Three major climate events are recognised: (1) LGCP beginning at ca. 28 000 cal. yr BP, ending at Termination I, ca. 18 000 cal. yr BP, and including a warmer and more variable phase between ca. 27 000 and 21 000 cal. yr BP, (2) LGIT between ca. 18 000 and 11 600 cal. yr BP, including a Lateglacial warm period from ca. 14 800 to 13 500 cal. yr BP and a Lateglacial climate reversal between ca. 13 500 and 11 600 cal. yr BP, and (3) Holocene interglacial conditions, with two phases of greatest warmth between ca. 11 600 and 10 800 cal. yr BP and from ca. 6 800 to 6 500 cal. yr BP. Some key boundaries coincide with volcanic tephras. Copyright © 2007 John Wiley & Sons, Ltd.