Solubility and Diffusivity of Water in Lunar Basalt, and Chemical Zonation in Olivine-Hosted Melt Inclusions

I report the solubility and diffusivity of water in lunar basalt and an iron-free basaltic analogue at 1 atm and 1350 °C. Such parameters are critical for understanding the degassing histories of lunar pyroclastic glasses. Solubility experiments have been conducted over a range of fO₂ conditions from three log units below to five log units above the iron-wüstite buffer (IW) and over a range of pH₂/pH₂O from 0.03 to 24. Quenched experimental glasses were analyzed by Fourier transform infrared spectroscopy (FTIR) and secondary ionization mass spectrometry (SIMS) and were found to contain up to ~420 ppm water. Results demonstrate that, under the conditions of our experiments: (1) hydroxyl is the only H-bearing species detected by FTIR; (2) the solubility of water is proportional to the square root of pH₂O in the furnace atmosphere and is independent of fO₂ and pH₂/pH₂O; (3) the solubility of water is very similar in both melt compositions; (4) the concentration of H₂ in our iron-free experiments is <3 ppm, even at oxygen fugacities as low as IW-2.3 and pH₂/pH₂O as high as 24; and (5) SIMS analyses of water in iron-rich glasses equilibrated under variable fO₂ conditions can be strongly influenced by matrix effects, even when the concentrations of water in the glasses are low. Our results can be used to constrain the entrapment pressure of the lunar melt inclusions of Hauri et al. (2011).

Diffusion experiments were conducted over a range of fO₂ conditions from IW-2.2 to IW+6.7 and over a range of pH₂/pH₂O from nominally zero to ~10. The water concentrations measured in our quenched experimental glasses by SIMS and FTIR vary from a few ppm to ~430 ppm. Water concentration gradients are well described by models in which the diffusivity of water (D^*_water) is assumed to be constant. The relationship between D^*_water and water concentration is well described by a modified speciation model (Ni et al. 2012) in which both molecular water and hydroxyl are allowed to diffuse. The success of this modified speciation model for describing our results suggests that we have resolved the diffusivity of hydroxyl in basaltic melt for the first time. Best-fit values of D^*_water for our experiments on lunar basalt vary within a factor of ~2 over a range of pH₂/pH₂O from 0.007 to 9.7, a range of fO₂ from IW-2.2 to IW+4.9, and a water concentration range from ~80 ppm to ~280 ppm. The relative insensitivity of our best-fit values of D^*_water to variations in pH₂ suggests that H₂ diffusion was not significant during degassing of the lunar glasses of Saal et al. (2008). D^*_water during dehydration and hydration in H₂/CO₂ gas mixtures are approximately the same, which supports an equilibrium boundary condition for these experiments. However, dehydration experiments into CO₂ and CO/CO₂ gas mixtures leave some scope for the importance of kinetics during dehydration into H-free environments. The value of D^*_water chosen by Saal et al. (2008) for modeling the diffusive degassing of the lunar volcanic glasses is within a factor of three of our measured value in our lunar basaltic melt at 1350 °C.

In Chapter 4 of this thesis, I document significant zonation in major, minor, trace, and volatile elements in naturally glassy olivine-hosted melt inclusions from the Siqueiros Fracture Zone and the Galapagos Islands. Components with a higher concentration in the host olivine than in the melt (MgO, FeO, Cr₂O₃, and MnO) are depleted at the edges of the zoned melt inclusions relative to their centers, whereas except for CaO, H₂O, and F, components with a lower concentration in the host olivine than in the melt (Al₂O₃, SiO₂, Na₂O, K₂O, TiO₂, S, and Cl) are enriched near the melt inclusion edges. This zonation is due to formation of an olivine-depleted boundary layer in the adjacent melt in response to cooling and crystallization of olivine on the walls of the melt inclusions concurrent with diffusive propagation of the boundary layer toward the inclusion center.

Concentration profiles of some components in the melt inclusions exhibit multicomponent diffusion effects such as uphill diffusion (CaO, FeO) or slowing of the diffusion of typically rapidly diffusing components (Na₂O, K₂O) by coupling to slow diffusing components such as SiO₂ and Al₂O₃. Concentrations of H2O and F decrease towards the edges of some of the Siqueiros melt inclusions, suggesting either that these components have been lost from the inclusions into the host olivine late in their cooling histories and/or that these components are exhibiting multicomponent diffusion effects.

A model has been developed of the time-dependent evolution of MgO concentration profiles in melt inclusions due to simultaneous depletion of MgO at the inclusion walls due to olivine growth and diffusion of MgO in the melt inclusions in response to this depletion. Observed concentration profiles were fit to this model to constrain their thermal histories. Cooling rates determined by a single-stage linear cooling model are 150–13,000 °C hr^-1 from the liquidus down to ~1000 °C, consistent with previously determined cooling rates for basaltic glasses; compositional trends with melt inclusion size observed in the Siqueiros melt inclusions are described well by this simple single-stage linear cooling model. Despite the overall success of the modeling of MgO concentration profiles using a single-stage cooling history, MgO concentration profiles in some melt inclusions are better fit by a two-stage cooling history with a slower-cooling first stage followed by a faster-cooling second stage; the inferred total duration of cooling from the liquidus down to ~1000 °C is 40 s to just over one hour.

Based on our observations and models, compositions of zoned melt inclusions (even if measured at the centers of the inclusions) will typically have been diffusively fractionated relative to the initially trapped melt; for such inclusions, the initial composition cannot be simply reconstructed based on olivine-addition calculations, so caution should be exercised in application of such reconstructions to correct for post-entrapment crystallization of olivine on inclusion walls. Off-center analyses of a melt inclusion can also give results significantly fractionated relative to simple olivine crystallization.

All melt inclusions from the Siqueiros and Galapagos sample suites exhibit zoning profiles, and this feature may be nearly universal in glassy, olivine-hosted inclusions. If so, zoning profiles in melt inclusions could be widely useful to constrain late-stage syneruptive processes and as natural diffusion experiments.

High frequency wave propagation in the earth: theory and observation

