Shear-wave velocity structure of the western part of the Mediterranean Sea from Rayleigh-wave analysis

The lithospheric structure of the western part of the Mediterranean Sea is shown by means of S-velocity maps, for depths ranging from 0 to 35 km, determined from Rayleigh-wave analysis. The traces of 55 earthquakes, which occurred from 2001 to 2003 in and around the study area have been used to obtain Rayleigh-wave dispersion. These earthquakes were registered by 10 broadband stations located on Iberia and the Balearic Islands. The dispersion curves were obtained for periods between 1 and 45 s, by digital filtering with a combination of MFT and TVF filtering techniques. After that, all seismic events were grouped in source zones to obtain a dispersion curve for each source-station path. These dispersion curves were regionalized and after inverted according to the generalized inversion theory, to obtain shear-wave velocity models for rectangular blocks with a size of 1° × 1°. The shear velocity structure obtained through this procedure is shown in the S-velocity maps plotted for several depths. These maps show the existence of lateral and vertical heterogeneity. In these maps is possible to distinguish several types of crust with an average S-wave velocity ranging from 2.6 to 3.9 km/s. The South Balearic Basin (SBB) is more characteristic of oceanic crust than the rest of the western Mediterranean region, as it is demonstrated by the crustal thickness. We also find a similar S-wave velocity (ranging from 2.6 km/s at the surface to 3.2 km/s at 10 km depth) for the Iberian Peninsula coast to Ibiza Island, the North Balearic Basin (NBB) and Mallorca Island. In the lower crust, the shear velocity reaches a value of 3.9 km/s. The base of the Moho is estimated from 15 to 20 km under Iberian Peninsula coast to Ibiza Island, continues towards NBB and increases to 20–25 km beneath Mallorca Island. While, the SBB is characterized by a thinner crust that ranges from 10 to 15 km, and a faster velocity. A gradual increase in velocity from the north to the south (especially in the upper 25 km) is obtained for the western part of the Mediterranean Sea. The base of the crust has a shear-wave velocity value around of 3.9 km/s for the western Mediterranean Sea area. This area is characterized by a thin crust in comparison with the crustal thickness of the eastern Mediterranean Sea area. This thin crust is related with the distensive tectonics that exists in this area. The low S-wave velocities obtained in the upper mantle might be an indication of a serpentinized mantle. The obtained results agree well with the geology and other geophysical results previously obtained. The shear velocity generally increases with depth for all paths analyzed in the study area.     

Structure of the crust and uppermost mantle of the Yanshan Belt and adjacent regions at the northeastern boundary of the North China Craton from Rayleigh Wave Dispersion Analysis

We have used Rayleigh wave dispersion analysis and inversion to produce a high resolution S-wave velocity imaging profile of the crust and uppermost mantle structure beneath the northeastern boundary regions of the North China Craton (NCC). Using waveform data from 45 broadband NCISP stations, Rayleigh wave phase velocities were measured at periods from 10 to 48 s and utilized in subsequent inversions to solve for the S-wave velocity structure from 15 km down to 120 km depth. The inverted lower crust and uppermost mantle velocities, about 3.75 km/s and 4.3 km/s on average, are low compared with the global average. The Moho was constrained in the depth range of 30–40 km, indicating a typical crustal thickness along the profile. However, a thin lithosphere of no more than 100 km was imaged under a large part of the profile, decreasing to only ~ 60 km under the Inner Mongolian Axis (IMA) where an abnormally slow anomaly was observed below 60 km depth. The overall structural features of the study region resemble those of typical continental rift zones and are probably associated with the lithospheric reactivation and tectonic extension widespread in the eastern NCC during Mesozoic–Cenozoic time. Distinctly high velocities, up to ~ 4.6 km/s, were found immediately to the south of the IMA beneath the northern Yanshan Belt (YSB), extending down to > 100-km depth. The anomalous velocities are interpreted as the cratonic lithospheric lid of the region, which may have not been affected by the Mesozoic–Cenozoic deformation process as strongly as other regions in the eastern NCC. Based on our S-wave velocity structural image and other geophysical observations, we propose a possible lithosphere–asthenosphere interaction scenario at the northeastern boundary of the NCC. We speculate that significant undulations of the base of the lithosphere, which might have resulted from the uneven Mesozoic–Cenozoic lithospheric thinning, may induce mantle flows concentrating beneath the weak IMA zone. The relatively thick lithospheric lid in the northern YSB may serve as a tectonic barrier separating the on-craton and off-craton regions into different upper mantle convection systems at the present time.     

