A two-layer model for studying 2D dissolved pollutant runoff over impermeable surfaces     

Taotao Zhang  Yang Xiao  Dongfang Liang  Hongwu Tang  Junzeng Xu  Saiyu Yuan  Nairu Wang  Bin Luan 《水文研究》2021,35(5):e14152

Dissolved pollutants in stormwater are a main contributor to water pollution in urban environments. However, many existing transport models are semi-empirical and only consider one-dimensional flows, which limit their predictive capacity. Combining the shallow water and the advection–diffusion equations, a two-dimensional physically based model is developed for dissolved pollutant transport by adopting the concept of a ‘control layer’. A series of laboratory experiments has been conducted to validate the proposed model, taking into account the effects of buildings and intermittent rainfalls. The predictions are found to be in good agreement with experimental observations, which supports the assumption that the depth of the control layer is constant. Based on the validated model, a parametric study is conducted, focusing on the characteristics of the pollutant distribution and transport rate over the depth. The hyetograph, including the intensity, duration and intermittency, of rainfall event has a significant influence on the pollutant transport rates. The depth of the control layer, rainfall intensity, surface roughness and area length are dominant factors that affect the dissolved pollutant transport. Finally, several perspectives of the new pollutant transport model are discussed. This study contributes to an in-depth understanding of the dissolved pollutant transport processes on impermeable surfaces and urban stormwater management.  相似文献   

Terrestrial LiDAR monitoring of coastal foredune evolution in managed and unmanaged systems     

I. Conery  K. Brodie  N. Spore  J. Walsh 《地球表面变化过程与地形》2020,45(4):877-892

Coastal dunes provide essential protection for infrastructure in developed regions, acting as the first line of defence against ocean-side flooding. Quantifying dune erosion, growth and recovery from storms is critical from management, resiliency and engineering with nature perspectives. This study utilizes 22 months of high-resolution terrestrial LiDAR (Riegl VZ-2000) observations to investigate the impact of management, anthropogenic modifications and four named storms on dune morphological evolution along ~100 m of an open-coast, recently nourished beach in Nags Head, NC. The influences of specific management strategies – such as fencing and plantings – were evaluated by comparing these to the morphologic response at an unmanaged control site at the USACE Field Research Facility (FRF) in Duck, NC (33 km to the north), which experienced similar environmental forcings. Various beach-dune morphological parameters were extracted (e.g. backshore-dune volume) and compared with aeolian and hydrodynamic forcing metrics between each survey interval. The results show that LiDAR is a useful tool for quantifying complex dune evolution over fine spatial and temporal scales. Under similar forcings, the managed dune grew 1.7 times faster than the unmanaged dune, due to a larger sediment supply and enhanced capture through fencing, plantings and walkovers. These factors at the managed site contributed to the welding of the incipient dune to the primary foredune over a short period of less than a year, which has been observed to take up to decades in natural systems. Storm events caused alongshore variable dune erosion primarily to the incipient dune, yet also caused significant accretion, particularly along the crest at the managed site, resulting in net dune growth. Traditional empirical Bagnold equations correlated with observed trends of backshore-dune growth but overpredicted magnitudes. This is likely because these formulations do not encompass supply-limiting factors and erosional processes. © 2019 John Wiley & Sons, Ltd.  相似文献   

Evaluating soil moisture estimation from ground-penetrating radar hyperbola fitting with respect to a systematic time-domain reflectometry data collection in a boreal podzolic agricultural field     

Chameera Illawathure  Gary Parkin  Sébastien Lambot  Lakshman Galagedara 《水文研究》2020,34(6):1428-1445

The upper 30 cm of the soil profile, which hosts the majority of the root biomass, can be considered as the shallow agricultural root zone of most temperate crops. The electromagnetic wave velocity in the soil obtained from reflection hyperbolas in ground-penetrating radar (GPR) data can be used to estimate soil moisture (SM). Finding shallow hyperbolas in a radargram and minimizing the subjective error associated with the hyperbola fitting are the main challenges in this approach. Nevertheless, we were motivated by the recent improvements of hyperbola fitting algorithms, which can reduce the subjective error and processing time. To overcome the difficulty of finding very shallow hyperbolas, we applied the hyperbola fitting method to reflections ranging from 27- to 50-cm depth using a 500-MHz centre-frequency GPR and compared the estimated moisture with vertically installed, 30-cm-long time-domain reflectometry (TDR) sensors. We also compared TDR and GPR sample areas in a 2-D plane using different GPR survey types and different hyperbola depths. SM measured with TDR and GPR were not significantly different according to Mann–Whitney's test. Our analyses showed that a root mean square error of 0.03 m³ m⁻³ was found between the two methods. In conclusion, the proposed method might be suitable to estimate SM with an acceptable accuracy within the root zone if the soil profile is fairly uniform within the application depth range.  相似文献   