The wave-theoretical analysis of acoustic and elastic waves refracted by a spherical boundary across which both velocity and density increase abruptly and thence either increase or decrease continuously with depth is formulated in terms of the general problem of waves generated at a steady point source and scattered by a radially heterogeneous spherical body. A displacement potential representation is used for the elastic problem that results in high frequency decoupling of P-SV motion in a spherically symmetric, radially heterogeneous medium. Through the application of an earth-flattening transformation on the radial solution and the Watson transform on the sum over eigenfunctions, the solution to the spherical problem for high frequencies is expressed as a Weyl integral for the corresponding half-space problem in which the effect of boundary curvature maps into an effective positive velocity gradient. The results of both analytical and numerical evaluation of this integral can be summarized as follows for body waves in the crust and upper mantle:

1) In the special case of a critical velocity gradient (a gradient equal and opposite to the effective curvature gradient), the critically refracted wave reduces to the classical head wave for flat, homogeneous layers.

2) For gradients more negative than critical, the amplitude of the critically refracted wave decays more rapidly with distance than the classical head wave.

3) For positive, null, and gradients less negative than critical, the amplitude of the critically refracted wave decays less rapidly with distance than the classical head wave, and at sufficiently large distances, the refracted wave can be adequately described in terms of ray-theoretical diving waves. At intermediate distances from the critical point, the spectral amplitude of the refracted wave is scalloped due to multiple diving wave interference.

These theoretical results applied to published amplitude data for P-waves refracted by the major crustal and upper mantle horizons (the Pg, P*, and Pn travel-time branches) suggest that the 'granitic' upper crust, the 'basaltic' lower crust, and the mantle lid all have negative or near-critical velocity gradients in the tectonically active western United States. On the other hand, the corresponding horizons in the stable eastern United States appear to have null or slightly positive velocity gradients. The distribution of negative and positive velocity gradients correlates closely with high heat flow in tectonic regions and normal heat flow in stable regions. The velocity gradients inferred from the amplitude data are generally consistent with those inferred from ultrasonic measurements of the effects of temperature and pressure on crustal and mantle rocks and probable geothermal gradients. A notable exception is the strong positive velocity gradient in the mantle lid beneath the eastern United States (2 x 10^-3 sec^-1), which appears to require a compositional gradient to counter the effect of even a small geothermal gradient.

New seismic-refraction data were recorded along a 800 km profile extending due south from the Canadian border across the Columbia Plateau into eastern Oregon. The source for the seismic waves was a series of 20 high-energy chemical explosions detonated by the Canadian government in Greenbush Lake, British Columbia. The first arrivals recorded along this profile are on the Pn travel-time branch. In northern Washington and central Oregon their travel time is described by T = Δ/8.0 + 7.7 sec, but in the Columbia Plateau the Pn arrivals are as much as 0.9 sec early with respect to this line. An interpretation of these Pn arrivals together with later crustal arrivals suggest that the crust under the Columbia Plateau is thinner by about 10 km and has a higher average P-wave velocity than the 35-km-thick, 62-km/sec crust under the granitic-metamorphic terrain of northern Washington. A tentative interpretation of later arrivals recorded beyond 500 km from the shots suggests that a thin 8.4-km/sec horizon may be present in the upper mantle beneath the Columbia Plateau and that this horizon may form the lid to a pronounced low-velocity zone extending to a depth of about 140 km.

Geologia da formação ferrífera do Serrote do Breu e de Alto das Pedras, Alagoas

A formação ferrífera do Serrote do Breu e de Alto das Pedras localiza-se no município de Campo Grande, Estado de Alagoas e está sendo pesquisada quanto ao seu potencial como minério de ferro. Ela está inserida em um domo de embasamento arqueano no interior da Faixa Sergipana, o Domo de Jirau de Ponciano. A área de estudo é caracterizada por dois altos topográficos denominados Serrote do Breu e Alto das Pedras, sustentados pela formação ferrífera, e que representam flancos opostos de um sinformal inclinado, com direção N60W e forte mergulho para sul, e extensão total de aproximadamente 2 km. A formação ferrífera ocorre em diversas camadas intercaladas em gnaisses quartzo-feldspáticos e em rochas metamáficas. Os primeiros foram agrupados na unidade de gnaisses quartzo-feldspáticos e as últimas na suíte intrusiva máfica-ultramáfica. Na porção interior do sinformal estão quartzitos e paragnaisses agrupados na unidade metassedimentar e cortando essas unidades há uma unidade de pegmatitos. A formação ferrífera é constituída por quartzo, hematita, anfibólio e magnetita. O anfibólio é em geral cummingtonita, mas riebeckita também ocorre subordinadamente. Os teores médios de SiO2, e Fe2O3t são 43,1% e 50,7%, respectivamente, e, assim como os demais elementos maiores, são compatíveis com outras formações ferríferas do mundo. Com base na petrografia e geoquímica de elementos terras raras os gnaisses quartzo-feldspáticos foram divididos em gnaisses bandados e gnaisses com titanita. Ambos apresentam composição riolítica e trend calcio-alcalino. Já as rochas metamáficas e metaultramáficas apresentam composição basáltica a andesítica e trend toleítico completamente dissociado daquele dos gnaisses. Acredita-se que os gnaisses quartzo-feldspáticos e as rochas metamáficas e metaultramáficas tenham se formado em ambientes tectônicos totalmente distintos, com as últimas tendo se formado provavelmente intrusivas nos primeiros.

Geochemical properties of the variolitic andesite from Susong in the Dabie orogenic belt, China