Crustal thickness map of Brazil: Data compilation and main features

We present a crustal thickness map of Brazil and adjacent areas based on a compilation of data published in the literature as well as new measurements. We used crustal thicknesses mainly derived from seismic datasets such as deep seismic refraction experiments, receiver function analyses, and surface-wave dispersion velocities. Crustal thicknesses derived from modelling gravity anomalies commonly depend on assumptions, such as constant density contrast across the Moho interface, which are not always easily verifiable and were considered only along the continental shelf to fill large gaps in the seismic data. Our compilation shows that the crust in the stable continental area onshore has an average thickness of 39 ± 5 km (1-σ deviation) and that no clear difference can be observed between low altitude, intracratonic sedimentary basins, NeoProterozoic foldbelts (except for the Borborema Province), and cratonic areas. The thinnest crust is found in the Borborema Province of NE Brazil (30–35 km) and along a narrow belt within Tocantins Province (∼35 km), roughly parallel to the Eastern border of the Amazon craton, while the thickest crust is found in the Amazon and São Francisco cratons (41 ± 4 km), and the Paraná Basin (42 ± 4 km). Both the Ponta Grossa and the Rio Grande Arches are areas of thinned crust, and the western border of the Brazilian platform, near the sub-Andean region, seems to be characterized by a crustal thickness of less than 40 km. Although sparse in data coverage, we expect the resulting crustal thickness map to be useful for future studies of isostasy, dynamic topography, and crustal evolution of the country.     

2-D seismic tomographic and ray tracing modelling of the crustal structure across the Sudetes Mountains basing on SUDETES 2003 experiment data

The seismic data obtained during SUDETES 2003 experiment are analysed, and detailed crustal structure for profiles S02, S03 and S06 is presented using three different 2-D techniques: (1) “smooth” tomography of refracted waves travel times, (2) ray tracing of reflected and refracted waves, and (3) joint velocity and depth of reflector tomographic inversion. In spite of different interpretation techniques used, the models of the crustal structure show common characteristic features. The low velocity (V_p < 4 km/s) sedimentary layer was documented in the northeastern part of the study area. The topmost basement has in general a velocity of 5.8–6.0 km/s, and velocities at ca. 20 km depth are 6.15–6.25 km/s. The strong reflecting boundaries were found at 20–23 and 25–28 km depth with a velocity contrast about 0.4 km/s, and the highest velocities in the lowermost crust are 6.8–7.2 km/s. In general, the crust of the Bohemian Massif is slightly thicker (33–35 km) than in the northern part of the area. Velocities beneath Moho are relatively low, of 7.95 km/s. On the basis of well recorded reflected waves, mantle reflectors were discovered in the depth interval ca. 40–70 km. Apart of new results for the geology and tectonics of the area, some conclusion could be made about different techniques used. In the 2-D case the “clasical” ray tracing method with using all correlated phases gives the most adequate model of the structure, because of full, manual control of the model creation. The “smooth” first arrival travel times tomography, although very fast, is not satisfactory enough to describe the complex structure. So, the best candidate in 3-D case seems to be travel time tomography for both refracted and reflected waves in multi-layers models.     

An analysis of the gravity field over north Singhbhum, Bihar, India

