Filtern
Erscheinungsjahr
- 2010 (90) (entfernen)
Dokumenttyp
- Wissenschaftlicher Artikel (66)
- Dissertation (21)
- Sonstiges (2)
- Monographie/Sammelband (1)
Schlagworte
- Arrayseismologie (2)
- Cosmogenic nuclides (2)
- Erdbeben (2)
- Erdbebenkatalog (2)
- Erdbebenschwarm 2008/09 (2)
- InSAR (2)
- Kosmogene Nuklide (2)
- Standorteffekte (2)
- Vogtland/West Bohemia (2)
- Vogtland/Westböhmen (2)
Institut
- Institut für Geowissenschaften (90) (entfernen)
Large-scale volcanic deformation recently detected by radar interferometry (InSAR) provides new information and thus new scientific challenges for understanding volcano-tectonic activity and magmatic systems. The destabilization of such a system at depth noticeably affects the surrounding environment through magma injection, ground displacement and volcanic eruptions. To determine the spatiotemporal evolution of the Lazufre volcanic area located in the central Andes, we combined short-term ground displacement acquired by InSAR with long-term geological observations. Ground displacement was first detected using InSAR in 1997. By 2008, this displacement affected 1800 km2 of the surface, an area comparable in size to the deformation observed at caldera systems. The original displacement was followed in 2000 by a second, small-scale, neighbouring deformation located on the Lastarria volcano. We performed a detailed analysis of the volcanic structures at Lazufre and found relationships with the volcano deformations observed with InSAR. We infer that these observations are both likely to be the surface expression of a long-lived magmatic system evolving at depth. It is not yet clear whether Lazufre may trigger larger unrest or volcanic eruptions; however, the second deformation detected at Lastarria and the clear increase of the large-scale deformation rate make this an area of particular interest for closer continuous monitoring.
African savannas are primarily used as pastures and are subject to changes in climate and management strategies. For sustainable management of these landscapes ecological knowledge on seasonal and long-term variability in plant community composition and the availability of green biomass is essential. In this study, we assessed the effects of dry and wet season on species richness and beta diversity for three sites along a gradient of increasing vegetation cover and precipitation in northwest Namibia. A hexagonal systematic sampling design was used to record floristic data. The Simple Matching, Soerensen, and multi-plot similarity coefficient and distance decay analyses were applied for examining beta diversity. Analyses were repeated while separating the plots according to the presence of woody vegetation. Species richness nearly doubled from dry to wet season; compositional similarity increased from dry to wet season and with increasing aridity of the study sites: distance decay was more pronounced in the dry season without any link to the precipitation gradient. Woody elements in the landscape, which occur along drainage lines or as tree islands, govern spatial and seasonal plant diversity fluctuations. Monitoring them is important for conservation strategies and for establishing grazing rules that ensure a sustainable use of savanna ecosystems.
The apparent isotope enrichment factor epsilon(macrophyte) of submerged plants (epsilon(macrophyte-DIC) = delta C-13(macrophyte) - delta C-13(DIC)) is indicative of dissolved inorganic carbon (DIC) supply in neutral to alkaline waters and is related to variations in aquatic productivity (Papadimitriou et al. in Limnol Oceanogr 50:1084-1095, 2005). This paper aims to evaluate the usage of epsilon(macrophyte) inferred from isotopic analyses of submerged plant fossils in addition to analyses of lake carbonate as a palaeolimnological proxy for former HCO3 (-) concentrations. Stable carbon isotopic analysis of modern Potamogeton pectinatus leaves and its host water DIC from the Tibetan Plateau and Central Yakutia (Russia) yielded values between -23.3 and +0.4aEuro degrees and between +14.0 and +6.5aEuro degrees, respectively. Values of epsilon (Potamogeton-DIC) (range -15.4 to +1.1aEuro degrees) from these lakes are significantly correlated with host water HCO3 (-) concentration (range 78-2,200 mg/l) (r = -0.86; P < 0.001), thus allowing for the development of a transfer function. Palaeo-epsilon (Potamogeton-ostracods) values from Luanhaizi Lake on the NE Tibetan Plateau, as inferred from the stable carbon isotope measurement of fossil Potamogeton pectinatus seeds (range -24 to +2.8aEuro degrees) and ostracods (range -7.8 to +7.5%) range between -14.8 and 1.6aEuro degrees. Phases of assumed disequilibrium between delta C-13(DIC) and delta C-13(ostracods) known to occur in charophyte swards (as indicated by the deposition of charophyte fossils) were excluded from the analysis of palaeo-epsilon. The application of the epsilon (Potamogeton-DIC)-HCO3 (-) transfer function yielded a median palaeo-HCO3 (-) -concentration of 290 mg/l. Variations in the dissolved organic carbon supply compare well with aquatic plant productivity changes and lake level variability as inferred from a multiproxy study of the same record including analyses of plant macrofossils, ostracods, carbonate and organic content.
Seismic wide-angle data were collected along a 40-km-long profile centered at the geothermal research well GrSk 3/90 in the Northeast German Basin. Tomographic inversion of travel time data provided a velocity and a vertical velocity gradient model, indicative of Cenozoic to Pre-Permian sediments. Wide-angle reflections are modeled and interpreted as top Zechstein and top Pre-Permian. Changes in velocity gradients are interpreted as the transition from mechanical to chemical compaction at 2-3 km depth, and localized salt structures are imaged, suggesting a previously unknown salt pillow in the southern part of the seismic profile. The Zechstein salt shows decreased velocities in the adjacent salt pillows compared to the salt lows, which is confirmed by sonic log data. This decrease in velocity could be explained by the mobilization of less dense salt, which moved and formed the salt pillows, whereas the denser salt remained in place at the salt lows. We interpret a narrow subvertical low-velocity zone under the salt pillow at GrSk 3/ 90 as a fault in the deep Permian to Pre-Permian. This WNW-ESE trending fault influenced the location of the salt tectonics and led to the formation of a fault-bounded graben in the Rotliegend sandstones with optimal mechanical conditions for geothermal production. Thermal modeling showed that salt pillows are related to chimney effects, a decrease in temperature, and increasing velocity. The assumed variations in salt lithology, density, and strain must thus be even higher to compensate for the temperature effect.
Three dimensional modelling of fractured and faulted reservoirs : framework and implementation
(2010)
Modelling of coupled physical processes in fractured and faulted media is a major challenge for the geoscience community. Due to the complexity related to the geometry of real fracture networks and fault systems, modelling studies have been mainly restricted either to two dimensional cases or to simplified orthogonal fracture systems consisting of vertical and horizontal fractures. An approach to generate three dimensional meshes for realistic fault geometries is presented. The method enables representation of faults in an arbitrary incline as two dimensional planes within a three dimensional, stratified porous matrix of a generic geometry. Based on a structural geological model, the method creates three dimensional unstructured tetrahedral meshes. These meshes can be used for finite element and finite volume numerical simulations. A simulation of a coupled fluid flow and heat transport problem for a two layered porous medium cut by two crossing faults is presented to test the reliability of the method.
Ecologists carry a well-stocked toolbox with a great variety of sampling methods, statistical analyses and modelling tools, and new methods are constantly appearing. Evaluation and optimisation of these methods is crucial to guide methodological choices. Simulating error-free data or taking high-quality data to qualify methods is common practice. Here, we emphasise the methodology of the 'virtual ecologist' (VE) approach where simulated data and observer models are used to mimic real species and how they are 'virtually' observed. This virtual data is then subjected to statistical analyses and modelling, and the results are evaluated against the 'true' simulated data. The VE approach is an intuitive and powerful evaluation framework that allows a quality assessment of sampling protocols, analyses and modelling tools. It works under controlled conditions as well as under consideration of confounding factors such as animal movement and biased observer behaviour. In this review, we promote the approach as a rigorous research tool, and demonstrate its capabilities and practical relevance. We explore past uses of VE in different ecological research fields, where it mainly has been used to test and improve sampling regimes as well as for testing and comparing models, for example species distribution models. We discuss its benefits as well as potential limitations, and provide some practical considerations for designing VE studies. Finally, research fields are identified for which the approach could be useful in the future. We conclude that VE could foster the integration of theoretical and empirical work and stimulate work that goes far beyond sampling methods, leading to new questions, theories, and better mechanistic understanding of ecological systems.
Global climate and the atmospheric partial pressure of carbon dioxide (p(CO2atm)) are correlated over recent glacial cycles, with lower p(CO2atm) during ice ages, but the causes of the p(CO2atm) changes are unknown. The modern Southern Ocean releases deeply sequestered CO2 to the atmosphere. Growing evidence suggests that the Southern Ocean CO2 'leak' was stemmed during ice ages, increasing ocean CO2 storage. Such a change would also have made the global ocean more alkaline, driving additional ocean CO2 uptake. This explanation for lower ice-age p(CO2atm), if correct, has much to teach us about the controls on current ocean processes.
Amphibole and mica Ar-40/Ar-39 ages as well as zircon, rutile and titanite U-Pb geochronology of eclogites and associated host rocks from the Higher Himalayan Crystalline Nappes (Indian Plate) in the Upper Kaghan Valley, Pakistan allow distinction of a multistage exhumation history. An Eocene age for peak-pressure metamorphism has been obtained by phengite Ar-40/Ar-39 (47.3 +/- 0.3 Ma) and zircon U-Pb (47.3 +/- 0.4 and 47.4 +/- 0.3 Ma) ages from cover and basement gneisses. A very short-lived metamorphic peak and rapid cooling is documented by an amphibole Ar-40/Ar-39 age of 46.6 +/- 0.5 Ma and a rutile U-Pb age of 44.1 +/- 1.3 Ma from eclogites. Phengite and biotite ages from cover and basement sequences metamorphosed during the Himalayan orogeny are 34.5 +/- 0.2 to 28.1 +/- 0.2 Ma whereas youngest biotites, yielding 23.6 +/- 0.1 and 21.7 +/- 0.2 Ma, probably reflect argon partial resetting. The amphibole age, together with those derived from phengite and zircon demonstrate a rate of initial exhumation of 86-143 mm/a i.e. an extremely rapid transport of the Indian Plate continental crust from ultra-high pressure (UHP) conditions back to crustal levels (47-46 Ma for transport from 140 to 40 km depth). Subsequent exhumation (46-41 Ma, 40-35 km) slowed to about 1 mm/a at the base of the continental crust but increased again later towards slightly higher exhumation rates of ca. 2 mm/a (41-34 Ma, 35- 20 km). This indicates a change from buoyancy-driven exhumation at mantle depths to compression forces related to continent-continent collision and accompanied crustal folding, thrusting and stacking that finally exposed the former deeply-buried rocks.
Temporal gravimeter observations, used in geodesy and geophysics to study variation of the Earth’s gravity field, are influenced by local water storage changes (WSC) and – from this perspective – add noise to the gravimeter signal records. At the same time, the part of the gravity signal caused by WSC may provide substantial information for hydrologists. Water storages are the fundamental state variable of hydrological systems, but comprehensive data on total WSC are practically inaccessible and their quantification is associated with a high level of uncertainty at the field scale. This study investigates the relationship between temporal gravity measurements and WSC in order to reduce the hydrological interfering signal from temporal gravity measurements and to explore the value of temporal gravity measurements for hydrology for the superconducting gravimeter (SG) of the Geodetic Observatory Wettzell, Germany. A 4D forward model with a spatially nested discretization domain was developed to simulate and calculate the local hydrological effect on the temporal gravity observations. An intensive measurement system was installed at the Geodetic Observatory Wettzell and WSC were measured in all relevant storage components, namely groundwater, saprolite, soil, top soil and snow storage. The monitoring system comprised also a suction-controlled, weighable, monolith-filled lysimeter, allowing an all time first comparison of a lysimeter and a gravimeter. Lysimeter data were used to estimate WSC at the field scale in combination with complementary observations and a hydrological 1D model. Total local WSC were derived, uncertainties were assessed and the hydrological gravity response was calculated from the WSC. A simple conceptual hydrological model was calibrated and evaluated against records of a superconducting gravimeter, soil moisture and groundwater time series. The model was evaluated by a split sample test and validated against independently estimated WSC from the lysimeter-based approach. A simulation of the hydrological gravity effect showed that WSC of one meter height along the topography caused a gravity response of 52 µGal, whereas, generally in geodesy, on flat terrain, the same water mass variation causes a gravity change of only 42 µGal (Bouguer approximation). The radius of influence of local water storage variations can be limited to 1000 m and 50 % to 80 % of the local hydro¬logical gravity signal is generated within a radius of 50 m around the gravimeter. At the Geodetic Observatory Wettzell, WSC in the snow pack, top soil, unsaturated saprolite and fractured aquifer are all important terms of the local water budget. With the exception of snow, all storage components have gravity responses of the same order of magnitude and are therefore relevant for gravity observations. The comparison of the total hydrological gravity response to the gravity residuals obtained from the SG, showed similarities in both short-term and seasonal dynamics. However, the results demonstrated the limitations of estimating total local WSC using hydrological point measurements. The results of the lysimeter-based approach showed that gravity residuals are caused to a larger extent by local WSC than previously estimated. A comparison of the results with other methods used in the past to correct temporal gravity observations for the local hydrological influence showed that the lysimeter measurements improved the independent estimation of WSC significantly and thus provided a better way of estimating the local hydrological gravity effect. In the context of hydrological noise reduction, at sites where temporal gravity observations are used for geophysical studies beyond local hydrology, the installation of a lysimeter in combination with complementary hydrological measurements is recommended. From the hydrological view point, using gravimeter data as a calibration constraint improved the model results in comparison to hydrological point measurements. Thanks to their capacity to integrate over different storage components and a larger area, gravimeters provide generalized information on total WSC at the field scale. Due to their integrative nature, gravity data must be interpreted with great care in hydrological studies. However, gravimeters can serve as a novel measurement instrument for hydrology and the application of gravimeters especially designed to study open research questions in hydrology is recommended.
Cities often have higher species diversity than the surrounding landscape. This diversity is important for both nature conservation and urban planning. The recreation of residents and the protection of species and habitats are simultaneous targets of maintaining urban green spaces. Data about the distribution and richness of species and their habitats have been compiled frequently; however, it is difficult and costly to measure the complete biodiversity of a region, necessitating useful surrogates. We tested species and habitat data in 27 protected areas in a Central German city and asked (1) whether the diversity of selected taxa acts as a surrogate for the diversity of other taxa and total investigated diversity, and (2) whether landscape structure and human impact explain species richness. Landscape structure metrics were based on soil and habitat types; human influence was measured as the degree of hemeroby. We tested and accounted for sample bias prior to analyses. (1) Vascular plant species richness explained total richness and single taxon richness best. (2) The size of a protected area was the most important predictor of species richness. After correcting for the effect of size, shape complexity, isolation, and matrix properties remained significant. Accordingly, the type of data frequently used for urban planning - collected over several years, by various persons for various purposes - is suitable regarding systematic conservation planning for species richness. The surrogate taxa concept applies in urban areas but with restrictions. Additionally, species richness should be examined in the context of both the city and its surrounding countryside.
The hummocky post-glacial soil landscapes with kettle holes as internal drainage systems are characterized by ponds that trap lateral fluxes in topographic depressions. A quantitative description is mostly limited by the unknown complexity of hydraulically relevant soil and sediment structures. This paper is focussing on a structure-based approach to identify relevant field-scale flow and transport processes. Illustrative examples demonstrate extreme variations in water table fluctuation for adjoining kettle holes. Explanations require a pedohydrologic concept of the arable soil landscape. Identification of structures is based on geophysical methods and soil hydraulic measurements. Electrical resistivity imaging yields 0.5 m-scale spatial structures that correspond with soil texture distributions. Electromagnetic induction provides larger-scale field maps that reflect major soil and sediment features. Results of both methods correspond within the limits of the different spatial resolutions. With geophysical exploration methods, colluvial areas with textural differences between upper and deeper soil layers, coarse-textured sediment lenses, and stony colluvial regions around kettle holes are identified as potentially relevant flow structures. The colluvial fringe around the pond seems to be a sensitive area with important lateral exchange fluxes. Tensiometer measurements perpendicular to this boundary indicate hydraulic gradients directed from the pond towards the partially saturated soil. The localized infiltration of trapped water in kettle holes can control large fractions of ground water recharge and may have implications for the fate of agricultural chemicals in post-glacial landscapes. While surface and subsurface hydraulic structures may be inferred using minimal-invasive techniques, better understanding of processes and properties governing lateral exchange fluxes between pond and surrounding soil are required.
To calibrate delta O-18 time-series from speleothems in the eastern Indian summer monsoon (ISM) region of India, and to understand the moisture regime over the northern Bay of Bengal (BoB) we analyze the delta O-18 and delta D of rainwater, collected in 2007 and 2008 near Cherrapunji, India. delta D values range from + 18.5 parts per thousand to 144.4 parts per thousand, while delta O-18 varies between +0.8 parts per thousand and 18.8 parts per thousand. The Local Meteoric Water Line (LMWL) is found to be indistinguishable from the Global Meteoric Water Line (GMWL). Late ISM (September-October) rainfall exhibits lowest delta O-18 and delta D values, with little relationship to the local precipitation amount. There is a trend to lighter isotope values over the course of the ISM, but it does not correlate with the patterns of temperature and rainfall amount delta O-18 and delta D time-series have to be interpreted with caution in terms of the 'amount effect' in this subtropical region. We find that the temporal trend in delta O-18 reflects increasing transport distance during the ISM, isotopic changes in the northern BoB surface waters during late ISM, and vapor re-equilibration with rain droplets. Using an isotope box model for surface ocean waters, we quantify the potential influence of river runoff on the isotopic composition of the seasonal freshwater plume in the northern BoB. Temporal variations in this source can contribute up to 25% of the observed changes in stable isotopes of precipitation in NE India. To delineate other moisture sources, we use backward trajectory computations and find a strong correlation between source region and isotopic composition. Palaeoclimatic stable isotope time-series from northeast Indian speleothems likely reflect changes in moisture source and transport pathway, as well as the isotopic composition of the BoB surface water, all of which in turn reflect ISM strength. Stalagmite records from the region can therefore be interpreted as integrated measures of the ISM strength.
The Lombok Island is part of the Lesser Sunda Islands (LSI) region – Indonesia, situated along the Sunda-Banda Arcs transition. It lies between zones characterized by the highest intensity geomagnetic anomalies of this region, remarkable as one of the eight most important features provided on the 1st edition of World Digital Magnetic Anomaly Map. The seismicity of this region during the last years is high, while the geological and tectonic structures of this region are still not known in detail. Some local magnetic surveys have been conducted previously during 2004–2005. However, due to the lower accuracy of the used equipment and a limited number of stations, the qualities of the previous measurements are questionable for more interpretations. Thus a more detailed study to better characterize the geomagnetic anomaly -spatially and temporally- over this region and to deeply explore the related regional geology, tectonic and seismicity is needed. The intriguing geomagnetic anomalies over this island region vis-à-vis the socio-cultural situations lead to a study with a special aim to contribute to the assessment of the potential of natural hazards (earthquakes) as well as a new natural resource of energy (geothermal potential).
This study is intended to discuss several crucial questions, including:
i. The real values and the general pattern of magnetic anomalies over the island, as well as their relation to the regional one.
ii. Any temporal changes of regional anomalies over the recent time.
iii. The relationships between the anomalies and the geology and tectonic of this region, especially new insights that can be gained from the geomagnetic observations.
iv. The relationships between the anomalies and the high seismicity of this region, especially some possible links between their variations to the earthquake occurrence.
First, all available geomagnetic data of this region and results of the previous measurements are evaluated. The new geomagnetic surveys carried out in 2006 and 2007/2008 are then presented in detail, followed by the general description of data processing and data quality evaluation. The new results show the general pattern of contiguous negative-positive anomalies, revealing an active arc related subduction region. They agree with earlier results obtained by satellite, aeromagnetic, and marine platforms; and provide a much more detailed picture of the strong anomalies on this island. The temporal characteristics of regional anomalies show a decreasing strength of the dipolar structure, where decreasing of the field intensities is faster than the regional secular variations as defined by the global model (the 10th generation of IGRF). However, some exceptions (increasing of anomalies) have to be noted and further analyzed for several locations.
Thereafter, simultaneous magnetic anomalies and gravity models are generated and interpreted in detail. Three profiles are investigated, providing new insights into the tectonics and geological evolution of the Lombok Island. Geological structure of this island can be divided as two main parts with different consecutive ages: an old part (from late Oligocene to late Miocene) in the South and a younger one (from Pliocene to Holocene) in the North. A new subduction in the back arc region (the Flores Thrust zone) is considered mature and active, showing a tendency of progressive subduction during 2005–2008. Geothermal potential in the northern part of this island can be mapped in more detail using these geomagnetic regional survey data. The earlier estimates of reservoir depth can be confirmed further to a depth of about 800 m. Evaluation of temporal changes of the anomalies gives some possible explanations related to the evolution of the back arc region, large stress accumulations over the LSI region, a specific electrical characteristic of the crust of the Lombok Island region, and a structural discontinuity over this island.
Based on the results, several possible advanced studies involving geomagnetic data and anomaly investigations over the Lombok Island region can be suggested for the future:
i. Monitoring the subduction activity of the back arc region (the Flores Thrust zone) and the accumulated stress over the LSI, that could contribute to middle term hazard assessment with a special attention to the earthquake occurrence in this region. Continuous geomagnetic field measurements from a geomagnetic observatory which can be established in the northern part of the Lombok Island and systematic measurements at several repeat stations can be useful in this regards.
ii. Investigating the specific electrical characteristic (high conductivity) of the crust, that is probably related to some aquifer layers or metal mineralization. It needs other complementary geophysical methods, such as magnetotelluric (MT) or preferably DC resistivity measurements.
iii. Determining the existence of an active structural fault over the Lombok Island, that could be related to long term hazard assessment over the LSI region. This needs an extension of geomagnetic investigations over the neighbouring islands (the Bali Island in the West and the Sumbawa Island in the East; probably also the Sumba and the Flores islands). This seems possible because the regional magnetic lineations might be used to delineate some structural discontinuities, based on the modelling of contrasts in crustal magnetizations.
Although the effects of grazing-induced savannah degradation on animal diversity are well documented, knowledge of how they affect space use or responding behaviour remains poor. In this study, we analysed space use of the spotted sand lizard (Pedioplanis l. lineoocellata) in degraded versus nondegraded habitats of southern Kalahari savannah habitats. Lizards were radio tracked, daily movement distances recorded and home range sizes calculated. In degraded Kalahari savannah habitats where plant diversity and perennial grass cover are low but shrub cover high, P. lineoocellata moves larger distances (40.88 +/- 6.42 m versus 27.43 +/- 5.08 m) and occupies larger home ranges (646.64 +/- 244.84 m(2) versus 209.15 +/- 109.84 m(2)) than in nondegraded habitats (high plant diversity, high perennial grass cover and low shrub cover). We assume that this increase in daily movement distances and home range sizes is a behavioural plasticity to limited food resources in degraded savannah habitats. Although P. lineoocellata is able to adjust to resource-poor savannah habitats, the increase in the lizard's movement activities is likely to result in a higher predation risk. This is supported by the lower availability of protective vegetation i.e. perennial grass cover. Hence, we conclude that despite behavioural plasticity of P. lineoocellata, overgrazing has a severe negative impact on the space use of P. lineoocellata.
Simulation of spatial sensor characteristics in the context of the EnMAP Hyperspectral mission
(2010)
The simulation of remote sensing images is a valuable tool for defining future Earth observation systems, optimizing instrument parameters, and developing and validating data-processing algorithms. A scene simulator for optical Earth observation data has been developed within the Environmental Mapping and Analysis Program (EnMAP) hyperspectral mission. It produces EnMAP-like data following a sequential processing approach consisting of five independent modules referred to as reflectance, atmospheric, spatial, spectral, and radiometric modules. From a modeling viewpoint, the spatial module is the most complex. The spatial simulation process considers the satellite-target geometry, which is adapted to the EnMAP orbit and operating characteristics, the instrument spatial response, and the sources of spatial nonuniformity (keystone, telescope distortion and smile, and detector coregistration). The spatial module of the EnMAP scene simulator is presented in this paper. The EnMAP spatial and geometric characteristics will be described, the simulation methodology will be presented in detail, and the capability of the EnMAP simulator will be shown by illustrative examples.
The hydrologic cycle of high mountainous catchments is frequently simulated with simple precipitation-discharge models representing the snow accumulation and ablation behavior of a very complex environment with a set of lumped equations accounting for altitudinal temperature and precipitation gradients. In this study, we present a methodology to include sparse snow depths measurements into the calibration process. Based on this methodology, we assess for a case study, the Rhonegletscher catchment (Switzerland), how much observed information we need to reliably calibrate the model, such that it reproduces the dominant system dynamics, discharge, as well as glacier mass balance. Here, we focus on the question whether observed discharge is sufficient as a calibration variable or whether we need annual or even seasonal glacier mass balance data. Introducing seasonally variable accumulation and ablation parameters is sufficient to enable the simple model to reproduce observed seasonal mass balances for the Rhonegletscher. Furthermore, our results suggest that calibrating the hydrological model exclusively on discharge can lead to wrong representations of the intra- annual accumulation and ablation processes and to a strong bias in long term glacier mass balance simulations. Adding only a few annual mass balance observations considerably reduces this bias. Calibrating exclusively on annual balance data can, in turn, lead to wrong seasonal mass balance simulations. Even if these results are case study specific, our conclusions provide valuable new insights into the benefit of different types of observations for calibrating hydrological models in glacier catchments. The presented multi-signal calibration framework and the simple method to calibrate a semi-lumped model on point observations has potential for application in other modeling contexts.
In the present study, we investigated the dispersion characteristics of medium-to-long period Rayleigh waves (2 s < T < 20 s) using both single-station techniques (multiple-filter analysis, and phase-match filter) and multichannel techniques (horizontal slowness [p] and angular frequency [omega] stack, and cross-correlation) to determine the velocity structure for the Mt. Etna volcano. We applied these techniques to a dataset of teleseisms, as regional and local earthquakes recorded by two broad-band seismic arrays installed at Mt. Etna in 2002 and 2005, during two seismic surveys organized by the Istituto Nazionale di Geofisica e Vulcanologia (INGV), sezione di Napoli. The dispersion curves obtained showed phase velocities ranging from 1.5 km/s to 4.0 km/s in the frequency band 0.05 Hz to 0.45 Hz. We inverted the average phase velocity dispersion curves using a non-linear approach, to obtain a set of shear-wave velocity models with maximum resolution depths of 25 km to 30 km. Moreover, the presence of lateral velocity contrasts was checked by dividing the whole array into seven triangular sub-arrays and inverting the dispersion curves relative to each triangle.
Applying metrics for hydrograph comparison is a central task in hydrological modelling, used both in model calibration and the evaluation of simulations or forecasts. Motivated by the shortcomings of standard objective metrics such as the Root Mean Square Error or the Mean Peak Time Error and the advantages of visual inspection as a powerful tool for simultaneous, case-specific and multi-criteria (yet subjective) evaluation, we propose a new objective metric termed Series Distance, which is in close accordance with visual evaluation. The Series Distance is an event-based method and consists of three parts, namely a Threat Score which evaluates overall agreement of event occurrence, and the overall distance of matching observed and simulated events with respect to amplitude and timing. The novelty of the latter two is the way in which matching point pairs on the observed and simulated hydrographs are identified, namely by the same relative position in matching segments (rise or recession) of matching events. Thus, amplitude and timing errors are calculated simultaneously but separately, from point pairs that also match visually, considering complete events rather than only individual points (which is for example the case with metrics related to Peak Time Errors). After presenting the Series Distance theory, we discuss its properties and compare it to those of standard metrics and visual inspection, both at the example of simple, artificial hydrographs and an ensemble of realistic forecasts. The results suggest that the Series Distance compares and evaluates hydrographs in a way comparable to visual inspection, but in an objective, reproducible way.