河南省内黄县河流地质考古研究

河流地质考古学是基于地层研究河流和考古遗址之间关系的学科。近年来我们在河南省内黄县开展的河流地质考古研究揭示了黄河复杂的演化历史,在此基础上进一步探讨了古代人类活动与周围环境的相互作用。本文主要介绍了2010~2016年我们在河南省内黄县3个全新世遗址（岸上、三杨庄和大张龙村）的地质考古工作中所取得的成果。研究区域内遗址的地层记录表明,许多考古遗址被深埋于地下,并可能影响了3000 a B.P.以来的河流沉积过程。我们在岸上遗址发掘了A、B、C、D共4处青铜时代的沟渠遗迹,这些沟渠的堆筑可能影响了后期的沉积过程并导致了遗址周边微地貌的改变;在三杨庄遗址识别出了多层不同时期的人为古土壤,包括新石器晚期、战国时期、汉代和唐代;在大张龙村发现了北宋时期黄河泛滥沉积物,其沉积过程可能受周边村落遗址的影响。根据测得的14C年代和沉积层厚度,本研究进一步对这3处遗址的沉积速率进行了估算,并与前人对华北平原沉积速率的相关研究进行了对比。结果表明,这3处遗址所显示的沉积速率自3000 a B.P.开始显著增加,与对早期历史时期黄河河道沉积速率的估算结果相吻合。因此,基于遗址的地质考古研究能够为探索人与环境的互动关系提供大量信息。未来的工作中,我们需要开展更多基于考古遗址的河流地质考古研究,以深入探讨华北平原的自然沉积过程与文明演进过程之间的关系。

     

南海张裂过程及其对晚中生代以来东南亚构造的启示——IODP建议书735-Full介绍

南海的形成揭示了大陆边缘张裂和盆地形成的复杂模式,尽管已经进行了广泛研究,但是关于基底岩石和深海盆沉积层的精确年代数据还很缺乏,这使得对南海张裂年代的估计存在很大的误差,对张裂机制和历史的各种假设没有得到验证.同时只有对南海的张裂过程有了精确地分析与刻画,才能更好地理解西太平洋边缘海盆地的形成以及它们在印支块体受印度-欧亚板块碰撞而向东南挤出、青藏高原隆升中可能起到的作用.2009年正式提交的国际综合大洋钻探计划(IODP)建议书735-Full建议在南海深海盆内的4个站位上实施钻探.这4个站位分布在南海盆地4个不同的次级构造单元上(南海东北部、西北次海盆、东部次海盆和西南次海盆),这样的站位设计会确保完成本建议书的整体研究目标,即揭示南海的张裂历史和它对晚中生代以来东南亚构造的启示.位于南海盆地最东北部的站位有助于确定该区域地壳的属性和验证古南海是否存在,位于西北次海盆的站住可能会提供南海的最早张裂年代,另外2个分别位于东部次海盆和西南次海盆的站位将重点确定2个次海盆的绝对年龄、基底矿物成分与磁化率以及2个次海盆的相对张裂次序.这些站位的水深大约在2 910～4 400 m,钻探深度预计到海底以下大约700～2 200 m,总的钻透深度为5 959 m,其中5 359 m穿透沉积层,另外600 m或400 m钻入基底.所有这些站位的位置是由已有的地球物理观测数据所确定,目前计划收集更多的地质与地球物理数据以满足IODP对井位调查数据的要求.     

Geophysical evidence for a large impact structure on the Falkland (Malvinas) Plateau

A large, roughly circular structural basin is recognised on the Falkland (Malvinas) Plateau to the NW of West Falkland (Gran Malvina) Island (S 51°00′, W 62°00′). The basin, seen in seismic‐reflection profiles and evident as a large negative gravity anomaly, has a diameter of ~250 km. The age of the basin is estimated to be Late Palaeozoic. It is completely buried by younger sediments and has no topographic expression on the sea floor. We propose that the basin and geophysical anomalies, especially the combination of a large circular negative gravity anomaly with a rim of positive anomalies, and a marked circular series of positive magnetic anomalies in the same area, may be best explained by the presence of a large buried impact structure.     

Is there a superior conceptual groundwater model structure for baseflow simulation?