The variolitic andesite from the Susong County in the Dabie Mountains implies that it was erupted in water. The mineralogy of the varioles is primarily radiate plagioclase (albite sind oligoclase), with little pyroxene, hornblende and quartz (derived from alteration). The pyroxene, hornblende and quartz are in the interstices between plagiocalse. The matrix consists of glass, hornblende, chlorite, epidote and zoisite. It is clearly subjected an extensive alteration. The andesite has an uncommon chemical composition. The SiO2 content is about 56.8%, TiO2 = 0.9%, MgO = 6.4%, Fe2O3 (tot) = 6.7%similar to 7.6%, 100Mg/(Mg + Fe) = 64.1 similar to 66.2. Mg-# is significantly high. The andesite has high abundances of large-lithophile trace elements (e.g. K, Ba. Sr, LREE), e.g. La/Nb = 5.56 similar to 6.07, low abundances of high-strength-field elements (HFSE e.g. Ta, Nb, P, Ti), particularly Ta and Nb strongly depleted. These are consistent with the characteristics of subduction-related magmas. In the spider diagram of trace elements, from Ce to right hand, the abundances of elements decrease quickly, showing a character of the continental margins. There has a strong punishment of light-rare-earth elements, with a significant diffraction of REEs (the mean value of (La/Yb)(N) is 32.84). No Eu anomaly, but there are anomaly high (La/Yb)(N) = 28.63 similar to 36.74, (La/Y)(N) = 70.33 similar to 82.4. The elements Y and Yb are depleted greatly, Y<20 g/g, Y-N = 2.74 similar to 2.84, Yb-N = 2.18 similar to 2.35. From the La-(La/Sm) diagram, the andesite is derived from partial melting. But the epsilone value of Nd is - 18.7 similar to -19.2, so that the material source may be the mantle materials affected by the crustal materials. The Nd model age is 1.9 Ga indicating that the variolitic basaltic andesite was resulted from the mantle wedge of North China block, which had the Nd model age of 2.5Ga, when the Yangze block which had the Nd model age of 1.7Ga subducted beneath it. So the variolitic andesite has characteristics of the island-are volconic rocks oil a continental basement in the vicinity of the destructive continental margin.

Major element, trace element, and Sr, Nd and Pb isotope studies of Cenozoic basalts from the South China Sea

The whole rock K-Ar ages of basalts from the South China Sea basin vary from 3.8 to 7.9 Ma, which suggest that intra-plate volcanism after the cessation of spreading of the South China Sea (SCS) is comparable to that in adjacent regions around the SCS, i.e., Leiqiong Peninsula, northern margin of the SCS, Indochina block, and so on. Based on detailed petrographic studies, we selected many fresh basaltic rocks and measured their major element, trace element, and Sr-Nd-Pb isotope compositions. Geochemical characteristics of major element and trace element show that these basaltic rocks belong to alkali basalt magma series, and are similar to OIB-type basalt. The extent of partial melting of mantle rock in source region is very low, and magma may experience crystallization differentiation and cumulation during the ascent to or storing in the high-level magma chamber. Sr-Nd-Pb isotopic data of these basaltic rocks imply an inhomogeneous mantle below the South China Sea. The nature of magma origin has a two end-member mixing model, one is EM2 (Enriched Mantle 2) which may be originated from mantle plume, the other is DMM (Depleted MORB Mantle). Pb isotopic characteristics show the Dupal anomaly in the South China Sea, and combined with newly found Dupal anomaly at Gakkel ridge in Arctic Ocean, this implies that Dupal anomaly is not only limited to South Hemisphere. In variation diagrams among Sr, Nd and Pb, the origin nature of mantle below the SCS is similar to those below Leiqiong peninsula, northern margin of the SCS and Indochina peninsula, and is different from those below north and northeast China. This study provides geochemical constraints on Hainan mantle plume.

Mineral provinces and material provenance of the surficial sediments near the Zhongsha Islands in the South China Sea

Characteristics and distribution patterns of elastic minerals (0.063 similar to 0.125 mm) in bottom sediments represent a significant indicator for the identification of the origin of sediment. One hundred and fourteen surface sediment samples, which were collected from the area near the Zhongsha Islands in the South China Sea, were analysed to identify the mineral suites and their distributions in the study area. The area can be divided into three mineral provinces: ( I) a province of biogenic minerals, which mainly originate from the Zhongsha Atoll; ( H) a province of volcanogenic minerals, which are mainly derived from local basaltic seamounts and small-scale volcanoes that are probably erupting, with some influences from the island-are volcanic region around the South China Sea; and (II) a mixed mineral province whose material source includes biogenic minerals, volcanogenic minerals and terrigenous minerals; the last province can be subdivided into a mixed mineral sub-province of the northeastern part of the study area, in which terrigenous minerals are mainly derived from China's Mainland and do not exceed 17 degrees N, and a mixed mineral sub-province of the southeastern part of the study area, in which terrigenous minerals are derived from Kalimantan and Indochina Peninsula and might be further transported into the deep sea basin through submarine canyons.

Geochemical and isotopic characteristics of volcanic rocks from the northern East China Sea shelf margin and the Okinawa Trough

Volcanic rocks both from the northern East China Sea (NECS) shelf margin and the northern Okinawa Trough are subalkaline less aluminous, and lower in High Field Strength Elements (HFSE). These rocks are higher in Large Ion Lithophile Elements (LILE), thorium and uranium contents, positive lead anomalies, negative Nb-Ta anomalies, and enrichment in Light Rare Earth Elements (LREE). Basalts from the NECS shelf margin are akin to Indian Ocean Mid-Ocean Ridge Basalt (MORB), and rhyolites from the northern Okinawa Trough have the highest Pb-207/Pb-208 and Pb-208/Pb-204 ratios. The NECS shelf margin basalts have lower Sr-87/Sr-86 ratios, epsilon(Nd) and sigma O-18 than the northern Okinawa Trough silicic rocks. According to K-40-Ar-40 isotopic ages of basalts from the NECS shelf margin, rifting of the Okinawa Trough may have been active since at least 3.65-3.86 Ma. The origin of the NECS shelf margin basalt can be explained by the interaction of melt derived from Indian Ocean MORB-like mantle with enriched subcontinental lithosphere. The basalts from both sides of the Okinawa Trough may have a similar origin during the initial rifting of the Okinawa Trough, and the formation of basaltic magmas closely relates to the thinning of continental crust. The source of the formation of the northern Okinawa Trough silicic rocks was different from that of the middle Okinawa Trough, which could have been generated by the interaction of basaltic melt with an enriched crustal component. From the Ryukyu island arc to East China, the Cenozoic basalts have apparently increasing trends of MgO contents and ratios of LREE to Heavy Rare Earth Elements (HREE), suggesting that the trace element variabilities of basalts may have been influenced by the subduction of the Philippine Sea plate, and that the effects of subduction of the Philippine Sea plate on the chemical composition of basaltic melts have had a decreasing effect from the Ryukyu island arc to East China.

