On characterizing the temporal dominance patterns of model parameters and processes

Diagnostic analyses of hydrological models intend to improve the understanding of how processes and their dynamics are represented in models. Temporal patterns of parameter dominance could be precisely characterized with a temporally resolved parameter sensitivity analysis. In this way, the discharge conditions are characterized, that lead to a parameter dominance in the model. To achieve this, the analysis of temporal dynamics in parameter sensitivity is enhanced by including additional information in a three‐tiered framework on different aggregation levels. Firstly, temporal dynamics of parameter sensitivity provide daily time series of their sensitivities to detect variations in the dominance of model parameters. Secondly, the daily sensitivities are related to the flow duration curve (FDC) to emphasize high sensitivities of model parameters in relation to specific discharge magnitudes. Thirdly, parameter sensitivities are monthly averaged separately for five segments of the FDC to detect typical patterns of parameter dominances for different discharge magnitudes. The three methodical steps are applied on two contrasting catchments (upland and lowland catchment) to demonstrate how the temporal patterns of parameter dynamics represent different hydrological regimes. The discharge dynamic in the lowland catchment is controlled by groundwater parameters for all discharge magnitudes. In contrast, different processes are relevant in the upland catchment, because the dominances of parameters from fast and slow runoff components in the upland catchment are changing over the year for the different discharge magnitudes. The joined interpretation of these three diagnostic steps provides deeper insights of how model parameters represent hydrological dynamics in models for different discharge magnitudes. Thus, this diagnostic framework leads to a better characterization of model parameters and their temporal dynamics and helps to understand the process behaviour in hydrological models. Copyright © 2015 John Wiley & Sons, Ltd.     

Late Glacial ice advances in southeast Tibet

The sensitivity of Tibetan glacial systems to North Atlantic climate forcing is a major issue in palaeoclimatology. In this study, we present surface exposure ages of erratic boulders from a valley system in the Hengduan Mountains, southeastern Tibet, showing evidence of an ice advance during Heinrich event 1. Cosmogenic nuclide analyses (¹⁰Be and ²¹Ne) revealed consistent exposure ages, indicating no major periods of burial or pre-exposure. Erosion-corrected (3 mm/ka) ¹⁰Be exposure ages range from 13.4 to 16.3 ka. This is in agreement with recalculated exposure ages from the same valley system by [Tschudi, S., Schäfer, J.M., Zhao, Z., Wu, X., Ivy-Ochs, S., Kubik, P.W., Schlüchter, C., 2003. Glacial advances in Tibet during the Younger Dryas? Evidence from cosmogenic ¹⁰Be, ²⁶Al, and ²¹Ne. Journal of Asian Earth Sciences 22, 301–306.]. Thus this indicates that local glaciers advanced in the investigated area as a response to Heinrich event 1 cooling and that periglacial surface adjustments during the Younger Dryas overprinted the glacial morphology, leading to deceptively young exposure ages of certain erratic boulders.     

A two-layer primitive equation model of an idealized Strait of Gibraltar connected to an eastern basin

We describe the solutions of a numerical two-layer primitive equation model of an idealized Strait of Gibraltar and an adjacent eastward basin. The quasi-steady circulation pattern in the basin features an anticyclonic gyre southward of the strait exit and an eastward boundary current attached to the southern boundary of the basin. A variation of the initial upper-layer depth and the target interface height at the eastern boundary causes minor changes of the upper-layer circulation in a basin with rectangular coastline and larger changes when a cape is added to the coastline of the southern boundary.     

The Use of Resonant Orbits in Satellite Geodesy: A Review

