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.  相似文献   

青藏高原及邻区未来地震活动性趋势数值分析          下载免费PDF全文

董培育  石耀霖  程惠红  乔学军 《地球物理学报》2020,63(3):1155-1169

区域中长期地震危险性数值分析研究,需要对其初始构造应力场有所了解,但目前以及未来一段时期内仍无法直接观测到深部孕震层区域的应力场状况.本文首先基于岩石库仑-摩尔破裂准则,利用青藏高原及邻区百年历史范围内的强震信息,来反演估算该区域的初始应力场.然后,考虑区域构造应力加载及强震造成的应力扰动共同作用,重现了历史强震的发展过程.然而对于初始应力场的反演估算,本文仅能给出区域其上下限的极限值,并不能唯一确定.因此,采用Monte Carlo随机法,进行大量独立的随机试验计算,生成数千种有差异的区域初始应力场模型,且保证每种模型都能令历史强震有序发生,但未来应力场演化过程不尽相同.最后,将数千种模型在未来时间段内的危险性预测结果集成为数理统计结果,据此给出了区域未来的地震危险性概率分布图.初步结果显示未来强震危险性概率较高地区集中在巴颜喀拉块体边界及鲜水河断裂带地区.  相似文献   

A nearly analytic symplectic partitioned Runge–Kutta method based on a locally one-dimensional technique for solving two-dimensional acoustic wave equations     

Chol Sim  Chunyou Sun  Nam Yun 《Geophysical Prospecting》2020,68(4):1253-1269

In this paper, we develop a new nearly analytic symplectic partitioned Runge–Kutta method based on locally one-dimensional technique for numerically solving two-dimensional acoustic wave equations. We first split two-dimensional acoustic wave equation into the local one-dimensional equations and transform each of the split equations into a Hamiltonian system. Then, we use both a nearly analytic discrete operator and a central difference operator to approximate the high-order spatial differential operators, which implies the symmetry of the discretized spatial differential operators, and we employ the partitioned second-order symplectic Runge–Kutta method to numerically solve the resulted semi-discrete Hamiltonian ordinary differential equations, which results in fully discretized scheme is symplectic unlike conventional nearly analytic symplectic partitioned Runge–Kutta methods. Theoretical analyses show that the nearly analytic symplectic partitioned Runge–Kutta method based on locally one-dimensional technique exhibits great higher stability limits and less numerical dispersion than the nearly analytic symplectic partitioned Runge–Kutta method. Numerical experiments are conducted to verify advantages of the nearly analytic symplectic partitioned Runge–Kutta method based on locally one-dimensional technique, such as their computational efficiency, stability, numerical dispersion and long-term calculation capability.  相似文献   

Fish assemblages respond to forest cover in small Amazonian basins     

《Limnologica》2020

The removal of native forest affects stream characteristics, processes, and organisms at the local scale. We compared the structure of fish assemblages between microbasins impacted by deforestation and those in pristine condition in the Amazonian Machado River basin, Brazil. Fish were collected with seine and dip nets along an 80-m stretch of 28 streams. At each site, we recorded physical, chemical, and land-water ecotone variables. We collected 6,586 specimens of 109 species, being 39 and 18 of them exclusively of forested and deforested streams, respectively. Non-significant differences were found for abundance and species richness between forested and deforested streams. A total of four main trophic groups were identified. Carnivores were more abundant in forested streams, whereas herbivores, omnivores and detritivores species were the most abundant in deforested streams. The deforested streams showed higher abundance and richness of algae and periphyton consumers, while forested streams presented higher abundance and richness of invertebrate consumers. Forested streams presented longer foodchains, higher occurrence and abundance of species that have more specialized habits and are intolerant to degraded environments, whereas generalist and tolerant species predominated in deforested streams. We conclude that species composition in Amazonian streams predictably responds to the degree of forest cover.  相似文献   

Coastal ridge constructive processes at a multi-decadal scale in Barreta Island (southern Portugal)     

Xabier Herrero  Susana Costas  Katerina Kombiadou 《地球表面变化过程与地形》2020,45(2):411-423

