major types and time space distribution of mesozoic ore deposits in south china and their geodynamic settings the ore deposits of the mesozoic age in south china can be divided into three groups each with different metal associations and spatial distributions and each related to major magmatic events the first event occurred in the late triassic 230 210 ma the second in the mid late jurassic 170 150 ma and the third in the early mid cretaceous 120 80 ma the late triassic magmatic event and associated mineralization is characterized by peraluminous granite related w sn nb ta mineral deposits the triassic ore deposits are considerably disturbed or overprinted by the later jurassic and cretaceous tectono thermal episodes the mid late jurassic magmatic and mineralization events consist of 170 160 ma porphyry skarn cu and pb zn ag vein deposits associated with i type granites and 160 150 ma metaluminous granite related polymetallic w sn deposits the late jurassic metaluminous granite related w sn deposits occur in a ne trending cluster in the interior of south china such as in the nanling area in the early mid cretaceous from about 120 to 80 ma but peaking at 100 90 ma subvolcanic related fe deposits developed and i type calc alkaline granitic intrusions formed porphyry cu mo and porphyry epithermal cu au ag mineral systems whereas s type peraluminous and or metaluminous granitic intrusions formed polymetallic sn deposits these cretaceous mineral deposits cluster in distinct areas and are controlled by pull apart basins along the south china continental margin based on mineral assemblage age and space time distribution of these mineral systems integrated with regional geological data and field observations we suggest that the three magmatic mineralization episodes are the result of distinct geodynamic regimes the triassic peraluminous granites and associated w sn nb ta mineralization formed during post collisional processes involving the south china block the north china craton and the indo china block mostly along the dabie sulu and songma sutures jurassic events were initially related to the shallow oblique subduction of the izanagi plate beneath the eurasian continent at about 175 ma but i type granitoids with porphyry cu and vein type pb zn ag deposits only began to form as a result of the breakup of the subducted plate at 170 160 ma along the nne trending qinzhou hangzhou belt also referred to as qin hang or shi hang belt which is the neoproterozoic suture that amalgamates the yangtze craton and cathaysia block a large subduction slab window is assumed to have formed in the nanling and adjacent areas in the interior of south china triggering the uprise of asthenospheric mantle into the upper crust and leading to the emplacement of metaluminous granitic magma and associated polymetallic w sn mineralization a relatively tectonically quiet period followed between 150 and 135 ma in south china from 135 ma onward the angle of convergence of the izanagi plate changed from oblique to parallel to the coastline resulting in continental extensional tectonics and reactivation of regional scale ne trending faults such as the tan lu fault this widespread extension also promoted the development of ne trending pull apart basins and metamorphic core complexes accompanied by volcanism and the formation of epithermal cu au deposits granite related polymetallic sn w deposits and hydrothermal u deposits between 120 and 80 ma with a peak activity at 100 90 ma

a late cretaceous tin metallogenic event in nanling w sn metallogenic province constraints from u pb ar ar geochronology at the jiepailing sn be f deposit hunan china abstract the jiepailing deposit located in southern hunan province china is a giant sn be f deposit in the nanling w sn province the sn be f mineralization is spatially associated with the jiepailing granite porphyry la mc icp ms zircon u pb dating of the jiepailing granite porphyry yielded a weighted mean 206 pb 238 u age of 90 5 0 9 ma mswd 0 32 which is interpreted as the emplacement age of the granite porphyry hydrothermal muscovite yields a plateau 40 ar 39 ar age of 92 1 0 7 ma mswd 0 9 which is well consistent with the zircon u pb age of the jiepailing granite porphyry responsible for the sn be f mineralization indicating a temporal link between the emplacement of the jiepailing granite porphyry and the sn be f mineralization our new high precise geochronological data suggest that the jiepailing giant sn be f deposit and related granite formed during the early late cretaceous 92 90 ma which provided convincing evidence for a previously unrecognized metallogenic event related to late cretaceous granitic magmatism in nanling region the occurrence of the late cretaceous sn metallogenic event identified in southern hunan province further highlights the importance of systematic metallogenic studies of the nanling w sn province integrated with high precise geochronological data obtained previously it is suggested that the nanling region mainly experienced three tectonomagmatic activities and lithospheric thinning events during mesozoic which are responsible for three sn w metallogenic events to form the world class nanling w sn metallogenic province

crustal concealing of small scale core field secular variation summary the earth s magnetic field is mainly produced within the earth s liquid and electrically conducting core as a result of a process known as the geodynamo many other sources also contribute to the magnetic signal accessible to observation at the earth s surface partly obscuring the main core magnetic field signal thanks to a series of very successful satellites and to advances in magnetic field modelling techniques considerable progress has however been made in the recent years toward better identifying the signal of each of these sources in particular temporal changes in the field of internal origin happen to be detectable now in spherical harmonic degrees up to perhaps 16 all of these changes are usually attributed to changes in the core field itself the secular variation on the ground that the lithospheric magnetization cannot produce such signals it has however been pointed out on empirical grounds that temporal changes in the field of internal origin produced by the induced part of the lithospheric magnetization could dominate the core field signal beyond degree 22 this short note revisits this issue by taking advantage of our improved knowledge of the small scale field changes and of the likely sources of the lithospheric field we rely on a simple extrapolation of the observed spatial spectrum of the field changes beyond degree 16 and use a forward approach based on a recent geological model of lithospheric magnetization this leads us to confirm that the main cause of the observed changes in the field of internal origin up to some critical degree n c is indeed likely to be the secular variation of the core field but that the signal produced by the time varying lithospheric field is bound to dominate and conceal the time varying core signal beyond that critical degree in very much the same way the permanent component of the lithospheric field dominates and conceals the permanent component of the core field beyond degree 14 all uncertainties taken into account we estimate n c to lie between 22 and 24 we however also note that in practice the main limitation to the observation of the core field small scale secular variation is not so much its concealing by the field of lithospheric origin but its fast changing nature and small magnitude this leads us to conclude that whereas cumulative small scale lithospheric field changes might be detected some day detection of core field secular variation beyond degree 18 is likely to remain a severe challenge for some time

absolute timing of sulfide and gold mineralization a comparison of re os molybdenite and ar ar mica methods from the tintina gold belt alaska new re os molybdenite dates from two lode gold deposits of the tintina gold belt alaska provide direct timing constraints for sulfide and gold mineralization at fort knox the re os molybdenite date is identical to the u pb zircon age for the host intrusion supporting an intrusive related origin for the deposit however 40ar 39ar dates from hydrothermal and igneous mica are considerably younger at the pogo deposit re os molybdenite dates are also much older than 40ar 39ar dates from hydrothermal mica but dissimilar to the age of local granites these age relationships indicate that the re os molybdenite method records the timing of sulfide and gold mineralization whereas much younger 40ar 39ar dates are affected by post ore thermal events slow cooling and or systemic analytical effects the results of this study complement a growing body of evidence to indicate that the re os chronometer in molybdenite can be an accurate and robust tool for establishing timing relations in ore systems

satellite magnetic anomalies related to seafloor spreading in the south atlantic ocean oceanic magnetic anomalies have been observed from satellite the strongest anomalies are the long wavelength components of the sea floor spreading signature unfortunately because of technical issues involving the treatment of satellite magnetic data these signals are obscured in the south atlantic ocean because they trend north south however a map does exist in which such features are observed essentially because of a better data processing technique further this map agrees with a physically motivated model based on non satellite magnetic input hence with properly treated data the magnetic anomaly maps should be useful for saying something about the geology rather than vice versa this situation will be considerably improved by ongoing advances in methods and new data sets the amplitude of the observations a factor of two larger than previous estimates confirms that the extrusive basaltic layer alone is inadequate to produce the signal and that deeper oceanic sources are required

scalar magnetic anomaly maps of earth derived from pogo and magsat data a new pogo scalar magnetic anomaly map at 400 km altitude is presented which consists of spherical harmonics of degree 15 60 on the basis of the common features of this map with two new magsat anomaly maps dawn and dusk two scalar magnetic anomaly maps of the earth are presented using two selection criteria with different levels of stringency these selection criteria suppress the noncrustal components of the original maps by different amounts the more stringent selection criteria seek to eliminate as much contamination as possible at the expense of suppressing some anomaly signal this map is represented by spherical harmonics of degree 15 60 the less stringent selection criteria seek to retain as much crustal signal as possible at the expense of also retaining some contaminating fields this map is represented by spherical harmonics of degree 15 65 the resulting two maps are highly correlated with degree correlation coefficients greater than 0 8

nontrivial decay of aftershock density with distance in southern california the decay of the aftershock density with distance plays an important role in the discussion of the dominant underlying cause of earthquake triggering here we provide evidence that its form is more complicated than typically assumed and that in particular a transition in the power law decay occurs at length scales comparable to the thickness of the crust this is supported by an analysis of a very recent high resolution catalog for southern california sc and surrogate catalogs generated by the epidemic type aftershock sequence etas model which take into account inhomogeneous background activity short term aftershock incompleteness anisotropic triggering and variations in the observational magnitude threshold our findings indicate specifically that the asymptotic decay in the aftershock density with distance is characterized by an exponent larger than 2 which is much bigger than the observed exponent of approximately 1 35 observed for shorter distances ranging from the main shock rupture length up to a length scale comparable to the thickness of the crust this has also important consequences for time dependent seismic hazard assessment based on the etas model

contribution of radar interferometry to a two step inversion of the kinematic process of the 1992 landers earthquake we study the rupture process of the 1992 landers earthquake to limit the trade off between slip amplitude and rupture time that affects solutions using only seismological data we adopt a two step approach we first constrain the slip distribution and its uncertainty by independent geodetic data to recover in the second step the temporal details of the rupture propagation the first step consists of an inversion of interferometric data and global positioning system measurements both independently and together to constrain slip distribution on a three segment fault model along both strike and dip direction we use a genetic algorithm to test the uniqueness of the solution and a least squares formulation to find the model which best fits the data we conclude from the results of these geodetic inversions that interferometric data are rich enough to access the slip distribution in the case of the landers earthquake since the surface deformations are more sensitive to shallow slip in our configuration the slip amplitude is better resolved near the surface than at depth the resulting slip distribution is in agreement with geological observations at the surface and confirms the heterogeneous nature of the landers earthquake most of the slip occurs at shallow depths on the homestead valley fault second segment with a maximum value of around 7 m another high slip zone is observed on the johnson valley fault first segment at 8 km depth in the second step we invert strong motion data with the a priori final slip amplitude and its uncertainty deduced from geodetic data to constrain the time history of the rupture process this second step emphasizes a strong variation of the temporal development of the earthquake fast rupture front velocities appear within high slip zones and the rupture slows when it encounters a resistance along the fault on average the rupture front propagates with velocities close to the s wave velocity and terminates about 20 s after initiation the large variations in both slip amplitude and rupture velocity suggest that the rupture process is better described by successively breaking asperities than by a pulse propagating with constant velocity