Development of the East Asian summer monsoon: Evidence from the sediment record in the South China Sea since 8.5 Ma

128 samples from Ocean Drilling Program (ODP) Site 1143 in the southern South China Sea were analyzed for grain size, clay minerals, biogenic opal content and quartz in order to reconstruct changes in East Asian monsoon climate since 8.5 Ma. An abrupt change of terrigenous mass accumulation rate (MAR), clay mineral assemblage, median grain size and biogenic opal MAR about 5.2 Ma suggests that between 8.5-5.2 Ma the source of terrigenous sediment was mainly in the region of surface uplift and basaltic volcanism in southern Vietnam. A simple model of East Asian summer monsoon evolution was based on the clay/feldspar ratio, kaolinite/chlorite ratio and biogenic opal MAR. The summer monsoon has two periods of maximum strength at 8.5-7.6 Ma and 7.1-6.2 Ma. Subsequently, there was a relatively stable period at 6.2-3.5 Ma, continued intensification about 3.5-2.5 Ma, and gradually weakening after 2.5 Ma. Since I Ma the monsoon has intensified, with remarkable high-frequency and amplitude variability. Simultaneous increase in sedimentation rates at ODP Sites 1143, 1146 and 1148, as well as in MAR of terrigenous materials, quartz, feldspar and clay minerals at ODP Site 1143 at 3.5-2.5 Ma, may be the erosional response to both global climatic deterioration and the strengthening of the East Asian summer monsoon after about 3-4 Ma. (c) 2006 Elsevier B.V. All rights reserved.

大兴安岭北段晚中生代火山岩地球化学特征及其地质意义

This thesis mainly concentrates on the geochronology, prtrology, elemental geochemistry and Sr-Nd-Pb-Hf isotopic geochemistry of the volcanic rocks in north Da’Hinggan Mountain. By analyzing the data obtained in this study and data from other people, this thesis explored the age distribution, petrology and mineralogy and geochemistry characteristics of the volcanic rocks in north Da’Hinggan Mountain. Furthermore, this thesis speculated upon the source characteristics of these volcanic rocks and their implications for the tectonic evolution and crust accretion. According to the twenty Ar-Ar ages, four zircon U-Pb SHRIMP ages and two Zircon U-Pb LA-ICP-MS ages, the duration of the eruption of the Late Mesozoic volcanic rocks in north Da’Hing Mountain was about 160Ma-106Ma. Most of these volcanic rocks belong to early Cretaceous and the late Jurassic volcanic rocks are only restricted in Manzhouli. The bulk of the late Mesozoic volcanic rocks are high-K calc-alkaline rocks. Only a small portion of these volcanic rocks are shoshonites. These rocks are mainly intermediate or acid and the basic rocks usually have higher alkaline contents. Rock types are very complex in this region. These volcanic rocks have a large TiO2 variation and the Al2O3 and alkaline contents are high. From the point of mineralogy, the plagioclases in these volcanic rocks are oligoclases, andesines and labradorites, and the labradorites are more common. Most pyroxenes in these volcanic rocks are augites which belong to clinopyroxene. The source of the Late Mesozoic volcanic rocks was an enriched lithospheric mantle. When the magma en route to the surface it was contaminated by crust material slightly and had some fractional crystallization. These rocks which mainly belong to high-K calc-alkaline series were one of the results of postorogenic tectonic-magmatic activities. The upwelling in late Mesozoic supplied heat to melt the enriched lithospheric mantle which was resulted from the subduction of paleo-Asian Ocean and/or Mengol-Okhotsk ocean. These late Mesozoic volcanic rocks are also important to the upper crustal accretion of north Da’Hinggan Mountain since the late Mesozoic. These volcanics and the contemporary emplacement of granites and the basaltic underplating in combination fulfilled the crust accretion history in north Da’Hinggan Mountain in Late Mesozoic.

华北克拉通新生代玄武岩携带的地幔橄榄岩Fe同位素地球化学及含金云母地幔岩岩石学研究

In recent years, thanks to the improvement of analytical methods and the use of MC-ICP-MS, Fe isotope can be measured precisely. Fe isotope shows considerable variation both in biological and inorganic processes (from low T to high T) in nature, Therefore, Fe isotope has become one of the exciting frontier sciences and has favorable prospects of the application to the geosciences and life sciences. Based on a comprehensive review of available references in the related field, this study focuses on the development of techniques for high-precision measurement of iron isotope using MC-ICP-MS, and application of the techniques developed to study the Fe isotopes as well as major and trace element compositions of minerals (Ol, Opx, Cpx and Sp) from spinel peridotitic xenoliths from Cenozoic alkaline basalts to investigate Fe isotopic features of the lithospheric mantle beneath the North China Craton. The minerals from these xenoliths are similar to those off-cratonic peridotites world-wide, but are remarkably different from those on-cratonic peridotites and clinopyroxenes from these spinel lherzolites exhibit two types of chondrite-normalized REE patterns i.e. LREE-depleted and flat or spoon-shaped. It is noted that total abundances of REE in clinopyroxenes from these peridotites show a broad negative correlation with Cr# numbers of Cpx and Sp. The Fe isotope results show that the spinel peridotitic xenoliths have small but distinguishable Fe isotopic variations in minerals (generally Ol < Opx < Cpx) and samples, and the isotopic range in spinel is relatively large. Positive linear relationship with the ε57Fecpx/ε57Feopx ratio close to one unit has been observed between Fe isotopes of coexistent Opx and Cpx, indicating that the Cpx and Opx have generally reached Fe isotopic equilibrium. However, Fe isotopes between the Ol and Sp show apparent disequilibrium. The broadly negative correlation between mineral Fe isotopes and oxygen fugacity (fo2), metasomatic indexes such as spinel Cr#, (La/Yb) N and (La/Sm) N ratios of clinopyroxenes suggest that Fe isotopic variations in different minerals and peridotites were probably produced by melt-peridotite interaction. This study further confirms the previous observation that the lithospheric mantle has distinguishable and heterogeneous Fe isotopic variations at a scale of xenoliths. Mantle metasomatism that induces the interaction of the lithospheric mantle peridotite with metasomatic agent is a most potential mechanism for the Fe isotope fractionation in mantle peridotites. Therefore, Fe isotope could be a new and powerful tool to probe the evolution of the lithospheric mantle. We also report mineral compositions, clinopyroxene trace element concentrations and Sr-Nd isotopes for newly-discovered phlogopite-bearing spinel lherzolite and olivine clinopyroxenite xenoliths from three different localities （Hannuoba, Hebei Province; Jining Sangyitang, Inner Mongolia; Hebi, Henan Province）of the North China Craton. Systematic comparisons with phlogopite-free spinel lherzolite xenolith from the same locality reveals that the phlogopite-bearing peridotitic xenoliths have relatively higher Al2O3, CaO, Na2O, K2O, TiO2 contents and lower MgO contents than those phogopite-free counterparts. The former also has higher LREE concentrations, but relatively less radiogenic Sr-Nd isotopic ratios. This demonstrates that mantle metasomatism can not only enrich the basaltic components and trace element concentrations, but also make a decrease in Mg# of the peridotites and olivines and a relative depletion in Sr-Nd isotopes. 87Rb/86Sr-87Sr/86Sr isochrons of the phlogopite-bearing xenoliths indicate that mantle metasomatism happened in the Mesozoic and/or Cenozoic time. The metasomatic agent was derived from the asthenosphere. The result also manifests that the widespread similarity of the geochemical features such as major and trace elements and isotopic compositions in the Cenozoic lithospheric mantle beneath the North China Craton to those “oceanic” lithospheric mantle could be as a result of the ubiquitous presence of the interaction between the old refractory peridotites and the infiltrated asthenospheric melt, rather than the actually newly-accreted lithospheric mantle.

