Dynamic river–groundwater exchange in the presence of a saline,semi‐confined aquifer

Understanding groundwater–surface water exchange in river banks is crucial for effective water management and a range of scientific disciplines. While there has been much research on bank storage, many studies assume idealized aquifer systems. This paper presents a field‐based study of the Tambo Catchment (southeast Australia) where the Tambo River interacts with both an unconfined aquifer containing relatively young and fresh groundwater (<500 μS/cm and <100 years old) and a semi‐confined artesian aquifer containing old and saline groundwater (electrical conductivity > 2500 μS/cm and >10 000 years old). Continuous groundwater elevation and electrical conductivity monitoring within the different aquifers and the river suggest that the degree of mixing between the two aquifers and the river varies significantly in response to changing hydrological conditions. Numerical modelling using MODFLOW and the solute transport package MT3DMS indicates that saline water in the river bank moves away from the river during flooding as hydraulic gradients reverse. This water then returns during flood recession as baseflow hydraulic gradients are re‐established. Modelling also indicates that the concentration of a simulated conservative groundwater solute can increase for up to ~34 days at distances of 20 and 40 m from the river in response to flood events approximately 10 m in height. For the same flood event, simulated solute concentrations within 10 m of the river increase for only ~15 days as the infiltrating low‐salinity river water drives groundwater dilution. Average groundwater fluxes to the river stretch estimated using Darcy's law were 7 m³/m/day compared with 26 and 3 m³/m/day for the same periods via mass balance using Radon (²²²Rn) and chloride (Cl), respectively. The study shows that by coupling numerical modelling with continuous groundwater–surface water monitoring, the transient nature of bank storage can be evaluated, leading to a better understanding of the hydrological system and better interpretation of hydrochemical data. Copyright © 2015 John Wiley & Sons, Ltd.     

The exponential decline in saturated hydraulic conductivity with depth: a novel method for exploring its effect on water flow paths and transit time distribution

The strong vertical gradient in soil and subsoil saturated hydraulic conductivity is characteristic feature of the hydrology of catchments. Despite the potential importance of these strong gradients, they have proven difficult to model using robust physically based schemes. This has hampered the testing of hypotheses about the implications of such vertical gradients for subsurface flow paths, residence times and transit time distribution. Here we present a general semi‐analytical solution for the simulation of 2D steady‐state saturated‐unsaturated flow in hillslopes with saturated hydraulic conductivity that declines exponentially with depth. The grid‐free solution satisfies mass balance exactly over the entire saturated and unsaturated zones. The new method provides continuous solutions for head, flow and velocity in both saturated and unsaturated zones without any interpolation process as is common in discrete numerical schemes. This solution efficiently generates flow pathlines and transit time distributions in hillslopes with the assumption of depth‐varying saturated hydraulic conductivity. The model outputs reveal the pronounced effect that changing the strength of the exponential decline in saturated hydraulic conductivity has on the flow pathlines, residence time and transit time distribution. This new steady‐state model may be useful to others for posing hypotheses about how different depth functions for hydraulic conductivity influence catchment hydrological response. Copyright © 2016 John Wiley & Sons, Ltd.     

A discontinuous finite element suspended sediment transport model for water quality assessments in river networks

Suspended sediment plays an important role in the distribution and transport of many pollutants (such as radionuclides) in rivers. Pollutants may adsorb on fine suspended particles (e.g. clay) and spread according to the suspended sediment movement. Hence, the simulation of the suspended sediment mechanism is indispensable for realistic transport modelling. This paper presents and tests a simple mathematical model for predicting the suspended sediment transport in river networks. The model is based on the van Rijn suspended load formula and the advection–diffusion equation with a source or sink term that represents the erosion or deposition fluxes. The transport equation is solved numerically with the discontinuous finite element method. The model evaluation was performed in two steps, first by comparing model simulations with the measured suspended sediment concentrations in the Grote Nete–Molse Nete River in Belgium, and second by a model intercomparison with the sediment transport model NST MIKE 11. The simulations reflect the measurements with a Nash‐Sutcliffe model efficiency of 0.6, while the efficiency between the proposed model and the NST MIKE 11 simulations is 0.96. Both evaluations indicate that the proposed sediment transport model, that is sufficiently simple to be practical, is providing realistic results.     

