杭州湾海域50年一遇波浪数值模拟研究

杭州湾外围海域岛礁众多,波浪传播机制复杂。为了了解该海域的波浪分布特征,采用MIKE 21 SW模块建立了杭州湾海域波浪数值模型,利用实测波浪资料对模型进行了验证,结果表明该模型适用于模拟杭州湾海域的波浪传播。利用1970—2002年嵊山海洋站实测极值波浪资料推算50年一遇波浪要素,将其作为模型边界条件,对杭州湾海域50年一遇的设计波浪进行了模拟,并对该海域的波浪传播特征进行了分析。结果表明,杭州湾内波高较之外海明显减小,外围众多岛礁起到了很好的遮挡保护作用。     

南海珠江口盆地东沙隆起区生物礁演化模式

南海是西太平洋最大的边缘海,其在形成演化过程中出现了有利于生物礁生长发育的环境。东沙隆起在中新世时期属于浅海环境,远离物源区,温度、盐度和水深适宜生物礁的生长和发育。通过地震资料的解释、测井、岩心资料的综合分析,并与相对海平面变化曲线对比,认为东沙隆起发育大量规模不等的生物礁,且其发育演化与相对海平面变化曲线和区域构造演化阶段可以对应起来。相对海平面的旋回性变化和复杂的构造演化过程结合,使得研究区生物礁具有较好的孔渗性,可以成为有潜力的油气储层。     

The role of silica in the hydrous metamorphism of chromite

Retrograde hydrous metamorphism has produced three types of microstructures in chromite grains from chromitites and enclosing rocks of the Tapo Ultramafic Massif (Central Peruvian Andes). In semi-massive chromitites (60–80 vol% chromite), (i) partly altered chromite with homogeneous cores surrounded by lower Al₂O₃ and MgO but higher Cr₂O₃ and FeO porous chromite with chlorite filling the pores. In serpentinites (ii) zoned chromite with homogeneous cores surrounded by extremely higher Fe₂O₃ non-porous chromite and magnetite rims, and (iii) non-porous chromite grains. The different patterns of zoning in chromite grains are the consequences of the infiltration of reducing and SiO₂-rich fluids and the subsequent heterogeneous interaction with more oxidizing and Fe-bearing fluids. During the first stage of alteration under reduced conditions magmatic chromite is dissolved meanwhile new metamorphogenic porous chromite crystallizes in equilibrium with chlorite. This reaction that involves dissolution and precipitation of minerals is here modeled thermodynamically for the first time. µSiO₂-µMgO pseudosection calculated for unaltered semi-massive chromitites at 2 kbar and 300 °C, the lowest P-T conditions inferred from the Tapo Ultramafic Massif and Marañón Complex, predicts that chromite + chlorite (i.e., partly altered chromite) is stable instead of chromite + chlorite + brucite at progressive higher µSiO₂ but lower µMgO. Our observation is twofold as it reveals that the important role of SiO₂ and MgO and the open-nature of this process. P-T-X diagrams computed using the different P-T pathways estimated for the enclosing Tapo Ultramafic Massif reproduce well the partial equilibrium sequence of mineral assemblages preserved in the chromitites. Nevertheless, it is restricted only to the P-T conditions of the metamorphic peak and that of the latest overprint. Our estimations reveal that a high fluid/rock ratio (1:40 ratio) is required to produce the microstructures and compositional changes observed in the chromitites from the Tapo Ultramafic Massif. The circulation of SiO₂-rich fluids and the mobilization of MgO from the chromitite bodies are linked with the formation of garnet amphibolites and carbonate-silica hydrothermalites (i.e., listwaenites and birbirites) in the ultramafic massif. The origin of these fluids is interpreted as a result of the dissolution of orthopyroxene and/or olivine from the metaharzburgites and metagabbros enclosed in the Tapo Ultramafic Massif.     

In situ Pb and bulk Sr isotope analysis of the Yinchanggou Pb-Zn deposit in Sichuan Province (SW China): Constraints on the origin and evolution of hydrothermal fluids