Previous work has shown that streamflow response during baseflow conditions is a function of storage, but also that this functional relationship varies among seasons and catchments. Traditionally, hydrological models incorporate conceptual groundwater models consisting of linear or non‐linear storage–outflow functions. Identification of the right model structure and model parameterization however is challenging. The aim of this paper is to systematically test different model structures in a set of catchments where different aquifer types govern baseflow generation processes. Nine different two‐parameter conceptual groundwater models are applied with multi‐objective calibration to transform two different groundwater recharge series derived from a soil‐atmosphere‐vegetation transfer model into baseflow separated from streamflow data. The relative performance differences of the model structures allow to systematically improve the understanding of baseflow generation processes and to identify most appropriate model structures for different aquifer types. We found more versatile and more aquifer‐specific optimal model structures and elucidate the role of interflow, flow paths, recharge regimes and partially contributing storages. Aquifer‐specific recommendations of storage models were found for fractured and karstic aquifers, whereas large storage capacities blur the identification of superior model structures for complex and porous aquifers. A model performance matrix is presented, which highlights the joint effects of different recharge inputs, calibration criteria, model structures and aquifer types. The matrix is a guidance to improve groundwater model structures towards their representation of the dominant baseflow generation processes of specific aquifer types. Copyright © 2014 John Wiley & Sons, Ltd.     

Landscape controls on long‐term runoff in subhumid heterogeneous Boreal Plains catchments

We compared median runoff (R) and precipitation (P) relationships over 25 years from 20 mesoscale (50 to 5,000 km²) catchments on the Boreal Plains, Alberta, Canada, to understand controls on water sink and source dynamics in water‐limited, low‐relief northern environments. Long‐term catchment R and runoff efficiency (RP^?1) were low and varied spatially by over an order of magnitude (3 to 119 mm/year, 1 to 27%). Intercatchment differences were not associated with small variations in climate. The partitioning of P into evapotranspiration (ET) and R instead reflected the interplay between underlying glacial deposit texture, overlying soil‐vegetation land cover, and regional slope. Correlation and principal component analyses results show that peatland‐swamp wetlands were the major source areas of water. The lowest estimates of median annual catchment ET (321 to 395 mm) and greatest R (60 to 119 mm, 13 to 27% of P) were observed in low‐relief, peatland‐swamp dominated catchments, within both fine‐textured clay‐plain and coarse‐textured glacial deposits. In contrast, open‐water wetlands and deciduous‐mixedwood forest land covers acted as water sinks, and less catchment R was observed with increases in proportional coverage of these land covers. In catchments dominated by hummocky moraines, long‐term runoff was restricted to 10 mm/year, or 2% of P. This reflects the poor surface‐drainage networks and slightly greater regional slope of the fine‐textured glacial deposit, coupled with the large soil‐water and depression storage and higher actual ET of associated shallow open‐water marsh wetland and deciduous‐forest land covers. This intercatchment study enhances current conceptual frameworks for predicting water yield in the Boreal Plains based on the sink and source functions of glacial landforms and soil‐vegetation land covers. It offers the capability within this hydro‐geoclimatic region to design reclaimed catchments with desired hydrological functionality and associated tolerances to climate or land‐use changes and inform land management decisions based on effective catchment‐scale conceptual understanding.     

A new magnitude category disaggregation approach for temporal high‐resolution rainfall intensities

Simulation of quick runoff components such as surface runoff and associated soil erosion requires temporal high‐resolution rainfall intensities. However, these data are often not available because such measurements are costly and time consuming. Current rainfall disaggregation methods have shortcomings, especially in generating the distribution of storm events. The objectives of this study were to improve point rainfall disaggregation using a new magnitude category rainfall disaggregation approach. The procedure is introduced using a coupled disaggregation approach (Hyetos and cascade) for multisite rainfall disaggregation. The new procedure was tested with ten long‐term precipitation data sets of central Germany using summer and winter precipitation to determine seasonal variability. Results showed that dividing the rainfall amount into four daily rainfall magnitude categories (1–10, 11–25, 26–50, >50 mm) improves the simulation of high rainfall intensity (convective rainfall). The Hyetos model category approach (Hyetos_Cat) with seasonal variation performs representative to observed hourly rainfall compared with without categories on each month. The mean absolute percentage accuracy of standard deviation for hourly rainfall is 89.7% in winter and 95.6% in summer. The proposed magnitude category method applied with the coupled Hyetos_Cat–cascade approach reproduces successfully the statistical behaviour of local 10‐min rainfall intensities in terms of intermittency as well as variability. The root mean square error performance statistics for disaggregated 10‐min rainfall depth ranges from 0.20 to 2.38 mm for summer and from 0.12 to 2.82 mm for the winter season in all categories. The coupled stochastic approach preserves the statistical self‐similarity and intermittency at each magnitude category with a relatively low computational burden. Copyright © 2014 John Wiley & Sons, Ltd.     

