Seasonal and event-scale variations in solute chemistry for four Sierra Nevada catchments

Hydrobiogeochemical processes controlling stream water chemistry were examined in four small (<5 km²) catchments having contrasting bedrock lithologies in the western Sierra Nevada foothills of California. The Mediterranean climate with its cool/wet and hot/dry cycle produces strong seasonal patterns in hydrological, biological and geochemical processes. Stream water solutes fall into three general groups according to seasonal fluctuation in concentration: strong, rainy season minimum–dry season maximum (Cl⁻, SO₄²⁻, base cations); weak, rainy season minimum–dry season maximum (Si); and rainy season maximum–dry season minimum (NO₃⁻ and K⁺). Solute dynamics in soil solutions and stream water suggest that mixing of drainage waters from bedrock and soil sources regulate stream water solute concentrations. Patterns are further altered by the leaching of solutes accumulated in the soil over the summer period of desiccation and the temporal discoupling of nutrient cycles that occurs due to differences in the timing between vegetation growth (late spring) and leaching (early winter). Solute concentrations are remarkably similar between watersheds with varying bedrock types, with the exception of nitrate, sulfate and bicarbonate. Three watersheds have nitrogen-bearing metasedimentary bedrock that contributes to elevated nitrate concentrations in stream waters. Watersheds whose bedrock includes mineralized veins of sulfide and carbonate minerals similarly have greater sulfate and bicarbonate concentrations in stream water. Hydrobiogeochemical processes are highly dynamic at the seasonal and storm-event temporal scales and spatially complex at the watershed scale making management of stream water chemical composition, such as nitrate concentrations, very challenging.     

Hydrogeochemical and isotopic evidence of groundwater salinization in a coastal aquifer: a case study in Jeju volcanic island, Korea

In order to identify the origin of saline groundwater in the eastern part of Jeju volcanic island, Korea, a hydrogeochemical and isotopic study has been carried out for 18 observation wells located in east and southeast coastal regions. The total dissolved solid contents of groundwaters are highly variable (77–21,782 mg/l). Oxygen, hydrogen, sulfur, and strontium isotopic data clearly show that the saline water results from mixing of groundwater with seawater. Strontium isotopic compositions and Br/Cl and I/Cl ratios strongly suggest that the source of salinity is modern seawater intrusion. Hydrogeochemical characteristics based on bivariate diagrams of major and minor ions show that changes in the chemical composition of groundwater are mainly controlled by the salinization process followed by cation-exchange reactions. The highly permeable aquifers at the east coastal region are characterized by low hydraulic gradient and discharge rate and high hydraulic conductivity as compared with other regions. These properties enhance the salinization of groundwater observed in the study area. Based on the Cl, Br, and δ¹⁸O data, seawater was determined to have intruded inland some 2.5 km from the coastline. Considering the poor correlation of sampling depth and Cl concentrations observed, the position of seawater-freshwater interface is not uniformly distributed in the study area, due to heterogeneities of the basaltic aquifers.     

Using lumped conceptual rainfall-runoff models to simulate daily isotope variability with fractionation in a nested mesoscale catchment

This paper presents 19 months of stable isotope (δ²H and δ¹⁸O) data to enhance understanding of water and solute transport at two spatial scales (2.3 km² and 122 km²) in the agricultural Lunan catchment, Scotland. Daily precipitation and stream isotope data, weekly lake and spring isotope data and monthly groundwater isotope data revealed important insights into flow pathways and mixing of water at both scales. In particular, a deeper groundwater flow path significantly contributes to total streamflow (25-50%). Upstream lake isotope dynamics, susceptible to evaporative fractionation, also appeared to have an important influence on the downstream isotope composition. This unique tracer data set facilitated the conceptualization of a lumped catchment-scale flow-tracer model. The incorporation of hydrological, mixing and fractionation processes based on these data improved simulations of the stream δ²H isotope response at the catchment outlet from 0.37 to 0.56 for the Nash-Sutcliffe statistic. The stable isotope data successfully aided model conceptualization and calibration in the quest for a simple water and solute transport model with improved representation of process dynamics.     