The Yinchanggou Pb-Zn deposit, located in southwestern Sichuan Province, western Yangtze Block, is stratigraphically controlled by late Ediacaran Dengying Formation and contains >0.3 Mt of metal reserves with 11 wt% Pb + Zn. A principal feature is that this deposit is structurally controlled by normal faults, whereas other typical deposits nearby (e.g. Maozu) are controlled by reverse faults. The origin of the Yinchanggou deposit is still controversial. Ore genetic models, based on conventional whole-rock isotope tracers, favor either sedimentary basin brine, magmatic water or metamorphic fluid sources. Here we use in situ Pb and bulk Sr isotope features of sulfide minerals to constrain the origin and evolution of hydrothermal fluids. The Pb isotope compositions of galena determined by femtosecond LA-MC-ICPMS are as follows: ²⁰⁶Pb/²⁰⁴Pb = 18.17–18.24, ²⁰⁷Pb/²⁰⁴Pb = 15.69–15.71, ²⁰⁸Pb/²⁰⁴Pb = 38.51–38.63. These in situ Pb isotope data overlap with bulk-chemistry Pb isotope compositions of sulfide minerals (²⁰⁶Pb/²⁰⁴Pb = 18.11–18.40, ²⁰⁷Pb/²⁰⁴Pb = 15.66–15.76, ²⁰⁸Pb/²⁰⁴Pb = 38.25–38.88), and both sets of data plotting above the Pb evolution curve of average upper continental crust. Such Pb isotope signatures suggest an upper crustal source of Pb. In addition, the coarse-grained galena in massive ore collected from the deep part has higher ²⁰⁶Pb/²⁰⁴Pb ratios (18.18–18.24) than the fine-grained galena in stockwork ore sampled from the shallow part (²⁰⁶Pb/²⁰⁴Pb = 18.17–18.19), whereas the latter has higher ²⁰⁸Pb/²⁰⁴Pb ratios (38.59–38.63) than the former (²⁰⁸Pb/²⁰⁴Pb = 38.51–38.59). However, both types of galena have the same ²⁰⁷Pb/²⁰⁴Pb ratios (15.69–15.71). This implies two independent Pb sources, and the metal Pb derived from the basement metamorphic rocks was dominant during the early phase of ore formation in the deep part, whereas the ore-hosting sedimentary rocks supplied the majority of metal Pb at the late phase in the shallow part. In addition, sphalerite separated from different levels has initial ⁸⁷Sr/⁸⁶Sr ratios ranging from 0.7101 to 0.7130, which are higher than the ore formation age-corrected ⁸⁷Sr/⁸⁶Sr ratios of country sedimentary rocks (⁸⁷Sr/⁸⁶Sr_200 Ma = 0.7083–0.7096), but are significantly lower than those of the ore formation age-corrected basement rocks (⁸⁷Sr/⁸⁶Sr_200 Ma = 0.7243–0.7288). Again, such Sr isotope signatures suggest that the above two Pb sources were involved in ore formation. Hence, the gradually mixing process of mineralizing elements and associated fluids plays a key role in the precipitation of sulfide minerals at the Yinchanggou ore district. Integrating all the evidence, we interpret the Yinchanggou deposit as a strata-bound, normal fault-controlled epigenetic deposit that formed during the late Indosinian. We also propose that the massive ore is formed earlier than the stockwork ore, and the temporal-spatial variations of Pb and Sr isotopes suggest a certain potential of ore prospecting in the deep mining area.     

Decoding Neoarchaean to Palaeoproterozoic tectonothermal events in the Rangnim Massif,North Korea: regional correlation and broader implications

Any knowledge about Archaean to Palaeoproterozoic magmatic and metamorphic events in North Korea has the potential to make a significant difference to our understanding of the early tectonic configuration and evolution of East Asia. This zircon U–Pb dating and Hf isotopic study documents multiple Neoarchaean to Palaeoproterozoic tectonothermal events from the meta-igneous complex in the Machollyong ‘Group’ of the Rangnim Massif. Two tonalitic-trondjemitic gneiss samples record a crystallization age of meta-igneous protoliths at ca. 2.56 Ga and multiple migmatization and metamorphism from 2.52 to 1.85 Ga. A meta-dolerite sample yields a magmatic emplacement age of ca. 1.83 Ga. In situ zircon Hf isotopic data indicate that most zircons from the gneiss samples have ε_Hf(t) values from –16.9 to + 3.1 and crustal model ages from 2.84 to 3.73 Ga, whereas magmatic zircons from the meta-dolerite dike record ε_Hf(t) values from –5.2 to + 5.2 and model ages of 2.05–2.44 Ga. The first-recognized Neoarchaean tonalitic-trondjemitic migmatite complex in the Rangnim Massif, together with previously identified tonalitic-trondhjemitic-granodioritic (TTG) rocks in the Rimjingang Belt and the coeval counterparts in western Gyeonggi massif, represents the oldest crustal nuclei in the Korean Peninsula. The multiple tectonothermal events in this study present reliable evidence not only for attesting to consanguinity of the basement between the Korean Peninsula and the North China Craton but also for defining the influence scope of the late Palaeoproterozoic orogeny in the Korean Peninsula.     

