Refine
Has Fulltext
- no (72) (remove)
Year of publication
- 2010 (72) (remove)
Document Type
- Article (66)
- Doctoral Thesis (5)
- Monograph/Edited Volume (1)
Is part of the Bibliography
- yes (72)
Keywords
- Lombok (1)
- geomagnetic (1)
- seismo-electromagnetic (1)
- tectonic (1)
Institute
- Institut für Geowissenschaften (72) (remove)
Barium, lanthanum, ytterbium, and yttrium partitioning experiments between fluid-saturated haplogranitic melts and aqueous solutions were conducted at 750 to 950 degrees C and 0.2 to 1 GPa to investigate the effects of melt and fluid composition, pressure, and temperature. Partition coefficients were determined using different experimental methods. On one hand quenched experiments were performed, and on the other hand, trace element contents in the aqueous fluid were determined directly using a hydrothermal diamond-anvil cell and synchrotron radiation X-ray fluorescence microanalysis of K-lines. The latter required a high excitation energy of 50 key due to the high energies necessary to excite the K-lines of the studied elements. The data from these two techniques showed good agreement for chloridic solutions, whereas quenching had a significant effect on results of the experiments with only water in the case of Ba. In Cl-free experiments, lanthanum and yttrium, trace element contents were even below detection limit in the quenched fluids, whereas small concentrations were detected in comparable in-situ experiments. This distinct difference is likely due to back reactions between fluid and melt upon cooling. The partitioning data of all elements show no dependence on the temperature and only small dependence on pressure. In contrast, the partitioning is strongly influenced by the composition of the starting fluid and melt. For chloridic fluids, there was a sharp increase in the Ba, La, Y and Yb partition coefficients with the alumina saturation index (ASI). The Ba partition coefficient increased from 0.002 at an ASI of 0.8 to 0.55 at an ASI of 1.07. At higher ASI, it decreased slightly to 0.2 at an ASI of similar to 1.3. Likewise, it was one to two orders of magnitude higher in chloridic fluids compared to those found in H2O experiments. Fluid-melt partition coefficients of La and Y increased from 0.002 at an ASI of similar to 0.8 to similar to 0.1 at an ASI of 1.2. In the same ASI range, the Yb partition coefficient increased to a maximum value of 0.02. Even at high salinities all elements fractionate into the melt. The compositional dependence of the partitioning data imply that both melt composition and fluid composition have a strong influence on trace element behavior and that complexation of Ba. REE and Y tin the fluid is not only controlled by the presence of Cl- in the fluid. Instead, interaction of these elements with major melt components dissolved in the fluid is very likely.
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.
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.
In piston cores from the open subarctic Pacific and the Okhotsk Sea, diatom-bound delta N-15 (delta N-15(db)), biogenic opal, calcium carbonate, and barium were measured from coretop to the previous glacial maximum (MIS 6). Glacial intervals are generally characterized by high delta N-15(db) (similar to 8 parts per thousand) and low productivity, whereas interglacial intervals have a lower delta N-15(db) (5.7-6.3 parts per thousand) and indicate high biogenic productivity. These data extend the regional swath of evidence for nearly complete surface nutrient utilization during glacial maxima, consistent with stronger upper water column stratification throughout the subarctic region during colder intervals. An early deglacial decline in delta N-15(db) of 2 parts per thousand at similar to 17.5 ka, previously observed in the Bering Sea, is found here in the open subarctic Pacific record and arguably also in the Okhotsk, and a case can be made that a similar decrease in delta N-15(db) occurred in both regions at the previous deglaciation as well. The early deglacial delta N-15(db) decrease, best explained by a decrease in surface nutrient utilization, appears synchronous with southern hemisphere-associated deglacial changes and with the Heinrich 1 event in the North Atlantic. This delta N-15(db) decrease may signal the initial deglacial weakening in subarctic North Pacific stratification and/or a deglacial increase in shallow subsurface nitrate concentration. If the former, it would be the North Pacific analogue to the increase in vertical exchange inferred for the Southern Ocean at the time of Heinrich Event 1. In either case, the lack of any clear change in paleoproductivity proxies during this interval would seem to require an early deglacial decrease in the iron-to-nitrate ratio of subsurface nutrient supply or the predominance of light limitation of phytoplankton growth during the deglaciation prior to Bolling-Allerod warming.
Understanding the interactions between the different processes that control the geothermal and fluid flow fields in sedimentary basins is crucial for exploitation of geothermal energy. Numerical models provide predictive and feasible information for a correct assessment of geothermal resources especially in areas where data acquisition is demanding. Here, we present results from numerical efforts to characterize the thermal structure and its interaction with the fluid system for the area of the North East German Basin (NEGB). The relative impact of the different (diffusive and advective) processes affecting the hydrothermal setting of the basin are investigated by means of three- dimensional numerical simulations. Lithospheric-scale numerical models are evaluated to understand the specific thermal signature of the relevant factors influencing the present-day conductive geothermal field in the NEGB. Shallow and deep structural controls on the thermal configuration of the basin are addressed and quantified. Interaction between the resulting thermal field and the active fluid system is investigated by means of three-dimensional simulations of coupled fluid flow and heat transport. Factors influencing stability and reliability of modeling predictions are discussed. The main effort is to build a physically consistent model for the basin which integrates the impacts of thermal gradients on the regional fluid regime and their coupling with the main geological units defining the basin.
Bedrock fracturing and rock strength are widely believed to influence landscape morphology and erosional resistance. Yet, understanding of the quantitative relationship between rock-mass strength and landscape evolution remains limited. Here we present a new application of seismic refraction surveys that uses variations in seismic velocity to interpret differences in bedrock fracture density with depth. We use a comparative study of Fiordland and the western Southern Alps of New Zealand to examine how differences in rock type and bedrock fracturing influence landscape morphology and landslide response to rock uplift. In both regions, slopes appear invariant with differential rock-uplift rates and slope distributions reveal modal hillslope angles of similar to 32 degrees. The majority of landslides initiate on slopes steeper than the modal hillslope angle, however, landslide magnitude-frequency distributions reveal order-of-magnitude differences between the regions, with Fiordland experiencing considerably smaller and less frequent landsliding events. Landslide-driven denudation rates of similar to 9 mm/yr in the western Southern Alps and between similar to 0.1 and 0.3 mm/yr in Fiordland approximate estimates of long-term rock-uplift rates for each region. The invariance of hillslope angles, near-normal slope distributions, predominance of landslide initiation on slopes steeper than modal values, and the apparent balance between rates of uplift and landslide-driven erosion suggest that hillslopes in both regions are at threshold angles. Their similar modal slopes further suggest that both ranges are characterized by equivalent rock-mass strength, despite striking differences in lithology. Additionally, our seismic analysis reveals nearly identical surface p-wave velocities. The unexpected equivalence of both modal slopes and surface velocities between these lithologically distinct ranges is attributed to contrasting degrees of surface fracturing that have differentially affected the intact rock properties, such that they now yield equivalent surface velocities and hillslope-scale strengths. Given that surface fractures help regulate threshold angles by modulating hillslope strength; we propose that shallow seismic velocities may provide a quantitative proxy for rock-mass strength. We define two contrasting fracture and landsliding environments. In Fiordland, dense geomorphic fracturing that is focused within the shallow subsurface appears to limit the depth and magnitude of most bedrock landslides. Conversely, in the western Southern Alps, tectonic forces produce pervasive fracturing with depth that results in larger, and deeper landslides. Our data suggest that bedrock fracturing at the Earth's surface modulates threshold hillslope angles, whereas the depth of bedrock fracturing influences the magnitude and frequency of landslide response to tectonic rock uplift.
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.
Shallowly situated evaporites in built-up areas are of relevance for urban and cultural development and hydrological regulation. The hazard of sinkholes, subrosion depressions and gypsum karst is often difficult to evaluate and may quickly change with anthropogenic influence. The geophysical exploration of evaporites in metropolitan areas is often not feasible with active industrial techniques. We collect and combine different passive geophysical data as microgravity, ambient vibrations, deformation and hydrological information to study the roof morphology of shallow evaporites beneath Hamburg, Northern Germany. The application of a novel gravity inversion technique leads to a 3-D depth model of the salt diapir under study. We compare the gravity-based depth model to pseudo-depths from H/V measurements and depth estimates from small-scale seismological array data. While the general range and trend of the diapir roof is consistent, a few anomalous regions are identified where H/V pseudo-depths indicate shallower structures not observed in gravity or array data. These are interpreted by shallow residual caprock floaters and zones of increased porosity. The shallow salt structure clearly correlates with a relative subsidence in the order of 2 mm yr(-1). The combined interpretation of roof morphology, yearly subsidence rates, chemical analyses of groundwater and of hydraulic head in aquifers indicates that the salt diapir beneath Hamburg is subject to significant ongoing dissolution that may possibly affect subrosion depressions, sinkhole distribution and land usage. The combined analysis of passive geophysical data may be exemplary for the study of shallow evaporites beneath other urban areas.
Phosphorus availability is commonly assumed to limit productivity in lowland tropical rainforests, yet there is relatively little information on the chemical forms of soil phosphorus in such ecosystems. We used the Hedley sequential fractionation scheme to assess phosphorus chemistry in five soils supporting tropical rainforest on Barro Colorado Island, Republic of Panama. The soils represented a range of orders (Inceptisols, Alfisols, and Oxisols) formed on contrasting geological substrates and topography, but under uniform climate and vegetation. Total phosphorus in surface horizons ranged between 315 and 1114 mg P kg(-1), being lowest on a soil derived from marine sediments and highest on soils derived from andesite. The majority of the phosphorus occurred in recalcitrant forms, although between 14% and 39% occurred as organic phosphorus. Readily-available phosphate, as extracted by anion-exchange membranes, occurred in small concentrations (4-13 mg P kg(-1)), although labile phosphorus, defined as phosphate extracted by anion-exchange membrane plus inorganic and organic phosphorus extracted by 0.5 M NaHCO3, accounted for between 4.7% and 11.4% of the total soil phosphorus. Our results indicate a strong control of geology and topography on soil phosphorus in tropical rainforests, which may have important implications for understanding the diversity and distribution of plant species in such ecosystems. Further, some of the most common soils on Barro Colorado Island, including those on the 50 ha forest dynamics plot, are rich in phosphorus despite their relatively advanced stage of pedogenesis.
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.
Ambient vibration measurements with small, temporary arrays that produce estimates of surface wave dispersion have become increasingly popular as a low-cost, non-invasive tool for site characterisation. An important requirement for these measurements to be meaningful, however, is the temporal consistency and repeatability of the resulting dispersion and spatial autocorrelation curve estimates. Data acquired within several European research projects (NERIES task JRA4, SESAME, and other multinational experiments) offer the chance to investigate the variability of the derived data products. The dataset analysed here consists of repeated array measurements, with several years of time elapsed between them. The measurements were conducted by different groups in different seasons, using different instrumentations and array layouts, at six sites in Greece and Italy. Ambient vibration amplitude spectra and locations of dominant sources vary between the two measurements at each location. Still, analysis indicates that this does not influence the derived dispersion information, which is stable in time and neither influenced by the instrumentation nor the analyst. The frequency range over which the dispersion curves and spatial autocorrelation curves can be reliably estimated depends on the array dimensions (minimum and maximum aperture) used in the specific deployment, though, and may accordingly vary between the repeated experiments. The relative contribution of Rayleigh and Love waves to the wavefield can likewise change between repeated measurements. The observed relative contribution of Rayleigh waves is generally at or below 50%, with especially low values for the rural sites. Besides, the visibility of higher modes depends on the noise wavefield conditions. The similarity of the dispersion and autocorrelation curves measured at each site indicates that the curves are stable, mainly determined by the sub-surface structure, and can thus be used to derive velocity information with depth. Differences between velocity models for the same site derived from independently determined dispersion and autocorrelation curves-as observed in other studies-are consequently not adequately explained by uncertainties in the measurement part.
We observed vertically displaced coastal and river markers after the 27 February 2010 Chilean earthquake [moment magnitude (Mw) 8.8]. Land-level changes range between 2.5 and -1 meters, evident along an ~500-kilometers- long segment identified here as the maximum length of coseismic rupture. A hinge line located 120 kilometers from the trench separates uplifted areas, to the west, from subsided regions. A simple elastic dislocation model fits these observations well; model parameters give a similar seismic moment to seismological estimates and suggest that most of the plate convergence since the 1835 great earthquake was elastically stored and then released during this event.
The closed-chamber method is the most common approach to determine CH4 fluxes in peatlands. The concentration change in the chamber is monitored over time, and the flux is usually calculated by the slope of a linear regression function. Theoretically, the gas exchange cannot be constant over time but has to decrease, when the concentration gradient between chamber headspace and soil air decreases. In this study, we test whether we can detect this non- linearity in the concentration change during the chamber closure with six air samples. We expect generally a low concentration gradient on dry sites (hummocks) and thus the occurrence of exponential concentration changes in the chamber due to a quick equilibrium of gas concentrations between peat and chamber headspace. On wet (flarks) and sedge- covered sites (lawns), we expect a high gradient and near-linear concentration changes in the chamber. To evaluate these model assumptions, we calculate both linear and exponential regressions for a test data set (n = 597) from a Finnish mire. We use the Akaike Information Criterion with small sample second order bias correction to select the best-fitted model. 13.6%, 19.2% and 9.8% of measurements on hummocks, lawns and flarks, respectively, were best fitted with an exponential regression model. A flux estimation derived from the slope of the exponential function at the beginning of the chamber closure can be significantly higher than using the slope of the linear regression function. Non-linear concentration-overtime curves occurred mostly during periods of changing water table. This could be due to either natural processes or chamber artefacts, e.g. initial pressure fluctuations during chamber deployment. To be able to exclude either natural processes or artefacts as cause of non-linearity, further information, e.g. CH4 concentration profile measurements in the peat, would be needed. If this is not available, the range of uncertainty can be substantial. We suggest to use the range between the slopes of the exponential regression at the beginning and at the end of the closure time as an estimate of the overall uncertainty.
In this paper, a unique set of samples from the uppermost crust down to the lithospheric mantle of Jordan is analyzed for composition and petrophysical properties (density. thermal conductivity, radiogenic heat production) These data, covering a vertical section of almost 65 km. are used in conjunction with surface heat flow to generate a detailed and comprehensive lithospheric thermal model that reflects the conditions of the Arabian Shield (AS) prior to the post- Oligocene onset of lithosphere thinning and Voluminous basaltic volcanism. The pre-Miocene model geotherms, based on conductive surface heat flows of 55 and 60 mW m(-2). (a) meet the range of lithosphere-asthenosphere boundary depths of 110-160 km known from seismology, (b) conform to results of thermomechanical models on the on.-in of the Dead Sea basin that started in Miocene time. and (c) are consistent with typical xenolith-derived geotherms for terranes of similar age and lithospheric thickness. Moho temperatures (at depths between 35 and 40 km) of the AS in pre-Miocene times were most likely in the order of 530-650 degrees C, with mantle heat flows averaging between 24 and 29 mW m(-2) Results contradict former views of the late Proterozoic/early Cambrian-stabilized AS being an anomalously cold terrane A "cold" thermal structure inferred from previously measured low surface heat flows (generally <= 45 mW m(-2)) is inconsistent with the thickness, composition, and petrophysical properties of the stable lithosphere of the shield.
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.
Investigation of transient soil moisture profiles yields valuable information of near- surface processes. A recently developed reconstruction algorithm based on the telegraph equation allows the inverse estimation of soil moisture profiles along coated, three rod TDR probes. Laboratory experiments were carried out to prove the results of the inversion and to understand the influence of probe rod deformation and solid objects close to the probe in heterogeneous media. Differences in rod geometry can lead to serious misinterpretations in the soil moisture profile, but have small influence on the average soil moisture along the probe. Solids in the integration volume have almost no effect on average soil moisture, but result in locally slightly decreased moisture values. Inverted profiles obtained in a loamy soil with a clay content of about 16% were in good agreement with independent measurements.
Investigation of transient soil moisture profiles yields valuable information of near- surface processes. A recently developed reconstruction algorithm based on the telegraph equation allows the inverse estimation of soil moisture profiles along coated, three rod TDR probes. Laboratory experiments were carried out to prove the results of the inversion and to understand the influence of probe rod deformation and solid objects close to the probe in heterogonous media. Differences in rod geometry can lead to serious misinterpretations in the soil moisture profile but have small influence on the average soil moisture along the probe. Solids in the integration volume have almost no effect on average soil moisture but result in locally slightly decreased moisture values. Inverted profiles obtained in a loamy soil with a clay content of about 16% were in good agreement with independent measurements.
A new reconstruction of Alpine Tethys combines plate-kinematic modelling with a wealth of geological data and seismic tomography to shed light on its evolution, from sea-floor spreading through subduction to collision in the Alps. Unlike previous models, which relate the fate of Alpine Tethys solely to relative motions of Africa, Iberia and Europe during opening of the Atlantic, our reconstruction additionally invokes independent microplates whose motions are constrained primarily by the geological record. The motions of these microplates (Adria, Iberia, Alcapia, Alkapecia, and Tiszia) relative to both Africa and Europe during Late Cretaceous to Cenozoic time involved the subduction of remnant Tethyan basins during the following three stages that are characterized by contrasting plate motions and driving forces: (1) 131-84 Ma intra-oceanic subduction of the Ligurian part of Alpine Tethys attached to Iberia coincided with Eo-alpine orogenesis in the Alcapia microplate, north of Africa. These events were triggered primarily by foundering of the older (170-131 Ma) Neotethyan subduction slab along the NE margin of the composite African-Adriatic plate; subduction was linked by a sinistral transform system to E-W opening of the Valais part of Alpine Tethys; (2) 84-35 Ma subduction of primarily the Piemont and Valais parts of Alpine Tethys which were then attached to the European plate beneath the overriding African and later Adriatic plates. NW translation of Adria with respect to Africa was accommodated primarily by slow widening of the Ionian Sea; (3) 35 Ma-Recent rollback subduction of the Ligurian part of Alpine Tethys coincided with Western Alpine orogenesis and involved the formation of the Gibraltar and Calabrian arcs. Rapid subduction and arc formation were driven primarily by the pull of the gravitationally unstable, retreating Adriatic and African slabs during slow convergence of Africa and Europe. The upper European-Iberian plate stretched to accommodate this slab retreat in a very mobile fashion, while the continental core of the Adriatic microplate acted as a rigid indenter within the Alpine collisional zone. The subducted lithosphere in this reconstruction can be correlated with slab material imaged by seismic tomography beneath the Alps and Apennines, as well as beneath parts of the Pannonian Basin, the Adriatic Sea, the Ligurian Sea, and the Western Mediterranean. The predicted amount of subducted lithosphere exceeds the estimated volume of slab material residing at depth by some 10-30%, indicating that parts of slabs may be superposed within the mantle transition zone and/or that some of this subducted lithosphere became seismically transparent.