白垩纪超静磁带地球磁场古强度测定

The long-term variations in the strength of the geomagnetic field provide important constraints on the chemical-physical processes of the Earth’s interior. Especially, the intensity of the geomagnetic field during the Cretaceous normal superchron (CNS) is crucial to understand the geodynamo. But a paucity of paleointensity further limits to obtain essential knowledge interior process in the deep earth. In order to improve the experimental efficiency, this study tried to apply two new rock magnetic methods (FORC diagram and low-temperature demagnetization technique) to determine the paleointensity. First, some problems in the theory and technique in paleointensity experiments were introduced. A combined palaeomagnetic and geochronologic study was further conducted on a basaltic lava sequence at Jianchang in Liaoning Province, northeastern China. Radiometric 40Ar/39Ar dating indicates that the volcanism occurred at about 119 Ma within the marine anomaly C34n in Cretaceous normal superchron (CNS). Rock magnetic investigations show that pseudo-single domain (PSD) titanium-poor titanomagnetite is dominant in the studied lava flows. Both stepwise thermal and alternating field demagnetizations isolate the well-defined normal characteristic remanent magnetization (ChRM) in three independent lava flows with a mean direction of D/I = 6.0/51.9 degree(a95 = 12.3degree). Palaeointensity was determined using the modified Thellier method with systematic partial thermoremanent magnetization (pTRM) checks on total 72 samples, but only 10 samples exhibit ideal linear behavior on the Arai-plot in the temperature interval of 300-560 C and yield an average paleointensity of (25.8+/-1.4)uT. In addition, slopes of the line defined by the initial and the final points on the Arai-plot for the other 18 samples with characteristic PSD features give an average paleointensity estimation value of (24.8+/-1.9)uT. The consistency of these two approaches confidently demonstrates the fidelity of our results. The overall mean field strength determined using both approaches are thus estimated to be (25.2+/-0.7 )uT. This value corresponds to the virtual dipole magnetic moments (VDM) of (4.5+/-0.1)E22 Am^2, which is about half of the value of present field. This finding suggests that palaeointensity just at the onset of the CNS is characterized by a weak magnetic field.

班公湖带多不杂超大型富金斑岩铜矿床的成岩成矿年代学、岩石学及高氧化岩浆¬流体-成矿作用

