Lower crustal melting and the role of open-system processes in the genesis of syn-orogenic quartz diorite–granite–leucogranite associations: constraints from Sr–Nd–O isotopes from the Bandombaai Complex, Namibia

The Bandombaai Complex (southern Kaoko Belt, Namibia) consists of three main intrusive rock types including metaluminous hornblende- and sphene-bearing quartz diorites, allanite-bearing granodiorites and granites, and peraluminous garnet- and muscovite-bearing leucogranites. Intrusion of the quartz diorites is constrained by a U–Pb zircon age of 540±3 Ma.

Quartz diorites, granodiorites and granites display heterogeneous initial Nd- and O isotope compositions (_{Nd (540 Ma)}=−6.3 to −19.8; δ¹⁸O=9.0–11.6‰) but rather low and uniform initial Sr isotope compositions (⁸⁷Sr/⁸⁶Sr_initial=0.70794–0.70982). Two leucogranites and one aplite have higher initial ⁸⁷Sr/⁸⁶Sr ratios (0.70828–0.71559), but similar initial _Nd (−11.9 to −15.8) and oxygen isotope values (10.5–12.9‰). The geochemical and isotopic characteristics of the Bandombaai Complex are distinct from other granitoids of the Kaoko Belt and the Central Zone of the Damara orogen. Our study suggests that the quartz diorites of the Bandombaai Complex are generated by melting of heterogeneous mafic lower crust. Based on a comparison with results from amphibolite-dehydration melting experiments, a lower crustal garnet- and amphibole-bearing metabasalt, probably enriched in K₂O, is a likely source rock for the quartz diorites. The granodiorites/granites show low Rb/Sr (<0.6) ratios and are probably generated by partial melting of meta-igneous (intermediate) lower crustal sources by amphibole-dehydration melting. Most of the leucogranites display higher Rb/Sr ratios (>1) and are most likely generated by biotite-dehydration melting of heterogeneous felsic lower crust. All segments of the lower crust underwent partial melting during the Pan-African orogeny at a time (540 Ma) when the middle crust of the central Damara orogen also underwent high T, medium P regional metamorphism and melting. Geochemical and isotope data from the Bandombaai Complex suggest that the Pan-African orogeny in this part of the orogen was not a major crust-forming episode. Instead, even the most primitive rock types of the region, the quartz diorites, represent recycled lower crustal material.     

Use of Geomorphological Information in Indirect Landslide Susceptibility Assessment

The objective of this paper is to evaluate the importance of geomorphological expert knowledge in the generation of landslide susceptibility maps, using GIS supported indirect bivariate statistical analysis. For a test area in the Alpago region in Italy a dataset was generated at scale 1:5,000. Detailed geomorphological maps were generated, with legends at different levels of complexity. Other factor maps, that were considered relevant for the assessment of landslide susceptibility, were also collected, such as lithology, structural geology, surficial materials, slope classes, land use, distance from streams, roads and houses. The weights of evidence method was used to generate statistically derived weights for all classes of the factor maps. On the basis of these weights, the most relevant maps were selected for the combination into landslide susceptibility maps. Six different combinations of factor maps were evaluated, with varying geomorphological input. Success rates were used to classify the weight maps into three qualitative landslide susceptibility classes. The resulting six maps were compared with a direct susceptibility map, which was made by direct assignment of susceptibility classes in the field. The analysis indicated that the use of detailed geomorphological information in the bivariate statistical analysis raised the overall accuracy of the final susceptibility map considerably. However, even with the use of a detailed geomorphological factor map, the difference with the separately prepared direct susceptibility map is still significant, due to the generalisations that are inherent to the bivariate statistical analysis technique.     

Timing of deformation in the Norseman-Wiluna Belt, Yilgarn Craton, Western Australia

Establishing relative and absolute time frameworks for the sedimentary, magmatic, tectonic and gold mineralisation events in the Norseman-Wiluna Belt of the Archean Yilgarn Craton of Western Australia, has long been the main aim of research efforts. Recently published constraints on the timing of sedimentation and absolute granite ages have emphasized the shortcomings of the established rationale used for interpreting the timing of deformation events. In this paper the assumptions underlying this rationale are scrutinized, and it is shown that they are the source of significant misinterpretations. A revised time chart for the deformation events of the belt is established. The first shortening phase to affect the belt, D₁, was preceded by an extensional event D_1e and accompanied by a change from volcanic-dominated to plutonic-dominated magmatism at approximately 2685–2675 Ma. Later extension (D_2e) controlled deposition of the ca 2655 Ma Kurrawang Sequence and was followed by D₂, a major shortening event, which folded this sequence. D₂ must therefore have started after 2655 Ma—at least 20 Ma later than previously thought and after the voluminous 2670–2655 Ma high-Ca granite intrusion. Younger transcurrent deformation, D₃–D₄, waned at around 2630 Ma, suggesting that the crustal shortening deformation cycle D₂–D₄ lasted approximately 20–30 Ma, contemporaneous with low-volume 2650–2630 Ma low-Ca granites and alkaline intrusions. Time constraints on gold deposits suggest a late mineralisation event between 2640–2630 Ma. Thus, D₂–D₄ deformation cycle and late felsic magmatism define a 20–30 Ma long tectonothermal event, which culminated with gold mineralisation. The finding that D₂ folding took place after voluminous high-Ca granite intrusion led to research into the role of competent bodies during folding by means of numerical models. Results suggest that buoyancy-driven doming of pre-tectonic competent bodies trigger growth of antiforms, whereas non-buoyant, competent granite bodies trigger growth of synforms. The conspicuous presence of pre-folding granites in the cores of anticlines may be a result from active buoyancy doming during folding.     