in situ la mc icp ms and id tims u pb geochronology of cassiterite in the giant furong tin deposit hunan province south china new constraints on the timing of tin polymetallic mineralization abstract the furong deposit located in southern hunan province china is a newly discovered giant tin deposit spatially associated with the a type qitianling granite the petrogenetic link between the giant furong tin deposit and the qitianling granite batholith remains controversial because of the lack of precise dating for the tin polymetallic mineralization here we report for the first time in situ u pb data on cassiterite obtained by la mc icp ms and the results are compared to id tims data the in situ la mc icp ms analyses of the cassiterites provide a reliable age constraint for tin polymetallic mineralization in the furong deposit yielding an isochron age of 159 9 1 9 ma at 95 confidence level mswd 18 which is indistinguishable from the id tims 206 pb 238 u weighted mean age of 158 2 0 4 ma the in situ u pb ages corroborate the ar ar dates published for the furong tin deposit which indicate that the mineralization is coeval with the emplacement of the qitianling granite batholith these results provide further evidence of a close temporal link between the pluton s emplacement and tin mineralization in the furong area the data obtained clearly demonstrate that in situ la mc icp ms u pb dating of cassiterite is a robust geochronometer for direct dating of tin polymetallic mineralization

the stress shadow problem in physics based aftershock forecasting does incorporation of secondary stress changes help main shocks are calculated to cast stress shadows across broad areas where aftershocks occur thus a key problem with stress based operational forecasts is that they can badly underestimate aftershock occurrence in the shadows we examine the performance of two physics based earthquake forecast models coulomb rate state crs based on coulomb stress changes and a rate and state friction law for their predictive power on the 1989 mw 6 9 loma prieta aftershock sequence the crs 1 model considers the stress perturbations associated with the main shock rupture only whereas crs 2 uses an updated stress field with stresses imparted by m 3 5 aftershocks including secondary triggering effects slightly improves predictability but physics based models still underestimate aftershock rates in locations of initial negative stress changes furthermore crs 2 does not explain aftershock occurrence where secondary stress changes enhance the initial stress shadow predicting earthquake occurrence in calculated stress shadow zones remains a challenge for stress based forecasts and additional triggering mechanisms must be invoked

geological modeling of the new champ magnetic anomaly maps using a geographical information system technique 1 reliable global crustal field anomaly maps produced from magnetic measurements of the champ satellite mission now allow for quantitative geological studies of crustal structure and composition we have developed a gis based forward modeling technique to model these anomaly maps on the basis of the geologic and tectonic maps of the world laboratory susceptibility values of the occurring rock types and the seismic thickness of the crust a vertically integrated susceptibility grid is generated in the gis system in addition a remanent magnetization grid is computed for the oceanic crust using a digital isochron map of the ocean floor and rotation models of the paleoplates combining the global vis and remanent magnetization grids the vertical magnetic field anomaly is computed at satellite altitude and compared with the corresponding champ magnetic anomaly map over the oceans induced and remanent magnetization explains well the prominent observed anomalies over the cretaceous quiet zones we also find a good agreement between predicted and observed anomalies over the continents remaining discrepancies between the predicted and observed anomalies can be used to adjust poorly known boundaries and the composition of the buried precambrian provinces until the recomputed anomalies fit the observed anomalies the feasibility of this approach is demonstrated on greenland the west african craton bangui in central africa and the kolyma omolon block in siberia we conclude that quantitative information on the lateral extent the composition and the thickness of the lower crust within a precambrian province can thus be inferred from the new satellite magnetic anomaly maps

new insights on intraplate volcanism in french polynesia from wavelet analysis of grace champ and sea surface data in this paper we discuss the origin of superswell volcanism on the basis of representation and analysis of recent gravity and magnetic satellite data with wavelets in spherical geometry we computed a refined gravity field in the south central pacific based on the grace satellite ggm02s global gravity field and the kms02 altimetric grid and a magnetic anomaly field based on champ data the magnetic anomalies are marked by the magnetic lineation of the seafloor spreading and by a strong anomaly in the tuamotu region which we interpret as evidence for crustal thickening we interpret our gravity field through a continuous wavelet analysis that allows to get a first idea of the internal density distribution we also compute the continuous wavelet analysis of the bathymetric contribution to discriminate between deep and superficial sources according to the gravity signature of the different chains as revealed by our analysis various processes are at the origin of the volcanism in french polynesia as evidence we show a large scale anomaly over the society islands that we interpret as the gravity signature of a deeply anchored mantle plume the gravity signature of the cook austral chain indicates a complex origin which may involve deep processes finally we discuss the particular location of the marquesas chain as suggesting that the origin of the volcanism may interfere with secondary convection rolls or may be controlled by lithospheric weakness due to the regional stress field or else related to the presence of the nearby tuamotu plateau

new empirical relationships among magnitude rupture length rupture width rupture area and surface displacement source parameters for historical earthquakes worldwide are compiled to develop a series of empirical relationships among moment magnitude m surface rupture length subsurface rupture length downdip rupture width rupture area and maximum and average displacement per event the resulting data base is a significant update of previous compilations and includes the additional source parameters of seismic moment moment magnitude subsurface rupture length downdip rupture width and average surface displacement each source parameter is classified as reliable or unreliable based on our evaluation of the accuracy of individual values only the reliable source parameters are used in the final analyses in comparing source parameters we note the following trends 1 generally the length of rupture at the surface is equal to 75 of the subsurface rupture length however the ratio of surface rupture length to subsurface rupture length increases with magnitude 2 the average surface displacement per event is about one half the maximum surface displacement per event and 3 the average subsurface displacement on the fault plane is less than the maximum surface displacement but more than the average surface displacement thus for most earthquakes in this data base slip on the fault plane at seismogenic depths is manifested by similar displacements at the surface log linear regressions between earthquake magnitude and surface rupture length subsurface rupture length and rupture area are especially well correlated showing standard deviations of 0 25 to 0 35 magnitude units most relationships are not statistically different at a 95 significance level as a function of the style of faulting thus we consider the regressions for all slip types to be appropriate for most applications regressions between magnitude and displacement magnitude and rupture width and between displacement and rupture length are less well correlated and have larger standard deviation than regressions between magnitude and length or area the large number of data points in most of these regressions and their statistical stability suggest that they are unlikely to change significantly in response to additional data separating the data according to extensional and compressional tectonic environments neither provides statistically different results nor improves the statistical significance of the regressions regressions for cases in which earthquake magnitude is either the independent or the dependent parameter can be used to estimate maximum earthquake magnitudes both for surface faults and for subsurface seismic sources such as blind faults and to estimate the expected surface displacement along a fault for a given size earthquake

coseismic displacements 1992 landers california earthquake the rupture associated with the 1992 landers earthquake ms 7 5 extended along the west edge of a 60 50 km trilateration network that had been surveyed two months before the earthquake and was resurveyed a month after it the coseismic displacement field within the network was calculated from the changes observed in the trilateration network supplemented by the displacements measured at four gps stations and other less timely trilateration measurements stations near the rupture were displaced by as much as 3 m the observed changes are reproduced reasonably well by dislocation modeling in which the landers rupture is represented by 8 vertical fault segments each extending from the surface to 9 km depth two additional segments represent the big bear and eureka peak faults the average modeled right lateral slip on the landers rupture is about 5 m however the residuals for the best fit slip distribution are about four times greater than expected from measurement error given that uncertainty a broad range of slip models furnish acceptable fits to the data the large residuals are attributed to inelastic deformation principally sympathetic slip on numerous unmodeled faults within the fault blocks assumed to be elastic in the model

molybdenite re os and muscovite 40ar 39ar dating of the xihuashan tungsten deposit central nanling district south china abstract the xihuashan tungsten deposit in the central nanling region south china is an important vein type ore deposit hosted in cambrian strata and mesozoic granitic intrusions wolframite and molybdenite are the principal ore minerals the gangue minerals are mainly quartz and muscovite wolframite and molybdenite are products of the first stage hydrothermal activity whereas muscovite formed dominantly at the second stage molybdenite re os and muscovite 40ar 39ar dating have been carried out to investigate the age of mineralization re os isotopic dating for molybdenite associated with wolframite yield a precise well constrained isochron age of 157 8 0 9 ma mswd 1 5 ar ar isotopic analyses of muscovite yield a plateau age of 152 8 1 6 ma in agreement with an inverse isochron age of 152 8 1 6 ma which is 5 mys younger than the re os age the molybdenite re os age is interpreted as the age of tungsten mineralization this age coincides well with the zircon u pb age of the host granitic intrusion reported previously the 5 mys difference between molybdenite re os and muscovite 40ar 39ar ages probably represents the duration of hydrothermal activity the results show that the xihuashan tungsten deposit is one of many important tungsten tin deposits formed during 150 to 160 ma associated with large scale lithospheric extension in south china

some problems concerning relationship between mesozoic cenozoic lithospheric extension and uranium metallogenesis in south china previous studies have shown that there existed large scale lithospheric extension events during the cretaceous tertiary period in south china and that large scale uranium mineralization occurred at the same time in this area which resulted in the formation of many important uranium deposits such as granite hosted volcanics hosted and carbonaceous siliceous pelitic sedimentary rock hosted ones this paper reviews previous researches on lithospheric extension and uranium metallogenesis in south china since mesozoic and preliminarily discusses their temporal and spatial associations in order to precisely reveal the relationship between lithospheric extension and uranium metallogenesis it is suggested that at least two fundamental problems should be solved the first is the accurate episodes of extension and uranium mineralization and the second is the mechanism governing the control of the tectonic regime over the formation of the uranium deposits to probe into these problems it is necessary to accurately date the uranium deposits and the mafic dikes intruded during the mesozoic extension and to determine the role and origin of gaseous components such as co2 in ore forming fluids