Multiple ridges across prograding coasts may display variable geometries, commonly expressed through varying elevations. Changes in ridge elevation have been traditionally related to the occurrence of fluctuating progradation rates, which might, in turn, be driven by shifting environmental conditions. Here, we explore the geometry and growth mechanisms of multiple ridges, generated at Barreta Island (Ria Formosa, southern Portugal), as a consequence of the rapid progradation of the island over the last 70 years, following the artificial fixation of the downdrift Faro-Olhão inlet with jetties in 1955. The variability in the morphology of these features was analysed in combination with available wind and wave data, in order to better distinguish growth mechanisms and understand the main parameters determining the final geometry of the observed ridges. The results suggest that (1) most of the identified ridges fall in the beach ridge classification, as they have been mostly built by marine processes, and (2) the parameters derived from, or closely related to wave climate variability (e.g. progradation rates, storm occurrence) can jointly explain most of the observed morphological changes, while aeolian processes played a secondary role. Indeed, ridge geometry appears mainly controlled by progradation rates, with higher ridges associated with lower progradation rates. Progradation rate, in turn, is mostly related to longshore wave power, storminess, and the occurrence storm groups. Yet, the final configuration of ridges can also be affected by runup levels and onshore winds. Therefore, establishing the relation between ridge geometry and wave climate is not a straightforward task, because of the complex processes and interactions that control coastal morphodynamics. © 2019 John Wiley & Sons, Ltd.  相似文献   

Simultaneous water and vapour transfer in an unsaturated mudstone in badlands terrain,southwestern Taiwan     

Kohei Higuchi  Masahiro Chigira  Der-Her Lee  Jian-Hong Wu 《地球表面变化过程与地形》2020,45(13):3323-3335

The simultaneous transfer of pore fluid and vapour was studied in the unsaturated shallow subsurface of a Plio-Pleistocene marine mudstone badland slope in southwestern Taiwan during the dry season using field monitoring data and numerical simulations. Data from field monitoring show mass-basis water contents of ~0.05 to ~0.10 that decrease towards the unsaturated ground surface and were invariant during the middle part of the dry season, except for daily fluctuations. In addition, the observed daily fluctuations in water content correlate with fluctuations in bedrock temperature, especially at depths of 2.5–5.0 cm. Periodic increases in water content occurred most notably during the day, when the bedrock temperature showed the greatest increase. Water contents then decreased to the previous state as bedrock temperature decreased during the night. Calculated vapour fluxes within the mudstone during the day increased up to 6 × 10⁻⁶–1 × 10⁻⁵ kg m⁻² s⁻¹, deriving a 0.01–0.02 increase in mass-basis water content at 2.5 cm depth for a 12-h period. This agrees with field monitoring data, suggesting that increases in water content occurred due to vapour intrusions into the bedrock. Pore water electrical conductivity (EC) showed periodic variations due to vapour intrusion, and gradually increased between the ground surface and depths of 2.5–5.0 cm. In contrast, pore water EC gradually decreased between 15 and 40 cm depth. Calculated water fluxes at depths of 2.5–40.0 cm varied from −4 × 10⁻⁶ to −2 × 10⁻⁹ kg m⁻² s⁻¹. These fluxes generated an increase in solute concentrations at the ground surface, with negative values of water flux indicating an upwards movement of water towards the surface. We show that the increase in solute content due to solute transfer from depth is highly dependent on variations in water flux with depth. © 2020 John Wiley & Sons, Ltd.  相似文献   

Characterization of sudden and sustained base flow jump hydrologic behaviour in the humid seasonal tropics of the Panama Canal Watershed     

Yanyan Cheng  Fred L. Ogden  Jianting Zhu 《水文研究》2020,34(3):569-582

Base flows are important for tropical regions with pronounced dry seasons, which are facing increasing water demands. Base flow generation, however, is one of the most challenging hydrological processes to characterize in the tropics. In many years during the May–December wet season in the Panama Canal Watershed (PCW), base flows in rivers abruptly increase. This increase persists until the start of the December–April dry season. Understanding this unusual base flow jump (BFJ) behaviour is critical to improve water provisioning in the seasonal tropics, especially during droughts and extended dry seasons. This study developed an integrated approach combining piecewise regression on cumulative average base flow and sensitivity analysis to calculate the timing and magnitude of BFJ. Rainfall, forest cover, mean land surface slope, catchment area, and estimated subsurface storage were tested as predictors for the occurrence and magnitude of the BFJs in seven subcatchments of the PCW. Sensitivity analysis on correlated predictors allowed ranking of predictor contributions due to isolated and cross-correlation effects. Correlations between observed BFJs and BFJs predicted by watershed and rainfall-related predictors were 0.92 and 0.65 for BFJ timing and magnitude, respectively. Forest cover was the second most significant predictor after cumulative rainfall for jump magnitude, owing to larger subsurface storage and groundwater recharge in forests than pastures. Catchments in the mountainous eastern PCW always generated larger jumps due to their higher rainfall and greater forest cover than the western PCW catchments. The cross-correlations between predictors contributed to more than 50% of the jump variances. The results demonstrate the importance of rainfall gradient and catchment characteristics in affecting the sudden and sustained BFJs, which can help inform land management decisions intended to enhance water supplies in the tropics. This study underscores the need for more research to further understand the hydrological processes involved in the BFJ phenomenon, including better BFJ models and field characterizations, to help improve tropical ecosystem services under a changing environment.  相似文献   