Modelling the short-term response of the Greenland ice-sheet to global warming

 A two-dimensional vertically integrated ice flow model has been developed to test the importance of various processes and concepts used for the prediction of the contribution of the Greenland ice-sheet to sea-level rise over the next 350 y (short-term response). The mass balance is modelled by the degree-day method and the energy-balance method. The lithosphere is considered to respond isostatically to a point load and the time evolution of the bedrock follows from a viscous asthenosphere. According to the IPCC-IS92a scenario (with constant aerosols after 1990) the Greenland ice-sheet is likely to cause a global sea level rise of 10.4 cm by 2100 AD. It is shown, however, that the result is sensitive to precise model formulations and that simplifications as used in the sea-level projection in the IPCC-96 report yield less accurate results. Our model results indicate that, on a time scale of a hundred years, including the dynamic response of the ice-sheet yields more mass loss than the fixed response in which changes in geometry are not incorporated. It appears to be important to consider sliding, as well as the fact that climate sensitivity increases for larger perturbations. Variations in predicted sea-level change on a time scale of hundred years depend mostly on the initial state of the ice-sheet. On a time scale of a few hundred years, however, the variability in the predicted melt is dominated by the variability in the climate scenarios. Received: 21 August 1996/Accepted: 12 May 1997     

A record of accelerated erosion in the recent sediments of Blelham Tarn in the English Lake district

Two frozen cores from Blelham Tarn were subsampled and measured using mineral magnetic, loss-on-ignition (LOI), radiometric, granulometric and diatom analyses. A detailed chronology was established using varves, radioisotopes and diatoms. This has enabled an accurately dated reconstruction of sedimentation over the past forty years. Despite a large increase in lake productivity, evidence suggests that the observed exponential increase in sedimentation rates can be attributed to erosion within the catchment. The predominant sediment source has been identified as surface soil. A comparison between the trend of accelerated sedimentation and the record of increased sheep stocking density for the area within which the most of the catchment lies, as well as observations of contemporary surface processes within the catchment, both suggest that much of the recent erosion is a direct response to increased pressure from sheep grazing.     

Flexible spectral methods for the generation of random fields with power-law semivariograms

Random field generators serve as a tool to model heterogeneous media for applications in hydrocarbon recovery and groundwater flow. Random fields with a power-law variogram structure, also termed fractional Brownian motion (fBm) fields, are of interest to study scale dependent heterogeneity effects on one-phase and two-phase flow. We show that such fields generated by the spectral method and the Inverse Fast Fourier Transform (IFFT) have an incorrect variogram structure and variance. To illustrate this we derive the prefactor of the fBm spectral density function, which is required to generate the fBm fields. We propose a new method to generate fBm fields that introduces weighting functions into the spectral method. It leads to a flexible and efficient algorithm. The flexibility permits an optimal choice of summation points (that is points in frequency space at which the weighting function is calculated) specific for the autocovariance structure of the field. As an illustration of the method, comparisons between estimated and expected statistics of fields with an exponential variogram and of fBm fields are presented. For power-law semivariograms, the proposed spectral method with a cylindrical distribution of the summation points gives optimal results.     

The identifiability of parameters in a water quality model of the Biebrza River, Poland

The identifiability of model parameters of a steady state water quality model of the Biebrza River and the resulting variation in model results was examined by applying the Monte Carlo method which combines calibration, identifiability analysis, uncertainty analysis, and sensitivity analysis. The water quality model simulates the steady state concentration profiles of chloride, phosphate, ammonium, and nitrate as a function of distance along a river. The water quality model with the best combination of parameter values simulates the observed concentrations very well. However, the range of possible modelled concentrations obtained for other more or less equally eligible combinations of parameter values is rather wide. This range in model outcomes reflects possible errors in the model parameters. Discrepancies between the range in model outcomes and the validation data set are only caused by errors in model structure, or (measurement) errors in boundary conditions or input variables. In this sense the validation procedure is a test of model capability, where the effects of calibration errors are filtered out. It is concluded that, despite some slight deviations between model outcome and observations, the model is successful in simulating the spatial pattern of nutrient concentrations in the Biebrza River.     

Influence of the Load Wavelength on the Permeability of a Viscosity Interface in the Mantle

Assuming a radially stratified Newtonian mantle in a steady-state approximation, we demonstrate that the permeability of a viscosity interface at 660-km depth strongly depends on the wavelength of buoyancy forces driving the flow. The flow induced by long-wavelength loads penetrates through the boundary freely even if the viscosity increases by two orders. In contrast, the boundary is practically impermeable for short-wavelength loads located in the upper mantle. Thus, a stepwise increase of viscosity is a significant obstacle for small descending features in the upper mantle, but huge upper mantle downwellings, or upwellings formed in the-lower mantle can overcome it easily. This indicates that certain care is necessary in interpreting the seismic structure of the mantle by means of flow models. The global tomographic image includes only the first few degrees of the harmonic series and, consequently, its interpretation in terms of a present-day flow field results in a predominantly whole-mantle circulation even for extreme viscosity contrasts.