first scalar magnetic anomaly map from champ satellite data indicates weak lithospheric field 1 satellite magnetic anomaly maps derived by different techniques from magsat pogo data vary by more than a factor of 2 in the deduced strength of the lithospheric magnetic field here we present a first anomaly map from new champ scalar magnetic field data after subtracting a recent orsted main and external field model we remove remaining unmodeled large scale external contributions from 120 track segments by subtracting a best fitting uniform field in order to preserve n s trending features the data are not filtered along track direct integration of the spherically gridded data yields the final degree 14 65 spherical harmonic expansion of the total intensity anomaly at the mean satellite altitude of 438 km apart from enhanced long wavelength features and a smoother general appearance our initial map is strikingly similar to one of the lower strength magsat pogo maps

in situ sims uraninite u pb dating and genesis of the xianshi granite hosted uranium deposit south china abstract the southeastern part of the nanling metallogenic province china is host to numerous granite hosted vein type hydrothermal uranium deposits the geology and geochemistry of these deposits have been extensively studied however accurate and precise ages for the uranium mineralization are scarce because the uranium minerals in these deposits are usually fine grained and may have formed in several stages therefore the ages previously obtained by the bulk dating techniques are possibly a mixed age the xianshi uranium deposit located in the southeastern part of the guidong granite complex is a major uranium deposit in south china the uranium mineralization from this deposit is mainly fine grained uraninite in quartz or calcite veins which are spatially associated with the cretaceous mantle derived mafic dykes micro raman spectroscopy and x ray diffraction analyses indicate that the dominant uranium mineral occurs as a rare form of uraninite u 3 o 7 three distinct generations of uranium minerals have been identified based on petrographic and field relations stage 1 uraninite has the lowest uo 2 and highest pbo contents whereas stage 3 uraninite has the highest uo 2 and lowest pbo contents uraninite from the xianshi deposit was dated using an in situ sims u pb dating technique the results show three distinct age groups 135 4 ma 113 2 ma and 104 2 ma which are in excellent agreement with the ages of three episodes of mantle derived mafic dykes therefore the xianshi uranium deposit has experienced at least three hydrothermal events that are responsible for the deposition of uranium ores which are genetically related to the emplacement of three sets of mafic dykes

mesozoic magmatism and granitic dome in the wugongshan massif jiangxi province and their genetical relationship to the tectonic events in southeast china abstract in se china a mesozoic granitic dome coeval to extensional tectonics is exposed in the wugongshan massif its center is occupied by mesozoic granitoids and granitic gneisses and the e w trending late paleozoic mesozoic pingxiang and anfu basins are located along its northern and southern sides respectively in this paper the petrological geochemical major incompatible and rare earth elements and isotopic characteristics of the granitoids and granitic gneisses of this dome allow us to constrain their geodynamic evolution the wugongshan mesozoic k feldspar enriched granitoids consist of six granitic bodies and lie within a paleozoic metamorphic basement the caledonian shanzhuang granodiorite containing some dioritic enclaves occurs in the study area and is a calc alkaline type granodiorite characterised by lower sio2 and composed of andesine biotite quartz hornblende magnetite and sphene the mesozoic granitoids either indosinian plutons or yanshanian bodies show a geochemical affinity of peraluminous s type granite having higher sio2 k2o contents al k na 2ca enriched in rb th and depleted in eu they also exhibit lree enriched patterns and marked negative eu anomalies the yanshanian plutons also contain characteristic minerals such as sillimanite garnet and monazite the wugongshan mesozoic granitoids have higher isr values 0 70981 0 72885 and lower end t 10 6 to 14 7 an average of 12 6 than those of the early paleozoic shanzhuang granodiorite the mesozoic granitoids were likely derived from ancient sedimentary rocks by partial melting the wugongshan mesozoic granitoids contain three petrological zones that is an early cretaceous massive granitic zone in the core zone i a foliated granitic zone zone ii and a triassic granitic gneissic zone in the outer rim zone iii the petrological and geochemical features of massive granitoids zone i are similar to those of foliated granites zone ii and of granitic gneisses zone iii in this mesozoic dome moreover geometric and kinematic features within the ductile parts units 2 and 3 are also similar to those within the brittle part unit 1 of the dome all these features show that the geologic evolution of the wugongshan doming is a continuous process lasting from triassic to early cretaceous time published 40ar 39ar and k ar ages in the wugongshan mesozoic granitic rocks fall into two major age fields corresponding to two geodynamic events of intracontinental deformation and magmatic activity the first tectonothermal event dated at 226 259 ma on the granitic gneisses in the outer rim of dome can be linked to indosinian collision between the south china plate and the north china plate during triassic time the second event took place during the early yanshanian dated at 180 120 ma on the granites of the core of dome and are responsible for final doming of the wugongshan plutons this event may be related to a westward subduction of the paleo pacific oceanic plate under the se china continent during late jurassic early cretaceous time

mechanism of mineralization in the changjiang uranium ore field south china evidence from fluid inclusions hydrothermal alteration and h o isotopes abstract the changjiang uranium ore field which contains 10 000 tonnes of recoverable u with a grade of 0 1 0 5 is hosted by triassic two mica and jurassic biotite granites and is one of the most important uranium ore fields in south china the minerals associated with alteration and mineralization can be divided into two stages namely syn ore and post ore the syn ore minerals are primarily quartz pitchblende hematite hydromica chlorite fluorite and pyrite the post ore minerals include quartz calcite fluorite pyrite and hematite the fluid inclusions of the early syn ore stage characteristically contain o 2 and those of the late syn ore and post ore stage contain h 2 and ch 4 the fluid inclusions in quartz of the syn ore stage include h 2 o h 2 o co 2 and co 2 types and they occur in clusters or along trails homogenization temperatures t h for the h 2 o co 2 and two phase h 2 o inclusions range from 106 c to 350 c and cluster in two distinct groups for each type salinities are lower than 10 wt nacl equiv the ore forming fluids underwent co 2 effervescence or phase separation at 250 c under a pressure of 1000 1100 bar the u th values of the altered granites are lowest close to the ore increase outwards but subsequently decrease close to unaltered granites from the unaltered granites to the ore the lowest fe 2 o 3 feo values become lower and the highest values higher the ree patterns of the altered granites and the ores are similar to each other the u contents of the ores show a positive correlation with total ree contents but a negative correlation with lree hree ratios indicating the pitchblende is ree bearing and selectively hree rich the eu values of the ore show a positive correlation with u contents indicating the early syn ore fluids were oxidizing the ce values show a negative correlation indicating the later mineralization environment became reducing the water rock interactions of the early syn ore stage resulted in oxidization of altered granites and reduction of the ore forming fluids and it was this reduction that led to the uranium mineralization during alteration in the early syn ore stage the oxidizing fluids leached uranium from granites close to faults and fe 2 o 3 feo ratios increased in the alteration zones the late syn ore and post ore alteration decreased the fe 2 o 3 feo ratios in the alteration zones the 18 o w smow values of the ore forming fluids range from 1 8 to 5 4 and the d w smow values range from 104 4 to 51 6 suggesting meteoric water the meteoric water underwent at least two stages of water rock interaction the first caused the fluids to become uranium bearing and the second stage which was primarily associated with ore bearing faults led to uranium deposition as pitchblende accompanied by co 2 effervescence

mesozoic metallogeny in east china and corresponding geodynamic settings an introduction to the special issue abstract the giant east china mesozoic metallogenic province hosts some of the world s largest resources of tungsten tin molybdenum antimony and bismuth ores of gold silver mercury lead zinc copper uranium and iron are also of major importance the province and its constituent metallogenic belts or regions south china middle lower yangtze river valley east qinling dabie interior of north china craton yan liao and north east china are the products of several pulses of igneous activity and mineralisation between 240 and 80 ma each successive stage has produced a distinctive suite of deposits that can be readily related to the geodynamic evolution of the region during the mesozoic this geodynamic evolution is linked to a complex series of tectonic events involving far field subduction plate collisions crustal thickening post collision collapse and rifting

resolution of direction of oceanic magnetic lineations by the sixth generation lithospheric magnetic field model from champ satellite magnetic measurements the champ satellite continues to provide highly accurate magnetic field measurements from decreasing orbital altitudes 350 km at solar minimum conditions using the latest 4 years 2004 2007 of readings from the champ fluxgate magnetometer including an improved scalar data product we have estimated the lithospheric magnetic field to spherical harmonic degree 120 corresponding to 333 km wavelength resolution the data were found to be sensitive to crustal field variations up to degree 150 down to 266 km wavelength but a clean separation of the lithospheric signal from ionospheric and magnetospheric noise sources was achieved only to degree 120 this new mf6 model is the first satellite based magnetic model to resolve the direction of oceanic magnetic lineations revealing the age structure of oceanic crust

the mackenzie river magnetic anomaly yukon and northwest territories canada evidence for early proterozoic magmatic arc crust at the edge of the north american craton abstract we characterize the nature of the source of the high amplitude long wavelength mackenzie river magnetic anomaly mra yukon and northwest territories canada based on magnetic field data collected at three different altitudes 300 m 3 5 km and 400 km the mra is the largest amplitude 13 nt satellite magnetic anomaly over canada within the extent of the mra source depth estimates 8 12 km from euler deconvolution of low altitude aeromagnetic data show coincidence with basement depths interpreted from reflection seismic data inversion of high altitude 3 5 km aeromagnetic data produces an average magnetization of 2 5 a m within a 15 to 35 km deep layer a value typical of magmatic arc complexes early proterozoic magmatic arc rocks have been sampled to the southeast of the mra within the fort simpson magnetic anomaly the mra is one of several broad scale magnetic highs that occur along the inboard margin of the cordillera in canada and alaska which are coincident with geometric changes in the thrust front transition from the mobile belt to stable cratonic north america the inferred early proterozoic magmatic arc complex along the western edge of the north american craton likely influenced later tectonic evolution by acting as a buttress along the inboard margin of the cordilleran fold and thrust belt

why no anomaly is visible over most of the continent ocean boundary in the global crustal magnetic field abstract it is generally believed that an anomaly should be expected over the continent ocean c o boundary in the global magnetic anomaly maps however such maps prepared at satellite altitude do not show an anomaly over most of the c o boundary we address this issue by a forward modelling technique we use a geographic information system gis technique to integrate information of the various geological units of continents and oceans the seismic crustal structure and standard susceptibility values of the rock types and compute a global vertically integrated susceptibility vis model in addition a remanent magnetisation model is used for the oceanic crust combining the induced and remanent magnetisation models the crustal field anomaly is predicted at a satellite altitude of 400 km for spherical harmonic degrees 16 80 the results show that an anomaly is indeed not predicted over the c o boundaries where the oceanic crust is flanked by continental regions phanerozoic in age we demonstrate that largely a c o anomaly over these regions is not even expected in the global magnetic anomaly map for degrees 1 80 confirming that it is not a long wavelength feature masked by the core field an anomaly however is predicted over some parts of the c o boundary where the flanking geological provinces over the continents are precambrian in age we argue that the vis values for the continents and oceans are largely comparable and therefore no anomaly is expected over most of the c o boundary