Tropical precipitation anomalies and d‐excess evolution during El Niño 2014‐16

The last 2014‐16 El Niño event was among the three strongest episodes on record. El Niño considerably changes annual and seasonal precipitation across the tropics. Here, we present a unique stable isotope data set of daily precipitation collected in Costa Rica prior to, during, and after El Niño 2014‐16, in combination with Lagrangian moisture source and precipitation anomaly diagnostics. δ²H composition ranged from ‐129.4 to +18.1 (‰) while δ¹⁸O ranged from ‐17.3 to +1.0 (‰). No significant difference was observed among δ¹⁸O (P=0.186) and δ²H (P=0.664) mean annual compositions. However, mean annual d‐excess showed a significant decreasing trend (from +13.3 to +8.7 ‰) (P<0.001) with values ranging from +26.6 to ‐13.9 ‰ prior to and during the El Niño evolution. The latter decrease in d‐excess can be partly explained by an enhanced moisture flux convergence across the southeastern Caribbean Sea coupled with moisture transport from northern South America by means of an increased Caribbean Low Level Jet regime. During 2014‐15, precipitation deficit across the Pacific domain averaged 46% resulting in a very severe drought; while a 94% precipitation surplus was observed in the Caribbean domain. Understanding these regional moisture transport mechanisms during a strong El Niño event may contribute to a) better understanding of precipitation anomalies in the tropics and b) re‐evaluate past stable isotope interpretations of ENSO events in paleoclimatic archives within the Central America region.     

相对湿度和PM2.5浓度对大气能见度的影响研究:基于小时资料的多站对比分析

为了全面分析浙江省不同区域能见度变化基本特征及影响机理,基于杭州、宁波、温州3个国家基本气象站2013-2014年逐时能见度观测资料,比较分析了3市能见度变化的基本特征。发现3市不同等级能见度出现频率基本一致,随着能见度等级的提高,出现频率逐渐降低;从能见度的日变化来看,07时（北京时）前后最低,之后缓慢上升,14-15时达到最高,随后逐渐下降;全年有两个能见度较低时段,分别出现在12月-次年2月和5-6月;总体而言,宁波能见度最优,杭州和温州大体相当。功率谱分析结果表明,3市能见度均有显著的日周期,高频波段呈现出多个显著谱峰,低频波段存在若干显著谱峰。进一步开展机理分析,发现相对湿度和PM_2.5浓度是调制大气能见度的关键因子,相对湿度增大、PM_2.5浓度升高导致能见度降低。在同一相对湿度等级下,初始阶段能见度随PM_2.5浓度的升高迅速降低,到达“拐点”之后降低速率趋于缓慢。在同一PM_2.5浓度水平下,相对湿度越大,能见度越低,说明水汽对能见度也有重要影响。基于相对湿度和PM_2.5浓度两个因子,采用非线性拟合方案构建了大气能见度定量统计模型,总体而言模型拟合效果较好。最后针对研究中存在的不足和未来值得进一步发掘的科学问题进行了讨论。     

中国中东部冬季霾日的形成与东亚冬季风和大气湿度的关系

利用1961-2013年中国地面台站长期观测资料和同期NCEP/NCAR再分析资料,以华北、江淮和华南为研究区,分析了中国中东部冬季霾日的形成与东亚冬季风以及大气湿度的关系。结果表明:(1)冬季霾日与东亚冬季风强度成显著的负相关。首先,东亚冬季风强度的减弱使得地面风速减小,进而导致冬季霾日增多。其中,华北7-8 m/s最大风速日数和江淮6-8 m/s最大风速日数的减少,及华南≤2 m/s最大风速日数的增多对各区冬季霾日的增多作用较大。其次,东亚冬季风减弱引起冬季气温的持续升高,易导致冬季霾日的增多,这在华北地区较之在江淮和华南更为明显。(2)由于气候变暖,冬季气温升高,使得近地面相对湿度减小。在江淮和华南地区,冬季霾日的增多与近地面相对湿度的减小显著相关,而在华北地区这种相关较弱。(3)冬季气温升高也有利于大气层结稳定度的增强,3个区域冬季霾日的增多均与大气层结稳定度的增强显著相关,特别是与对流层中低层(850-500 hPa)大气饱和度的降低显著相关。(4)冬季霾日数变化与区域水汽输送关系密切。其中,华北地区的冬季霾日数与水汽总收入成显著正相关,江淮地区与纬向水汽收入成显著正相关,与经向水汽收入成显著负相关,华南地区与经向水汽收入成显著负相关。     

