Upper Oligocene-Lower Miocene shallowing-upwardlacustrine sequences controlled by periodic and non-periodic processes(Ebro Basin, northeastern Spain)

The lacustrine facies from two sections (Candasnos and Fraga) ofthe Oligocene-Miocene Torrente de Cinca lithostratigraphic Unit in thecentral part of the Ebro Basin (Spain) have been analysed to determine theinfluence of orbital parameters in lacustrine sedimentation. The unit ispredominantly composed of limestones and marls, and represents a shallowlacustrine freshwater system. The sedimentological features of the faciesstudied demonstrate that the lower part of the Candasnos section representsoffshore lacustrine subenvironments whereas the upper part, and the whole ofthe Fraga section, characterise marginal lacustrine areas. Series of stratalthickness variations of limestone, marl, and limestone/marl couplets fromboth sections have been analysed using spectral analysis. This shows thatinformation corresponding to periodic cycles only appears in the offshorefacies, that is to say, in the lower part of Candasnos section, and disappearsin the marginal facies where non-periodic cycles exist. Furthermore, thespectral analysis of the offshore facies highlights the existence of a peak inthe power spectrum with a period of around 7 (6.8 to 7.8) that can berecognised in the field as shallowing-upward lacustrine sequences.Magnetostratigraphic data from the Candasnos section allow us to establish atime span of 2,808 years for the limestone/marl couplet from the lower partof this section, and between 19,000 and 22,000 years for the periodic cycleidentified, thus representing the climatic precession cycle. Shallowingsequences from marginal areas do not correspond with any periodiccycle.     

Use of Noise to Augment Training Data: A Neural Network Method of Mineral–Potential Mapping in Regions of Limited Known Deposit Examples

One of the main factors that affects the performance of MLP neural networks trained using the backpropagation algorithm in mineral-potential mapping isthe paucity of deposit relative to barren training patterns. To overcome this problem, random noise is added to the original training patterns in order to create additional synthetic deposit training data. Experiments on the effect of the number of deposits available for training in the Kalgoorlie Terrane orogenic gold province show that both the classification performance of a trained network and the quality of the resultant prospectivity map increasesignificantly with increased numbers of deposit patterns. Experiments are conducted to determine the optimum amount of noise using both uniform and normally distributed random noise. Through the addition of noise to the original deposit training data, the number of deposit training patterns is increased from approximately 50 to 1000. The percentage of correct classifications significantly improves for the independent test set as well as for deposit patterns in the test set. For example, using ±40% uniform random noise, the test-set classification performance increases from 67.9% and 68.0% to 72.8% and 77.1% (for test-set overall and test-set deposit patterns, respectively). Indices for the quality of the resultant prospectivity map, (i.e. D/A, D × (D/A), where D is the percentage of deposits and A is the percentage of the total area for the highest prospectivity map-class, and area under an ROC curve) also increase from 8.2, 105, 0.79 to 17.9, 226, 0.87, respectively. Increasing the size of the training-stop data set results in a further increase in classification performance to 73.5%, 77.4%, 14.7, 296, 0.87 for test-set overall and test-set deposit patterns, D/A, D × (D/A), and area under the ROC curve, respectively.     

Use of Fuzzy Membership Input Layers to Combine Subjective Geological Knowledge and Empirical Data in a Neural Network Method for Mineral-Potential Mapping

Use of GIS layers, in which the cell values represent fuzzy membership variables, is an effective method of combining subjective geological knowledge with empirical data in a neural network approach to mineral-prospectivity mapping. In this study, multilayer perceptron (MLP), neural networks are used to combine up to 17 regional exploration variables to predict the potential for orogenic gold deposits in the form of prospectivity maps in the Archean Kalgoorlie Terrane of Western Australia. Two types of fuzzy membership layers are used. In the first type of layer, the statistical relationships between known gold deposits and variables in the GIS thematic layer are used to determine fuzzy membership values. For example, GIS layers depicting solid geology and rock-type combinations of categorical data at the nearest lithological boundary for each cell are converted to fuzzy membership layers representing favorable lithologies and favorable lithological boundaries, respectively. This type of fuzzy-membership input is a useful alternative to the 1-of-N coding used for categorical inputs, particularly if there are a large number of classes. Rheological contrast at lithological boundaries is modeled using a second type of fuzzy membership layer, in which the assignment of fuzzy membership value, although based on geological field data, is subjective. The methods used here could be applied to a large range of subjective data (e.g., favorability of tectonic environment, host stratigraphy, or reactivation along major faults) currently used in regional exploration programs, but which normally would not be included as inputs in an empirical neural network approach.     

Optimal Model for Geoid Determination from Airborne Gravity

Two different approaches for transformation of airborne gravity disturbances, derived from gravity observations at low-elevation flying platforms, into geoidal undulations are formulated, tested and discussed in this contribution. Their mathematical models are based on Green's integral equations. They are in these two approaches defined at two different levels and also applied in a mutually reversed order. While one of these approaches corresponds to the classical method commonly applied in processing of ground gravity data, the other approach represents a new method for processing of gravity data in geoid determination that is unique to airborne gravimetry. Although theoretically equivalent in the continuous sense, both approaches are tested numerically for possible numerical advantages, especially due to the inverse of discretized Fredholm's integral equation of the first kind applied on different data. High-frequency synthetic gravity data burdened by the 2-mGal random noise, that are expected from current airborne gravity systems, are used for numerical testing. The results show that both approaches can deliver for the given data a comparable cm-level accuracy of the geoidal undulations. The new approach has, however, significantly higher computational efficiency. It would be thus recommended for real life geoid computations. Additional errors related to regularization of gravity data and the geoid, and to accuracy of the reference field, that would further deteriorate the quality of estimated geoidal undulations, are not considered in this study.     