About seven hundred gravity stations were established 2–3 miles apart over the Precambrian terrain of Singhbhum that lies between latitude 22° 15′ to 23°°15′N and longitude 85° to 87°E. Bouguer anomalies ranging from +4 to −62 mGal are found in the area. The observed Bouguer anomaly map was analyzed into regional and residual components. The residual anomaly map shows an excellent correlation with geology. The Singhbhum granite batholith is associated with several gravity lows. The residual anomaly map outlines nine plutonic granitic masses within the Singhbhum batholith. Negative residuals are also observed over some intrusive granites outside the batholith. Residual gravity highs are noted over the Dalma hill as well as over the Dhanjori lava complex on the eastern part of the Singhbhum batholith.Two-dimensional models suggestive of subsurface configuration of several major geologic units in the area are presented. These indicate that some of the plutonic granites within the Singhbhum batholith are of relatively large dimensions. The basin containing the Iron Ore Group of rocks to the west of the batholith, as well as the basin containing Singhbhum Group of rocks outside the Copper Belt thrust, may have sedimentary thicknesses of the order of 6–7 km. The Dalma lavas attain their maximum thickness of about 2.5 km in the form of a syncline, underneath which the Singhbhum Group of rocks is also found to be the thickest. The Copper Belt thrust, a major Precambrian fracture around the Singhbhum batholith, is moderately north-dipping near the surface but possibly attains a steeper slope at depth. The thrust appears to be quite deep seated. A threedimensional computer-based model for the Dhanjori lava—gabbro complex on the eastern part of the Singhbhum batholith has been deduced. Maximum thickness of these basic rocks is found to exist underneath a thin cap of granophyre. The geological implication of these results is discussed.Variation in the regional anomalies seems to be attributable to a mass deficiency under the Singhbhum batholith. The batholith may extend subsurfacially towards the north across the Copper Belt thrust. The northern tip of the batholith probably became dissected along the line of intersection of the two orogenic trends in the area and subsided. Over this subsided part, the Singhbhum Group of rocks was deposited at a later stage. Gravity data suggest a fairly large amount of subsidence in the area.     

Palaeozoic structural development along the Tornquist Zone, Kattegat area, Denmark

The interpretation of newly released commercial 2D reflection seismic data in the Kattegat area, Denmark, has provided us with a better understanding of the Palaeozoic tectonic processes along the Tornquist Fault Zone. A Base Palaeozoic time structure map, a Lower Palaeozoic TWT isopach map, a “true” Lower Palaeozoic TWT isopach map, an Upper Carboniferous/Lower Permian syn-rift TWT isopach map, a Top pre-Zechstein time structure map and a Zechstein combined TWT isopach and Palaeogeography map have been generated. The uniform Lower Palaeozoic sequence thickness in the Kattegat, both inside and outside the Tornquist Zone indicates only minor lateral movements if any, whereas the extensive Upper Silurian sequence, increasing in thickness to the north, indicates a relatively fast regional subsidence. The Base Palaeozoic time structure map and the Late Palaeozoic syn-rift isopach map show a clear Late Palaeozoic extension in the area. The syn-rift isopach map, in combination with the time-equivalent opening of the Skagerrak graben at right angles to the Tornquist Zone in the Kattegat, indicates that this extensional tectonic event had a dextral slip component. Measurements on internal extensional faults in the Tornquist Zone, give a minimum right-lateral displacement of 10.4 km. The footwall blocks were deeply eroded during the Early Permian rifting, and at Zechstein times the area became a peneplane. The Tornquist Zone was later exposed to several tectonic phases, where dextral slip played a role, indicated by the “push up” and “pull down” structures caused by restraining and releasing bends of the Børglum Fault. The dextral displacement along the Børglum Fault since the beginning of the Permian is in the order of 5–7 km based on the displacement of a Lower Palaeozoic local depocentre. Early Permian depocentres and faults, which gives a total amount of right-lateral displacement since the Early Palaeozoic in the order of 15–20 km. The continuously repeated tectonic episodes along the Tornquist Zone throughout most of the Phanerozoic, show that the zone was easily reactivated, implying deep-seated basement faults. The Tornquist Zone can be seen as a “buffer zone”, between continental blocks, whenever changes in the regional stress field are induced.     

The fine structure of the lower lithosphere underneath the Blue Road profile in northern Scandinavia

On the basis of DSS data from the recordings of compressional and shear waves made during the Blue Road shootings in 1972, velocity—depth models for the lower lithosphere are presented. A lid with a velocity of 8.32 km/s was found underneath the Mohorovic̆ić discontinuity. An upper-mantle fine structure consisting of several high- and low-velocity zones was derived. At a depth of 96 km, a P-velocity of 8.5 km/s was reached.A minimum shear-wave velocity of 4.70 km/s has been determined for the uppermost part of the mantle underneath the Mohorovičić discontinuity. This velocity allows a shear low-velocity zone of a few tens of kilometers in the uppermost post of the mantle.     

The mechanics of a magma chamber-fault system in trans-tension with application to Coso

