Single-grain zircon dating of the metamorphic and granitic rocks from the Biscayarhalvøya-Holtedahlfonna zone, north-west Spitsbergen

The Biscayarhalvøya-Holtedahlfonna zone (BHZ) in north-western Spitsbergen is a north-south trending, narrow horst, with crystalline basement rocks exposed under a Devonian unconformity. Previous K-Ar, ⁴⁰Ar/³⁹Ar and Rb-Sr analyses have confirmed the occurrence of Caledonian thermal events, and Grenvillian ages have been obtained by conventional zircon U-Pb and single-zircon Pb evaporation methods. A total of 55 zircon grains from three samples (an augen metagranite, a micaceous schist and a granitic neosome of migmatite) have been analysed by the single-zircon Pb evaporation method. The grains with the age range of ca. 950-1100 My (million years) are the major component in all three samples, suggesting tectono-thermal activity in that period. The detrital versus resorption orgin of the rounded shapes of these grains from the granitic neosome is not clear yet. Therefore, the ages of the migmatization and of the sedimentary protoliths are not concluded. The youngest presumed detrital grain from the granitic neosome is 1060 My old. The metagranite, cutting the Richarddalen unit, yielded grains with an age of ca. 950 Mya. A granite dyke with an age range of 955-968 My cuts the Biscayarhuken unit in the northern Liefdefjorden area. These indicate the sedimentary protoliths of the Richarddalen and Biscayarhuken units are pre-Neoproterozoic. The youngest detrital zircon ages of ca. 940 My indicate Neoproterozoic sedimentary protoliths of the Solanderfjellet micaceous schists. A significant population of zircon grains with an age range of 1600-1900 My in all three samples suggests a wide exposure of these rocks in the source areas during Meso- and Neoproterozoic times. Several Archean ages have also been obtained. The results are generally conformable with those obtained from north-western Spitsbergen.     

Wind erosion from military training lands in the Mojave Desert, California, U.S.A.

Military training activities reduce vegetation cover, disturb crusts, and degrade soil aggregates, making the land more vulnerable to wind erosion. The objective of this study was to quantify wind erosion rates for typical conditions at the Marine Corps Air Ground Combat Center, Twentynine Palms, CA, U.S.A. Five Big Spring Number Eight (BSNE) sampler stations were installed at each of five sites. Each BSNE station consisted of five BSNE samplers with the lowest sampler at 0·05 m and the highest sampler at 1·0 m above the soil surface. Once a month, sediment was collected from the samplers for analysis. Occurrence of saltating soil aggregates was recorded every hour using Sensits, one at each site. The site with the most erosion had a sediment discharge of 311 kg m⁻¹ over a period of 17 months. Other sites eroded much less because of significant rock cover or the presence of a crust. Hourly sediment discharge was estimated combining hourly Sensit count and monthly sediment discharge measured using BSNE samplers. More simultaneously measured data are needed to better characterize the relationship between these two and reconstruct a detailed time-series of wind erosion. This measured time-series can then be used for comparison with simulation results from process-based wind erosion models such as the Wind Erosion Prediction System (WEPS), once it has been adapted to the unique aspects of military lands.     

Towards non-linear inversion for characterization of time-lapse phenomena through numerical modelling

We compare two geophysical survey measurements of the same type made at different times in order to characterize the change in the geological medium during the elapsed time. The aim of this study is to develop a strategy using a full non-linear inversion algorithm as the interpretation tool. In this way, not only the location and the form of the changes are recovered, but also the changes in the material parameters of the geological medium can be estimated. In order to solve this fully non-linear problem, the so-called ‘multiplicative regularized contrast source inversion’ (MR-CSI) method is employed. The unique property of this iterative method is that it does not solve the forward problem at each iterative step. This makes it possible to use the non-linear inversion algorithm for large-scale computation problems. The numerical results show that by taking into account the non-linear nature of the problem, interpretation of the time-lapse data can be significantly improved, compared with that obtained using linear inversion.     

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.