Analytical methods that use natural heat as a tracer to quantify surface water–groundwater exchange, evaluated using field temperature records

Two methods applying natural heat as a tracer to quantify surface water–groundwater exchange were evaluated using field data. Arrays capable of monitoring and recording the streambed response to diurnal temperature variations in the surface water were deployed for a 2-month period in three locations in perennial pools at Maules Creek, New South Wales, Australia. Multi-level array design, field deployment and parameter estimation are discussed. The applicability of analytical solutions derived from the heat transport equation to the streambed environments was analysed using the recorded temperature time series. The stream was found to lose water to the aquifer, which was supported by simultaneously recorded hydraulic gradients. However, the one-dimensional (1D) analytical solutions did not adequately describe the observed streambed thermal response at two locations. The resulting artefacts in the estimated flow velocities are discussed. It was hypothesised that the artefacts originate from model limitation due to streambed heterogeneity and application of 1D solutions to multi-dimensional and dynamic streambed flow. This consequently imposes limitations on the field applicability of the methods. Nevertheless, in combination with time series of surface water and streambed water levels, the use of heat as a tracer provided a powerful tool for better understanding the shallow hydrogeological system.     

Mean Flow Near Edges and Within Cavities Situated Inside Dense Canopies

A streamfunction-vorticity formulation is used to explore the extent to which turbulent and turbulently inviscid solutions to the mean momentum balance explain the mean flow across forest edges and within cavities situated inside dense forested canopies. The turbulent solution is based on the mean momentum balance where first-order closure principles are used to model turbulent stresses. The turbulently inviscid solution retains all the key terms in the mean momentum balance but for the turbulent stress gradients. Both exit and entry versions of the forest edge problem are explored. The turbulent solution is found to describe sufficiently the bulk spatial patterns of the mean flow near the edge including signatures of different length scales reported in canopy transition studies. Next, the ‘clearing inside canopy’ or the so-called ‘cavity’ problem is solved for the inviscid and turbulent solutions and then compared against flume experiments. The inviscid solution describes the bulk flow dynamics in much of the zones within the cavity. In particular, the solution can capture the correct position of the bulk recirculation zone within the cavity, although with a weaker magnitude. The inviscid solution cannot capture the large vertical heterogeneity in the mean velocity above the canopy, as expected. These features are better captured via the first-order closure representation of the turbulent solution. Given the ability of this vorticity formulation to capture the mean pressure variations and the mean advective acceleration terms, it is ideal for exploring the distributions of scalars and roughness-induced flow adjustments on complex topography.     

Estimation of the Variability of Tailings Dams Properties in Order to Perform Probabilistic Assessment

Dams are often used for storing mine tailings. These structures present failure mechanisms that can lead to major risks for surrounding populations. In view to managing these risks, it is necessary to take a probabilistic approach in predicting their behaviour. The application of such approaches is limited by the difficulty of obtaining experimental data to estimate the variability of the parameters and conditioned by the relevance of the probability models chosen to represent this variability. This article proposes (1) a probabilistic modelling of the index properties of the mine tailings constituting these dams based on statistical analyses and (2) a method using dynamic penetration tests to estimate on site the mine tailings friction angle and its variability. This method, applied to chilean tailings dams proposes a single model for all tailings dams in order to associate a probability law to the effective friction angle (ϕ′). The procedure is illustrated on the probabilistic study of slope stability carried out at the global scale of a dam and also at the local scale of each of its constituent layers.     

Spatial distribution and seasonal variation in 18O/16O of modern precipitation and river water across the conterminous USA

We report a quantitative analysis of regional differences in the the oxygen isotope composition of river water and precipitation across the USA because data are now available to undertake a more geographically and temporally extensive analysis than was formerly possible. Maps of modern, mean annual δ¹⁸O values for both precipitation (δ¹⁸O_PPT) and river water (δ¹⁸O_RIV) across the 48 contiguous states of the USA have been generated using latitude and elevation as the primary predictors of stable isotope composition while also incorporating regional and local deviations based on available isotopic data. The difference between these two maps was calculated to determine regions where δ¹⁸O_RIV is significantly offset from local δ¹⁸O_PPT. Additional maps depicting seasonal and extreme values for δ¹⁸O_RIV and δ¹⁸O_PPT were also constructed. This exercise confirms the presence of regions characterized by differences in δ¹⁸O_RIV and δ¹⁸O_PPT and specifically identifies the magnitude and regional extent of these offsets. In particular, the Great Plains has δ¹⁸O_RIV values that are more positive than precipitation, while much of the western USA is characterized by significantly lower δ¹⁸O_RIV values in comparison with local δ¹⁸O_PPT. The most salient feature that emerged from this comparison is the ‘catchment effect’ for the rivers. Because river water is largely derived from precipitation that fell upstream of the sample locality (i.e. at higher elevations) δ¹⁸O_RIV values are often lower than local δ¹⁸O_PPT values, particularly in catchments with high‐elevation gradients. Seasonal patterns in the isotopic data substantiate the generally accepted notion that amplitudes of δ¹⁸O variation are greatly dampened in river water relative to those of local precipitation. Copyright © 2005 John Wiley & Sons, Ltd.     