comparing crustal and mantle fabric from the north american craton using magnetics and seismic anisotropy a central target in earth sciences is the study of deformation at various depth levels within the earth seismology has offered a remarkable tool for doing this via seismic anisotropy it is however not always clear how to interpret those observations a question of interest is to understand the relation between the deformation of the mantle and the crust and in studying the relation between the two mantle deformation is expressed in seismic anisotropy in this paper we seek an objective way of extracting information about crustal fabric as well to be able to compare with seismic anisotropy the magnetization of crustal rocks offers an attractive possibility for doing this we thus explore the use of magnetic data and we compare magnetic crustal fabric orientation with mantle fabric observations from seismic anisotropy we apply our technique to the north american craton for which we have an excellent magnetic dataset and we show that there is a clear relation between crustal and mantle fabric for the cratonic region this has important implications for crustal formation and for interpreting seismic anisotropy observations

a new estimation of the decay of aftershock density with distance to the mainshock we investigate how aftershocks are spatially distributed relative to the mainshock compared to previous studies ours focuses on earthquakes causally related to the mainshock rather than on aftershocks of previous aftershocks we show that this distinction can be made objectively but becomes uncertain at long time scales and large distances analyzing a regional earthquake data set it is found that at time t following a mainshock of magnitude m the probability of finding an aftershock at distance r relative to the mainshock fault decays as r where is typically between 1 7 and 2 1 for 3 m 6 and is independent of m for r less than 10 to 20 km and t less than 1 day uncertainties on this probability at larger r and t do not allow for a correct estimation of this spatial decay we further show that a static stress model coupled with a rate and state friction model predicts a similar decay with an exponent 2 2 in the same space and time intervals this suggests that static stress changes could explain the repartition of aftershocks around the mainshock even at distances larger than 10 times the rupture length

sulfides associated with podiform bodies of chromite at tsangli eretria greece small bodies of pyrrhotite chalcopyrite minor pentlandite and magnetite occur at the peripheries of podiform bodies of chromite in ultramafic ophiolitic rocks at tsangli eretria central greece banding of magnetite and sulfide within the bodies is reminiscent of magmatic banding a magmatic origin has been proposed for similar sulfide masses in the troodos ophiolite panayiotou 1980 the compositions of the host rocks chromite and of the sulfides have been investigated on average the sulfide mineralization recalculated to metal content in 100 sulfide contains 0 55 ni 5 15 cu 0 29 co 9 ppb pd 179 ppb pt 16 ppb rh 112 ppb ru 31 ppb ir 58 ppb os and 212 ppb au these metal contents particularly the high cu cu ni ratio of 0 78 and the pt pd pt ratio of 0 95 are inconsistent with the sulfides having reached equilibrium with their ni rich host rocks at magmatic temperatures and accordingly it is concluded that they are not of magmatic origin the average 34s value of the sulfide bodies is 2 while that of a sample of pyrite from country rock schist is 15 6 these values are inconclusive as to the origin of the sulfur it is suggested that the sulfides have been precipitated by hydrothermal fluids possibly those responsible for the serpentinization of the host rocks the source of the metals may have been the host rocks themselves

aftershock modeling based on uncertain stress calculations 1 we discuss the impact of uncertainties in computed coseismic stress perturbations on the seismicity rate changes forecasted through a rate and state dependent frictional model we aim to understand how the variability of coulomb stress changes affects the correlation between predicted and observed changes in the rate of earthquake production we use the aftershock activity following the 1992 m7 3 landers california earthquake as a case study to accomplish these tasks we first analyze the variability of stress changes resulting from the use of different published slip distributions we find that the standard deviation of the uncertainty is of the same size as the absolute stress change and that their ratio the coefficient of variation cv is approximately constant in space this uncertainty has a strong impact on the forecasted aftershock activity if a rate and state frictional model is considered we use the early aftershocks to invert for friction parameters and the coefficient of variation by means of the maximum likelihood method we show that when the uncertainties are properly taken into account the inversion yields stable results which fit the spatiotemporal aftershock sequence the analysis of the 1992 landers sequence demonstrates that accounting for realistic uncertainties in stress changes strongly improves the correlation between modeled and observed seismicity rate changes for this sequence we measure a friction parameter asn 0 017 mpa and a coefficient of stress variation cv 0 95

a search for evidence of secondary static stress triggering during the 1992 mw7 3 landers california earthquake sequence secondary triggering of aftershocks is widely observed and often ascribed to secondary static stress transfer however small to moderate earthquakes are generally disregarded in estimates of coulomb stress changes cfs either because of source parameter uncertainties or a perceived lack of importance we use recently published high quality focal mechanisms and hypocenters to reassess the role of small to moderate earthquakes for static stress triggering of aftershocks during the 1992 mw7 3 landers california earthquake sequence we compare the cfs imparted by aftershocks 2 m 6 onto subsequent aftershocks with the total cfs induced by the m 6 main shocks we find that incremental stress changes between aftershock pairs are potentially more often positive than expected over intermediate distances cumulative aftershock stress changes are not reliable for receivers with nearby aftershock stress sources because we exclude unrealistic aftershock stress shadows that result from uniform slip models nonetheless 27 of aftershocks receive greater positive stress from aftershocks than from the main shocks overall 85 of aftershocks are encouraged by the main shocks while adding secondary stress encourages only 79 we infer that source parameter uncertainties of small aftershocks remain too large to convincingly demonstrate or rule out that secondary stress transfer induces aftershocks an important exception concerns aftershocks in main shock stress shadows well resolved secondary stress from detected aftershocks rarely compensates negative main shock stress these aftershocks require a different triggering mechanism

new global vector magnetic anomaly maps derived from magsat data new global maps of the vector magnetic anomaly field of earth are derived from magsat data improvement on earlier published vector and scalar maps is achieved with the following procedure 1 the noise level of vector measurements associated with the discontinuities of the in flight sensor attitude determination is reduced by selecting the data 2 a monthly model of the main field is computed to account for secular variation using only quiet day data 3 the time dependence of the ring current effect is taken into account using hourly values of the dst index based on ground obseivations 4 sectorial indices based on ground observations are introduced to select profiles corresponding to low magnetic activity these indices provide a better discrimination of quiet data than the global kp index usually used for that purpose 5 a time and longitude dependent dynamic model of the equatorial electrojet is derived and used to correct magnetic measurements at low latitudes anomaly maps of x and z components are obtained using selected measurements within the reduced altitude range 350 425 km these maps are compared with earlier magsat anomaly maps and their respective resolutions are discussed these maps evidence large anomalies of lithospheric origin within the oceanic domain in particular the z map displays a characteristic pattern of long wavelength anomalies in these areas it is first shown that some of these anomalies particularly in the north atlantic and indian oceans may be correlated with the age of the lithosphere then a significant similarity between magnetic and geoid anomalies above ocean basins is evidenced this suggests that long wavelength oceanic magnetic anomalies may also be attributed to variations of the thickness of the magnetized layer associated with ocean floor topography

mineral chemistry trace elements and sr nd hf isotope geochemistry and petrogenesis of cailing and furong granites and mafic enclaves from the qitianling batholith in the shi hang zone south china abstract the shi hang zone is an important ne trending mesozoic magmatic belt composed of granites with relative high end t values and young tdm model ages in south china however the petrogenesis and the tectonic environment for the shi hang zone magmatic rocks remain controversial we report here mineral chemistry geochemical and sr nd hf isotopic data for the cailing and furong granites and mafic microgranular enclaves mmes from the qitianling granite batholith in southern hunan province south china the qitianling granite batholith is a multi staged composite pluton with three phases cailing furong and huangtangling according to their ages and petrography the cailing 163 160 ma and furong 157 153 ma phases are mainly composed of porphyritic amphibole biotite monzogranite and they share similar geochemical and isotopic characteristics both of them show similar sio2 contents from 66 50 to 70 28 and metaluminous a cnk values of 0 80 to 0 98 the granites are characterized by high contents of large ion lithosphile elements lile such as rb th u pb high field strength elements hfse such as nb ta zr hf and zr nb ce y contents 350 ppm and high 10 000 ga al ratios 2 6 chondrite normalized ree patterns show relative enrichment of light rare earth elements lrees and significant negative eu anomalies mineralogical and geochemical features suggest that the cailing and furong granites are a type which can be further classified as a2 subtype they have relatively lower 87sr 86sr i ratios 0 7091 0 7132 higher end t values 5 5 to 7 6 and younger nd isotopic model ages 1 48 1 56 ga than those common s type granites in south china zircon ehf t values vary from 8 1 to 3 7 the mmes in the cailing phase show similar trace element and sr nd isotopic characteristics with the host granites but zircons from the mmes show different ehf t values 6 4 2 6 with those from the host granites 8 1 to 3 7 this indicates that the mmes and host granites were crystallized from different sources of magmas providing direct evidence for mafic felsic magma mixing processes the isotope data indicate that the cailing and furong granites from the qitianling batholith were derived from a hybrid magma consisting of about 80 felsic magma derived from old crust and about 20 mantle derived mafic magma the strong magma mixing at about 160 155 ma caused by intra arc rifting or back arc extension related to subduction of the paleo pacific plate is favored to explain the petrogenesis of the cailing and furong granites as well as the shi hang zone

earth s lithospheric magnetic field determined to spherical harmonic degree 90 from champ satellite measurements summary the champ magnetic field mission is providing highly reliable measurements from which the global lithospheric magnetic field can be determined in unprecedented resolution and accuracy using almost 5 yr of data we derive our fourth generation lithospheric field model termed mf4 which is expanded to spherical harmonic degree and order 90 after subtracting from the full magnetic field observations predicted fields from an internal field model up to degree 15 an external field model up to degree two and the predicted magnetic field signatures for the eight dominant ocean tidal constituents we fit and remove remaining external fields and polar electrojet signatures in a track by track scheme from a subset of least disturbed tracks we estimate the mf4 model by least squares damping ill determined coefficients by regularization the resulting mf4 model provides a good representation of the lithospheric field down to an altitude of about 50 km at lower latitudes with reduced accuracy in the polar regions crustal features come out significantly sharper than in previous models in particular bands of magnetic anomalies along subduction zones become visible by satellite for the first time