Drying and dying of a subducted slab: Coupled Li and B isotope variations in Western Anatolia Cenozoic Volcanism

In lavas spanning ~ 10Ma of subduction-related volcanism in Western Anatolia, we observe remarkably similar patterns of δ⁷Li and δ¹¹B variation. In this setting, magmatism records a transition from calc–alkaline to ultrapotassic character, consistent with overall lower mean extents of melting, and a changing mantle source that reflects a fractionating, higher temperature slab input consistent with the gradual cessation of subduction. Subsequent rift-related intraplate magmatism record δ⁷Li signatures within the range observed for MORBs and OIBs, indicating an abrupt transition to a mantle source unmodified by subduction.     

Fractal variations of the Transcarpathians,West Ukraine,seismicity and their potential relation to changing phases of local seismic cycles

The evolution of spatial arrangement of the Transcarpathian earthquakes is studied by means of fractal analysis. Spatial correlation dimension is evaluated for a subset of fixed number of catalogue events consecutively shifted in time by one event forward. Distinctive patterns that may be related to different phases of local seismic cycle/cycles have been identified, corresponding to extreme values of the dimension. Their potential relation to evolution of the seismogenic stress field is discussed. The time it takes for the capacity dimension to reach the so-called percolation threshold is estimated for a subset consecutively extended by adding next catalogue events. Temporal and spatial frames of percolation analysis are discussed in relation to parameters of the future earthquake.     

Wind-wave affected phosphate loading variations and their relationship to redox condition in Lake Taihu

Variation of wind speed and the physico-chemical parameters, such as dissolved phosphate, ferrous and manganese in lake water were observed on site. Together with the chemistry analysis and simulated experiment in lab, the change of phosphate concentration in lake water was analyzed. The variation of ferrous/phosphate ratio explained that along with the enhancement of wind-wave effect and the oxidation ability of lake water, the effects of co-precipitation and removal of dissolved phosphate and iron in the lake water were reinforced. The ferrous/phosphate ratio in pore water was less than 2.0, demonstrating that the dissolved phosphate can be released into the overlying water. But, in the lake water, the stability of phosphate was controlled by the water dynamics. The phosphate release experiment showed that molecular release was only part of the whole and the direct discharge of phosphate in the pore water was also a part. The mineralization and biological process of suspended particulates in the water may be another important reason for the whole phosphate loadings.     

Relationships between deformation and mixing processes in lava flows: a case study from Salina (Aeolian Islands,Tyrrhenian Sea)

 The massive unit of a lava flow from Porri volcano (Salina, Aeolian Islands) displays many unusual structures related to the physical interaction between two different magmas. The magma A represents approximately 80% of the exposed lava surface; it has a crystal content of 51 vol.% and a dacitic glass composition (SiO₂=63–64 wt.%). The magma B has a basaltic-andesite glass composition (SiO₂=54–55 wt.%) and a crystal content of approximately 18 vol.%. It occurs as pillow-like enclaves, banding, boudin-like and rolling structures which are hosted in magma A. Structural analysis suggests that banding and boudin-like structures are the result of the deformation of enclaves at different shear strain. The linear correlation between strain and stratigraphic height of the measured elements indicates a single mode of deformation. We deduce that the component B deformed according to a simple shear model. Glass analyses of the A–B boundary indicate that A and B liquids mix together at high shear strain, whereas only mingling occurs at low shear strain. This suggests that the amount of deformation (i.e. forced convection) plays an important role in the formation of hybrid magmas. High shear strain may induce stretching, shearing and rolling of fluids which promote both forced convection and dynamical diffusion processes. These processes allow mixing of magmas with large differences in their physical properties. Received: 15 July 1995 / Accepted: 30 May 1996     

Soil moisture variation in relation to topography and land use in a hillslope catchment of the Loess Plateau, China

