Robust representation of dry cells in single-layer MODFLOW models

Dewatered or "dry" grid cells in the USGS ground water modeling software MODFLOW may cause nonphysical artifacts, trigger convergence failures, or interfere with parameter estimation. These difficulties can be avoided in two dimensions by modifying the spatial differencing scheme and the iterative procedure used to resolve nonlinearities. Specifically, the spatial differencing scheme is modified to use the water level on the upstream side of a pair of adjacent cells to calculate the saturated thickness and hence intercell conductance for the pair. This makes it possible to explicitly constrain the water level in a cell to be at or above the cell bottom elevation without introducing nonphysical artifacts. Thus constrained, all initially active cells will remain active throughout the simulation. It was necessary to replace MODFLOW's Picard iteration method with the Newton-Raphson method to achieve convergence in demanding applications involving many dry cells. Tests using a MODFLOW variant based on the new method produced results nearly identical to conventional MODFLOW in situations where conventional MODFLOW converges. The new method is extremely robust and converged in scenarios where conventional MODFLOW failed to converge, such as when almost all cells dewatered. An example application to the Edwards Aquifer in south-central Texas further demonstrates the utility of the new method.     

河西走廊中部砾质戈壁蛛形纲节肢动物群落特征

戈壁中蕴含了丰富的动植物资源,蛛形纲节肢动物是戈壁中重要的节肢动物类群,在戈壁生态系统食物网中扮演重要角色.以河西走廊中部典型砾质戈壁为研究区,2018年和2019年逐月监测了蛛形纲节肢动物种类及数量.结果表明:砾质戈壁有22个蛛形纲节肢动物种或属,Steatoda albomaculata、Xysticus spp....     

2005年南京市PM_(10)浓度时间序列分析

根据南京市2005年的PM10日均浓度和气象观测资料,利用Maflab小波分析工具,对PM10浓度的年时间序列进行分析,得出了该地区PM10日均浓度的变化规律：2—5月较大,6-9月较小,10—12月较大;并且全年共有5个突变点。结合HYSPLIT4污染扩散模式,对污染源进行追踪,结果表明：造成南京市大气污染主要有局地污染物的排放和外来污染物的输送两种类型;在垂直方向上,污染物的长距离输送一直稳定在低层。     

一个降水预报的时间序列线性平稳预测模型

本文利用兰州１９４４～１９９７年的月平均降水资料,建立了线性平稳序列的降水预测模型,该模型使用了功谱密度函数中的最大熵法（或叫全极模型）,并将特征多项式模大于１的根反射到单位圆内,再返回修正后的线性预测的系数。并对１９８６～１９９７年１１年的月降水做了预测试验,试验结果表明,该模型具有一定的预报能力,其中取１５阶预报效果较好。此方法在短期气候预测业务中,可作为台站月、季、年降水预测走向客观化、定量化方法的一种初步尝试     

A mechanism for Atlantic multidecadal variability in the Kiel Climate Model

Atlantic Multidecadal Variability (AMV) is investigated in a millennial control simulation with the Kiel Climate Model (KCM), a coupled atmosphere–ocean–sea ice model. An oscillatory mode with approximately 60 years period and characteristics similar to observations is identified with the aid of three-dimensional temperature and salinity joint empirical orthogonal function analysis. The mode explains 30 % of variability on centennial and shorter timescales in the upper 2,000 m of the North Atlantic. It is associated with changes in the Atlantic Meridional Overturning Circulation (AMOC) of ±1–2 Sv and Atlantic Sea Surface Temperature (SST) of ±0.2 °C. AMV in KCM results from an out-of-phase interaction between horizontal and vertical ocean circulation, coupled through Irminger Sea convection. Wintertime convection in this region is mainly controlled by salinity anomalies transported by the Subpolar Gyre (SPG). Increased (decreased) dense water formation in this region leads to a stronger (weaker) AMOC after 15 years, and this in turn leads to a weaker (stronger) SPG after another 15 years. The key role of salinity variations in the subpolar North Atlantic for AMV is confirmed in a 1,000 year long simulation with salinity restored to model climatology: No low frequency variations in convection are simulated, and the 60 year mode of variability is absent.     

西藏麦地卡湿地的浮游植物——2.功能群特征及其与环境因子的关系

为揭示青藏高原湿地浮游植物功能群时空分布特征与环境因子的关系,本文于2019年7月(丰水期)和2019年10月(枯水期)在麦地卡湿地调查浮游植物功能群组成和生物量,并运用Spearman相关性分析、冗余分析(RDA)和典范对应分析(CCA)等方法,对麦地卡湿地浮游植物功能群时空分布特征及其与环境因子的关系进行分析.结果表明:1)根据功能群分类方法,麦地卡湿地浮游植物可划分为27个功能群,分别为A、B、C、D、E、F、G、H1、J、K、L_M、L_O、M、MP、N、P、S1、S2、T、T_B、T_C、W1、W2、X1、X2、X3和Z;在时间维度上,丰水期有26个功能群,枯水期有23个功能群;在空间维度上,河流湿地有24个功能群,湖泊湿地有23个功能群,沼泽湿地有12个功能群.2)丰水期,以功能群A、D、L_M、L_O、MP、N、P、T、T_B和X3为主,枯水期以功能群A、D、J、L_M、L_O、MP、N、P、T和T_B为主,其中A和T功能群为2个水情期的共有绝对优势功能群,功能群J仅为枯水期的绝对优势功能群.3)生态状态指数(Q)范围为1.80～4.86,平均值为3.80,表明麦地卡湿地水质状态为“好”.4)相关性分析表明,浊度、COD_Cr、总磷和硝态氮与部分优势功能群生物量具有显著相关性;RDA分析表明,盐度和总氮是影响丰水期优势功能群生物量的主要环境因子.CCA分析表明,总氮和硝态氮是影响枯水期优势功能群生物量的主要环境因子.     