Dynamic resonance, arising from commensurate (orbital or rotational) periods of satellites or planets with each other, has been a strong force in the development of the solar system. The repetition of conditions over the commensurate periods can result in amplified long-term changes in the positions of the bodies involved. Such resonant phenomena driven by the commensurability between the mean motion of certain artificial Earth satellites and the Earth’s rotation originally contributed to the evaluation and assessment of the Stokes parameters (harmonic geopotential coefficients) that specify the Earth’s gravitational field. The technique constrains linear combinations of the harmonic coefficients that are of relevant resonant order (lumped coefficients). The attraction of the method eventually dwindled, but the very accurate orbits of CHAMP and GRACE have recently led to more general insights for commensurate orbits applied to satellite geodesy involving the best resolution for all coefficients, not just resonant ones. From the GRACE mission, we learnt how to explain and predict temporary decreases in the resolution and accuracy of derived geopotential parameters, due to passages through low-order commensurabilities, which lead to low-density ground-track patterns. For GOCE we suggest how to change a repeat orbit height slightly, to achieve the best feasible recovery of the field parameters derived from on-board gradiometric measurements by direct inversion from the measurements to the harmonic geopotential coefficients, not by the way of lumped coefficients. For orbiters of Mars, we have suggestions which orbits should be avoided. The slow rotation of Venus results in dense ground-tracks and excellent gravitational recovery for almost all orbiters.     

Pressure-specific and multiple pressure response of fish assemblages in European running waters

We classified homogenous river types across Europe and searched for fish metrics qualified to show responses to specific pressures (hydromorphological pressures or water quality pressures) vs. multiple pressures in these river types. We analysed fish taxa lists from 3105 sites in 16 ecoregions and 14 countries. Sites were pre-classified for 15 selected pressures to separate unimpacted from impacted sites. Hierarchical cluster analysis was used to split unimpacted sites into four homogenous river types based on species composition and geographical location. Classification trees were employed to predict associated river types for impacted sites with four environmental variables. We defined a set of 129 candidate fish metrics to select the best reacting metrics for each river type. The candidate metrics represented tolerances/intolerances of species associated with six metric types: habitat, migration, water quality sensitivity, reproduction, trophic level and biodiversity. The results showed that 17 uncorrelated metrics reacted to pressures in the four river types. Metrics responded specifically to water quality pressures and hydromorphological pressures in three river types and to multiple pressures in all river types. Four metrics associated with water quality sensitivity showed a significant reaction in up to three river types, whereas 13 metrics were specific to individual river types. Our results contribute to the better understanding of fish assemblage response to human pressures at a pan-European scale. The results are especially important for European river management and restoration, as it is necessary to uncover underlying processes and effects of human pressures on aquatic communities.     

Fusion of multi‐resolution surface (terrestrial laser scanning) and subsurface geodata (ERT,SRT) for karst landform investigation and geomorphometric quantification

A multi‐method research design based on terrestrial laser scanning, GIS, geophysical prospecting (electrical resistivity tomography, refraction seismics) and sedimentology is applied for the first time to investigate enclosed karst depressions in an integrated way. Fusing multi‐resolution surface and subsurface geodata provides profound insights into the formation, geometry and geomorphologic processes of dolines. The studied landforms, which are located in the Dikti Mountains of East Crete, are shown to be filled by loose sediments of thicknesses of up to 30 m that mainly consist of fine‐grained material overlying solid bedrock at depths below 35 to 45 m. By combining subsurface observations with geomorphometric calculations, local doline genesis can be traced back to initial collapse of fractured bedrock followed by subsequent infilling with colluvials. In order to define crucial methodological requirements and guidelines for data fusion, both the impact of different elevation models and the influence of data resolution are assessed. Surface volumes of depressions derived by the digital surface model are 7–21% higher than the results obtained from the terrain model due to vegetation. Similarly, estimates of infill volume calculated on the basis of geophysical outcomes and elevation data differ by up to 13%. Calculations of the landforms' current volumes (i.e. total surface and subsurface volume), however, are fairly insensitive to raster resolution. Hence, the distinct geomorphologic properties of landforms (e.g. shape, terrain roughness, slope inclination) substantially determine the geomorphometric analysis of both surface and subsurface data. As shown by the findings, data fusion to integrate digital terrain, geophysical and sedimentological datasets of varied resolutions benefits geomorphologic studies and helps provide a comprehensive image of landforms. Copyright © 2013 John Wiley & Sons, Ltd.     