The profile characteristics and the temporal dynamics of soil moisture variation were studied at 26 locations in Da Nangou catchment (3.5 km²) in the loess area of China. Soil moisture measurements were performed biweekly at five depths in the soil profile (0–5, 10–15, 20–25, 40–45 and 70–75 cm) from May to October 1998 using Delta-T theta probe. Soil moisture profile type and temporal variation type and their relationship to topography and land use were identified by detrended canonical correspondence analysis (DCCA) and correlation analysis. The profile distribution of time-averaged soil moisture content can be classified into three types i.e. decreasing-type, waving-type and increasing-type. The profile features of soil moisture (e.g. profile gradient and profile variability) are influenced by different environmental factors. The profile type of soil moisture is only attributed to land use while profile gradient and profile variability of soil moisture is mainly related to land use and topography (e.g. landform type and slope). The temporal dynamics of layer-averaged soil moisture content is grouped into three types including three-peak type, synchro-four-peak type and lagged-four-peak type. These types are controlled by topography rather than by land use. The temporal dynamic type of soil moisture shows significant correlation with relative elevation, slope, aspect, while temporal variance displays significant relation with slope shape. The mean soil moisture is related to both the profile and dynamics features of soil moisture and is controlled by both land use and topography (e.g. aspect, position, slope and relative elevation). The spatial variability of soil moisture across landscape varies with both soil depths and temporal evolution.     

Short term fluctuations in dissolved organic matter concentrations in streamflow draining a forested watershed and their relation to the catchment budget

Dissolved organic matter (DOM) concentrations in throughfall, throughflow, and runoff were recorded in a small (95 ha) woodland catchment in North Warwickshire for a period of eight weeks following a prolonged dry spell. DOM is shown to be positively related to stream discharge over the sampling period, although distinctive flushing effects were identified. The net contribution of DOM to total dissolved solids load carried in the river was only 2.4 per cent, and represented significantly less than published carbon losses by mineralization of soil organic matter. Throughfall inputs were some 100 times higher than streamflow outputs from the catchment.     

Oxygen and carbon isotopic characteristics of rainwater, drip water and present speleothems in a cave in Guilin area, and their environmental meanings

The studies of the oxgen and carbon isotopes of the rainwater in Guilin area, the drip water and the present carbonate deposit in Panlong cave of Guilin show that: (i) as to the general characteristics of the oxygen isotopes of the rainwater within a year and between years, the δ18O values decrease with an increase of air temperature and the rainfall, and the correlation between δ18O values and the mean monthly air temperature is much better than that between δ18O values and the rainfall, and the δ18O values of the rainwater during the summer monsoon are much smaller than those during winter monsoon; (ii) δ18O values of the drip water have a quite good correlation with the δ18O values of the rainwater in the same period; (iii) when the conditions are appropriate, δ13C can be used as an environmental proxy, that is, the smaller δ13C of speleothems is, the larger the proportion of C3 plants is and the more plentiful the rainfall is. On the contrary, C4 plants may be prevailing or the environment may be a stony desert caused by climate changes or human activity.     

Using stable isotopes to estimate travel times in a data‐sparse Arctic catchment: Challenges and possible solutions

Use of isotopes to quantify the temporal dynamics of the transformation of precipitation into run‐off has revealed fundamental new insights into catchment flow paths and mixing processes that influence biogeochemical transport. However, catchments underlain by permafrost have received little attention in isotope‐based studies, despite their global importance in terms of rapid environmental change. These high‐latitude regions offer limited access for data collection during critical periods (e.g., early phases of snowmelt). Additionally, spatio‐temporal variable freeze–thaw cycles, together with the development of an active layer, have a time variant influence on catchment hydrology. All of these characteristics make the application of traditional transit time estimation approaches challenging. We describe an isotope‐based study undertaken to provide a preliminary assessment of travel times at Siksik Creek in the western Canadian Arctic. We adopted a model–data fusion approach to estimate the volumes and isotopic characteristics of snowpack and meltwater. Using samples collected in the spring/summer, we characterize the isotopic composition of summer rainfall, melt from snow, soil water, and stream water. In addition, soil moisture dynamics and the temporal evolution of the active layer profile were monitored. First approximations of transit times were estimated for soil and streamwater compositions using lumped convolution integral models and temporally variable inputs including snowmelt, ice thaw, and summer rainfall. Comparing transit time estimates using a variety of inputs revealed that transit time was best estimated using all available inflows (i.e., snowmelt, soil ice thaw, and rainfall). Early spring transit times were short, dominated by snowmelt and soil ice thaw and limited catchment storage when soils are predominantly frozen. However, significant and increasing mixing with water in the active layer during the summer resulted in more damped steam water variation and longer mean travel times (~1.5 years). The study has also highlighted key data needs to better constrain travel time estimates in permafrost catchments.     