Duobuza copper deposit, newly discovered typical gold-rich porphyry copper deposit with superlarge potential, is located in the Tiegelong Mesozoic tectonic -magmatic arc of the southern edge of Qiangtang block and the northern margin of Bangonghu-Nujiang suture. Quartz diorite porphyrite and grandiorite porphyry, occurred in stock, are the main ore-bearing porphyries. As the emplacement of porphyry stock, a wide range of hydrothermal alteration has developed. Within the framework of the ore district, abundant hydrothermal magnetite developed, and the relationship between precipitation of copper and gold and hydrothermal magnetite seems much close. Correspondingly, a series of veinlets and network veinlets occurred in all alteration zones. Therefore, systematic research on such a superlarge high-grade Duobuza gold-rich porphyry copper deposit can fully revealed the metallogenic characteristics of gold-rich porphyry copper deposits in this region, establish metallogenetic model and prospecting criteria, and has important practical significance on the promotion of regional exploration. In addition, this research on it can enrich metallogenic theory of strong oxidation magma-fluid to gold-rich porphyry copper deposit, and will be helpful to understand the metallogenic characteristics in early of subduction of Gangdese arc stages and its entire evolution history of the Qinghai-Tibet Plateau, the temporal and spatial distribution of ore deposits and their geodynamics settings. Northern ore body of Duobuza copper deposit have been controlled with width (north-south) about 100 ~ 400 m, length (east-west) about 1400 m, dip of 200 °, angle of dip 65 °~ 80 °. And controlled resource amount is of 2.7 million tons Cu with grade 0.94% and 13 tons Au with 0.21g/tAu. Overall features of ore body are large scale, higher grade copper, gold-rich. Ore occurred in the body of granodiotite porphyry and quartz diorite porphyrite and its contact zone with wall rock. Through the detailed mapping and field work studies, some typies of alteration are identificated as follows: albitization, biotititation, sericitization, silication, epidotization, chloritization, carbonatization, illitization, kaolinization and so on. The range of alteration is more than 10km2. Wall alteration zone can be divided into potassic alteration, moderate argillization alteration, argillization, illite-hydromuscovite or propylitization from ore-bearing porphyry center outwards, but phyllic alteration has not well developed and only sericite-quartz veins occurred in local area. Moreover, micro-fracture is development in ore district , and correspondingly a series of veinlets are development as follows: biotite vein (EB type), K-feldspar-biotite-chalcopyrite-quartz vein, magnetite-antinolite-K-feldspar vein, quartz-chalcopyrite-magnetite veins (A-type), quartz-magnetite-biotite-K-feldspar vein, chalcopyrite veinlets in potassic alteration zone; (2) chalcopyrite occurring in the center vein–quartz vein (B type), chalcopyrite veinlets, chalcopyrite-gypsum vein in intermediate argillization alteration; (3) chalcopyrite- pyrite-quartz vein, pyrite-quartz vein, chalcopyrite-gypsum veins, quartz-gypsum- molybdenite-chalcopyrite vein in argillization alteration; (4) gypsum veins, quartz-(molybdenite)-chalcopyrite vein, quartz-pyrite vein, gypsum- chalcopyrite vein, potassium feldspar veinlets, Carbonate veins, quartz-magnetite veins in the wall rock. In short, various veins are very abundant within the framework of the ore district. The results of electronic probe microscopy analysis (EMPA) indicate that Albite (Ab 91.5~99.7%) occurred along the rim of plagioclase phenocryst and fracture, and respresents the earliest stages of alteration. K-feldspar (Or 75.1~96.9%) altered plagioclase phenocryst and matrix or formed secondary potassium feldspar veinlets. Secondary biotite occurred mainly in phenocryst, matrix and veinlets, belong to magnesium-rich biotite formed under the conditions of high-oxidation magma- hydrothermal. Chloritization developed in all alteration zones and alterd iron- magnesium minerals such as biotite and hornblende and then formed chlorite veinlets. As the temperature rises, Si in the tetrahedral site of chlorite decreased, and chlorite component evolved from diabantite to ripiolite. The consistent 280℃~360℃ of formation temperature hinted that chlorite formed on the same temperature range in all alteration zones. However, formation temperature range of chlorite from the gypsum-carbonate-chlorite vein was 190℃~220℃, and it may be the product of the latest stage of hydrothermal activity. The closely relationship between biotite and rutile indicate that most of rutiles are precipitated in the process of biotite alteration and recrystallization. In addition, the V2O3 concentration of rutile from ore body in Duobuza gold-rich porphyry copper deposit is >0.4％, indicate that V concentration in rutile has important significance on marking main ore body of porphyry copper deposit. Apatites from Duobuza deposit all are F-rich. And apatite in the wall rock contained low MnO content and relatively high FeO content, which may due to the basaltic composition of the wall rocks. The MnO in apatite from altered porphyry show a strong positive correlation with FeO. In addition, Cl/F ratio of apatite from wall rock was highest, followed by the potassic alteration zone and potassic alteration zone overprinted by moderate argillization alteration was the lowest. SO2 in Apatite are in the scope of 0 to 0.66%, biotite in the apatite has the highest SO2, followed by the potassic alteration zone, potassic alteration zone overprinted by moderate argillization alteration, and the lowest in the surrounding rocks, which may be caused by the decrease of oxygen fugacity of hydrothermal fluid and S exhaust by sulfide precipitation in potassic alteration. Magnetite in the wall rock have higher Cr2O3 and lower Al2O3 features compared with altered porphyry, this may be due to basalt wall rock generally has high Cr content. And magnetites have higher TiO2 content in potassic alteration than moderate argillization alteration overprinted by potassic alteration, argillization and wall rock, suggested that its formation temperature in potassic alteration was the highest among them. The ore minerals mainly are chalcopyrite and bornite, and Au contents of chalcopyrite, bornite, and pyrite are similar with chalcopyrite slightly higher. The Eu* negative anomaly of disseminated chalcopyrite was relatively lower than chalcopyrite in veinlets. Within a drill hole, the Eu* negative anomaly of disseminated chalcopyrite was gradually larger from bottom to top. Magnetite has the same distribution model, with obvious negative Eu* abnormal, and ΣREE in great changes. The gypsum has the highest ΣREE content and the obvious negative anomaly, and biotite obviously has the Eu* abnormal. Based on the petrographic and geochemical characteristics, five series of magmatic rocks can be broadly classified; they are volcanic rocks of the normal island arc, high-Nb basaltic rocks, adakites, altered porphyry and diorite. The Sr, Nd, Hf isotopes and geochemistry of various series of magmatic rock show that they may be the result of mixing between basic magma and various degrees of acid magma coming from lower crust melted by high temperature basic underplating from partial melting of the subduction sediment melt metasomatic mantle wedge. Furthermore S isotope and Pb isotope of the sulfide, ore-bearing porphyries and volcanic rocks indicated ore-forming source is the mantle wedge metasomatied by subduction sediment melt. Oxygen fugacity of magma estimated by Fe2O3/FeO of whole rock and zircon Ce4+/Ce3+ indicated that the oxidation of basalt-andesitic rocks is higher than ore-forming porphyry, and might imply high-oxidation characteristics of underplated basic magma. Its high oxidative mechanism is likely mantle sources metasomatied by subduction sediment magma, including water and Fe3+. And such high oxidation of basaltic magma is conducive to the mantle of sulfides in the effective access to melt. And the An component of dark part within plagioclase phenocryst zoning belong to bytownite (An 74%), and its may be a result of magma composition changes refreshment by basaltic magma injection. SHRIMP zircon U-Pb and LA-ICP-MS zircon U-Pb geochronology study showed that the intrusions and volcanic rocks from Duobuza porphyry copper deposit belong to early Cretaceous magma series (126~105Ma). The magma evolution series are as follows: the earliest diorite and diorite porphyrite → ore-bearing porphyry and barren grandiorite porphyry →basaltic andesite → diorite porphyrite → andesite → basaltic andesite, and magma component shows a evolution trend from intermediate to intermediate-acid to basic. Based on the field evidences, the formation age of high-Nb basalt may be the latest. The Ar-Ar geochronology of altered secondary biotite, K-feldspar and sericite shows that the main mineralization lasting a interval of about 4 Ma, the duration limit of whole magma-hydrothermal evolution of about 6 Ma, and possibly such a long duration limit may result in the formation of Duobuza super-large copper deposit. Moreover, tectonic diagram and trace element geochemistry of volcanic rocks and diorite from Duobuza porphyry copper deposit confirm that it formed in a continental margin arc environment. Zircon U-Pb age of volcanic rocks and porphyry fall in the range of 105~121Ma, and Duobuza porphyry copper deposit locating in the north of the Bangonghu- Nujiang suture zone, suggested that Neo-Tethys ocean still subducted northward at least early Cretaceous, and its closure time should be later than 105 Ma. Three major inclusion types and ten subtypes are distinguished from quartz phenocrysts and various quartz veins. Vapor generally coexisting with brine inclusions, suggest that fluid boiling may be the main ore-forming mechanism. Raman spectrums of fluid inclusions display that the content of vapor and liquid inclusion mainly contain water, and vapor occasionally contain a little CO2. In addition, the component of liquid inclusions mainly include Cl-, SO42-, Na+, K+, a small amount of Ca2+, F-; and Cl- and Na+ show good correlation. Vapor mainly contains water, a small amount of CO2, CH4 and C2H6 and so on. The daughter minerals identified by Laman spectroscopy and SEM include gypsum, chalcopyrite, halite, sylvite, rutile, potassium feldspar, Fe-Mn-chloride and other minerals, and ore-forming fluid belong to a complex hydrothermal system containing H2O-NaCl-KClFeCl2CaCl2. H and O isotopic analysis of quartz phenocryst, vein quartz, magnetite, chlorite and gypsum from all alteration zones show that the ore-forming fluid of Duobuza gold-rich porphyry copper deposit consisted mainly of magmatic water, without addition of meteric water. Duobuza gold-rich porphyry copper deposit formed by the primary magmatic fluid (600-950C), which has high oxidation, ultra-high salinity and metallogenic element-rich, exsolution direct from the magma, and it is representative of the typical orthomagmatic end member of the porphyry continuum. Moreover, the fluid evolution model of Duobuza gold-rich porphyry copper deposit has been established. Furthermore, two key factors for formation of large Au-rich porphyry copper deposit have been summed up, which are ore-forming fluids earlier separated from magma and high oxidation magma-mineralization fluid system.

