Quantitative evaluation of different hydrological modelling approaches in a partly glacierized Swiss watershed

Mountain water resources management often requires hydrological models that need to handle both snow and ice melt. In this study, we compared two different model types for a partly glacierized watershed in central Switzerland: (1) an energy‐balance model primarily designed for snow simulations; and (2) a temperature‐index model developed for glacier simulations. The models were forced with data extrapolated from long‐term measurement records to mimic the typical input data situation for climate change assessments. By using different methods to distribute precipitation, we also assessed how various snow cover patterns influenced the modelled runoff. The energy‐balance model provided accurate discharge estimations during periods dominated by snow melt, but dropped in performance during the glacier ablation season. The glacier melt rates were sensitive to the modelled snow cover patterns and to the parameterization of turbulent heat fluxes. In contrast, the temperature‐index model poorly reproduced snow melt runoff, but provided accurate discharge estimations during the periods dominated by glacier ablation, almost independently of the method used to distribute precipitation. Apparently, the calibration of this model compensated for the inaccurate precipitation input with biased parameters. Our results show that accurate estimates of snow cover patterns are needed either to correctly constrain the melt parameters of the temperature‐index model or to ensure appropriate glacier surface albedos required by the energy‐balance model. Thus, particularly when only distant meteorological stations are available, carefully selected input data and efficient extrapolation methods of meteorological variables improve the reliability of runoff simulations in high alpine watersheds. Copyright © 2011 John Wiley & Sons, Ltd.     

Stereotomography assisted by migration of attributes

Depth velocity model building remains a difficult step within the seismic depth imaging sequence. Stereotomography provides an efficient solution to this problem but was limited until now to a picking of seismic data in the prestack time un-migrated domain. We propose here a method for stereotomographic data picking in the depth migrated domain. Picking in the depth migrated domain exhibits the advantage of a better signal-to-noise ratio and of a more regular distribution of picked events in the model, leading to a better constrained tomographic inverse problem. Moreover, any improvement on the velocity model will improve the migrated results, again leading to improved picking. Our strategy for obtaining a stereotomographic dataset from a prestack depth migration is based on migration of attributes (and not on a kinematic demigration approach!). For any locally coherent event in the migrated image, migration of attributes allows one to compute ray parameter attributes corresponding to the specular reflection angle and dip. For application to stereotomography, the necessary attributes are the source/receiver locations, the traveltime and the data slopes. For the data slope, when the migration velocity model is erroneous, some additional corrections have to be applied to the result of migration of the attributes. Applying these corrections, our picking method is theoretically valid whatever the quality of the migration velocity model. We first present the theoretical aspects of the method and then validate it on 2D synthetic and real seismic reflection data sets.     

Exploring debris‐flow history and process dynamics using an integrative approach on a dolomitic cone in western Austria

The evolution of a debris‐flow cone depends on a multitude of factors in the hydrogeomorphic system. Investigations of debris‐flow history and cone dynamics in highly active catchments therefore require an integrative approach with a temporal and spatial resolution appropriate for the goals of the study. We present the use of an orthophoto time series to augment standard dendrogeomorphic techniques to describe the spatio‐temporal dynamics of debris flows on a highly active cone in the western Austrian Alps. Analysis of seven orthophotos since 1951 revealed a migration of active deposition areas with a resulting severe loss of forest cover (> 80%) and a mean tree loss per year of 10·4 (range 1·3–16·6 trees per year). Analysis of 193 Pinus mugo ssp. uncinata trees allowed the identification of 161 growth disturbances corresponding to 16 debris flows since 1839 and an average decadal frequency of 0·9 events. As a result of the severe loss of forest cover, we speculate that < 20% of the more recent events were actually captured in the tree‐ring record, giving a decadal return interval of ~7·5 events for a period of 60 years. Based on three annual field observations, it is evident that this catchment (the Bärenrüfe) produces very frequent (< 1 yr), small (in the order of a few 10 to 100 m³) debris flows with minor material relocation. The specific challenges of tree‐ring analysis in this tree species and in highly active environments are explicitly addressed in the discussion and underline the necessity of employing complementary methods of analysis in an integrative manner. Copyright © 2012 John Wiley & Sons, Ltd.     

Seasonal Water Storage Variations as Impacted by Water Abstractions: Comparing the Output of a Global Hydrological Model with GRACE and GPS Observations