genesis of the yuanlingzhai porphyry molybdenum deposit jiangxi province south china constraints from petrochemistry and geochronology abstract the newly discovered yuanlingzhai porphyry molybdenum mo deposit in southern jiangxi province belongs to the group of mo only deposits in the nanling region the mineralization developed at contact zones between the yuanlingzhai granite porphyry and neoproterozoic metamorphic rocks of the xunwu formation precise la mc icpms zircon u pb dating of the yuanlingzhai porphyry as well as the adjacent western keshubei and eastern keshubei granites yielded ages of 165 49 0 59 ma 159 68 0 43 ma and 185 13 0 52 195 14 0 63 ma respectively molybdenite re os isochron ages of the ores are 160 1 162 7 1 1 ma which is consistent with the age of large scale w sn deposits in south china the yuanlingzhai porphyry is characterized by high k2o p2o5 and a cnk 1 33 1 59 and low cao and na2o the rock shows relatively enriched lree without significant eu anomalies eu eu 0 80 0 90 geochemical and mineralogical characteristics indicate that the ore hosting porphyry is a typical s type granite generated from the partial melting of crustal material with only minor mantle contribution both harker and evolutionary discrimination diagrams indicate that the yuanlangzhai and western keshubei granites are not products of co magmatic evolution the keshubei granites and xunwu formation were not significant sources for the components in the porphyry mineralization but the yuanlangzhai granite may have supplied some ore forming material however the main ore forming material was carried by fluids from deep sources as demonstrated by fluid inclusion and stable isotope data from the molybdenum deposit the mo porphyry deposit formed in an extensional setting and was possibly associated with jurassic subduction of the izanagi plate

geopotential field anomaly continuation with multi altitude observations abstract conventional gravity and magnetic anomaly continuation invokes the standard poisson boundary condition of a zero anomaly at an infinite vertical distance from the observation surface this simple continuation is limited however where multiple altitude slices of the anomaly field have been observed increasingly areas are becoming available which are constrained by multiple boundary conditions from surface airborne and satellite surveys this paper describes the implementation of continuation with multi altitude boundary conditions in cartesian and spherical coordinates and investigates the advantages and limitations of these applications continuations by eps equivalent point source inversion and the ft fourier transform as well as by scha spherical cap harmonic analysis are considered these methods were selected because they are especially well suited for analyzing multi altitude data over finite patches of the earth such as those covered by the admap database in general continuations constrained by multi altitude data surfaces are invariably superior to those constrained by a single altitude data surface due to anomaly measurement errors and the non uniqueness of continuation

noon time equatorial electrojet its spatial features as determined by the champ satellite 1 new observations obtained by the satellite champ reveal a detailed picture of the noon time equatorial electrojet eej the low orbit of champ and its high precision magnetometers reveal the spatial structure of the eej with unprecedented accuracy data from more than two and a half years have been used to investigate average features and also the global characteristics of the eej rather than interpreting the magnetic signatures we determined the horizontal current distribution by using a very general current model series of line currents this makes the results independent of satellite altitude and ambient field geometry the procedure for determining the structure of the electrojet is fully automated giving an objective response some of the spatial features of the noon time eej are as follows the electrojet current peaks right at the dip equator there is no deviation from it either on a seasonal basis or with longitude the width of the eej 4 in latitude at half the peak value of the current density profile is for a given longitude fairly constant independent of the amplitude return currents north and south of the eastward current are a common feature of the eej they peak at latitudes some 5 away from the dip equator the intensity of the eej varies strongly from day to day the average peak current density exhibits a clear dependence on longitude peaks show up over south america and indonesia the average current density follows closely the monthly mean of the solar flux index f10 7 the total eej eastward current is about three times as strong as the return current the total current and the peak current density are related to each other by a power law this suggests that the longitude dependence of the eej intensity can be explained by the varying cross sectional area of the cowling channel

evidence for short spatial correlation lengths of the noontime equatorial electrojet inferred from a comparison of satellite and ground magnetic data 1 the current density of the noontime equatorial electrojet eej as determined from champ data is highly variable between successive passes of the satellite which are separated by 23 in distance and 93 min in time an open question is to which extent this variability is caused by temporal or spatial variations in the ionosphere another important question is the connection between eej and global solar quiet sq current systems we try to answer these questions by comparing the eej current density estimated from highquality scalar magnetic field measurements of the champ satellite with the magnetic horizontal intensity variations at six equatorial observatory pairs distributed across the globe data taken during the period 2000 2002 were used for the present study we apply corrections for the effect of local time lt and sq fields by estimating the correlation coefficients between the ground and satellite data as a function of distances between measurements new insights into the spatial structure of the eej have been obtained the high correlation when champ passes directly over an observatory decays quickly in eastern and western directions typically within 15 of longitudinal separation between satellite and observatory the correlation falls well below the statistical significance level this observation holds for all longitude sectors interestingly the correlation between champ inferred eej strength and observatory differences breaks down for the observatory pairs outside of a 4 latitudinal band this implies that the eej and sq variations are uncorrelated for periods up to 1 hour additionally it was found that monitoring of the eej can be performed best if the reference observatory is 4 to 5 apart from the dip equator

genetic relationship between silver lead zinc mineralization in the wutong deposit guangxi province and mesozoic granitic magmatism in the nanling belt southeast china more than 50 of the world s total reserves of tungsten are in china and most tungsten deposits are located in the nanling range in southeast china this study explores the potential genetic relationship between tungsten tin w sn mineralization and shallower ag pb zn deposits in the nanling range based on data from the wutong deposit guangxi province the lead oxygen carbon sulfur and strontium isotopic compositions of minerals at wutong indicate that a single crustal derived fluid was responsible for mineralization wutong likely formed at relatively low temperatures 200 300 c and low pressures as indicated by the similarity between homogenization temperatures of fluid inclusions and those estimated from s isotopic compositions of minerals the hubnerite age 92 3 104 4 ma indicates that the wutong mineralization is likely related to nearby late yanshanian cretaceous s type granites derived from proterozoic crust this mineralization event coincides with the last w sn mineralization event and the cretaceous peak of mineralization in the nanling range

isotopic and paleomagnetic constraints on the mesozoic tectonic evolution of south china in order to better constrain the paleogeographic evolution of south china we measured sm nd and rb sr isotopic compositions for 23 mesozoic granites that crop out throughout the area tightly grouped neodymium depleted mantle model ages 1 4 0 3 ga suggest the region is underlain by relatively homogeneous proterozoic crust and fail to define crustal provinces neither the isotopic nor geologic data suggest that a mesozoic suture exists however granites possessing anomalously high sm 8 ppm and nd 45 ppm concentrations relatively high initial epsilon neodymium 4 to 8 and high but variable initial 87sr 86sr 0 759 to 0 713 form a northeast trending zone that coincides with two prominent mesozoic basins southeast of the zone lie the majority of mesozoic intrusives and upper triassic to lower cretaceous extensional basins found in south china mesozoic paleomagnetic poles are well clustered northwest of the zone pre cretaceous poles southeast of it are discordant with respect to those from the northwest the only recognized tectonostratigraphic terrane in south china lies southeast of the zone the terrane is bordered by a northeast trending sinistral fault that was active in the mesozoic other faults in south china have similar attitudes ages and sense of shear together the observations suggest that the mesozoic tectonic regime in south china consisted of strike slip activity plus concomitant rifting as terranes or fragments of similar crust were transported north along sinistral faults the zone defined by the granites enriched in nd and sm demarcates displaced terranes to the southeast from relatively stable land to the northwest

extensional tectonics within a subduction type orogen the case study of the wugongshan dome jiangxi province southeastern china abstract in southeastern china the wugongshan massif consists of metamorphic rocks and granites generally considered as basement rocks formed during an early paleozoic orogeny our structural and radiometric reappraisal of this massif shows that it results from a mesozoic doming in the metamorphic rocks the foliation defines an e w elliptical shape and bears a n s stretching lineation the axial zone of the dome is occupied by oriented plutons the outer part of the dome consists of devonian quartzite carboniferous permian limestone and triassic jurassic sandstone ductile deformation observed in the permian carbonates and older units develops contrasting kinematics the northern and southern flanks present top to the north and top to the south motions respectively the same divergent pattern of motion is recorded from brittle structures up to the jurassic sandstone the late cretaceous red beds are undeformed in the axial zone coaxial flow indicated by symmetric pressure shadows and quartz fabrics predominates along the n and s flanks of the dome post folial recumbent folds are symmetrically overturned to the n and s respectively but at the two extremities of the dome the fold axes are parallel to the stretching lineation 40 ar 39 ar ages of biotite and muscovite from mylonite and granite suggest a late triassic age for the ductile deformation and an early cretaceous age for the final doming the wugongshan dome is similar to a metamorphic core complex formed in an extensional tectonic regime a model which emphasizes plutonism is discussed from permian to early cretaceous southeastern china was a subduction type orogen the crust of which was thickened by accretion of mantle derived magmas the emplacement of a large amount of mesozoic granitoids in southeastern china promoted thermal softening and gravitational instability that allowed the crust to extend the formation of the wugongshan dome and other metamorphic core complexes appears to be a the direct consequence of plutonic activity and crust rheology

can electromagnetic disturbances related to the recent great earthquakes be detected by satellite magnetometers abstract on december 26 2004 the world s fourth largest earthquake since 1900 and the largest since the 1964 prince william sound alaska earthquake occurred off the west coast of northern sumatra with a magnitude of 9 3 on march 28 2005 another event of magnitude 8 7 took place in the same region the december 26 2004 earthquake has prompted scientists to investigate possible electromagnetic signatures of this event using ground magnetic observations iyemori et al iyemori t et al 2005 geomagnetic pulsations caused by the sumatra earthquake on december 26 2004 geophys res lett 32 l20807 doi 10 1029 2005gl024083 have suggested that a 3 6 min long geomagnetic pulsation observed shortly after this event was generated by the earthquake they have speculated that a 30 s magnetic pulsation was also caused by the earthquake here for the first time champ satellite magnetic and electron density data are examined to find out if electromagnetic signatures which are possibly related to these recent megathrust earthquakes are observed in satellite magnetic data we have shown that some specific features are observed after the two earthquakes with periods of about 16 and 30 s our results favor an external source origin for the 30 s pulsation moreover after more than 1 h champ magnetic data indicate the existence of a feature characterized by the same parameters duration amplitude and frequency content which could be associated with each earthquake respectively further investigations are required in order to answer the question of whether these signals can be associated with earthquakes and to assign their possible usefulness with respect to earthquake development