北山及邻区古生代-早中生代增生与碰撞大地构造格局

The Beishan orogenic collage locates at the triple-joint among Xinjiang, Gansu, and Inner Mongolia Provinces, at which the Siberian, Tarim and North China plates join together. It also occupies the central segment of the southern Central Asian Orogenic Belt (CAOB). The main study area in the present suty focused on the southwest part of the Beishan Mountain, which can be subdivided into four units southernward, the Mazhongshan continental block, Huaniushan Arc, Liuyuan suture zone and Shibanshan-Daqishan Arc. 1. The Huaniushan Arc was formed by northernward dipping subduction from the Orcovician to Permian, in which volcanic rocks ranging from basic to acidic with island arc affinity were widely developed. The granitiod intrusions become smaller and younger southward, whichs indicates a southward rollback of slab. The granitiod intrusions are mainly composed of I type granites, and their geochemical compositions suggest that they have affinities of island arc settings. In the early Paleozoic（440Ma-390Ma）. The Shibanshan-Daqishan Arc, however, were produced in the southernward dipping subduction system from Carboniferous to Permian. Volcanic rocks from basic to acidic rocks are typical calcic-alkaline rocks. The granitiod intrusions become smaller and younger northernward, indicating subdution with a northernward rollback. The granitiod intrusions mainly consist of I-type granites, of which geochemical data support they belong to island arc granite. 2. Two series of adakite intrusions and eruptive rocks have been discovered in the southern margin of the Huaniushan Island Arc. The older series formed during Silurian (441.7±2.5Ma) are gneiss granitoid. These adakite granites intruded the early Paleozoic Liuyuan accretionary complex, and have the same age as most of the granite intrusions in the Huanniushan Arc. Their geochemical compostions demonstrate that they were derived from partial melting of the subudcted oceanic slab. These characteristics indicate a young oceanic crust subduction in the early Paleozoic. The late stage adakites with compositons of dacites associate with Nb-enriched basalts, and island arc basalts and dacites. Their geochemistries demonstrate that the adakites are the products of subducted slab melts, whereas the Nb-enriched basalt is products of the mantle wedge which have metasomatized by adakite melts. Such a association indicates the existences of a young ocean slab subduction. 3. The Liuyuan suture zone is composed of late Paleozoic ophiolites and two series of accretionary complexes with age of early Paleozoic. The early Paleozoic accretionary complex extensively intruded by early Palozioc granites is composed of metamorphic clastics, marble, flysch, various metamorphic igneous rocks (ultramafic, mafic and dacite), and eclogite blocks, which are connected by faults. The original compositions of the rocks in this complex are highly varied, including MORB, E-MORB, arc rocks. Geochronological study indicates that they were formed during the Silurian (420.9±2.5Ma and 421.1±4.3Ma). Large-scale granitiods intruded in the accretionary complex suggest a fast growth effect at the south margin of the Huaniushan arc. During late Paleozoic, island arc were developed on this accretionary complex. The late Paleozoic ophiolite has an age of early Permian (285.7±2.2Ma), in which the rock assemblage includes ultra-mafic, gabbros, gabbros veins, massive basalts, pillow basalt, basaltic clastic breccias, and thin layer tuff, with chert on the top.These igneous rocks have both arc and MORB affinities, indicating their belonging to SSZ type ophiolite. Therefore, oceanic basins area were still existed in the Liuyuan area in the early Permian. 4. The mafic-ultramafic complexes are distributed along major faults, and composed of zoned cumulate rocks, in which peridotites are surrounded by pyroxenite, hornblendites, gabbros norite and diorite outward. They have island-arc affinities and are consistent with typical Alaska-type mafic-ultramafic complexes. The geochronological results indicate that they were formed in the early Permian. 5. The Liuyuan A-type granite were formed under post-collisional settings during the late Triassic (230.9±2.5Ma), indicating the persistence of orogenic process till the late Triassic in the study area. Geochronological results suggested that A-type granites become younger southward from the Wulungu A-type granite belt to Liuyuan A-type granite belt, which is in good agreement with the accretionary direction of the CAOB in this area, which indicate that the Liuyuan suture is the final sture of the Paleo-Asin Ocean. 6. Structural geological evidence demonstrate the W-E spreading of main tectonic terrenes. These terrenes had mainly underwent through S-N direction contraction and NE strike-faulting. The study area had experienced a S-N direction compression after the Permian, indicating a collisional event after the Permian. Based on the evidene from sedimentary geology, paleontology, and geomagnetism, our studies indicate that the orogenic process can be subdivided into five stages: (1) the pre-orogenic stage occurred before the Ordovicain; (2) the subduction orogenic stage occurred from the Orcovician to the Permian; (3) the collisional orogenic stage occurred from the late Permian to the late Triassic; (4) the post-collision stage occurred after the Triassic. The Liuyuan areas have a long and complex tectonic evolutional history, and the Liuyuan suture zone is one of the most important sutures. It is the finally suture zone of the paleo-Asian ocean in the Beishan area.