Better quantification of continental water storage variations is expected to improve our understanding of water flows, including evapotranspiration, runoff and river discharge as well as human water abstractions. For the first time, total water storage (TWS) on the land area of the globe as computed by the global water model WaterGAP (Water Global Assessment and Prognosis) was compared to both gravity recovery and climate experiment (GRACE) and global positioning system (GPS) observations. The GRACE satellites sense the effect of TWS on the dynamic gravity field of the Earth. GPS reference points are displaced due to crustal deformation caused by time-varying TWS. Unfortunately, the worldwide coverage of the GPS tracking network is irregular, while GRACE provides global coverage albeit with low spatial resolution. Detrended TWS time series were analyzed by determining scaling factors for mean annual amplitude (f _GRACE) and time series of monthly TWS (f _GPS). Both GRACE and GPS indicate that WaterGAP underestimates seasonal variations of TWS on most of the land area of the globe. In addition, seasonal maximum TWS occurs 1 month earlier according to WaterGAP than according to GRACE on most land areas. While WaterGAP TWS is sensitive to the applied climate input data, none of the two data sets result in a clearly better fit to the observations. Due to the low number of GPS sites, GPS observations are less useful for validating global hydrological models than GRACE observations, but they serve to support the validity of GRACE TWS as observational target for hydrological modeling. For unknown reasons, WaterGAP appears to fit better to GPS than to GRACE. Both GPS and GRACE data, however, are rather uncertain due to a number of reasons, in particular in dry regions. It is not possible to benefit from either GPS or GRACE observations to monitor and quantify human water abstractions if only detrended (seasonal) TWS variations are considered. Regarding GRACE, this is mainly caused by the attenuation of the TWS differences between water abstraction variants due to the filtering required for GRACE TWS. Regarding GPS, station density is too low. Only if water abstractions lead to long-term changes in TWS by depletion or restoration of water storage in groundwater or large surface water bodies, GRACE may be used to support the quantification of human water abstractions.     

The critical period of high mortality of larvae fish —A discussion based on current research

Most fishes suffer high mortality during early development stages. The existence of a “critical period” was hypothesized by Hjort (1914). However, several decades of research have resulted in considerable controversy over the hypothesis. This paper attempts to clarify the critical period concept, demonstrate its usefulness, and resolve some of the controversy surrounding it. Based on our own studies and the literature, we find that within the “aquaculture school” (those concerned with culturing fishes, particulary larval stages), there is a great deal of evidence for the critical period concept. But we redefine the concept of “critical period” of mortality of fish larvae in this way. The deepest, sharpest, and fastest physiological and ecological changes during the development of the larvae are produced with the shift from endogenous nutrition to exogenous. At this stage, the old functions are quickly being replaced by new functions inside the body. Because of this shifting of function and imperfection of organs, the larvae are most sensitive at this time to environmental factors, particularly food supply. A higher specific mortality rate often occurs immediately folowing the period of strictly endogenous yolk feeding, and during the period of first exogenous feeding. Then it may be followed by a lower specific mortality rate. The significance of a “critical period” for fish populations is discussed. A good understanding of the “critical period” would provide a means for predicting the strength of year classes in natural fish populations, and also suggests approaches for increasing survival in aquaculture systems.     

Stable isotope study of cave percolation waters in subtropical Brazil: Implications for paleoclimate inferences from speleothems

We analyze the interannual monthly variability of oxygen isotope ratios in data from IAEA stations along the Atlantic coast of South America between 23° and 34° S to evaluate the influence of parameters such as temperature, rainfall amount and moisture source contribution on meteoric water recharging two karst systems in subtropical Brazil. In addition, a 2 year monitoring program performed on soil and cave drip and rimstone pool waters from sampling sites with contrasting discharge values and located at 100 and 300 m below the surface in the Santana Cave System (24°31′ S; 48°43′ W), is used to test the influence of hydrologic and geologic features on the temporal variations of seepage water δ¹⁸O.

Interannual monthly variations in δ¹⁸O of rainfall reflect primarily regional changes in moisture source contribution related to seasonal shifts in atmospheric circulation from a more monsoonal regime in summer (negative values of δ¹⁸O) to a more extratropical regime in winter (positive values of δ¹⁸O). Variations in groundwater δ¹⁸O indicate that the climatic signal of recent rainfall events is rapidly transmitted through the relatively deep karst aquifer to the cave drip waters, regardless of location of collection in the cave. In addition, the data also suggest that water replenishment in the system is triggered by the increase in hydraulic head during periods when recharge exceeds the storage capacity of the soil and epikarst reservoirs. Significant perturbations in the groundwater composition, characterized by more positive values of δ¹⁸O, are probably connected to an increased Atlantic moisture contribution associated with extratropical precipitation. This implies that the δ¹⁸O of speleothems from caves in this region may be a suitable proxy for studying tropical–extratropical interactions over South America, a feature that is intrinsically related to the global atmospheric circulation.     

The significance of flow in the matrix of the Chalk unsaturated zone