Modelling climate and land-use change impacts with SWIM: lessons learnt from multiple applications

Abstract

The Soil and Water Integrated Model (SWIM) is a continuous-time semi-distributed ecohydrological model, integrating hydrological processes, vegetation, nutrients and erosion. It was developed for impact assessment at the river basin scale. SWIM is coupled to GIS and has modest data requirements. During the last decade SWIM was extensively tested in mesoscale and large catchments for hydrological processes (discharge, groundwater), nutrients, extreme events (floods and low flows), crop yield and erosion. Several modules were developed further (wetlands and snow dynamics) or introduced (glaciers, reservoirs). After validation, SWIM can be applied for impact assessment. Four exemplary studies are presented here, and several questions important to the impact modelling community are discussed. For which processes and areas can the model be used? Where are the limits in model application? How to apply the model in data-poor situations or in ungauged basins? How to use the model in basins subject to strong anthropogenic pressure?

Editor D. Koutsoyiannis; Associate editor C. Perrin     

Estimation of the temporal autocorrelation structure by the collocation technique with an emphasis on soil moisture studies

Abstract

The collocation technique has become a popular tool in oceanography and hydrology for estimating the error variances of different data sources such as in situ sensors, models and remote sensing products. It is also possible to determine calibration constants, for example to account for an off-set between the data sources. So far, the temporal autocorrelation structure of the errors has not been studied, although it is known that it has detrimental effects on the results of the collocation technique, in particular when calibration constants are also determined. This paper shows how the (triple) collocation estimators can be adapted to retrieve the autocovariance functions; the statistical properties as well as the structural deficencies are described. The coupling between the autocorrelation of the error and the estimation of calibration constants is studied in detail, due to its importance for analysing temporal changes. In soil moisture applications, such time variations can be induced, for example, by seasonal changes in the vegetation cover, which affect both models and remote sensing products. The limitations of the proposed technique associated with these considerations are analysed using remote sensing and in situ soil moisture data. The variability of the inter-sensor calibration and the autocovariance are shown to be closely related to temporal patterns of the data.

Editor D. Koutsoyiannis

Citation Zwieback, S., Dorigo, W., and Wagner, W., 2013. Estimation of the temporal autocorrelation structure by the collocation technique with an emphasis on soil moisture studies. Hydrological Sciences Journal, 58 (8), 1729–1747.     

Lake Floor Morphology and Sediment Architecture of Lake Torneträsk,Northern Sweden

Here we present datasets from a hydroacoustic survey in July 2011 at Lake Torneträsk, northern Sweden. Our hydroacoustic data exhibit lake floor morphologies formed by glacial erosion and accumulation processes, insights into lacustrine sediment accumulation since the beginning of deglaciation, and information on seismic activity along the Pärvie Fault. Features of glacial scouring with a high‐energy relief, steep slopes, and relative reliefs of more than 50 m are observed in the large W‐basin. The remainder of the lacustrine subsurface appears to host a broad variety of well preserved formations from glacial accumulation related to the last retreat of the Fennoscandian ice sheet. Deposition of glaciolacustrine and lacustrine sediments is focused in areas situated in proximity to major inlets. Sediment accumulation in distal areas of the lake seldom exceeds 2 m or is not observable. We assume that lack of sediment deposition in the lake is a result of different factors, including low rates of erosion in the catchment, a previously high lake level leading to deposition of sediments in higher elevated paleodeltas, tributaries carrying low suspension loads as a result of sedimentation in upstream lakes, and an overall low productivity in the lake. A clear off‐shore trace of the Pärvie Fault could not be detected from our hydroacoustic data. However, an absence of sediment disturbance in close proximity to the presumed fault trace implies minimal seismic activity since deposition of the glaciolacustrine and lacustrine sediments.