Quantifying net ecosystem exchange by multilevel ecophysiological and turbulent transport models

To quantify the interplay between scalar sources and sinks (S_c) and net ecosystem exchange (NEE), “forward” and “inverse” approaches have been proposed. The canonical form of forward approaches is a one-dimensional ecophysiological-radiative transfer scheme coupled to turbulent transport theory. In contrast, inverse approaches strictly rely on turbulent transport theory and mean scalar concentration as their primary input to infer S_c and NEE. While the formulation of both approaches have evolved over the past decade, no systematic comparison between them was undertaken for the same data set, and over a wide range of atmospheric conditions. Our objective is to compare the predicted S_c and NEE from these two approaches with eddy-covariance measurements. The results show that the forward method outperformed all three inverse methods for unstable and neutral conditions on short time scales (30 min) but yielded comparable results at longer time scales. Poor agreement was obtained under stable conditions for all models. Hence, for modeling event-based flux variations, forward models are preferred. Since the forward method requires detailed knowledge of ecophysiological, drag, radiative transfer and other canopy attributes, all of which are difficult to obtain on a routine basis, a symbiotic use of forward and inverse approaches is most advantageous.     

Validation of IPCC AR4 models over the Iberian Peninsula

This paper reports analysis of the ability of 24 coupled global climate models that were used in the Fourth Assessment Report (AR4) of the Intergovernmental Panel on Climate Change to simulate the current monthly seasonal cycle of sea level pressure, surface air temperature and precipitation over the Iberian Peninsula in the last two decades of the twentieth century. The period investigated runs from 1979 to 1998. In order to assess the performance of the models, averaged seasonal cycles and probability density functions (PDFs) calculated from model simulations are compared with the corresponding seasonal cycles, whilst PDFs are also obtained using the data from the ERA40 reanalysis and the Global Precipitation Climatology Project. We found that simulated PDFs generally provided a better fit to actual PDFs than seasonal cycles do. This conclusion indicates that when evaluating model performance, the climate variability as measured by means of PDFs is not the only climatic element that should be tested. Regarding the comparison based on the seasonal cycle, results also show that the root mean square skill score is more useful than the r skill score. To rank the AR4 models, sea level pressure, surface air temperature and precipitation variables were selected and a group of five AR4 models were identified as the models which best reproduce current climate in the area: MIROC3.2-HIRES, MPI-ECHAM5, GFDL-CM2.1, BCCR-BCM2.0 and UKMO-HADGEM1. The rank obtained should not be understood in a hierarchical manner because there is a certain degree of internal variability in the model ensembles. Finally, it should be noted that these results are in good agreement with other classifications found in the scientific literature.     

Double-difference Relocation of the Aftershocks of the Tecom��n, Colima, Mexico Earthquake of 22 January 2003

On 22 January 2003, the M _w?=?7.6 Tecomán earthquake struck offshore of the state of Colima, Mexico, near the diffuse triple junction between the Cocos, Rivera, and North American plates. Three-hundred and fifty aftershocks of the Tecomán earthquake with magnitudes between 2.6 and 5.8, each recorded by at least 7 stations, are relocated using the double difference method. Initial locations are determined using P and S readings from the Red Sismológica Telemétrica del Estado de Colima (RESCO) and a 1-D velocity model. Because only eight RESCO stations were operating immediately following the Tecomán earthquake, uncertainties in the initial locations and depths are fairly large, with average uncertainties of 8.0?km in depth and 1.4?km in the north?Csouth and east?Cwest directions. Events occurring between 24 January and 31 January were located using not only RESCO phase readings but also additional P and S readings from 11 temporary stations. Average uncertainties decrease to 0.8?km in depth, 0.3?km in the east?Cwest direction, and 0.7?km in the north?Csouth direction for events occurring while the temporary stations were deployed. While some preliminary studies of the early aftershocks suggested that they were dominated by shallow events above the plate interface, our results place the majority of aftershocks along the plate interface, for a slab dipping between approximately 20° and 30°. This is consistent with the slab positions inferred from geodetic studies. We do see some upper plate aftershocks that may correspond to forearc fault zones, and faults inland in the upper plate, particularly among events occurring more than 3?months after the mainshock.     

Connections between hydrodynamics, benthic landscape and associated fauna in the Balearic Islands, western Mediterranean

