Contrasting responses of dimictic and polymictic lakes to environmental change: a spatial and temporal study

Although comparative analyses between dimictic and polymictic lakes have noted the key role of mixing regime in governing water quality, limnologists have historically focused on dimictic lakes, leaving polymictic lakes relatively understudied. In this study, we investigated whether the effects of agricultural development on water quality differed between dimictic and polymictic lakes in a naturally nutrient-rich region of Alberta, Canada. Through a spatial limnological analysis of 36 sites, we found that the relationship between surface water total phosphorus concentration and the percent of agriculture (% Agr) in the catchments differed between dimictic and polymictic lakes, where the proportion of variance explained was much more pronounced in the dimictic (79% explained) than in the polymictic systems (7% explained). Paleolimnological analyses of subfossil chironomids in surface sediment samples (0–1 cm) from 18 of the 36 study lakes, and in sediment core profiles from the dimictic and polymictic basins of a eutrophic lake, showed that water quality differed between mixis groups. We found that the surface sediment chironomid assemblages differed significantly between polymictic and dimictic lakes. Detailed analyses of the sediment cores demonstrated that the two basin types differed in their responses to land-use change through time, as only the dimictic basin showed a parallel increase in anoxia with increasing agricultural development. We suggest that in naturally-fertile landscapes, external nutrient loading exerts a more notable effect on dimictic lakes, whereas internal nutrient loading is more important in polymictic systems.     

On the combination of global and local data in collocation theory

Due to the successful operation of dedicated satellite gravity missions, nowadays high-accuracy global gravity field models have become available. This triggers the challenge to optimally combine this long to medium wavelength gravity field information derived from space-borne data with high-resolution terrestrial gravity data. In this paper, the least squares collocation concept is revised with the attempt to consistently unify the combination procedure in such a way that the full information contained in both data sets is merged. For example, in local or regional geoid determination the remove-restore method is usually applied only partially taking into account the accuracy of the global model coefficients used for the long-wavelength reduction. The key advantage of the extended formulation is the fact that it automatically accounts for the error covariance of all data types involved. The applicability, feasibility and performance of the proposed method is investigated in the frame of numerical closed-loop simulations. The two main fields of application, i.e., the improvement of a global gravity field model by terrestrial gravity field data, and, vice versa, the support to a regional geoid solution by the incorporation of a global gravity field model, have been analyzed and assessed. Although applied under simplified conditions, it could be shown that the method works and is practically applicable.     

Climate changes during the last glacial termination inferred from diatom-based temperatures and pollen in a sediment core from Längsee (Austria)

A sediment core section from Längsee, a small meromictic lake in the southern Alpine lowland (Carinthia, Austria) close to the Würmian ice margin, was investigated by means of diatoms and pollen. The main aims of the study were to reconstruct water temperature as a signal of climate change during the last glacial termination, compare the aquatic and terrestrial response to the changing climate, and place our findings into a climatic frame on the northern hemispheric scale. A calibration data set (ALPS06) of 116 lakes was constructed using data from newly studied lakes and from two previously published data sets and we established a transfer function for predicting summer epilimnetic water temperatures (SEWT). A locally weighted weighted average regression and calibration model (R _jack ²  = 0.89; RSMEP = 1.82°C) was applied to the fossil diatom assemblages in order to reconstruct SEWT. Three major sections were distinguished in the time window of approximately 19–13 cal ka BP, which fitted well with the oxygen isotope curve and the isotope-event stratigraphy from the Greenland ice-core GRIP. The first section was a warming period (SEWT range from 11.6 to 18.0°C; average 15.8°C = ca. 6°C below present) called the Längsee oscillation, which probably correlates with the warmer sub-section (GS-2b) of the Greenland Stadial 2. The subsequent section represents a climate cooling, called the Längsee cold period (SEWT range between 10.6 and 15.9°C; average 12.9°C), which probably corresponds with the sub-section GS-2a of the Greenland Stadial 2, the Heinrich 1 cold event of the North Atlantic, and partially the Gschnitz Stadial in the Alps. The Längsee cold period shows a tri-partition: Two colder phases are separated by a warmer inter-phase. The passive ordination of the core sample scores along maximum water depth indicated that the Längsee cold period was drier than the Längsee oscillation. Strong short-term fluctuations during the Längsee oscillation and the Längsee cold period indicate climate instability. The third section represented climate warming during the Längsee late glacial interstadial (=Greenland Interstadial 1, GI-1) with an average SEWT of 17.5°C. From the minor climatic fluctuations during this interstadial, mainly indicated by pollen, the fluctuation most likely related to the Gerzensee oscillation showed a SEWT decline. During the early immigration and expansion period of shrubs and trees, aquatic and terrestrial records showed distinct discrepancies that might have arose because of time lags in response and differences in sensitivity.     