东天山后碰撞镁铁-超镁铁岩镍铜硫化物矿床PGE地球化学示踪与富集规律研究

Eastern Tianshan area, a Paleozoic complex trench-arc-basin system, experienced multi-period sudbuction and collision in geological history. A large number of Early Permian mafic-ultramafic intrusions emplaced along deep faults in post-collision extension tectonic stage and hosted a series of magmatic Cu-Ni sulfide deposits. This paper sets newly-discovered Tianyu magmatic Cu-Ni sulfide deposit related to small intrusion as an example. On basis of the study of ore-forming mechanism of Cu-Ni deposit, we compared PGE content and discussed enrichment mechanism and regularity of mafic-ultramafic rocks and ores in Jueluotage tectonic belt and Central Tianshan Massif. PGE and Cu, Ni, S contents correlate with each other. PGE is mainly controlled by S content.Samilar PGE distribution patterns of mafic-ultramafic show that complex originated from the same parental magma; Pd/Ir and Ni/Cu ratios indicate that high-Mg basaltic magma and deep sulfide segregation happened during magma evolution process. PGE and Cu-Ni ores are enriched in liquid sulfide and only individual samples completely control by monosulfide solid solution. Comparison of all control factors, early segregation of sulfide and quality of supply of magma may be the key factors leading to the Eastern Tianshan magmatic copper-nickel sulfide deposits don’t reach PGE grade, but we do not rule out the possibility of occurrence continuous mafic-ultramafic rocks and PGE-rich layer in deep.

乍得西南部中新生代火山作用与盆地演化

Located in the Central and West African, Chad, which is not well geological explored, is characterized by Mesozoic- Cenozoic intra-continental rift basins. The boreholes exposed that, during Mesozoic-Cenozoic times, volcanic activities were intense in these basins, but study on volcanic rocks is very weak, especially on those embedded in rift basins, and so far systematic and detailed work has still no carried out. Based on the project of China National Oil and Gas Exploration and Development Corporation, “The analysis of reservoir condition and the evaluation of exploration targets of seven basins in block H in Chad”, and the cooperative project between Institute of Geology and Geophysics, CAS and CNPC International (Chad) Co. Ltd., “Chronology and geochemistry studies on Mesozoic-Cenozoic volcanic rocks from southwestern Chad Basins”, systematic geochronology, geochemistry and Sr-Nd-Pb isotopic geochemistry studies on volcanic rocks from southwestern Chad basins have been done in the thesis for the first time. Detailed geochronological study using whole-rock K-Ar and Ar-Ar methods shows the mainly eruption ages of these volcanic rocks are Late Cretaceous- Paleogene. Volcanic rocks in the well Nere-1 and Figuier-1 from Doba basin are products of the Late Cretaceous which majority of the K-Ar (Ar-Ar) ages fall in the interval 95-75 Ma, whereas volcanic rocks in the well Ronier-1 from Bongor Basin and the Well Acacia-1 from Lake Chad Basin formed in the Paleogene which the ages concentrated in 66-52Ma. Two main periods of volcanic activity can be recognized in the study area, namely, the Late Cretaceous period and the Paleogene period. Volcanic activities have a general trend of south to north migration, but this may be only a local expression, and farther future studies should be carried on. Petrology study exhibits these volcanic rocks from southwestern Chad basins are mainly tholeiitic basalt. Major- and trace elements as well as Sr-Nd-Pb isotopic geochemistry studies show that the late Cretaceous and the Paleogene basalts have a definitely genetic relationship, and magmas which the basalts in southwestern Chad basins derived from were produced by fractional crystallization of olivine and clinopyroxene and had not do suffered from crustal contamination. These basalts are prominently enriched light rare earth elements (LREE), large-ion lithophile elements (LILE) and high field strength elements (HFSE) and depleted compatible elements. They have positive Ba, Pb, Sr, Nb, Ta, Zr, Hf anomalies and negative Th, U, P,Y anomalies. It is possible that the basalts from southwestern Chad basins mainly formed by mixing of depleted mantle (DM) and enriched mantle (EMⅡ) sources. The late Cretaceous basalts have higher (87Sr/86Sr)i ratios than the Paleogene basalts’, whereas have lower (143Nd/144Nd)i ratios than the latter, showing a significant temporal evolution. The mantle sources of the Late Cretaceous basalts may have more enriched mantle(EMⅡ) compositions, whereas those of the Paleogene basalts are relatively more asthenospheric mantle (DM) components. The mantle components with temporal change observed in basalts from Chad basins were probably correlated with the asthenospheric mantle upwelling and lithospheric thinning in Central and Western Africa since Mesozoic. Mesozoic- Cenozoic Volcanism in Chad basins probably is a product of intra- plate extensional stress regime, corresponded to the tectonic setting of the whole West and Central African during Cretaceous. Volcanism is closely correlated with rifting. As time passed from early period to late, the basaltic magma of Chad basins, characterized with shallower genetic depth, higher density and smaller viscosity, probably indicates the gradual strengthening evolution of the rifting. In the initial rife stage, volcanic activities are absent in the study area. Volcanic activities are basiccally corresponded with the strong extensional period of Chad basins, and the eruption of basalts was slightly lagged behind the extensional period. In the post-rift stage (30-0Ma), these basins shifted to the thermal sag phase, volcanic activities in the study area significantly decreased and then terminated.