西南印度洋中脊三维地震探测中炮点与海底地震仪的位置校正

在海上实施三维地震探测过程中,人工震源枪阵中心与船上GPS的距离及地震探测作业中的船行方向造成炮点实际位置与预设位置有一定偏差;自由落体投放的OBS由于海流的影响会偏离原定设计位置(投放点),因此,炮点与海底地震仪(OBS)的位置校正是三维地震结构研究中的基本环节.本文利用艏向信息校正了炮点位置;采用蒙特卡洛和最小二乘法方法对海底地震仪的位置进行了校正,并探讨了直达水波曲线特征.结果表明 OBS位置一般偏离设计点1 km左右,其误差范围在20 m以内,校正后的OBS记录剖面展示了真实的记录情况.该研究结果为下一步西南印度洋的三维层析成像研究提供了坚实数据基础,同时为今后南海的三维深部地壳结构探测提供经验与借鉴.     

High resolution WRF climatic simulations for the Iberian Peninsula: Model validation

A high resolution atmospheric modelling study was done for a 20-year recent historical period. The dynamic downscaling approach adopted used the Max Planck Institute Earth System Model (MPI-ESM) to drive the WRF running in climate mode. Three online nested domains were used covering part of the North Atlantic and Europe, with a resolution 81 km, and reaching 9 km in the innermost domain which covers the Iberian Peninsula.This paper presents the validation of the WRF configuration, which is based on historic simulations between 1986 and 2005 and observational datasets of near surface temperature and precipitation for the same period. The validation was done in terms of comparison of probability distributions between model results and observations, as daily climatologies, spatially averaged inside subdomains obtained with cluster analysis of the observations, for each of the four seasons. In addition, Taylor diagrams are presented for each of the seasons and subdomains. This validation approach was repeated with the results of a new WRF simulation with the same parameterisations but forced by the ERA-Interim reanalysis. The capacity of the MPI-ESM driven WRF configuration to compare with observations and in a manner similar to the ERA-Interim driven WRF, ensures the capacity of the configuration for climate and climate change studies.Considering the difficulty to simulate extremes in long term simulations, the results showed a comfortable comparison of both models (forced by climate model and reanalysis results) with observations. This provides us confidence on the continuity of using the MPI-ESM driven WRF configuration for climate studies.     

Robust inferences on climate change patterns of precipitation extremes in the Iberian Peninsula

This work presents a methodology to make statistical significant and robust inferences on climate change from an ensemble of model simulations. This methodology is used to assess climate change projections of the Iberian daily-total precipitation for a near-future (2021–2050) and a distant-future (2069–2098) climates, relatively to a reference past climate (1961–1990).Climate changes of precipitation spatial patterns are estimated for annual and seasonal values of: (i) total amount of precipitation (PRCTOT), (ii) maximum number of consecutive dry days (CDD), (iii) maximum of total amount of 5-consecutive wet days (Rx5day), and (iv) percentage of total precipitation occurred in days with precipitation above the 95th percentile of the reference climate (R95T). Daily-total data were obtained from the multi-model ensemble of fifteen Regional Climate Model simulations provided by the European project ENSEMBLES. These regional models were driven by boundary conditions imposed by Global Climate Models that ran under the 20C3M conditions from 1961 to 2000, and under the A1B scenario, from 2001 to 2100, defined by the Special Report on Emission Scenarios of the Intergovernmental Panel on Climate Change.Non-parametric statistical methods are used for significant climate change detection: linear trends for the entire period (1961–2098) estimated by the Theil-Sen method with a statistical significance given by the Mann-Kendall test, and climate-median differences between the two future climates and the past climate with a statistical significance given by the Mann-Whitney test. Significant inferences of climate change spatial patterns are made after these non-parametric statistics of the multi-model ensemble median, while the associated uncertainties are quantified by the spread of these statistics across the multi-model ensemble. Significant and robust climate change inferences of the spatial patterns are then obtained by building the climate change patterns using only the grid points where a significant climate change is found with a predefined low uncertainty.Results highlight the importance of taking into account the spread across an ensemble of climate simulations when making inferences on climate change from the ensemble-mean or ensemble-median. This is specially true for climate projections of extreme indices such CDD and R95T. For PRCTOT, a decrease in annual precipitation over the entire peninsula is projected, specially in the north and northwest where it can decrease down to 400 mm by the middle of the 21st century. This decrease is expected to occur throughout the year except in winter. Annual CDD is projected to increase till the middle of the 21st century overall the peninsula, reaching more than three weeks in the southwest. This increase is projected to occur in summer and spring. For Rx5day, a decrease is projected to occur during spring and autumn in the major part of the peninsula, and during summer in northern Iberia. Finally, R95T is projected to decrease around 20% in northern Iberia in summer, and around 15% in the south-southwest in autumn.     