origin of the carbonaceous debris and its implication for mineralization within the qianjiadian uranium deposit southern songliao basin abstract appearance and coal petrology of the carbonaceous debris cd from the qianjiadian uranium deposit southern songliao basin are illustrated in order to reveal the relationship between the cd and uranium mineralization response of the cd and uranium mineralization to diabase dikes and origin of the cd three types of cd have been identified in the yaojia formation cd serving as trapped sediments cd serving as laminae of cross bedding and cd occurring in lumpy sandstones with an average content of 91 2 in the cd organic macerals are primarily composed of homoeollinite desmocollinite and inertinite inorganic macerals of the cd mainly include pyrite with an average content of 8 8 average contents of the desmocollinite inertinite and pyrite from the cd in the ore bearing section are higher than that in the non ore bearing section respectively the study demonstrates that the cd rich in inertinite and desmocollinite facilitates uranium mineralization uranium is frequently adsorbed in the desmocollinite and inertinite and coffinite is embedded in the inertinite inferring that roles of the cd in uranium mineralization are both adsorption and reduction vitrinite reflectance vr of the cd ranges from 0 24 to 0 513 with an average value of 0 335 the vr of the cd in the ore bearing section is generally larger than that in the non ore bearing section uranium mineralization and enlargement of the vr are both attributed to the diabase dikes the cd originated from the late cretaceous wildfire related woody plants and it has been transported into fluvial channels and preserved in sandstones

impact of earthquake rupture extensions on parameter estimations of point process models abstract stochastic point processes are widely applied to model spatiotemporal earthquake occurrence in particular the epidemic type aftershock sequence etas model has been shown to successfully reproduce the short term clustering of earthquakes an important parameter of the model is the value describing the scaling of the aftershock productivity with magnitude of the triggering earthquake according to 10 m fitting of the space dependent etas model to empirical data yields values that are typically much smaller than the scaling inverted from more simple stacking of aftershock sequences we show by means of synthetic simulations that this is likely to result from assuming spatial isotropy of aftershock occurrence that in fact aligns along the mainshock rupture we fit the space dependent and space independent etas models to simulations where each earthquake is a line source with an empirical magnitude length relation although the space time model describes past activity quite well it overestimates the forecasted earthquake rate on the other hand the application of the space independent etas model predicts future seismicity well and can therefore be applied for forecasting purposes our test for the observed aftershock sequence following the 1992 m 7 3 landers earthquake supports these results

precise molybdenite re os and mica ar ar dating of the mesozoic yaogangxian tungsten deposit central nanling district south china the yaogangxian deposit in the central nanling region south china consists of vein type ore bodies hosted in cambrian to jurassic strata and mesozoic granitic intrusions wolframite and molybdenite are the dominant ore minerals intergrown with gangue minerals of quartz feldspar phlogopite and muscovite we have carried out molybdenite re os and phlogopite and muscovite 40ar 39ar dating to better understand the timing and genesis of mineralization re os dating of eight molybdenite samples yielded model ages ranging from 152 0 3 5 to 161 1 4 5 ma with an average of 156 0 ma the re os analyses give a well defined 187re 187os isochron with an age of 154 9 2 6 ma mswd 2 4 hydrothermal phlogopite and muscovite display extremely flat 40ar 39ar age spectra phlogopite yields a 40ar 39ar plateau age of 153 0 1 1 ma whereas muscovite yields a plateau age of 155 1 1 1 ma both 40ar 39ar ages are in good agreement with the re os ages placing the timing of tungsten mineralization at about 154 ma this age is consistent with the field relationships our new data when combined with published geochronological results from other major deposits in this region suggest that large scale w sn mineralization occurred throughout the central nanling region in the late jurassic

improving empirical aftershock modeling based on additional source information 1 several mechanisms are proposed to underlie earthquake triggering including static stress interactions and dynamic stress transfer significant differences of these mechanisms are particularly expected in the spatial distribution of aftershocks however testing the different hypotheses is challenging because it requires the consideration of the large uncertainties involved in stress calculations as well as the appropriate consideration of secondary aftershock triggering which is related to stress changes induced by smaller pre and aftershocks in order to evaluate the forecast capability of different mechanisms we take the effect of smaller magnitude earthquakes into account by using the epidemic type aftershock sequence etas model where the spatial probability distribution of direct aftershocks if available is correlated to alternative source information and mechanisms we test surface shaking rupture geometry and slip distributions as an approximation of the shaking level we use shakemap data which are available in near real time after a main shock and thus could be used for first order forecasts of the spatial aftershock distribution alternatively we test the use of empirical decay laws related to minimum fault distance and coulomb stress change calculations based on published and random slip models for comparison we analyze the likelihood values of the different model combinations in the case of several well known aftershock sequences 1992 landers 1999 hector mine 2004 parkfield our test shows that the fault geometry is the most valuable information for improving aftershock forecasts furthermore we find that static stress maps can additionally improve the forecasts of off fault aftershock locations while the integration of ground shaking data could not upgrade the results significantly

mineralization and fluid inclusion study of the shizhuyuan w sn bi mo f skarn deposit hunan province china the shizhuyuan deposit china is a world class w sn bi mo f skarn deposit hosted by devonian limestone in the thermal aureole of the qianlishan granite complex the qianlishan complex comprises five separate intrusions including fine grained porphyritic biotite granite 182 187 ma granite 1 medium grained biotite k feldspar granite 158 163 ma granite 2 fine grained biotite and k feldspar granite granite 3 granitic porphyry 144 146 ma granite 4 and diabase 142 ma the upper part of the granite complex is characterized by extensive greisen alteration skarn zones distributed around the intrusions are mainly calcic consisting of garnet garnet pyroxene vesuvianite garnet and wollastonite vesuvianite progressively outward from the intrusion following the primary skarn formation some of the skarn zones underwent retrograde alteration the high garnet pyroxene ratio and the diopside rich and andradite rich compositions of the pyroxene and garnet indicate that the skarn belongs to the oxidized type mineralization consists of sn be veinlet ore type i in marble and porphyry massive w bi mo sn skarn ore type ii stockwork w sn bi mo f ore type iii and w sn mo bi greisen ore type iv mainly associated with granite 2 emplacement of granite 2 was accompanied by late intense fracturing characterized by stockwork mineralization type iii which was superimposed on massive skarn and greisen zones the stockwork ore consists mainly of greisen and skarn veins and veinlets with scheelite wolframite molybdenite cassiterite bismuthinite and fluorite the sm nd method has been used to analyze the sm nd concentration and data for pyroxene and garnet in the massive skarn type ii ore associated with granite 2 yields an isochron age of 157 ma this age corresponds closely to the age of granite 2 fluid inclusion studies of four types of ore samples reveal four types of inclusions aqueous two phase inclusions h2o co2 inclusions gas rich inclusions and daughter mineral bearing inclusions containing halite or calcite homogenization temperatures of fluid inclusions in skarn minerals range from 350 to 535 c the homogenization temperatures of fluid inclusions in greisen and stockwork are lower than that in skarn and range from 200 to 360 c the fluid inclusion data indicate that there are two types of fluids associated with the massive skarn and greisen having salinities from 26 to 41 wt percent nacl equiv and from 1 to 21 wt percent nacl equiv respectively the high salinity inclusions are daughter mineral bearing whereas the low salinity inclusions are dominantly aqueous with a few gas rich inclusions the results suggest the evolution of the ore forming fluids either by immiscibility or by mixing between a high temperature high salinity magmatic water and a low temperature low salinity fluid such as meteoric water consistent with previously published isotopic studies the shizhuyuan deposit is a multi element skarn highly enriched in w sn bi mo and f this is interpreted to reflect multiple sources and stages of mineralization which in turn may account for the large size of the shizhuyuan deposit

spatial and temporal distribution of slip for the 1992 landers california earthquake we have determined a source rupture model for the 1992 landers earthquake m w 7 2 compatible with multiple data sets spanning a frequency range from zero to 0 5 hz geodetic survey displacements near field and regional strong motions broadband teleseismic waveforms and surface offset measurements have been used explicitly to constrain both the spatial and temporal slip variations along the model fault surface our fault parameterization involves a variable slip multiple segment finite fault model which treats the diverse data sets in a self consistent manner allowing them to be inverted both independently and in unison the high quality data available for the landers earthquake provide an unprecedented opportunity for direct comparison of rupture models determined from independent data sets that sample both a wide frequency range and a diverse spatial station orientation with respect to the earthquake slip and radiation pattern in all models consistent features include the following 1 similar overall dislocation patterns and amplitudes with seismic moments of 7 to 8 10 26 dyne cm seismic potency of 2 3 to 2 7 km 3 2 very heterogeneous unilateral strike slip distributed over a fault length of 65 km and over a width of at least 15 km though slip is limited to shallower regions in some areas 3 a total rupture duration of 24 sec and an average rupture velocity of 2 7 km sec and 4 substantial variations of slip with depth relative to measured surface offsets the extended rupture length and duration of the landers earthquake also allowed imaging of the propagating rupture front with better resolution than for those of prior shorter duration strike slip events our imaging allows visualization of the rupture evolution including local differences in slip durations and variations in rupture velocity rupture velocity decreases markedly at shallow depths as well as near regions of slip transfer from one fault segment to the next as rupture propagates northwestward along the multiply segmented fault length the rupture front slows as it reaches the northern limit of the johnson valley landers faults where slip is transferred to the southern homestead valley fault an abrupt acceleration is apparent following the transfer this process is repeated and is more pronounced as slip is again passed from the northern homestead valley fault to the emerson fault although the largest surface offsets were observed at the northern end of the rupture our modeling indicates that substantial rupture was also relatively shallow less than 10 km in this region

weakly magnetic crust in the canadian cordillera abstract current models of continental crust favor an increase in magnetization with depth here we report a counter example from the canadian cordillera where almost a full thickness of non magnetic continental crust is suggested by joint interpretation of magnetic and seismic data the magnetic field over the cordillera is characterized by complex short wavelength 100 km portion of the cordilleran field is subdued and mainly featureless and suggests a lack of magnetic sources at greater depths seismic reflection and refraction data from three major transects in the yukon and british columbia canada support this interpretation and indicate that sedimentary like formations make up the majority of the crust the dominance of shallow upper crustal magnetization in the canadian cordillera contrasts with the generally held view that the lower continental crust is the primary source for long wavelength magnetic anomalies sources for these anomalies are often assumed to be located in the lower crust when surface magnetizations are insufficient to produce such anomalies or no correlation exists between the magnetic field and the mapped surface geology the canadian cordillera appears to be an example of a non magnetic lower crust overlain by a more magnetic upper crust that is however not magnetized strongly enough to produce significant long wavelength magnetic anomalies