The mechanical interaction between an elliptically shaped magma chamber and a fault subject to transtension is investigated with particular reference to the Coso geothermal field. The geologic setting of the Coso field is interpreted as a releasing bend step-over structure formed by the Airport Lake and Owens Valley dextral strike-slip fault system. The role of the Coso volcano-magmatic center in the development of the “over-step” structure is examined by treating the magma chamber as a liquid inclusion in a viscoelastic crust containing a fault (Airport Lake). The problem is numerically solved using a 2D viscoelastic finite element model with thermally activated viscosity to account for thermal weakening of the rock. The temperature distribution around the magma body is calculated based on a 3D steady-state approach and using the mesh-less numerical method. The fault is modeled as a frictionless contact. The simulated distributions of stress and strain around the inclusion display a rotation caused by the shearing component of the applied transtension. The results indicate that the fault tends to overstep the chamber in a geometric pattern similar to a step-over. There is good correspondence between the computed distributions of the maximum shear stress in the vicinity of the magma chamber and the map of earthquake epicenters at a depth of 7–10 km in Coso.     

Deep structure under Yellowstone National Park U.S.A.: A continental “hot spot”

In order to understand the origin of long-lived loci of volcanism (sometimes called “hot spots”) and their possible role in global tectonic processes, it is essential to know their deep structure. Even though some work has been done on the crustal, upper-mantle, and deep-mantle structure under some of these “hot spots”, the picture is far from clear. In an attempt to study the structure under the Yellowstone National Park U.S.A., which is considered to be such a “hot spot”, we recorded teleseisms using 26 telemetered seismic stations and three groups of portable stations. The network was operated within a 150 km radius centered on the Yellowstone caldera, the major, Quaternary volcanic feature of the Yellowstone region. Teleseismic delays of about 1.5 sec are found inside the caldera, and the delays remain high over a 100 km wide area around the caldera. The spatial distribution and magnitude of the delays indicate the presence of a large body of low-velocity material with horizontal dimensions corresponding approximately to the caldera size (40 km × 80 km) near the surface and extending to a depth of 200–250 km under the caldera. Using ray-tracing and inversion techniques, it is estimated that the compressional velocity inside the anomalous body is lower than in the surrounding rock by about 15% in the upper crust and by 5% in the lower crust and upper mantle. It is postulated that the body is partly composed of molten rock with a high degree of partial melting at shallow depths and is responsible for the observed Yellowstone volcanism. The large size of the partially molten body, taken together with its location at the head of a 350 km zone of volcanic propagation along the axis of the Snake River Plain, indicates that the volcanism associated with Yellowstone has its origin below the lithosphere and is relatively stationary with respect to plate motion. Using our techniques, we are unable to detect any measurable velocity contrast in the mantle beneath the low-velocity body, and, hence, we are unable to determine whether the Yellowstone melting anomaly is associated with a deep heat source or with any deep phenomenon such as a convection plume, chemical plume, or gravitational anchor.     

Crustal structure from deep seismic soundings along the Koyna II (Kelsi-Loni) profile in the Deccan Trap area, India

The crustal depth section obtained from deep seismic soundings along the Koyna II (Kelsi-Loni) profile, which lies near latitude 18°N roughly in the east-west direction in that part of the Deccan Trap Maharashtra State, India, shows a number of reflection segments below the Deccan Traps down to the Moho discontinuity. A deep fault below the Deccan Traps 13 km east of Mahad divides the entire cross-section including the Moho boundary into two crustal blocks. The reflection segments show updip towards the west coast in the western block. The Moho discontinuity which is at a depth of 39 km near the deep fault starts rising towards the coast, reaching a depth of 31.5 km at the west coast. The eastern block is thrown up by 1.5 km with respect to the western block along the deep fault. A structural contour map of the Moho discontinuity for the Koyna reservoir area has been prepared from the present results and the crustal information obtained along the Koyna I profile (Kaila et al., 1979a), shows that the deep fault in the Koyna area is aligned in the NNW-SSE direction.Refraction seismic data analysis by the wave front method reveals that the thickness of the Deccan Trap increases towards the west coast. The Deccan Trap is 600–700 m thick in the eastern region between Nira (SP 130) and Loni (SP 200) and attains a thickness of 1500 m at 10 km east of the west coast. The longitudinal wave velocity in the Deccan Traps along the profile varies from 4.8 to 5.0 km/sec and in the crystalline basement from 6.0 to 6.15 km/sec. A tentative isopach contour map of the Deccan Traps and a tentative structural contour map of the Pre-Deccan Trap contact have been prepared for the Koyna reservoir area from the results along the Koyna II and Koyna I profiles. A flexure aligned in a NNW-SSE direction, in the Pre-Deccan Trap contact, which is an expression of the deep fault into the basement, has been clearly brought out. The flexure coincides in general with the orientation of the Deccan volcanic scarp in this area.     