大震前地震活动平静现象一种可能的机制:断层的蠕滑-匀阻化

许多大地震前地震活动存在明显的平静现象,由于其普遍性,地震平静现象被认为是具有物理意义的地震前兆并被作为地震预报的重要依据.关于大震前地震平静现象的物理机制,国内外学者从不同角度进行了研究,其中包括用破裂强度的双峰式分布、扩容硬化、震前滑动造成的应力松弛、震源体内外介质性质的差异和相互作用等机制解释强震孕育过程中短期地震活动的平静现象.尽管在一些模型中把岩石力学实验结果作为证据之一,但来自与地震活动最相似的岩石变形声发射的实验证据并不多.因为以往大量的实验结果表明,在恒定的加载条件下(等位移速率、等应力速率等),岩石破坏失稳前声发射活动多表现为加速特征,尽管一些岩石破坏前声发射活动加速现象不明显,但也难以观测到平静现象.然而,在我们开展的一系列断层结构和介质不同的中尺度标本的声发射实验中,一些标本破坏失稳前声发射活动表现出明显的相对平静现象,据此我们对大地震前地震活动平静现象的物理机制进行了讨论.     

粤北热水岩体SHRIMP锆石U-Pb年龄及地球化学特征研究

热水岩体位于贵东岩体中部南侧,SHRIMP锆石U-Pb定年结果为162.8±5.8 Ma(MSWD=4.2),属于燕山早期岩浆活动产物。该岩体具有富硅(SiO2平均为76.32%)、富铝(A/CNK值平均为1.14)、钾大于钠(K2O/Na2O值平均为1.32)和高的P2O5含量(平均为0.41%);大离子元素富集,Ba、Sr、Ti、Nb、Ta亏损,具有高的Rb/Sr(平均为28.5)和Rb/Nb值(平均为19.8);轻、重稀土之间分馏不明显,配分模式呈水平型,Eu亏损明显(δEu值平均为0.14);εNd(t)值低(平均为-10.2),Nd模式年龄古老(平均为1777 Ma);岩浆源区物质由泥质岩组成。这些特征表明,热水岩体属于典型的S型或壳源型花岗岩。热水岩体可能是在燕山早期太平洋板块俯冲导致先存的东西向大东山-贵东断裂带重新复活的构造背景下,由古—中元古代泥质变质岩低程度部分熔融的方式形成。     

粤北铀矿定位场

粤北铀矿区是南岭铀成矿带核心地段,也是我国花岗岩型铀矿首选突破、集中发育地区。我们用１：２０万铀矿地质编图方法、地质异常区和构造岩相域理论,从新的角度总结研究了粤北铀矿系列时空结构。得出大型地质异常区,多期次地热活动中心,岩石脆性变形产生的剪切脆裂岩相域,巨型结构和放射性元素异常场共同构成铀矿定位场,据此进行了成矿预测。     

Crystallochemical effects of heat treatment on Fe-dominant tourmalines from Dolní Bory (Czech Republic) and Vlachovo (Slovakia)

Heat treatment was performed on selected Fe-dominant tourmalines to establish the nature of any change in optical properties. Two tourmaline samples from Dolní Bory, Czech Republic (TDB) and Vlachovo, Slovakia (TVL) were heated at 450, 700 and 900°C at 0.1 mPa and ambient oxidation conditions for 8 h. EMPA study shows that tourmaline from Vlachovo has schorlitic composition and tourmaline from Dolní Bory is alkali-depleted schorl to foitite. Although the black colour remained unchanged after heating at 450°C, it changed to brown at 700°C and reddish brown at 900°C. No significant changes of chemical composition were observed during heating. X-ray diffraction, infrared and Mössbauer study showed negligible oxidation of tourmaline heated at 450°C, but a significant change in iron valency state and deprotonization at 700°C. The oxidation of Fe is the main cause of tourmaline colour change, and the substitution vector for oxidation of Fe is Fe³⁺OFe _?1 ²⁺ (OH)_?1. The predicted deprotonization of OH was confirmed by infrared spectroscopy, which documented a decrease in OH groups in both samples, mainly at the V site. The oxidation of Fe is mostly significant in the Y site as documented on the compression of the Y-site octahedra and subsequent decrease in the a lattice parameter. This feature is consistent with lattice dimensions in the transition from schorl and foitite dimensions to those consistent with fluor-buergerite. The Z-site octahedra did not compressed and were not affected by heating-induced Fe oxidation, which indicates only negligible content of ^Z Fe²⁺ in original samples. After heating at 900°C, the tourmaline structure collapsed likely due to the thermally induced weakening of bonds in Y and Z octahedra, which results in amorphization of tourmaline. Subsequently, breakdown products including Fe-oxides and mullite replaced alkali-depleted amorphized tourmaline.