Extracting flow structures from tracer data

A method of visualizing structures in closed chaotic flows out of homogenous particle distributions is presented in the example of models of a meandering jet. To this end, the system will be leaked or opened up by defining a region of the flow, so that a particle is considered to be escaped if it leaves this region. By applying this method to an ensemble of nonescaped tracers, we are able to characterize mixing processes by visualizing the converging and stretching filamentations (stable and unstable manifolds) in the flow without using additional mathematical tools. The possibility of applying the algorithm to analyze buoy data, and a comparison with the finite time manifolds are discussed.     

Coupled inverse modelling of groundwater flow and mass transport and the worth of concentration data

This paper presents the extension of the self-calibrating method to the coupled inverse modelling of groundwater flow and mass transport. The method generates equally likely solutions to the inverse problem that display the variability as observed in the field and are not affected by a linearisation of the state equations. Conditioning to the state variables is measured by an objective function including, among others, the mismatch between the simulated and measured concentrations. Conditioning is achieved by minimising the objective function by gradient-based methods. The gradient contains the partial derivatives of the objective function with respect to: log conductivities, log storativities, prescribed heads at boundaries, retardation coefficients and mass sources. The derivatives of the objective function with respect to log conductivity are the most cumbersome and need the most CPU-time to be evaluated. For this reason, to compute this derivative only advective transport is considered. The gradient is calculated by the adjoint-state method. The method is demonstrated in a controlled, synthetic study, in which the worth of concentration data is analysed. It is shown that concentration data are essential to improve transport predictions and also help to improve aquifer characterisation and flow predictions, especially in the upstream part of the aquifer, even in the case that a considerable amount of other experimental data like conductivities and heads are available. Besides, conditioning to concentration data reduces the ensemble variances of estimated transmissivity, hydraulic head and concentration.     

Digital model of corrected accumulated flow for peak discharge data acquisition and drainage system design

The design of a drainage system for a roofing slate quarry was implemented by the enhancement of discharge peak estimation, and the uncertainty inevitably associated with the engineering model was reduced.

The development of a topographical, geological, and vegetation cover database developed from a Geographical Information System (GIS) allowed for the definition of the drainage network for a hydraulic system, along with the calculation of the runoff coefficient. This is applied to the digital model of accumulated flow (DMF) as a weight correction coefficient, using a matrix-based model at 5×5 m resolution. The new digital model of corrected accumulated flow (DMCF) is the result of combining the thematic maps with the map of slope <3%, which was previously created from the slope model. It is demonstrated that this new model allows to apply the “Rational Method” on cartographic units defined by the GIS.

The DMCF is compared with other traditional applications of the Rational Method based on the calculation of the discharge peak considering: (1) the drainage basin as a single watershed or (2) defining an average runoff coefficient in each sub-watershed. Both approaches have bigger discharge peaks than those obtained by the DMCF since the slope, lithology, and vegetation cover have average values, and the runoff coefficient is poorly defined, increasing the uncertainty in the discharge peak.     

Lower crustal granulite xenoliths from the Pannonian Basin,Hungary, Part 2: Sr–Nd–Pb–Hf and O isotope evidence for formation of continental lower crust by tectonic emplacement of oceanic crust

Mafic granulite xenoliths from the lower crust of the Pannonian Basin are dominated by LREE-depleted bulk-rock compositions. Many of these have MORB-like ¹⁴³Nd/¹⁴⁴Nd but ⁸⁷Sr/⁸⁶Sr is elevated relative to most MORBs. Their '¹⁸O values cover a wide range from +3.8 to +9.5‰. A group of LREE-enriched mafic granulites have higher ⁸⁷Sr/⁸⁶Sr (0.704-0.708) and lower ¹⁴³Nd/¹⁴⁴Nd (0.5128-0.5124), with higher '¹⁸O values on average (+7.8 to +10.6‰) than the LREE-depleted granulites. The LREE-enriched granulites are, however, isotopically similar to newly discovered metasedimentary granulite xenoliths. A sublinear correlation in )Hf-)Nd isotope space has a shallower slope than the crust-mantle array, with the metasedimentary rocks forming the low )Hf end member; the radiogenic end is restricted to the LREE-depleted granulites and these overlap the field of MORB. Pb isotopes for the LREE-depleted samples are less radiogenic on average than those of the LREE-enriched and metasedimentary xenoliths, and metasedimentary granulites have consistently higher ²⁰⁸Pb/²⁰⁴Pb. The wide range in '¹⁸O over a restricted range in Nd and Sr isotope values, in combination with the predominance of LREE-depleted trace-element compositions, is consistent with an origin as a package of hydrothermally altered oceanic crust. The existence of '¹⁸O values lower than average MORB and/or mantle peridotite requires that at least some of these rocks were hydrothermally altered at high temperature, presumably in the oceanic lower crust. The low ¹⁴³Nd/¹⁴⁴Nd of the LREE-enriched mafic granulites cannot be explained by simple mixing between a LREE-depleted melt and an enriched component, represented by the recovered metasediments. Instead, we interpret these rocks as the metamorphic equivalent of the shallowest levels of the ocean crust where pillow basalts are intimately intercalated with oceanic sediments. A possible model is accretion of oceanic crustal slices during subduction and convergence, followed by high-grade metamorphism during the Alpine orogeny.