A model for the deep structure of the East African rift system from simultaneous inversion of teleseismic data

Simultaneous inversion of teleseismic data results in a model for the western branch of the East African Rift system between 1°S and 10°N characterized by a 35 km crust and thin high-velocity lid, overlying a channel possessing both low S and low P velocities (4.47 and 7.69 km/sec, respectively). A strong reflector at 140 km depth marks the boundary to high-velocity material (4.67 and 8.44 km/sec, respectively).The eastern branch has a crustal thickness of 40 km and is characterized by low S velocities, 4.43 km/sec in the lid to a depth of 78 km, 4.09–4.21 km/sec in a channel which extends to 161 km depth. The S velocities remain relatively low at greater depths, but cannot be precisely determined because of the limited resolution.We have interpreted these and other geophysical data as due to a diapiric intrusion of material with a high pyroxene content and possibly with some low-density eclogite, coming from the mesosphere. A mushroom of this rock has stalled below the western branch at about 55 km below the surface and is cooling. The diapirism in the eastern branch has broken through to the surface. The mantle of the African shield to the west of the Rift shows seismic velocities which indicate that it is depleted in basaltic components.     

Application of high-resolution seismic techniques in the evaluation of earthquake site response, Ottawa Valley, Canada

High-resolution seismic surveys, including P- and S-wave studies, have been conducted in an area of the Ottawa River valley located 80 km east of Ottawa (Canada). Based on dating of paleolandslides, the existence of paleoearthquake activity has been postulated in this area. The target zone for the seismic survey is characterized by surface disturbance and sediment deformation. P-wave seismic imaging was used to map the overburden–bedrock interface as well as to indicate reflecting boundaries within the overburden. The area of surface disturbance was found to overlie a buried bedrock basin, 8 km in diameter, infilled with a maximum thickness of 180 m of unconsolidated Quaternary sediments. Preliminary results of core logging show the presence of sand overlain by deformed fine sediments within the disturbed area. Shear-refraction studies reveal differences in the velocity–depth profiles between the disturbed area and the surrounding undisturbed areas. The shear-wave reflection method was used to produce a fundamental resonant period map for the area. Surface sediment disturbance was probably due to a combination of ground-motion amplification due to the basin (thick soft sediments) and the presence of water-saturated sand at depth.     

The state of the upper mantle beneath southern Africa

We present a new upper mantle seismic model for southern Africa based on the fitting of a large (3622 waveforms) multi-mode surface wave data set with propagation paths significantly shorter (≤ 6000 km) than those in globally-derived surface wave models. The seismic lithosphere beneath the cratonic region of southern Africa in this model is about 175 ± 25 km thick, consistent with other recent surface wave models, but significantly thinner than indicated by teleseismic body-wave tomography. We determine the in situ geotherm from kimberlite nodules from beneath the same region and find the thermal lithosphere model that best fits the nodule data has a mechanical boundary layer thickness of 186 km and a thermal lithosphere thickness of 204 km, in very good agreement with the seismic measurement. The shear wave velocity determined from analyzes of the kimberlite nodule compositions agree with the seismic shear wave velocity to a depth of 150 km. However, the shear wave velocity decrease at the base of the lid seen in the seismic model does not correspond to a change in mineralogy. Recent experimental studies of the shear wave velocity in olivine as a function of temperature and period of oscillation demonstrate that this wave speed decrease can result from grain boundary relaxation at high temperatures at the period of seismic waves. This decrease in velocity occurs where the mantle temperature is close to the melting temperature (within 100 °C).     

Crustal seismic velocity structure near Faial and Pico Islands (AZORES), from local earthquake tomography