fifth generation lithospheric magnetic field model from champ satellite measurements 1 six years of low orbit champ satellite magnetic measurements have provided an exceptionally high quality data resource for lithospheric magnetic field modeling and interpretation here we describe the fifth generation satellite only magnetic field model mf5 the model extends to spherical harmonic degree 100 as a result of careful data selection extensive corrections filtering and line leveling the model has low noise levels even if evaluated at the earth s surface the model is particularly suited for inferring large scale structure and composition of the lithosphere it is also meant to serve as the long wavelength part of continental and global scale marine and aeromagnetic anomaly maps

onto what planes should coulomb stress perturbations be resolved 1 coulomb stress maps are produced by computing the tensorial stress perturbation due to an earthquake rupture and resolving this tensor onto planes of a particular orientation it is often assumed that aftershock fault planes are optimally oriented in other words the regional stress and coseismic stress change are used to compute the orientation of planes most likely to fail and the coseismic stress is resolved onto these orientations this practice assumes that faults capable of sustaining aftershocks exist at all orientations an assumption contradicted by the observation that aftershock focal mechanisms have strong preferred orientations consistent with mapped structural trends here we systematically investigate the best planes onto which stress should be resolved for the landers hector mine loma prieta and northridge earthquakes by quantitatively comparing observed aftershock distributions with stress maps based on optimally oriented planes two and three dimensional main shock orientation and regional structural trend we find that the best model differs between different tectonic regions but that in all cases models that incorporate the regional stress field tend to produce stress maps that best fit the observed aftershock distributions although not all such models do so equally well our results suggest that when the regional stress field is poorly defined or in structurally complex areas the best model may be to fix the strike of the planes upon which the stress is to be resolved to that of the main shock but allow the dip and rake to vary

slip distribution of the 1992 landers earthquake and its implications for earthquake source mechanics we use near source 10 to 164 km displacement seismograms to image the slip distribution and rupture history of the 28 june 1992 mw 7 3 landers california earthquake aftershock seismograms from similar distances are modeled to find the velocity model and frequency range 0 05 to 0 25 hz over which theoretical green s functions are most accurate and the measure of fit is used as an upper bound on theoretical error in the mainshock inversion we represent the rupture surface with three planar segments divided into 3 by 3 km elements extending from the surface to 18 km depth and solve for the slip distribution and the rupture time model that minimizes misfit to both recorded seismograms and mapped surface displacement in a least squares sense we investigate two faulting models one with a uniformly short 3 sec rise time everywhere on the fault surface and the other with a variable rise time 2 to 6 sec we perform sensitivity tests using synthetic data from a hypothetical earthquake to evaluate model resolution and solution stability in the sensitivity tests both fault models recover similar slip and rupture features but neither is capable of imaging details of the rise time for the landers earthquake we find an average rupture velocity of 2 5 km sec and use this average for a starting model in a linearized inversion for rupture time in our solutions of slip amplitude distribution the southernmost johnson valley fault segment has 20 of the total seismic moment 6 to 8 1019 n m with small displacements near the hypocenter the homestead valley segment contributes half of the moment with the largest slip amplitudes 25 to 35 km northwest of the hypocenter at 4 to 12 km depth and the camp rock emerson segment contributes the remaining moment with the largest slip amplitude 35 to 50 km northwest of the hypocenter in the shallow crust 9 km there is some evidence that the rupture front is delayed as it encounters high slip regions suggesting that prior to the mainshock these areas were further from failure owing either to greater strength or lower prestress

frequency domain inversion of strong motions application to the 1992 landers earthquake we present a frequency domain inversion in which the observed earthquake strong ground motions are used to constrain the space time dependence of slip on a fault green s functions are numerically evaluated and the parameters describing the rupture are the local slip rupture time and rise time these parameters are simultaneously evaluated without additional constraints this procedure allows for large variations in the local rupture velocity the june 28 1992 landers earthquake mw 7 3 offers an exceptional opportunity to apply this technique to a major strike slip event we model the rupture evolution including local differences in slip durations and variations in rupture velocity our final results are in good agreement with other inversion studies geodetic and surface observations the main discrepancies occurred at depth and at the end of the johnson valley fault we show that a relatively low resolution could be an explanation for these differences rupture velocity and slip are extremely heterogeneous both along strike and with depth a moment of 0 90 1020 n m was found the slip distribution shows that this event consists of a series of regions of high slip subevents separated by regions of relative low slip approximately 50 of the moment was released on the homestead valley fault in this region of large slip the rupture velocity inferred by our inversion is well constrained and is equal to 3 0 km s at depth and 2 5 km s near the surface our inversion favors the hypothesis that the duration of the slip at each point of the fault is of the order of the duration of rupture of each subevent

the satellite along track analysis in planetary magnetism summary many physical characteristics of planets such as their topography magnetic and gravitational fields are routinely detected and measured by spacecrafts at satellite altitudes even if little is known about the measured signal it is possible to separate the large scale components from other contributions by a spectral analysis carried out along the spacecraft orbit this procedure which dates back to the early age of the satellite era is routinely applied with spherical harmonics analysis for the study of large scale planetary magnetic fields particularly those that vary rapidly in time in this paper we review the approximations of this procedure for investigating internal and external magnetic fields of planets we show that the magnetic field analyses along the orbits are limited by a finite frequency resolution that is known to introduce a spectral leakage this leakage may lead to artificial spatio temporal variations of the magnetic field such as an apparent internal field secular variation an asymmetry of the magnetospheric field and a regional distortion of the lithospheric field structures we quantify the errors for the earth s lithospheric field and display its distribution in the spatial and the spectral domains we also discuss these limits for the moon and mars lithospheric magnetic fields our analytical results illustrate the pros and cons of the orbit by orbit analyses and should allow us to avoid the pitfall of geophysical overinterpretation of some artificial magnetic field variations

cretaceous tectonic evolution of south china a preliminary synthesis abstract the cretaceous tectonic evolution of south china is characterized by widespread extensional basin and dome generation voluminous magma intrusion eruption and associated polymetallic mineralization all of which are of world wide interests that have stimulated the attention of geologists for more than half a century due to the lack of a comprehensive understanding of regional tectonic evolution and geodynamics of south china many controversies regarding the nature and origin of these features remain this paper attempts to make a review by synthesizing existing structural petrological geochronological and geochemical data of the cretaceous structures and magmas which guides us to propose a three stage tectono thermal evolutionary history of south china during the cretaceous period the earliest cretaceous 145 137 ma tectonic stage was characterized by syn orogenic shortening deformation and metamorphism under a nw se compressional setting which generated voluminous porphyry cu au ore bearing adakitic rocks in the lower yangtze river belt and gneissic granites in the coastal area its geodynamic origin was interpreted as combined effects of the pacific izanagi ridge subduction beneath the lower yangtze river belt and the pacific subduction beneath the cathaysia block this syn orogenic shortening stage was followed by a post orogenic stage 136 86 ma comprising two episodes of alternate extensional and shortening events the extension 136 118 ma in the earlier episode was dominated by a nw se extensional regime it led to significant taphrogenesis manifested by large scale extensional basins and voluminous magma intrusion eruption this extensional event was associated with a combination of the slab window opening during the ridge subduction and the rollback of the subducted pacific slab the subsequent nw se transpressional event led to cessation of active marginal magmatism and resulted in tectonic inversion of previous rift basins this transpression lasted from 117 ma to 108 ma and its driving mechanism was attributed to the collision between the eastern asian margin and the west philippine block the extension 107 86 ma in the later episode was dominated by a wnw ese extensional regime which led to a second phase of basin subsidence and produced numerous a and i type granites and bimodal volcanoes the tectonic regime then changed at ca 85 ma to wnw ese oriented transpression causing inversion of the late cretaceous rift basins and cessation of extension related magmatism the change of stress field from wnw ese extension to wnw ese transpression was related to the variations of the subducted slab dynamics i e the transition from either the ese ward retreat to wnw ward subduction of the pacific slab or the alternate stress during slab break off during the latest cretaceous a drastic change of regional tectonic stress orientation from wnw ese to n s occurred in south china the third stage predominated by n s extension was developed which caused a new phase of regional crustal subsidence along the youngest e w trending extensional structures

a spatially continuous magnetization model for mars 1 using a three component magnetic field data set at over 100 000 satellite points previously compiled for spherical harmonic analysis we have produced a continuously varying magnetization model for mars the magnetized layer was assumed to be 40 km thick an average value based on previous studies of the topography and gravity field the severe nonuniqueness in magnetization modeling is addressed by seeking the model with minimum root mean square rms magnetization for a given fit to the data with the trade off between rms magnetization and fit controlled by a damping parameter our preferred model has magnetization amplitudes up to 20 a m it is expressed as a linear combination of the green s functions relating each observation to magnetization at the point of interest within the crust leading to a linear system of equations of dimension the number of data points although this is impractically large for direct solution most of the matrix elements relating data to model parameters are negligibly small we therefore apply methods applicable to sparse systems allowing us to preserve the resolution of the original data set thus we produce more detailed models than any previously published although they share many similarities we find that tectonism in the valles marineris region has a magnetic signature and we show that volcanism south of the dichotomy boundary has both a magnetic and gravity signature the method can also be used to downward continue magnetic data and a comparison with other leveling techniques at mars surface is favorable

an improved space time etas model for inverting the rupture geometry from seismicity triggering this study incorporates the rupture geometry of big earthquakes in the formulation of the epidemic type aftershock sequence etas model which is a point process model widely applied in the study of spatiotemporal seismicity rather than regarding every earthquake occurring at a point in space and time we apply the new model to the catalog from sichuan province china between 1990 and 2013 during which the wenchuan mw7 9 earthquake occurred in may 2008 our results show that the modified model has better performance in both data fitting and aftershock simulation confirming that the elliptic aftershock zone is caused by the superposition of the isotropic triggering effect from each patch of the rupture zone moreover using the technique of stochastic reconstruction we inverted the fault geometry and verified that direct aftershocks of the main shock more likely occur in the transitive parts from high slip parts to low median slip parts of the main shock fault area