A modified numerical model for moisture-salt transport in unsaturated sandy soil under evaporation     

Wei Wen  YuanMing Lai  ZheMin You 《寒旱区科学》2020,12(3):125-133

Soil salinization, caused by salt migration and accumulation underneath the soil surface, will corrode structures. To analyze the moisture-salt migration and salt precipitation in soil under evaporation conditions, a mathematical model consisting of a series of theoretical equations is briefly presented. The filling effect of precipitated salts on tortuosity factor and evaporation rate are taken into account in relevant equations. Besides, a transition equation to link the solute transport equation before and after salt precipitation is proposed. Meanwhile, a new relative humidity equation deduced from Pitzer ions model is used to modify the vapor transport flux equation. The results show that the calculated values are in good agreement with the published experimental data, especially for the simulation of volume water content and evaporation rate of Toyoura sand, which confirm the reliability and applicability of the proposed model.  相似文献   

Upper ocean high resolution regional modeling of the Arabian Sea and Bay of Bengal     

Suneet Dwivedi  Alok Kumar Mishr  Atul Srivastava 《海洋学报(英文版)》2019,38(5):32-50

In this paper, effort is made to demonstrate the quality of high-resolution regional ocean circulation model in realistically simulating the circulation and variability properties of the northern Indian Ocean(10°S–25°N,45°–100°E) covering the Arabian Sea(AS) and Bay of Bengal(BoB). The model run using the open boundary conditions is carried out at 10 km horizontal resolution and highest vertical resolution of 2 m in the upper ocean.The surface and sub-surface structure of hydrographic variables(temperature and salinity) and currents is compared against the observations during 1998–2014(17 years). In particular, the seasonal variability of the sea surface temperature, sea surface salinity, and surface currents over the model domain is studied. The highresolution model's ability in correct estimation of the spatio-temporal mixed layer depth(MLD) variability of the AS and BoB is also shown. The lowest MLD values are observed during spring(March-April-May) and highest during winter(December-January-February) seasons. The maximum MLD in the AS(BoB) during December to February reaches 150 m (67 m). On the other hand, the minimum MLD in these regions during March-April-May becomes as low as 11–12 m. The influence of wind stress, net heat flux and freshwater flux on the seasonal variability of the MLD is discussed. The physical processes controlling the seasonal cycle of sea surface temperature are investigated by carrying out mixed layer heat budget analysis. It is found that air-sea fluxes play a dominant role in the seasonal evolution of sea surface temperature of the northern Indian Ocean and the contribution of horizontal advection, vertical entrainment and diffusion processes is small. The upper ocean zonal and meridional volume transport across different sections in the AS and BoB is also computed. The seasonal variability of the transports is studied in the context of monsoonal currents.  相似文献   

An extended finite volume method and fixed-stress approach for modeling fluid injection–induced tensile opening in fractured reservoirs     

Rajdeep Deb  Patrick Jenny 《国际地质力学数值与分析法杂志》2020,44(8):1128-1144

Fluid injection–induced tensile opening is modeled using an extended finite volume method (XFVM). An embedded fracture strategy is used for the flow problem, that is, the fractures are discretized using finite volume segments without resolving the grid around them. Further, the discontinuities across fractures are modeled using special basis functions. The fracture openings due to enhanced fluid pressure and the associated shear slip due to traction free boundary condition on the fracture segments are both modeled using these special discontinuity basis functions. Mass transfer between fractures and matrix is modeled using the pressure difference. The enhancement of fracture storativity due to tensile opening leads to stronger coupling between flow and mechanics. An iterative scheme relying on the fixed-stress approach for fractures, which conserves the stress dependent terms over each iteration of the flow problem, has been introduced. Tensile opening has been simulated for single fractures embedded in two- and three-dimensional matrices. The convergence criterion for sequentially implicit fixed-stress scheme for fractures embedded in elastic media is established and has been validated numerically. Further, for 2D simulations, the effect of the matrix permeability for fracture propagation due to tensile opening has been studied.  相似文献   

Topographical controls on small‐volume pyroclastic flows     

Joan Martí  Domenico M. Doronzo  Dario Pedrazzi  Ferran Colombo 《Sedimentology》2019,66(6):2297-2317