The Azores archipelago (Portugal) is located on an oceanic plateau, in a geodynamic environment prone to intense seismic and volcanic activity. In order to investigate the crustal structure in this region, we have conducted a local earthquake tomography study in the area of the islands of Faial, Pico and S. Jorge using data recorded in July 1998. The July 9th 1998 earthquake, near Faial Island, triggered an aftershock sequence of thousands of events that lasted for several months and were recorded by a total of 14 stations located on the three islands surrounding the epicentral area. In the upper crustal layers, consistency is seen between the tomographic results and the islands' surface volcanic units. Beneath the Faial central volcano a low Vp (< 6.0 km/s) anomaly roughly located at 3–7 km depth, suggests a connection to the plumbing system, possibly the presence of a magma chamber. In NE Faial, a high Vp (> 6.3 km/s) body was found located at mid-lower crust, most likely an intrusion of gabbroic composition, that is bordered by the registered seismic activity; its shape suggesting a tectonic controlled mechanism. The relocated hypocenters, together with the overall analysis of the Tomographic model, suggest a tectonic segmentation of Faial Island. The crustal thickness under the islands volcanic buildings of the Faial–Pico area was estimated at around 14 km.     

Crustal structure of mainland China from deep seismic sounding data

Since 1958, about ninety seismic refraction/wide angle reflection profiles, with a cumulative length of more than sixty thousand kilometers, have been completed in mainland China. We summarize the results in the form of (1) a new contour map of crustal thickness, (2) fourteen representative crustal seismic velocity–depth columns for various tectonic units, and, (3) a Pn velocity map. We found a north–south-trending belt with a strong lateral gradient in crustal thickness in central China. This belt divides China into an eastern region, with a crustal thickness of 30–45 km, and a western region, with a thickness of 45–75 km. The crust in these two regions has experienced different evolutionary processes, and currently lies within distinct tectonic stress fields. Our compilation finds that there is a high-velocity (7.1–7.4 km/s) layer in the lower crust of the stable Tarim basin and Ordos plateau. However, in young orogenic belts, including parts of eastern China, the Tianshan and the Tibetan plateau, this layer is often absent. One exception is southern Tibet, where the presence of a high-velocity layer is related to the northward injection of the cold Indian plate. This high-velocity layer is absent in northern Tibet. In orogenic belts, there usually is a low-velocity layer (LVL) in the crust, but in stable regions this layer seldom exists. The Pn velocities in eastern China generally range from 7.9 to 8.1 km/s and tend to be isotropic. Pn velocities in western China are more variable, ranging from 7.7 to 8.2 km/s, and may display azimuthal anisotropy.     

Landslide susceptibility mapping of the slopes in the residual soils of the Mengen region (Turkey) by deterministic stability analyses and image processing techniques

The aim of the present study is to prepare a landslide susceptibility map of a region of about 120 km², between Gökcesu and Pazarköy (around Mengen, NW Turkey) at approximately 10 km north of the North Anatolian Fault Zone, where frequent landslides occur. For this purpose, mechanisms of the landslides were studied by two-dimensional stability analyses together with field observations, and the parameters controlling the development of such slides were identified. Field observations indicated that the failures generally developed within the unconsolidated and/or semiconsolidated soil units in forms of rotational, successive shallow landslides within the weathered zone in Mengen, Cukurca and Sazlar formations. Although consisting of residual soils, Capak and Gökdag formations do not exhibit landslides as the natural slopes formed on these, do not exceed the critical slope angles. Statistical evaluations and distribution of the landslides on the topographical map showed that such parameters as cohesion, angle of internal friction, slope, relative height, orientation of slopes, proximity to drainage pattern, vegetation cover and proximity to major faults were the common features on the landslides. Digital images were obtained to represent all these parameters on gray scale on the SPOT image and on the digital elevation model (DEM) of the area using image processing techniques. Soil mechanics tests were carried out on 36 representative samples collected from different units, and parameters were determined for two-dimensional stability analyses basing on “sensitivity approach” and for the preparation of digital shear strength map. In order to determine the critical slope angle values for the residual soils, a series of sensitivity analyses were realized by using two-dimensional deterministic slope stability analyses techniques for varying values of cohesion, angle of internal friction and slope height along with varying saturation conditions. According to the results of the sensitivity analyses, the Mengen formation was found to be most susceptible unit to landslides, covering about 33.5% of the region studied in terms of surface area. The distribution of the critical slopes were determined by superimposing the critical slope values from sensitivity analyses on slope map of the study area. On the other hand, iso-cohesion and iso-friction maps were produced by locating the values of cohesion and internal friction angles in a geographic coordinate system such that they coincide with sample locations on the DEM and by further interpolation of the values concerned. The pixel values were evaluated in gray scale from 0 to 255, 0 representing the lowest pixel value and 255 representing the highest. Sensitivity analyses on cohesion and angle of internal friction to investigate the effects of these parameters only on stability, revealed that cohesion was effective at a rate of 70% by itself, while angle of internal friction alone controlled the stability by a rate of 30%. The iso-cohesion and iso-friction maps previously obtained were digitally combined in these rates and a “shear strength map” was prepared. The geographic setting of the study area is such that northern slopes usually receive dense precipitation. In relation to this fact, about 42% of the landslides are due north. Thus, a slope orientation map was prepared using the DEM, and slopes facing north were evaluated as being more susceptible to sliding. Proximity to the drainage pattern was another important factor in the evaluation, as streams could adversely affect the stability by either eroding the toe or saturating the slope, or both. When considered together, in conjunction with the field observations, faults and landslides showed a close association. In the area, about 88% of the landslides were detected within an area closer than 250 m to major faults, therefore, a main discontinuity map was produced using the SPOT image of the region, and “proximity to major faults” was evaluated as a parameter as most of the landslides developed in areas where the vegetation was rather sparse. A vegetation cover map was therefore obtained from the SPOT image, and the areas with denser vegetation were considered to be less susceptible to sliding with respect to the areas with less or no vegetation. Having prepared the maps accounting for the distribution of critical slopes, shear strength properties, relative height, slope angle, orientation of the slopes, vegetation cover, proximity to the drainage pattern, geographic corrections were carried on each of these, and a potential failure map was obtained for the residual soils by superimposing all these maps. Next, a classification was performed on the final map and five relative zones of susceptibility were defined. When compared with this map, all of the landslides identified in the field were found to be located in the most susceptible zone. The performance of the method used in processing the images appears to be quite high, the zones determined on the map being the zones of relative susceptibility.     