comparison of short term and time independent earthquake forecast models for southern california we have initially developed a time independent forecast for southern california by smoothing the locations of magnitude 2 and larger earthquakes we show that using small m 2 earthquakes gives a reasonably good prediction of m 5 earthquakes our forecast outperforms other time independent models kagan and jackson 1994 frankel et al 1997 mostly because it has higher spatial resolution we have then developed a method to estimate daily earthquake probabilities in south ern california by using the epidemic type earthquake sequence model kagan and knopoff 1987 ogata 1988 kagan and jackson 2000 the forecasted seismicity rate is the sum of a constant background seismicity proportional to our time independent model and of the aftershocks of all past earthquakes each earthquake triggers aftershocks with a rate that increases exponentially with its magnitude and decreases with time following omori s law we use an isotropic kernel to model the spatial distribution of aftershocks for small m 5 5 mainshocks for larger events we smooth the density of early aftershocks to model the density of future aftershocks the model also assumes that all earthquake magnitudes follow the gutenberg richter law with a uniform b value we use a maximum likelihood method to estimate the model parameters and test the short term and time independent forecasts a retro spective test using a daily update of the forecasts between 1 january 1985 and 10 march 2004 shows that the short term model increases the average probability of an earthquake occurrence by a factor 11 5 compared with the time independent forecast

ocean tidal signals in observatory and satellite magnetic measurements 1 ocean flow moves sea water through the earth s magnetic field inducing electric fields currents and secondary magnetic fields these motionally induced magnetic fields have a potential for the remote sensing of ocean flow variability a first goal must be to gain a better understanding of magnetic field generation by tidal ocean flow we predict the motionally induced magnetic fields for the six major tidal constituents and compare their amplitudes with the spectra of night time observatory and satellite magnetic measurements for the indian ocean the magnetic variations at the solar s2 kl and p1 periods turn out to be dominated by unrelated external fields in contrast observed lunar m2 and n2 tidal signals are in fair agreement with predictions from motional induction the lunar diurnal o1 signal visible at some observatories could be caused by ocean flow but disagrees in amplitude with our predictions

a constitutive law for rate of earthquake production and its application to earthquake clustering seismicity is modeled as a sequence of earthquake nucleation events in which the distribution of initial conditions over the population of nucleation sources and stress ing history control the timing of earthquakes the model is implemented using solutions for nucleation of unstable fault slip on faults with experimentally derived rate and state dependent fault properties this yields a general state variable constitutive formulation for rate of earthquake production resulting from an applied stressing history to illustrate and test the model some characteristics of seismicity following a stress step have been explored it is proposed that various features of earthquake clustering arise from sensi tivity of nucleation times to the stress changes induced by prior earthquakes the model gives the characteristic omori aftershock decay law and interprets aftershock parameters in terms of stress change and stressing rate earthquake data appear to support a model prediction that aftershock duration defined as the time for rates to return to the back ground seismicity rate is proportional to mainshock recurrence time observed spatial and temporal clustering of earthquake pairs arises as a consequence of the spatial depen dence of stress changes of the first event of the pair and stress sensitive time dependent nucleation applications of the constitutive formulation are not restricted to the simple stress step models investigated here it may be applied to stressing histories of arbitrary complexity the apparent success at modeling clustering phenomena suggests the possibil ity of using the formulation to estimate short to intermediate term earthquake probabil ities following occurrence of other earthquakes and for inversion of temporal variations of earthquake rates for changes in driving stress

geology and genesis of kafang cu sn deposit gejiu district sw china kafang is one of the main ore deposits in the world class gejiu polymetallic tin district sw china there are three main mineralization types in the kafang deposit i e skarn cu sn ores stratiform cu ores hosted by basalt and stratiform cu sn ores hosted by carbonate the skarn mainly consists of garnet and pyroxene and retrograde altered rocks these retrograde altered rocks are superimposed on the skarn and are composed of actinolite chlorite epidote and phlogopite major ore minerals are chalcopyrite pyrrhotite cassiterite pyrite and scheelite sulfur and pb isotopic components hint that the sources of different types of mineralization are distinctive and indicate that the skarn ore mainly originated from granitic magma whereas the basalt hosted cu ores mainly derived from basalt microthermometry results of fluid inclusions display a gradual change during the ore forming process the homogenization temperature of different types of inclusions continuously decreases from early to late mineralization stages the salinities and freezing temperatures exhibit similar evolutionary tendencies with the t homogenizaion while the densities of the different types keep constant the majority being less than 1 oxygen and hydrogen isotopic values o and d of the hydrothermal fluids fall within ranges of 3 1 to 7 7 with an average of 6 15 calculated at the corresponding homogenization temperature and 73 and 98 with an average of 86 5 respectively microthermometry data and h o isotopes indicate that the ore forming fluid of the kafang deposit is mainly derived from magma in the early stage and a mixture of meteoric and magmatic water in late stage molybdenite re os age of the skarn type mineralization is 83 4 2 1 ma and the stratiform ores hosted by basalt is 84 2 7 3 ma which are consistent with the la icp ms zircon age of the xinshan granite intrusion 83 1 0 4 ma the evidence listed above reflects the fact that different ore styles in the kafang deposit belong to the same mineralization system

multiple mesozoic mineralization events in south china an introduction to the thematic issue mesozoic mineral deposits in south china include world class deposits of w sn and sb and those that provide the major sources of ta cu hg as tl pb zn au and ag for the entire country these deposits can be classified into polymetallic hydrothermal systems closely related to felsic intrusive rocks sn w mo granites cu porphyries polymetallic and fe skarns and polymetallic vein deposits and low temperature hydrothermal systems with no direct connection to igneous activities mvt deposits epithermal au and sb deposits recent studies have shown that they formed in the triassic indosinian jurassic cretaceous early yanshanian and cretaceous late yanshanian stages indosinian deposits include major mvt pb zn ag deposits and granite related w sn deposits early yanshanian deposits are low temperature sb au and high temperature w sn and cu porphyry types many late yanshanian deposits are low temperature au as sb hg and u deposits and also include high temperature w sn polymetallic deposits the formation of these deposits is linked with a specific tectonothermal evolution and igneous activities this special issue brings together some of the latest information in eight papers that deal with the origins and tectonic environments of mineral deposits formed in these stages we anticipate that this issue will stimulate more interests in these ore deposits in south china

national geophysical data center candidate for the world digital magnetic anomaly map marine and airborne magnetic anomaly data have been collected for more than half a century providing global coverage of the earth furthermore the german champ satellite is providing increasingly accurate information on large scale magnetic anomalies the world digital magnetic anomaly map project is an international effort to integrate all available near surface and satellite magnetic anomaly data into a global map database teams of researchers were invited to produce candidate maps using a common pool of data sets here we present the national geophysical data center ngdc candidate to produce a homogeneous map the near surface data were first line leveled and then merged by least squares collocation long wavelengths were found to agree surprisingly well with independent satellite information this validates our final processing step of merging the short wavelength part of the near surface data with long wavelength satellite magnetic anomalies

he pb and s isotopic constraints on the relationship between the a type qitianling granite and the furong tin deposit hunan province china abstract the furong tin deposit located in southern hunan province china is a large newly discovered deposit with an estimated sn reserve of about 700 000 tons the deposit is spatially and temporally associated with the a type qitianling granite the 3he 4he ratios of fluid inclusions trapped in sulfides from the deposit range from 0 13 to 2 95 ra indicating a mixed crust mantle source similar to that of the qitianling granite complex lead isotopes in feldspars of the granite mostly in the range of 206pb 204pb 18 547 19 180 207pb 204pb 15 598 15 825 208pb 204pb 37 912 39 068 and in sulfide minerals 206pb 204pb 18 467 18 836 207pb 204pb 15 503 15 772 208pb 204pb 38 607 39 099 are essentially the same indicating that the lead in both was derived mainly from the upper crust with a small mantle component the sulfur isotopic values of the ore minerals vary widely ranging from 26 1 to 10 4 but are mainly within the range of 0 2 to 10 4 implying that the s in the hydrothermal fluids was derived mainly from magmatic fluids with minor contributions from the mantle sedimentary rocks and biogenic material the he pb and s isotopes all confirm that the ore forming fluids of the furong deposit were magmatic in origin derived from the melts that produced the qitianling a type granite the tin mineralization occurs within or along the margins of the qitianling granite complex and the timing of the mineralization coincides with the intrusive age of the granite the available data suggest that the emplacement of the granites and the mineralization was related to mantle upwelling and extension of the lithosphere in south china during the mesozoic

extending earthquakes reach through cascading earthquakes whatever their size can trigger other earthquakes mainshocks cause aftershocks to occur which in turn activate their own local aftershock sequences resulting in a cascade of triggering that extends the reach of the initial mainshock a long lasting difficulty is to determine which earthquakes are connected either directly or indirectly here we show that this causal structure can be found probabilistically with no a priori model nor parameterization large regional earthquakes are found to have a short direct influence in comparison to the overall aftershock sequence duration relative to these large mainshocks small earthquakes collectively have a greater effect on triggering hence cascade triggering is a key component in earthquake interactions

in situ elemental and isotopic compositions of apatite and zircon from the shuikoushan and xihuashan granitic plutons implication for jurassic granitoid related cu pb zn and w mineralization in the nanling range south china abstract the nanling range south china was world famous for hosting abundant granitoid related copper lead zinc cu pb zn polymetallic ore deposits and quartz vein type wolframite w ore deposits both of which have roughly similar formation ages predominately between 150 and 160 ma with a peak of ca 156 ma in this study accessory minerals from typical cu pb zn bearing and w bearing granitic rocks are compared in terms of elemental and isotopic compositions apatite from w bearing granites in the xihuashan pluton has lower end t value 11 9 to 8 6 than that from cu pb zn bearing granitic rocks in shuikoushan with end t value of 8 7 to 4 2 and initial 87sr 86sr ratio of 0 7097 0 7109 zircon grains in xihuashan yielded ehf t value from 14 9 to 11 4 18o value from 8 6 to 10 4 and highly variable and negative 7li value 45 8 to 3 8 whereas those in shuikoushan have relatively higher ehf t value 10 6 to 8 1 lower 18o value 8 4 9 7 and highly variable 7li values from 12 7 to 17 6 in situ hf o li isotopic compositions of zircon and sr nd isotopes of apatite suggesting that the shuikoushan granitic pluton was likely generated from dehydration melting of amphibolite from a metal fertile mafic source in the middle to lower crust whereas the xihuashan granitic pluton could be derived from partial melting of metapelite with minor amphibolite in the middle to upper crust the geochemical records in accessory minerals fingerprint that the shuikoushan granitic magma was characterized by high cl content 0 11 1 44 wt and logfo2 value fmq 1 whereas the xihuashan granitic magma have elevated f 3 51 4 80 wt and li 3 49 42 4 ppm contents with low logfo2 value