Using hysteretic behaviour and hydrograph classification to identify hydrological function across the “hillslope–depression–stream” continuum in a karst catchment     

Runrun Zhang  Xi Chen  Zhicai Zhang  Chris Soulsby 《水文研究》2020,34(16):3464-3480

In cockpit karst landscapes, fluxes from upland areas contribute large volumes of water to low-lying depressions and stream flow. Hydrograph hysteresis and similarity between monitoring sites is important for understanding the space–time variability of hydrologic responses across the “hillslope–depression–stream” continuum. In this study, the hysteretic feature of hydrographs was assessed by characterizing the loop-like relationships between responses at upstream sites relative to subsurface discharge at the outlet of a small karst catchment. A classification of hydrograph responses based on the multi-scale smoothing Kernel -derived distance classifies the hydrograph responses on the basis of similarities between hillslope and depression sites, and those at the catchment outlet. Results demonstrate that the temporal and spatial variability of hydrograph hysteresis and similarity between hillslope flow and outlet stream flow can be explained by the local heterogeneity of depression aquifer. Large depression storage deficits emerging in the highly heterogeneous aquifer produce strong hysteresis and multiple relationships of upstream hydrographs relative to the outlet subsurface discharge. In contrast, when depression storage deficits are filled during consecutive rainfall events, depression hydrographs at the high permeability sites are almost synchronous or exhibit a monotonous function with the hydrographs at the outlet. This reduced hydrograph hysteresis enhances preferential flow paths in fractured rocks and conduits that can accelerate the hillslope flow to the outlet. Therefore, classification of hydrograph similarities between any upstream sites and the catchment outlet can help to identify the dominant hydrological functions in the heterogeneous karst catchment.  相似文献   

Characterizing spatial and temporal variation in 18O and 2H content of New Zealand river water for better understanding of hydrologic processes     

Jing Yang  Bruce D. Dudley  Kelsey Montgomery  Will Hodgetts 《水文研究》2020,34(26):5474-5488

Time series of hydrogen and oxygen stable isotope ratios (δ²H and δ¹⁸O) in rivers can be used to quantify groundwater contributions to streamflow, and timescales of catchment storage. However, these isotope hydrology techniques rely on distinct spatial or temporal patterns of δ²H and δ¹⁸O within the hydrologic cycle. In New Zealand, lack of understanding of spatial and temporal patterns of δ²H and δ¹⁸O of river water hinders development of regional and national-scale hydrological models. We measured δ²H and δ¹⁸O monthly, together with river flow rates at 58 locations across New Zealand over a two-year period. Results show: (a) general patterns of decreasing δ²H and δ¹⁸O with increasing latitude were altered by New Zealand's major mountain ranges; δ²H and δ¹⁸O were distinctly lower in rivers fed from higher elevation catchments, and in eastern rain-shadow areas of both islands; (b) river water δ²H and δ¹⁸O values were partly controlled by local catchment characteristics (catchment slope, PET, catchment elevation, and upstream lake area) that influence evaporation processes; (c) regional differences in evaporation caused the slope of the river water line (i.e., the relationship between δ²H and δ¹⁸O in river water) for the (warmer) North Island to be lower than that of the (cooler, mountain-dominated) South Island; (d) δ²H seasonal offsets (i.e., the difference between seasonal peak and mean values) for individual sites ranged from 0.50‰ to 5.07‰. Peak values of δ¹⁸O and δ²H were in late summer, but values peaked 1 month later at the South Island sites, likely due to greater snow-melt contributions to streamflow. Strong spatial differences in river water δ²H and δ¹⁸O caused by orographic rainfall effects and evaporation may inform studies of water mixing across landscapes. Generally distinct seasonal isotope cycles, despite the large catchment sizes of rivers studied, are encouraging for transit time analysis applications.  相似文献   

基岩倾角对盆地地震响应的影响分析          下载免费PDF全文

卓发成 《地震工程学报》2020,42(4):1035-1042

为了研究盆地基岩倾角对盆地地表动力反应的影响,选取了一个跨度2 km,深度500 m的二维成层盆地为研究对象进行有限元动力反应分析。以持时0.25 s的狄拉克脉冲作为输入地震动,探讨了P波和SV波垂直入射时,盆地基岩倾角从10°\,20°\,30°\,40°依次增加的情况下盆地地震反应。结果表明:入射波在倾斜基岩处发生波型转化,同时产生面波,在盆地中心区域发生汇聚;随着基岩倾角的增加地表中心处的卓越频率有减小的趋势。  相似文献