Some observations on the geochemistry of gold mineralization in the weathered zone at Norseman, Western Australia

At Norseman Au mineralization is hosted by parallel easterly dipping quartz reefs in a westerly dipping sequence of Archaean basaltic flows and As is strongly correlated with this mineralization. Near-surface exploration is hampered by an extensive cover of deeply weathered soils and a veneer of silty calcareous material, probably of aeolian origin, which, because of its low Fe content of < 5%, has a very low (< 10 ppm) and often indiscernible As signature over mineralization. At a test traverse over the southern end of the Mararoa quartz reef, drill sampling to 10 m depth in highly weathered in situ material beneath the calcareous veneer produced As anomalies which are strongly correlated with the extrapolation to surface of a shear zone hosting Au mineralization which is probably blind. The As anomalies can be extended into the surface horizon by manipulating the data to compensate for the low Fe content of the calcareous layer.The Au values in the drill samples were very low and poorly correlated with the extrapolation to surface of the Mararoa shear. However, the absence of a Au halo (at a 10 ppb detection limit) may be due to the shear being barren near surface, or may be due to the samples being taken from the strongly leached zones of a laterite profile. In the area of the test traverse at the southern Mararoa area, the upper horizons of the original laterite profile have been eroded away and the remaining horizons covered by a veneer of calcareous material. There is a better correlation of the Au values of the shallow drill samples with the projection to surface of a porphyry, which hosts minor low-grade Au mineralization.The results for Au dispersion around the Mararoa shear contrast with dispersion at the nearby Hinemoa Mine where gold mineralization crops out. Channel sampling in costeans showed that an Au anomaly (threshold 40 ppb) extends 50 to 150 m into oxidized wall rocks. This well defined halo may be due to the secondary dispersion of Au from the quartz reef into the ferruginous zone of the laterite profile.Gold in the leaves and twigs of Melaleuca pauperiflora F. Muell. from the traverse at the southern end of the Mararoa shear was determined by neutron activation analysis and found to be broadly correlated with the projected Au mineralization. However, the lack of correlation with the soil Au values to 10 m may be due to the plants taking up most of their water from below the depth of sampling.     

5: Late Variscan mineralizing systems related to orogenic processes: The French Massif Central