Late quaternary sea‐level changes,crustal movements and coastal evolution in Northumberland,UK

Northumberland lies in the transition between Holocene emergence and submergence and is thus a critical zone for testing models of isostatic rebound. We have collected data from this area to reconstruct relative sea‐level changes and lateral coastline movements for the last 14000 y. These are deposits from tidal marsh, back‐barrier wetland and terrestrial environments producing 47 sea‐level index points from 12 sites. There is no unequivocal evidence for Late Devensian sea levels above present and the reliable sea‐level index points are restricted between −6 m and +2.5 m relative to present and 9.0–2.5 kyr cal. BP. Analysis of these quantifies differential responses to glacio‐ and hydroisostatic rebound, with the northern sites recording a mid‐Holocene sea‐level maximum ca. 2.5 m above present, whereas the southern sites show a maximum ca. 0.5 m above present. These observations show a reasonable fit with the predictions from quantitative models of glacio and hydroisostatic rebound, but there is currently no unique solution of Earth and ice model parameters that will explain all the sea‐level observations. Copyright © 2000 John Wiley & Sons, Ltd.     

Bayesian chemistry-assisted hydrograph separation (BACH) and nutrient load partitioning from monthly stream phosphorus and nitrogen concentrations

A Bayesian chemistry-assisted hydrograph separation (BACH) approach was developed, based on calibration of a three-component recursive digital filter, that requires monthly water quality data only. This enables BACH to be applied to the large number of rural catchments for which continuous flow records and monthly water chemistry time series exist from ‘state of environment’ monitoring programmes, but little supplementary data required for more sophisticated analysis techniques. As well as estimating fast, medium, and slow flow components, BACH also estimates for each flow component a time-invariant concentration of the chemical tracers chosen, allowing flow path-specific loads to be calculated. The method was demonstrated using 15 years of total phosphorus (TP) and total nitrogen (TN) data from eight mesoscale catchments in the Waikato region of New Zealand’s North Island. Calibration was done separately for three 5-year data periods, and validated against data from the following 5-year period. Flow path separation and concentration predictions were consistent between data periods, indicating that the TP–TN combination contained sufficient information to reliably identify three flow paths in each catchment; an event-response near-surface flow path with high concentrations of both phosphorus and nitrogen, a seasonal shallow groundwater flow path with lower concentrations of TP but high concentrations of TN, and a deeper slower groundwater flow path characterised by generally low concentrations of both TP and TN. Based on this analysis, the catchments were able to be grouped in three hydro-types. This shows that commonly available water quality data can support robust, objective flow separation and nutrient load apportionment, even in the absence of other supporting data, provided appropriate modelling methods are used.     

Fault slip and identification of the second fault plane in the Varzeghan earthquake doublet

An intraplate earthquake doublet, with 11-min delay between the events, devastated the city of Varzeghan in northwestern Iran on August 11, 2012. The first Mw 6.5 strike-slip earthquake, which occurred after more than 200 years of low seismicity, was followed by an Mw 6.4 oblique thrust event at an epicentral separation of about 6 km. While the first event can be associated with a distinct surface rupture, the absence of a surface fault trace and no clear aftershock signature makes it challenging to identify the fault plane of the second event. We use teleseismic body wave inversion to deduce the slip distribution in the first event. Using both P and SH waves stabilize the inversion and we further constrain the result with the surface rupture extent and the aftershock distribution. The obtained slip pattern shows two distinct slip patches with dissimilar slip directions where aftershocks avoid high-slip areas. Using the estimated slip for the first event, we calculate the induced Coulomb stress change on the nodal planes of the second event and find a preference for higher Coulomb stress on the N-S nodal plane. Assuming a simple slip model for the second event, we estimate the combined Coulomb stress changes from the two events on the focal planes of the largest aftershocks. We find that 90% of the aftershocks show increased Coulomb stress on one of their nodal planes when the N-S plane of the second event is assumed to be the correct fault plane.     

