Seismic Hazard Estimate at the Iberian Peninsula

—?Seismic hazard at the Iberian Peninsula has been evaluated by using a methodology which combines both zonified and non-zonified probabilistic methods. Seismic sources are used when considering zones where certain calculation parameters may be considered homogeneous, as in zonified methods, while, on the other hand, earthquakes are considered wherever it has taken place, as in non-zonified methods. The methodology which is applied in this paper has been originally used to calculate the seismic hazard maps in the United States. In our case, it has been necessary to adapt the method to the specific features of the seismicity in the Iberian Peninsula and its geographical surroundings, not only with respect to its distribution and characteristics, but also with respect to the properties of the seismic catalog used.¶Geographically, the main feature of the result is the fact that it reflects both historical seismicity and current seismic clusters of the region. Despite the smoothing, maps show marked differences between several seismic zones; these differences becoming more noticeable as exposure time increases. Maximum seismic hazard is found to be in the southwestern region of the Peninsula, especially in the area of the Cape St. Vicent, and around Lisbon. The uncertainty of the results, without considering that due to the attenuation laws, as deduced from the other evaluation parameters, is quite stable, being more sensitive to the parameters b and m _max of the Gutenberg-Richter relation.     

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.     

On the circulation in the East Frisian Wadden Sea: numerical modeling and data analysis

 In this paper we use a combination of numerical modeling and data analysis to gain a better understanding of the major characteristics of the circulation in the East Frisian Wadden Sea. In particular, we concentrate on the asymmetry of the tidal wave and its modulation in the coastal area, which results in a complex pattern of responses to the sea-level forcing from the North Sea. The numerical simulations are based on the 3-D primitive equation General Estuarine Transport Model (GETM) with a horizontal resolution of 200 m and terrain-following vertical coordinates. The model is forced at its open boundaries with sea-level data from an operational model for the German Bight (German Hydrographic Office). The validation data for our model simulations include time series of tidal gauge data and surface currents measured at a pile in the back-barrier basin of the Island Langeoog, as well as several ADCP transects in the Accumer Ee tidal inlet. Circulation and turbulence characteristics are investigated for typical situations driven by spring and neap tides, and the analysis is focused on dominating temporal and spatial patterns. By investigating the response of five back-barrier basins with rather different morphologies to external forcing, an attempt is made to elucidate the dominating physical balances controlling the circulation in the individual sub-basins. It is demonstrated that the friction at the seabed tends to slow down the tidal signal in the shallow water. This leads to the establishment of flood dominance in the shallow sea north of the barrier islands. South of the islands, where the water volume of the channels at low tide is smaller than the tidal prism, the asymmetry of the tidal signal is shifted towards ebb dominance, a feature which is particularly pronounced at spring tide. At the northern open boundary, the tidal wave propagating from west to east generates a sea-level difference of ∼1 m along the boundary, and thereby triggers vigorous alongshore currents. The frictional control in the model is located in the inlets, as well as along the northern boundary. The correlation between velocity and turbulent kinetic energy tends to the establishment of a net southward transport, giving theoretical support to the observed accumulation of sediments on the intertidal flats. Weak turbulence along the northern shores of the barrier islands and the small magnitude of the residual currents there promote accumulation of suspended matter in these areas, although wave action will generally counteract this effect. Received: 29 May 2002 / Accepted: 26 September 2002 Responsible Editor: Jean-Marie Beckers Acknowledgements We are indebted to S. Dick for providing the data from the operational model of BSH and to B. Flemming for the useful discussions. The topography data and Fig. 1 have been prepared in cooperation with F. Meyer. Figure 2 has been prepared by G. Brink-Spalink. We also thank for the comments from an anonymous reviewer which helped to improve our paper.     

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.