What drives high flow events in the Swiss Alps? Recent developments in wavelet spectral analysis and their application to hydrology

Extreme hydrological events are often triggered by exceptional co-variations of the relevant hydrometeorological processes and in particular by exceptional co-oscillations at various temporal scales. Wavelet and cross wavelet spectral analysis offers promising time-scale resolved analysis methods to detect and analyze such exceptional co-oscillations. This paper presents the state-of-the-art methods of wavelet spectral analysis, discusses related subtleties, potential pitfalls and recently developed solutions to overcome them and shows how wavelet spectral analysis, if combined to a rigorous significance test, can lead to reliable new insights into hydrometeorological processes for real-world applications. The presented methods are applied to detect potentially flood triggering situations in a high Alpine catchment for which a recent re-estimation of design floods encountered significant problems simulating the observed high flows. For this case study, wavelet spectral analysis of precipitation, temperature and discharge offers a powerful tool to help detecting potentially flood producing meteorological situations and to distinguish between different types of floods with respect to the prevailing critical hydrometeorological conditions. This opens very new perspectives for the analysis of model performances focusing on the occurrence and non-occurrence of different types of high flow events. Based on the obtained results, the paper summarizes important recommendations for future applications of wavelet spectral analysis in hydrology.     

K/Ar dates of hydrothermal clays from core hole VC-2B, Valles caldera, New Mexico and their relation to alteration in a large hydrothermal system

Seventeen K/Ar dates were obtained on illitic clays within Valles caldera (1.13 Ma) to investigate the impact of hydrothermal alteration on Quaternary to Precambrian intracaldera and pre-caldera rocks in a large, long-lived hydrothermal system ( 1.0 Ma to present). Clay samples came from scientific core hole VC-2B (295°C at 1762 m) which was spudded in the Sulphur Springs thermal area and drilled into the boundary between the central resurgent dome and the western ring-fracture zone. Six illitic clays within Quaternary caldera-fill debris flow, tuffaceous sediment, and ash-flow tuff (48 to 587 m depth) yield ages from 0.35 to 1.09 Ma. Illite from Miocene pre-caldera sandstone (765 m) gives an age of 6.74 Ma. Two dates on illite from sandstones in Permian red beds (1008 and 1187 m) are 4.33 and 4.07 Ma, respectively. Surprisingly, three dates on illites from altered andesite pebbles within the red beds (1010–1014 m) are 0.95 to 1.06 Ma. Four illite dates on variably altered Precambrian quartz monzonite (1615–1762 m) range from 2.90 to 276 Ma.Post-Valles age illite is not correlated with alteration style (argillic to propylitic). Rather, post-Valles ages are uniformly obtained from illites in highly fractured, intensely altered, caldera-fill rocks and the Permian volcanic clasts. Generally, finer clay fractions from identical samples yield younger ages. Plots of ⁴⁰Ar/³⁶Ar versus ⁴⁰K/³⁶Ar and ⁴⁰Ar* versus ⁴⁰K for the illites in caldera-fill rocks lie close to a 1-Ma isochron. Most illite dates older than Valles caldera are difficult to interpret because they correspond to the ages of pre-Valles volcanic and hydrothermal episodes in the Jemez volcanic field ( 13 Ma). In addition, older dates may be caused by co-mingling of different illites during sample preparation, or by inherited argon or lost argon in illites from rocks with potentially complex hydrothermal histories. However, the range of ages obtained from illites in Permian sands and pebbles and from Precambrian crystalline rocks indicates that Valles hydrothermal activity is overwhelming illite produced by earlier geologic events.     