The French Massif Central constitutes an exceptional study area due to the diversity of its metallic deposits, its internal position in the Variscan belt, and the abundance of available geological, geophysical and metallogenic data obtained within the GeoFrance 3D programme. The deposits, formed towards the end of the orogenic evolution, represent the economic products of two distinct mineralizing systems, a Au ± Sb hydrothermal system and a W ± Sn and rare-metals magmatic–hydrothermal system, which were simultaneously active during a short time span between ca. 310 and 300 Ma.Two types of gold deposit can be distinguished on the basis of their depth of emplacement: “deep-seated” gold deposits developed under lithostatic to hydrostatic pressure during rapid exhumation, and “shallow” gold deposits emplaced under hydrostatic pressure with no significant uplift.Deposits of W ± Sn and rare-metals were emplaced in the upper crust during final crystallization of specialized magmas after their rapid ascent, perhaps enhanced by simultaneous regional uplift. The gold-bearing systems are associated with a complex network of re-activated crustal-scale faults initially active during the period between 335 and 315 Ma. Normal motion along the faults, coeval with 335 to 315 Ma granite–migmatite domes, played a major role in the 3D distribution of the hydrothermal plumbing system. Gold and related metals were carried within huge hydrothermal cells, which reached ca. 100 km by 10 km in area, and 30 km in depth. In contrast, granites rich in magmatophile elements (W, Sn, rare-metals) generated smaller hydrothermal cells (10 km by 10 km in area, and < 6 km deep). Extraction of metals, by both deep-seated fluids and specialized magmas, occurred during granulitization of the lower crust at 300 ± 15 Ma. In the French Massif Central, the genesis of the two late Carboniferous mineralizing systems coincided with the end of syn-collisional extension and ended just before post-collisional extension.     

The upper crust low velocity layer; a Rayleigh (Rg) phase velocity study from SE Norway

A recent surge of interest in short-period fundamental mode Rayleigh waves (Rg) propagation is motivated by their usefulness for structural mapping and estimation of attenuation (Q) properties of the uppermost crust, their importance in scattering studies and as a nuisance (false alarms) in detection seismology. We have investigated Rg phase velocity dispersion characteristics, using recordings from 16 out of 22 subarrays of the NORSAR array of aperture 100 km. The observed phase velocity in combination with a one-layer and two-layers model over halfspace models were inverted using a maximum likelihood scheme. In view of the smoothness and simplicity of the dispersion curves, the number of unknowns were limited to three parameters, namely shear velocities and thickness, using a Poisson's ratio of 0.25 for P-S velocity coupling and density fixed. In the general array siting area the derived low velocity layer has an average thickness of 1.02 km, and the associated shear velocity is 2.82 km s^-1. The underlying model half space appears to be uniform, and its derived shear velocity is 3.55 km s^-1, which is typical of the upper crust. In contrast to other similar studies, no obvious correlation was found with local geology, which mainly consists of Precambrian crystalline rock sequences. Group velocities were calculated from the structural models and compared to the observational ones. The agreement between calculated and observed group velocities was reasonable.     

Structure of the upper mantle in the Northeastern Atlantic close to the azores-gibraltar ridge from surface-wave and body-wave observations

Phase velocities of teleseismic Rayleigh waves have been measured in the central North Atlantic on both sides of the Azores-Gibraltar Ridge (AGR) by means of a specially designed long-period station network. The dispersion data obtained were regionalized and then subjected to a “hedgehog” inversion, which gives a set of upper mantle models compatible with the observational data within specified error bounds.Reasonable model solutions were selected by using regional body-wave observations, such as P_n- and S_n-wave velocities determined from earthquakes along the AGR. The S_(itn) velocities measured indicate that the shear-wave velocity in the mantle part of the lithosphere is much higher on the northern side of the AGR. Strongly negative P-wave residuals in this area indicate faster seismic propagation than implied by the Jeffreys-Bullen travel-time tables, while propagation is much slower in the Gulf of Cadiz area. Furthermore the residuals show a clear difference for paths through oceanic and continental domains and suggest that the transition between these two domains extends much further into the ocean on the southern side of the AGR than on the northern side.The proposed model for the structure of the upper mantle in that region shows that there exists a pronounced velocity contrast across the AGR. Thickening of the lithospheric plate with increasing plate age is indicated to the south of the ridge. The greatest thickness is reached close to the continental margin within a zone about 500 km wide, whose velocity close to the Canary Islands and Madeira is significantly lower, probably due to the well-known volcanic activity there. These observations together with the travel time residuals reveal that this zone seems to be of a transitional nature somewhere between a continental and oceanic structure.