Recent trends of extreme precipitation indices in the Iberian Peninsula using observations and WRF model results

Spatial and temporal distributions of the trends of extreme precipitation indices were analysed between 1986 and 2005, over the Iberian Peninsula (IP). The knowledge of the patterns of extreme precipitation is important for impacts assessment, development of adaptation and mitigation strategies. As such, there is a growing need for a more detailed knowledge of precipitation climate change.This analysis was performed for Portuguese and Spanish observational datasets and results performed by the Weather Research and Forecast (WRF) model forced by the ERA-Interim reanalysis. Extreme precipitation indices recommended by the Expert Team for Climate Change Detection Monitoring and Indices were computed, by year and season. Then, annual and seasonal trends of the indices were estimated by Theil-Sen method and their significance was tested by the Mann-Kendal test. Additionally, a second simulation forced by the Max Planck Institute Earth System Model (MPI-ESM), was considered. This second modelling configuration was created in order to assess its performance when simulating extremes of precipitation.The annual trends estimated for the 1986–2005, from the observational datasets and from the ERA-driven simulation reveal: 1) negative statistically significant trends of the CWD index in the Galicia and in the centre of the IP; 2) positive statistically significant trends of the CDD index over the south of the IP and negative statistically significant trends in Galicia, north and centre of Portugal; 3) positive statistically significant trends of the R75p index in some regions of the north of the IP; 4) positive statistically significant trends in the R95pTOT index in the Central Mountains Chain, Leon Mountains and in the north of Portugal.Seasonally, negative statistically significant trends of the CWD index were found in Galicia, in winter and in the south of the IP, in summer. Positive statistically significant trends of the CWD index were identified in the Leon Mountains, in spring, and in Galicia, in autumn. For the CDD index, negative statistically significant trends were seen in Valencia, in the spring, and, in Galicia and Portugal (north and centre), in summer. Positive statistically significant trends of the CDD index were found: in the east of the IP, in the winter; in the Cantabrian Mountain, in the spring; and, in the south of the IP, in summer. Regarding to the R75p index, negative statistically significant trends were found in Galicia, in winter and positive statistically significant trends in the north of Portugal, in spring and in the Central Mountains Chain and north of Portugal, in autumn. For the R95pTOT index, negative statistically significant trends were found over the Sierra Cuenca and Sierra Cazorla, in winter and positive statistically significant trends were found over the Sierra Cebollera, in winter and in Castile-la Mancha region, in spring.The results of the annual and seasonal trends of the extreme precipitation indices performed for observational datasets and the simulation forced by ERA-Interim, are similar. The results obtained for the simulation forced by MPI-ESM are not satisfactory, and can be a source of criticism for the use of simulation forced by MPI-ESM in this type of climate change studies. Even for the relatively short period used, the WRF model, when properly forced is a useful tool due to the similar results of Portuguese and Spanish observational datasets and the simulation forced by ERA-Interim.     

Mineralization of pegmatites in parts of the Oban Massif, Southeastern Nigeria： A preliminary analysis

The Oban Basement Massif of southeastern Nigeria is composed of metamorphosed rocks including phyllites, schists, gneisses and amphibolites cut by pegmatitic dykes of varying length and thickness, which intruded the metamorphic rocks. Preliminary geochemical study and analysis of these pegmatites from western Oban Massif at Uyanga, Akwa Ibami, Iwuru I, Iwuru Ⅱ and Igbofia showed that the pegmatites are highly albitized. This is incon-sistent with earlier postulations that the pegmatites in this part of Nige...