The development of a rapid screening technique to measure antibiotic activity in effluents and surface water samples

Use of antibiotics in medicine, farming and aquaculture results in a continual supply of these pharmaceuticals and their breakdown products reaching the environment. A method has been developed to use for screening environmental samples to establish whether they contain antibiotics. The assay is based on a commercially available test kit and uses some modifications of the test procedure and SPE techniques to deliver a fast, semi-quantifiable result. Development of the assay and results from spiked and environmental samples are discussed.     

Habitat suitability of the invasive water hyacinth and its relation to water quality and macroinvertebrate diversity in a tropical reservoir

In this study, we assessed the relationship between the occurrence of the invasive water hyacinth (Eichhornia crassipes) and water quality properties as well as macroinvertebrate diversity in a tropical reservoir, situated in western Ecuador. Macroinvertebrates and physico-chemical water quality variables were sampled at 32 locations (during the dry season of 2013) in both sites covered and non-covered by water hyacinth in the Daule-Peripa reservoir. The results indicated that, in terms of water quality, only turbidity was significantly different between sampling sites with and without water hyacinth (Mann–Whitney U-test, p < 0.01). The habitat suitability model showed that water hyacinth was present at sites with a low turbidity. The percentage water hyacinth cover increased with decreasing turbidity. The Biological Monitoring Working Party-Colombia score and the Margalef diversity index were significantly higher (Mann–Whitney U-test, p < 0.01) at sampling sites where water hyacinth was present compared to water hyacinth absent sites. However, there were no significant differences in the Shannon–Wiener index, Evenness index and Simpson index between the sampling sites with and without water hyacinth. Our results suggest that water hyacinth cover was an important variable affecting the diversity of macroinvertebrates in the Daule-Peripa reservoir, with intermediate levels of water hyacinth cover having a positive effect on the diversity of macroinvertebrates. Information on the habitat suitability of water hyacinth and its effect on the physico-chemical water quality and the macroinvertebrate community are essential to develop conservation and management programs for large tropical reservoirs such as the Daule-Peripa reservoir and the Guayas river basin, where water resources are being at high risk due to expansion of agricultural and industrial development activities.     

Dust pollen distribution on a continental scale and its relation to present-day vegetation along north-south transects in east China

A series of dust pollen samples was collected along N-S transects in east China (18°N to 53°N latitudes). Sample sites extend from the cold-temperate zone in the north to the tropical region in the south. Pollen taxa characterize each region and reflect the natural and devastated vegetation as well as corre- sponding climatic zones. The quantitative pollen data can be used to estimate the spatial distribution of planted and introduced species. Valuable information of human disturbance of the natural forest is evaluated by quantitative comparison between dust pollen and in-situ pollen of protected forest. In addition, percentages of grass pollen vary regularly from north to south that is consistent with spatial distribution of net primary productivity in east China. Among all grasses, Artemisia and the Gramineae carry the clearest signal: their ratio increases northwards and therefore represent a suitable and con- venient tool for palaeoclimate reconstructions.     

Significant variations in the productivity of green macroalgae in a mesotidal estuary: implications to the nutrient loading of the system and the adjacent coastal area

A spatially dynamic model for the productivity of spores and adults of green macroalgae (Enteromorpha sp.) was developed for a mesotidal estuary (Mondego estuary, Portugal). Many of the algal processes and parameters included in the model were experimentally obtained. Model predictions were compared to a real time series (1993-1997) of macroalgal biomass variation and the two sets show a good agreement (ANOVA, P<0.001). Results suggest that algal growth is highly sensitive to small changes in depth and exhibits different patterns of variation in different seasons. On a yearly basis, global calculations for the south channel of the estuary (137 ha) suggest that during bloom years, macroalgal biomass may reach about 21,205 ton DW compared to 240 ton DW in regular years. On a seasonal basis, the difference may be even more significant. The consequences of such variations on the nitrogen and phosphorus loading of the system and the adjacent coastal area are discussed.     