Linking the abiotic and biotic traits of ecosystems is a critical step towards understanding their structure and functioning. Here we attempt to determine the connections between the hydrodynamics, benthic landscape and the associated fish communities on the coastal continental shelf off the Balearic Islands (western Mediterranean). Specifically we investigate the role the hydrodynamics play in shaping the benthic landscape, and whether the hydrodynamics affect the composition and structure of demersal fish communities. A realistic numerical model was used to establish the hydrodynamic characteristics of the area. The study area showed high hydrodynamic variability on a medium spatial scale (tens of km) in terms of mean water velocity (ū). Principal component analysis was used to determine the main gradients of macro-epibenthic variability. Redundancy analysis (RDA) was used to model the effect of the hydrodynamics on macro-epibenthic species. RDA was also used to model the effect of the hydrodynamics and macro-epibenthos on the abundance of the associated fish fauna, and on its biomass at a community level using biomass spectra classes. The results showed that the hydrodynamics had a significant influence on the distribution of both macro-epibenthic species and the associated fish species. The latter was also influenced by the macro-epibenthos. Fish size appeared to be a key attribute for the distribution of species across gradients of ū and macro-epibenthic change. Our findings can be applied in ecosystem-based fisheries management, as they show that it is necessary to take into account both the biotic and abiotic traits of the habitats when the habitat use and requirements of the associated species are defined.     

Black shales and massive sulfide deposits: causal or casual relationships? Insights from Rammelsberg, Tharsis, and Draa Sfar

Black shales and massive sulfides represent reduced lithofacies that require isolation from oxic environments to be preserved. This, together with the sedimentary affinity of both lithofacies, can explain their common concurrence in the geologic record. The present study is based on the comparison of Rammelsberg in Germany, Tharsis in Spain, and Draa Sfar in Morocco, three massive sulfide deposits closely associated with black shales that are distributed along the European and North African Variscan orogen. The study entails geochemical, biostratigraphic, and stratigraphic analyses of the black shale sequences hosting the three deposits and mineralogical and textural analyses of the sulfides. All three deposits were formed in immature, tectonically unstable basins within an active continental margin or continental magmatic arc. Their stratigraphic records consist of a sequence of black shales enclosing massive sulfides and variable proportions of bimodal volcanic and subvolcanic rocks. The major differences among the three deposits concern the size, composition, and mineralogy. Regarding age, they are diachronous and younger southward: Rammelsberg is middle Eifelian, Tharsis latest Famennian, and Draa Sfar late Viséan. The study of redox conditions of the paleoenvironment using organic and inorganic proxies highlights similarities and significant differences among the three ore-hosting basins during massive sulfide and black shale deposition. The black shales generally display low C_org and high S_tot contents. At Rammelsberg, the S_tot/C_tot ratios provide values typical for normal Middle Devonian marine environments, which suggests that the original reactive organic C is now fixed in carbonates. At Tharsis, most of the samples have C_org >1 and S_tot/C_org values equivalent to those of Devonian?CCarboniferous normal marine sediments. However, some pyritic hanging-wall samples have C_org <1 and S_tot up to 5?wt.%, suggesting the epigenetic addition of HS^?. The S_tot/C_org ratio for the Draa Sfar samples resembles that of Middle Carboniferous normal marine environments. Geochemical inorganic proxies used to define the environmental conditions include the enrichment factors of U (U_EF) and Mo (Mo_EF) together with V/Cr and V/(V?+?Ni) ratios. Footwall shales at Filón Norte (Tharsis) show positive and eventually elevated U_EF and Mo_EF values, which suggests anoxic conditions, whereas at Rammelsberg and Draa Sfar oxic bottom water is indicated. The relations V/Cr and V/(V + Ni) in all three cases point to a redox boundary near the sediment?Cwater interface, although at Tharsis some samples indicate anoxic/euxinic conditions (i.e., V/(V + Ni) >0.9). Regarding the environmental conditions of the source areas, feldspar illitization and selective depletion in Na and Ca occurred at the three studies sites. Available sulfur isotopic data from the Rammelsberg and Tharsis sulfide ore indicate that biogenic reduction of marine sulfate was a major sulfur source during massive sulfide generation. Nevertheless, a hydrothermal sulfur source has also been detected. At Rammelsberg, this is indicated from the polymetallic sulfides that replace sedimentary and diagenetic pyrite. At Tharsis, the bacteriogenic sulfur signature is also restricted to sulfide with less evolved textures, whereas a hydrothermal source is more evident in sulfides showing evidence of recrystallization. Both geochemical and isotopic data suggest that the bacteriogenic reduction process was inhibited by rapid burial. The sedimentation rates calculated for Rammelsberg, Tharsis, and Draa Sfar were in the range 7?C13, 8?C14, and 19?C27?cm/ka, respectively. Continuous sedimentation of black shale favored the isolation of the massive sulfides and organic material from bottom waters and hence favored their preservation. Accordingly, the relationships between black shales and massive sulfides are considered to be casual. Nevertheless, the tectono-sedimentary evolution of each basin controlled the deposition of both black shales and massive sulfides and the parameters that favored their coeval deposition.