复杂地形城市多尺度气流对大气污染影响的研究进展

复杂地形城市一般建立在山地、丘陵、沿海地带。复杂地形下的大气污染传输、扩散机制是一个复杂的问题。本文依据城市地形地貌将复杂地形归类为河谷地形、三面环山临海地形、盆地地形、马蹄型地形和峡口地形。结合国内外对这五类地形下城市污染物传输扩散及污染形成机制的研究成果,介绍了不同地理位置、不同复杂地形城市多尺度气流相互或交替作用的特点及其对污染传输扩散的影响,期望能够为其它复杂地形城市污染形成机制研究和大气污染防治提供一些借鉴。     

无锡梅雨期湿沉降综合分析

利用2008—2014年梅雨期间酸雨观测资料及2014年6月16—27日降水个例加密采样资料,结合大气污染物资料分析了近7 a无锡梅雨期酸雨特征,研究降水过程中空气污染物、p H值、电导率的变化及降水对污染物的清除作用。结果表明:无锡市梅雨期酸雨年平均p H值呈现逐年递增趋势。降水过程中,颗粒物质量浓度显著降低;气体浓度变化受其自身日变化及排放源影响大于雨水的清除作用;样品的p H值、K值每个过程变化并不一致,K值变化与颗粒物质量浓度变化大致保持一致。降水、风对颗粒物质量浓度影响大于对气体浓度的影响。长时间连续降水时,降水对颗粒污染物的清除存在极限。小时雨量在0~0.5 mm时,降水对颗粒物浓度做负清除,其值反而略有增加;小时雨量在0.6~5.0 mm时,降水对颗粒物质量浓度做正清除;小时雨量达到5.1 mm及以上时,对PM_(2.5)和PM_(2.5-10)做正清除,对PM_(10)做负清除。降水对SO_2有稀释清除作用;对NO_2的稀释作用取决于其开始浓度值;对CO、O_3的清除作用不显著。     

同化FY3B/IRAS资料对GRAPES分析场影响诊断

GRAPES(Globe and Regional Assimilation and Prediction System)变分系统能够同化常规资料和非常规卫星资料,这些被同化的资料究竟对同化系统得到的分析场有何影响,目前国内外尚未见相关的研究文献。为此,首次采用基于信息熵信号自由度思想,诊断风云三号B星(FY3B)红外分光计(Infrared Atmospheric Sounder,IRAS)资料对GRAPES分析场的影响。诊断过程中,采用数值逼近方法,统计2012年12月24日18时到2013年1月22日00时共114个时次IRAS资料对GRAPES分析场影响,结果表明,IRAS中高层通道亮温资料对GRAPES分析场影响比地表通道20观测亮温的影响大,地表通道8和9观测亮温对分析场影响较大。前24个GRAPES变分同化时次每个时次IRAS通道亮温对分析场影响的贡献率分析结果显示,高层通道和H_2O通道贡献率较大。个例分析结果表明,在同化探空资料基础上加入IRAS资料后,温度和湿度增量场变化幅度较大,表明IRAS资料对分析场有降温和增湿作用。     

AAC算法的大气校正复合改进算法

利用模拟数据,评价Autonomous Atmospheric Compensation(AAC)算法的抗噪性,认为AAC算法的抗噪性较弱。基于TASI实测数据,利用AAC算法开展反演计算时,计算结果呈现出多样性问题。结合In-scene Atmospheric Compensation(ISAC)算法中黑体像元的标定方法,提出了一种复合改进算法。首先,利用ISAC算法反演的大气透过率和路径辐射,重新计算AAC算法中大气透过率之比(Tr)和相邻两强弱吸收通道的路径辐射之差(Pd),再次,运用经验公式获得稳定的大气反演结果(大气透过率和路径辐射),有效解决了计算结果多样性的问题。利用复合改进算法,开展的温度与发射率分离实验,证明反演得到的发射率波谱更接近野外实测波谱。