Low‐P and high‐T metamorphism of basalts: Insights from the Sudbury impact melt sheet aureole and thermodynamic modelling

Low‐pressure and high‐temperature (LP–HT) metamorphism of basaltic rocks, which occurs globally and throughout geological time, is rarely constrained by forward phase equilibrium modelling, yet such calculations provide valuable supplementary thermometric information and constraints on anatexis that are not possible to obtain from conventional thermometry. Metabasalts along the southern margin of the Sudbury Igneous Complex (SIC) record evidence of high‐grade contact metamorphism involving partial melting and melt segregation. Peak metamorphic temperatures reached at least ~925°C at ~1–3 kbar near the SIC contact. Preservation of the peak mineral assemblage indicates that most of the generated melt escaped from these rocks leaving a residuum characterized by a plagioclase–orthopyroxene–clinopyroxene–ilmenite‐magnetite±melt assemblage. Peak temperatures reached ~875°C up to 500 m from the SIC lower contact, which marks the transition to metabasalts that only experienced incipient partial melting without melt loss. Metabasalts ~500 to 750 m from the SIC contact are characterized by a similar two‐pyroxene mineral assemblage, but typically contain abundant hornblende that overgrew clino‐ and orthopyroxene along an isobaric cooling path. Metabasalts ~750 to 1,000 m from the SIC contact are characterized by a hornblende–plagioclase–quartz–ilmenite assemblage indicating temperatures up to ~680°C. Mass balance and phase equilibria calculations indicate that anatexis resulted in 10–20% melt generation in the inner ~500 m of the aureole, with even higher degrees of melting towards the contact. Comparison of multiple models, experiments, and natural samples indicates that modelling in the Na₂O–CaO–FeO–MgO–Al₂O₃–SiO₂–H₂O–TiO₂–O₂ (NCFMASHTO) system results in the most reliable predictions for the temperature of the solidus. Incorporation of K₂O in the most recent amphibole solution model now successfully predicts dehydration melting by the coexistence of high‐Ca amphibole and silicate melt at relatively low pressures (~1.5 kbar). However, inclusion of K₂O as a system component results in prediction of the solidus at too low a temperature. Although there are discrepancies between modelling predictions and experimental results, this study demonstrates that the pseudosection approach to mafic rocks is an invaluable tool to constrain metamorphic processes at LP–HT conditions.     

High‐temperature gold‐copper extraction with chloride flux in lava tubes of Tolbachik volcano (Kamchatka)

Subvolcanic environments in supra‐subduction zones are renowned for hosting epithermal deposits that often contain electrum and native gold, including bonanza examples. This study examined mineral assemblages and processes occurring in shallow‐crust volcanic settings using recent eruption (2012–2013) of the basaltic Tolbachik volcano in the Kamchatka arc. The Tolbachik eruptive system is characterized by an extensive system of lava tubes. After cessation of magma input, the tubes maintained the flow of hot oxidized gases that episodically interacted with the lava surfaces and sulphate‐chloride precipitates from volcanic gases on these surfaces. The gas‐rock interaction had strong pyrometamorphic effects that resulted in the formation of molten salt, oxidized (tenorite, hematite, Cu‐rich magnesioferrite) and skarn‐like silicate mineral assemblages. By analogy with experimental studies, we propose that a combination of these processes was responsible for extraction of metals from the basaltic wall rocks and deposition of Cu‐, Fe‐ and Cu‐Fe‐oxides and native gold.