The significance of flow in the matrix of the Chalk unsaturated zone, in comparison with flow in fractures, has been the subject of much debate. In this article, important elements of the literature are discussed in detail and several simple modelling analyses based on steady-state flow are presented. A study of the sensitivity of solute spreading to fracture spacing in models that ignore matrix flow shows that this latter assumption is generally incompatible with observed solute profiles, unless unrealistically small fracture spacings are assumed. The effect of air phase continuities (e.g. bedding planes) on matrix flow has also been examined. These discontinuities are frequently interrupted by points of connectivity between matrix blocks. An issue therefore is the relationship between connectivity and its effect on inter-block conductance. A simple analysis of the Laplace equation shows that just 1% connectivity represents an effective pathway equivalent to 18% of the local rock hydraulic conductivity. Obviously, when there is no fracture flow, solute spreading is significantly reduced. However, dual permeability model simulations show that matrix flow reduces solute spreading in the presence of persistent fracture flow. All of the above studies suggest that flow in the matrix of the Chalk unsaturated zone is significant and that ignoring it may result in a serious misunderstanding of the system.     

Kyanite eclogites from the Rietfontein kimberlite pipe,mier coloured reserve,Gordonia, Cape Province,South Africa

The kimberlite pipe at Rietfontein has been found to contain a suite of ultrabasic xenoliths, amongst which are kyanite eclogites of unusual freshness. The kyanite crystals are of varying shades of blue. Physical properties of some of the primary minerals together with the chemical composition of four kyanite eclogites and of a light and a dark blue kyanite crystal are presented.The xenoliths found in the Rietfontein Pipe are compared with those found in the Zagadochnaya Pipe and certain similarities and differences noted.     

Homogeneous reprocessing of GPS,GLONASS and SLR observations

The International GNSS Service (IGS) provides operational products for the GPS and GLONASS constellation. Homogeneously processed time series of parameters from the IGS are only available for GPS. Reprocessed GLONASS series are provided only by individual Analysis Centers (i. e. CODE and ESA), making it difficult to fully include the GLONASS system into a rigorous GNSS analysis. In view of the increasing number of active GLONASS satellites and a steadily growing number of GPS+GLONASS-tracking stations available over the past few years, Technische Universität Dresden, Technische Universität München, Universität Bern and Eidgenössische Technische Hochschule Zürich performed a combined reprocessing of GPS and GLONASS observations. Also, SLR observations to GPS and GLONASS are included in this reprocessing effort. Here, we show only SLR results from a GNSS orbit validation. In total, 18 years of data (1994–2011) have been processed from altogether 340 GNSS and 70 SLR stations. The use of GLONASS observations in addition to GPS has no impact on the estimated linear terrestrial reference frame parameters. However, daily station positions show an RMS reduction of 0.3 mm on average for the height component when additional GLONASS observations can be used for the time series determination. Analyzing satellite orbit overlaps, the rigorous combination of GPS and GLONASS neither improves nor degrades the GPS orbit precision. For GLONASS, however, the quality of the microwave-derived GLONASS orbits improves due to the combination. These findings are confirmed using independent SLR observations for a GNSS orbit validation. In comparison to previous studies, mean SLR biases for satellites GPS-35 and GPS-36 could be reduced in magnitude from \(-35\) and \(-38\)  mm to \(-12\) and \(-13\)  mm, respectively. Our results show that remaining SLR biases depend on the satellite type and the use of coated or uncoated retro-reflectors. For Earth rotation parameters, the increasing number of GLONASS satellites and tracking stations over the past few years leads to differences between GPS-only and GPS+GLONASS combined solutions which are most pronounced in the pole rate estimates with maximum 0.2 mas/day in magnitude. At the same time, the difference between GLONASS-only and combined solutions decreases. Derived GNSS orbits are used to estimate combined GPS+GLONASS satellite clocks, with first results presented in this paper. Phase observation residuals from a precise point positioning are at the level of 2 mm and particularly reveal poorly modeled yaw maneuver periods.     

Time‐geographic Derivation of Feasible Co‐presence Opportunities from Network‐constrained Episodic Movement Data

Certain datasets on moving objects are episodic in nature – that is, the data is characterized by time gaps during which the position of the object is unknown. In this article, a model is developed to study the sparsely sampled network‐constrained movement of several objects by calculating both potential and feasible (i.e. more likely) co‐presence opportunities over time. The approach is applied to the context of a static sensor network, where the location of an object is only registered when passing a sensor location along a road network. Feasibility is incorporated based on the deviation from the shortest path. As an illustration, the model is applied to a large Bluetooth tracking dataset gathered at a mass event. The model output consists of maps showing the temporal evolution of the distribution of feasible co‐presence opportunities of tracked visitors over the network (i.e. the number of visitors that could have been present together). We demonstrate the model's usefulness in studying the movement and distribution of a crowd over a study area with relatively few sampling locations. Finally, we discuss the results with a special emphasis on the distinction between feasible and actual presence, the need for further validation and calibration, and the performance of the implementation.