Emplacement of small‐volume (<0·1 km³) pyroclastic flows is significantly influenced by topography. The Arico ignimbrite on Tenerife (Canary Islands) is a characteristic small‐volume pyroclastic flow deposit emplaced on high relief topography. The pyroclastic flow flowed down pre‐existing valleys on the southern slopes of the island. In proximal areas deep (up to 100 m) valleys acted as efficient conduits for the pyroclastic flow, which was mostly channelled; in this particular area the ignimbrite corresponds to a homogeneous, moderately welded deposit, consisting of flattened pumices in an abundant ashy matrix with a relatively low lithic fragment content. In intermediate zones significant changes occur in the steepness of the slope and, although still channelled, here the pyroclastic flow was influenced by hydraulic jumps. In this area, two different units can be clearly distinguished in the ignimbrite: the lower unit is composed of a lithic‐rich ground‐layer deposit that formed at the turbulent, highly concentrated head of the flow; the upper unit consists of a well welded pumice‐rich deposit that occasionally reveals a basal layer formed by shearing with the lower part. This division into two units is maintained as far as distal areas near the present‐day coastline, where the slope is very gentle or null and the ignimbrite is not channelled. The ground layer is not found in distal areas. The ignimbrite here only consists of the upper unit, which is occasionally repeated due to a surging process provoked by the lower flow speed, as the pyroclastic flow spread out of the channelled zone. A theoretical model on how topography controlled the deposition of the Arico ignimbrite is derived by interpreting the observed lithological and sedimentological variations in terms of changes in topography and bedrock morphology. This new model is of general applicability and will help to explain other deposits of similar characteristics.  相似文献   

“一带一路”背景下中国国际航空网络的空间格局及演变     

杜方叶  王姣娥  谢家昊  杜德林 《地理科学进展》2019,38(7):963-972

2013年“一带一路”倡议提出将设施联通作为合作战略重点之一,航空网络作为设施联通的重要组成部分,其地位提升到新的高度。为评估“一带一路”倡议的提出对中国国际航空网络的影响,论文选取了2013年和2018年数据,从国际航线、国际航班以及通航城市3个方面分析“一带一路”倡议提出以来中国国际航空网络的空间格局及其演变特征。研究结果表明：① 整体上来看,“一带一路”倡议并未从根本上改变中国国际航空网络的空间格局,仍以东亚、东南亚、美国以及俄罗斯等地区和国家为主,5 a来其变化主要集中在进一步加深与已有城市的联系;② 从沿线国家来看,一带一路”倡议增加了与中国直接通航的沿线国家数量,并提升了沿线国家在中国国际航空网络中的地位;③ 越来越多的境内城市参与到中国对外开放格局中,且航线和航班在北京、上海、广州的聚集程度下降,而其他城市的航线和航班占比逐渐增加。  相似文献   

频率域海洋可控源电磁垂直各向异性三维反演   总被引：1，自引：0，他引：1       下载免费PDF全文

彭荣华  胡祥云  李建慧  刘亚军 《地球物理学报》2019,62(6):2165-2175

地层宏观电性各向异性会对可控源电磁响应产生重要影响.由于海底地层电性结构常表现为电导率各向异性，若仅对海洋可控源电磁（MCSEM）数据进行常规各向同性反演，有可能无法获得准确的反演解释结果，从而削弱MCSEM技术的可靠性.本文实现了电导率垂直各向异性（VTI）条件下频率域海洋可控源电磁数据三维反演算法.其中，三维正演采用基于二次场控制方程的交错网格有限体积法，并利用直接矩阵分解技术来求解离散所得的大型线性方程组，有利于快速计算多场源的响应.反演采用具有近似二次收敛性的高斯牛顿算法对目标函数进行最优化.最后，对具有VTI电性各向异性特征的盐丘构造模型的MCSEM合成数据分别进行了电导率各向同性和垂直各向异性三维反演，结果表明：各向同性三维反演算法无法对受VTI介质影响的MCSEM数据进行正确的反演解释，而垂直各向异性三维反演能够获得更为可靠的地下电阻率结构和异常体分布，展现出对海底电性各向异性结构更为优良的反演解释能力.  相似文献   

A Holistic Approach to Upwelling and Downwelling along the South-West Coast of India     

Phiros Shah  R. Sajeev  K. J. Thara  Grinson George  Muhammad Shafeeque  S. Akash 《Marine Geodesy》2019,42(1):64-84

An attempt has been made to develop a holistic understanding of upwelling and downwelling along the south-west coast of India. The main objective was to elucidate the roles of different forcings involved in the vertical motion along this coast. The south-west coast of India was characterized by upwelling during the south-west monsoon (May to September) and by downwelling during the north-east monsoon and winter (November to February). The average vertical velocity calculated along the south-west coast from the vertical shift of the 26?°C isotherm is 0.57?m/day during upwelling and 0.698?m/day during downwelling. It was concluded that upwelling along the south-west coast of India is driven by offshore Ekman transport due to the alongshore wind, Ekman pumping, horizontal divergence of currents and by the propagation of coastally trapped waves. Whereas downwelling along the coast is driven only by convergence of currents and the propagation of coastally trapped Kelvin waves. Along the west coast of India, the downwelling-favorable Kelvin waves come from the equator and upwelling-favorable waves come from the Gulf of Mannar region.  相似文献