Data requirements for a 5-mm quasi-geoid in the Netherlands

We assess the surface gravity data requirements for a 5-mm quasi-geoid model for the Netherlands mainland and continental shelf in terms of omission and commission errors. The omission error critically depends on the roughness of the topography and bathymetry. For the Netherlands continental shelf, Central and Northern Netherlands, the omission error is well described by the model 0.32d mm, where d is the data spacing in km. For the more hilly Southern Netherlands, the omission error model is 0.92d mm. The commission error depends on the kernel modification, the data spacing, and the data accuracy. When using the spheroidal Stokes kernel, it is well described by 0.277 d σ_Δg mm, where σ_Δg is the noise standard deviation of surface gravity data in mGal. An upper bound of the commission error of the state-of-the-art satellite-only gravity model GOCO05S over the Netherlands is e^{0.03676L–11.419} mm, where L is the maximum degree up to which this model is used. Only if this model is truncated at a sufficiently low degree, e.g., at degree 100, its contribution to the total commission error can be neglected. We determine the total error as the sum of commission and omission errors. Hence, to realize a 5-mm quasi-geoid model for the Netherlands mainland and continental shelf, a data spacing of 3.5 km is needed when assuming a noise standard deviation of 1.5 mGal for surface gravity data. The currently available land-based gravity data fulfill this requirement. This does not apply to the situation at sea, where the density of the shipboard gravity data and the accuracy of the radar altimeter-derived data do not allow the realization of a 5-mm quasi-geoid model.     

Chemical Characterisation of NIST Silicate Glass Certified Reference Material SRM 610 by ICP-MS, TIMS, LIMS, SSMS, INAA, AAS and PIXE

National Institute of Science and Technology (NIST) silicate glass SRM 610 is widely used as a certified reference material for various micro-analytical techniques such as SIMS or laser ablation ICP-MS. SRM 610 has been nominally doped with sixty one trace elements at the 500 μg g⁻¹ level, but certified concentration data exist for only a few of these elements. This study reports concentration data for fifty nine trace elements obtained by ICP-MS, SSMS, LIMS, TIMS, INAA, AAS, and PIXE analyses of two different SRM 610 wafers. Most elements fall within a 10% band around a median value of about 440 μg g⁻¹. The REE concentrations are shown to be constant to 3% (1 σ), thus emphasizing the value of SRM 610 as a reference material for REE analyses.
Comparison of our values with published data suggests that different SRM 610 wafers are, within errors, chemically identical for most elements. Exceptions to this general rule appear to be restricted to elements which were partly lost during the production of the glass, e.g. Ag and Br. On the basis of six independent determinations of Rb concentrations, which are systematically lower by a few percent than the reported NIST value, we argue that the certified Rb concentration may not be representative for all distributed SRM 610 wafers.     

Groundwater—Surface waters interactions at slope and catchment scales: implications for landsliding in clay‐rich slopes

Understanding water infiltration and transfer in soft‐clay shales slopes is an important scientific issue, especially for landsliding. Geochemical investigations are carried out at the Super‐Sauze and Draix‐Laval landslides, both developed in the Callovo‐Oxfordian black marls, with the objective to define the origin of the groundwater. In situ investigations, soil leaching experiments and geochemical modeling are combined to identify the boundaries of the hydrological systems. At Super‐Sauze, the observations indicate that an external water flow occurs in the upper part of the landslide at the contact between the weathered black marls and the overlying formations, or at the landslide basement through a fault network. Such external origin of water is not observed at the local scale of the Draix‐Laval landslide but is detected at the catchment scale with the influence of deep waters in the streamwater quality of low river flows. Hydrogeological conceptual models are proposed emphasizing the role of the interactions between local (slope) and regional (catchment) flow systems. The observations suggest that this situation is a common case in the Alpine area. Expected consequences of the regional flows on slope stability are discussed in term of rise of pore water pressures and physicochemical weathering of the clay shales.