Diatom thanatocoenosis in a middle Galician Ría: Spatial patterns and their relationship to the seasonal diatom cycle in the water column and hydrographic conditions

Modern diatom distribution patterns in the surface sediment of the Ferrol Ría and their relationship to the hydrography and diatom patterns in the water column were studied to determine the hydrographic influence on the record of these biogenic components. Diatom abundance in the water column was assessed for different oceanographic periods and compared with the biosiliceous sedimentary record. Very low abundances were found in the water column during the winter, whereas in spring and summer, diatoms proliferated. Chaetoceros spp. formed the bulk of the water column community during spring and summer, followed by Thalassionema nitzschioides and Rhizosolenia spp. Nitzschia longissima represented a significant portion of the winter assemblage, together with Paralia sulcata and benthic taxa. Leptocylindrus danicus, N. longissima and Skeletonema costatum characterized the autumn campaign, when stratification of the waters occurs, with L. danicus being especially abundant in the outer ría.Seasonal hydrographic and associated productivity patterns govern the abundance and assemblage of the diatoms preserved in surface sediments. Samples located in the inner ría area and its margins exhibited the highest abundances of diatoms, and were primarily dominated by benthic species. The freshwater group, crysophycean cysts and phytoliths were present in the landward stations influenced by river runoff. The middle ría was characterized by P. sulcata and Thalassiosira spp., with minor occurrences of the benthic and freshwater group. Chaetoceros R.S., L. danicus R.S. and T. nitzschioides typified the outer ría, an assemblage that corresponds to nutrient-rich coastal areas of high productivity influenced by oceanic waters, demonstrating the impact of oceanic waters flowing into the embayment due to enhanced tidal mixing through the narrow channel. Therefore, sediment diatom assemblages reflect diatom production patterns in the water column of the Ría. However, we must proceed with caution when interpreting the paleorecord in the inner area due to the high contribution of allochthonous taxa, which is indicative of low water depths. This paper contributes to a better understanding of diatom thanatocoenosis in the Galician Rías, where very few studies of this kind have been done to date.     

Simulated changes due to global warming in the variability of precipitation, and their interpretation using a gamma-distributed stochastic model

The interannual variability of monthly mean January and July precipitation and its possible change due to global warming are assessed using a five-member ensemble of climate for the period 1871–2100, simulated by the CSIRO Mark 2 global coupled atmosphere–ocean model. In the 1961–1990 climate, for much of the middle to high latitudes the standard deviation of precipitation for both months is roughly proportional to the mean, with the coefficient of variation (C) typically 0.3–0.5. The variability there is shown to be largely consistent with that from a first-order Markov chain model of the daily rainfall occurrence, with the distribution of wet-day amounts approximated by a gamma distribution. Global distributions of Mark 2-based parameters of this stochastic model, commonly used in weather generators, are presented. In low latitudes, however, the variability from the coupled model is typically double that anticipated by the stochastic model, as quantified by an ‘overdispersion ratio’. C often exceeds one at subtropical locations, where rain is less frequent, but sometimes relatively heavy.The standard deviation of monthly mean precipitation S generally increases as the global model warms, with the global mean S in 2071–2100 in January (July) being 9.0% (11.5%) larger than in 1961–1990. Decreases in some subtropical locations occur, particularly where mean precipitation decreases. The global pattern of overdispersion is largely unchanged, however, and the changes in S can be related to those in the stochastic model parameters. Much of the increase in S is associated with increases in the scale parameter of the gamma distribution of wet-day amounts. Changes in C, which is unaffected by this parameter, are generally small. Increases in C in several subtropical bands and over northern midlatitude land in July are related to a decreased frequency of precipitation, and (to a lesser degree) changes in the gamma shape parameter. Some potential applications of the results to downscaling are discussed, and illustrated using observed rainfall from southeast Australia.