首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
Anthropogenic modifications to the landscape, with agricultural activities being a primary driver, have resulted in significant alterations to the hydrologic cycle. Artificial drainage, including surface and subsurface drainage (tile drains), is one of the most extensive manipulations in agricultural landscapes and thus is expected to provide a distinct signature of anthropogenic modification. This study adopts a data synthesis approach in an effort to characterize the signature of artificial subsurface drainage. Daily discharge data from 24 basins across the state of Iowa, which encapsulate a range of anthropogenic modifications, are assessed using a variety of flow metrics. Results indicate that the presence of artificial subsurface drainage leads to a homogenization of landscape hydrologic response. Non‐tiled watersheds exhibit a decrease in the area‐normalized peak discharge and an increase in the baseflow ratio (baseflow/streamflow) with increases in the spatial scale, while scale invariance is apparent in tiled basins. Within‐basin variability in hydrograph recession coefficients also appears to decrease with increases in the proportion of the catchment that is artificially drained. Finally, the differences between tiled and non‐tiled landscapes disappear at scales greater than approximately 2200 km2, indicating that this may be a threshold scale for studying the effects of tile drainage. This decrease in within‐basin variability and the scale invariance of hydrologic metrics in artificially drained watersheds are attributed to the creation of a bypass flow hydrologic pathway that bypasses the complexity of the catchment travel paths. Spatial homogeneity in responses implies that it may be possible to develop more parsimonious hydrologic models for these regions. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

2.
Distributed, continuous hydrologic models promote better understanding of hydrology and enable integrated hydrologic analyses by providing a more detailed picture of water transport processes across the varying landscape. However, such models are not widely used in routine modelling practices, due in part to the extensive data input requirements, computational demands, and complexity of routing algorithms. We developed a two‐dimensional continuous hydrologic model, HYSTAR, using a time‐area method within a grid‐based spatial data model with the goal of providing an alternative way to simulate spatiotemporally varied watershed‐scale hydrologic processes. The model calculates the direct runoff hydrograph by coupling a time‐area routing scheme with a dynamic rainfall excess sub‐model implemented here using a modified curve number method with an hourly time step, explicitly considering downstream ‘reinfiltration’ of routed surface runoff. Soil moisture content is determined at each time interval based on a water balance equation, and overland and channel runoff is routed on time‐area maps, representing spatial variation in hydraulic characteristics for each time interval in a storm event. Simulating runoff hydrographs does not depend on unit hydrograph theory or on solution of the Saint Venant equation, yet retains the simplicity of a unit hydrograph approach and the capability of explicitly simulating two‐dimensional flow routing. The model provided acceptable performance in predicting daily and monthly runoff for a 6‐year period for a watershed in Virginia (USA) using readily available geographic information about the watershed landscape. Spatial and temporal variability in simulated effective runoff depth and time area maps dynamically show the areas of the watershed contributing to the direct runoff hydrograph at the outlet over time, consistent with the variable source area overland flow generation mechanism. The model offers a way to simulate watershed processes and runoff hydrographs using the time‐area method, providing a simple, efficient, and sound framework that explicitly represents mechanisms of spatially and temporally varied hydrologic processes. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

3.
The objective of this paper is to investigate the variation of geomorphology and runoff characteristics in saturated areas under different partial contributing area (PCA) conditions. Geomorphologic information and hydrologic records from two mid‐size watersheds in northern Taiwan were selected for analysis. The PCA ratio in the watershed during a storm was assumed equal to the ratio of the surface‐flow volume to the direct runoff volume from measured hydrologic data. The extents of PCA regions were then determined by using a topographic‐index threshold. Consequently, the geomorphologic factors in saturated and unsaturated areas could be calculated using a digital elevation model, and these factors could then be linked to a geomorphology‐based IUH model for runoff simulation, which can consider both the surface‐ and subsurface‐flow processes in saturated and unsaturated areas, respectively. The results show that geomorphologic characteristics in the saturated areas vary significantly with different PCA ratios especially for higher order streams. A large PCA ratio results in a sharp hydrograph because the quick surface flow dominates the runoff process, whereas the hydrologic response in a low PCA case is dominated by the delayed subsurface flow. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

4.
Spatially distributed hydrologic models can be effectively utilized for flood event simulation over basins where a complex system of reservoirs affecting the natural flow regime is present. Flood peak attenuation through mountain reservoirs can, in fact, mitigate the impact of major floods in flood‐prone areas of the lower river valley. Assessment of this effect for a complex reservoir system is performed with a spatially distributed hydrologic model where the surface runoff formation and the hydraulic routing through each reservoir and the river system are performed at a fine spatial and time resolution. The Toce River basin is presented as a case study, because of the presence of 14 active hydroelectric dams that affect the natural flow regime. A recent extreme flood event is simulated using a multi‐realization kriging method for modelling the spatial distribution of rainfall. A sensitivity analysis of the key elements of the distributed hydrologic model is also performed. The flood hydrograph attenuation is assessed. Several possible reservoir storage conditions are used to characterize the initial condition of each reservoir. The results demonstrate how a distributed hydrologic model can contribute to defining strategies for reservoir management in flood mitigation. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

5.
Terrestrial sediment yield – often measured as suspended sediment load in stream channels – commonly scales with drainage area within homogeneous land surface regions. But the effect of drainage area has not usually been recognized in comparative sediment yield analyses, rendering most comparisons of sediment yield from disparate source areas invalid. The procedure to discount scale differences for comparative purposes is presented. Mathematical scaling varies according to landscape condition and provides a physical interpretation of that condition. The results open the way for rational construction of a ‘sediment delivery ratio’. Copyright © 2017 John Wiley & Sons, Ltd.  相似文献   

6.
This study demonstrates the importance of the including and appropriately parameterizing peatlands and forestlands for basin‐scale integrated surface–subsurface models in the northern boreal forest, with particular emphasis on the Athabasca River Basin (ARB). With a long‐term water balance approach to the ARB, we investigate reasons why downstream mean annual stream flow rates are consistently higher than upstream, despite the subhumid water deficit conditions in the downstream regimes. A high‐resolution 3D variably saturated subsurface and surface water flow and evapotranspiration model of the ARB is constructed based on the bedrock and surficial geology and the spatial distribution of peatlands and their corresponding eco‐regions. Historical climate data were used to drive the model for calibration against 40‐year long‐term average surface flow and groundwater observations during the historic instrumental period. The simulation results demonstrate that at the basin‐scale, peatlands and forestlands can have a strong influence on the surface–subsurface hydrologic systems. In particular, peatlands in the midstream and downstream regimes of the ARB increase the water availability to the surface–subsurface water systems by reducing water loss through evapotranspiration. Based on the comparison of forestland evapotranspiration between observation and simulation, the overall spatial average evapotranspiration in downstream forestlands is larger than that in peatlands and thus the water contribution to the stream flow in downstream areas is relatively minor. Therefore, appropriate representation of peatlands and forestlands within the basin‐scale hydrologic model is critical to reproduce the water balance of the ARB.  相似文献   

7.
The urban environment modifies the hydrologic cycle resulting in increased runoff rates, volumes, and peak flows. Green infrastructure, which uses best management practices (BMPs), is a natural system approach used to mitigate the impacts of urbanization onto stormwater runoff. Patterns of stormwater runoff from urban environments are complex, and it is unclear how efficiently green infrastructure will improve the urban water cycle. These challenges arise from issues of scale, the merits of BMPs depend on changes to small‐scale hydrologic processes aggregated up from the neighborhood to the urban watershed. Here, we use a hyper‐resolution (1 m), physically based hydrologic model of the urban hydrologic cycle with explicit inclusion of the built environment. This model represents the changes to hydrology at the BMP scale (~1 m) and represents each individual BMP explicitly to represent response over the urban watershed. Our study varies both the percentage of BMP emplacement and their spatial location for storm events of increasing intensity in an urban watershed. We develop a metric of effectiveness that indicates a nonlinear relationship that is seen between percent BMP emplacement and storm intensity. Results indicate that BMP effectiveness varies with spatial location and that type and emplacement within the urban watershed may be more important than overall percent.  相似文献   

8.
D. Raje  P. Priya  R. Krishnan 《水文研究》2014,28(4):1874-1889
In climate‐change studies, a macroscale hydrologic model (MHM) operating over large scales can be an important tool in developing consistent hydrological variability estimates over large basins. MHMs, which can operate at coarse grid resolutions of about 1° latitude by longitude, have been used previously to study climate change impacts on the hydrology of continental scale or global river basins. They can provide a connection between global atmospheric models and water resource systems on large spatial scales and long timescales. In this study, the variable infiltration capacity (VIC) MHM is used to study large scale hydrologic impacts of climate change for Indian river basins. Large‐scale changes in runoff, evapotranspiration and soil moisture for India, as well as station‐scale changes in discharges for three major river basins with distinct climatic and geographic characteristics are examined in this study. Climate model projections for meteorological variables (precipitation, temperature and wind speed) from three general circulation models (GCMs) and three emissions scenarios are used to drive the VIC MHM. GCM projections are first interpolated to a 1° by 1° hydrologic model grid and then bias‐corrected using a quantile–quantile mapping. The VIC model is able to reproduce observed statistics for discharges in the Ganga, Narmada and Krishna basins reasonably well, even at the coarse grid resolution employed using a calibration period for years 1965–1970 and testing period from 1971–1973/1974. An increasing trend is projected for summer monsoon surface runoff, evapotranspiration and soil moisture in most central Indian river basins, whereas a decrease in runoff and soil moisture is projected for some regions in southern India, with important differences arising from GCM and scenario variability. Discharge statistics show increases in mid‐flow and low flow at Farakka station on Ganga River, increased high flows at Jamtara station upstream of Narmada, and increased high, mid‐flow and low flow for Vijayawada station on Krishna River in the future. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

9.
The objective of this study was to test the practicability of defining hydrologic response units as combinations of soil, land use and topography for modelling infiltration at the hillslope and catchment scales. In an experimental catchment in the East African Highlands (Kwalei, Tanzania), three methods of measuring infiltration were compared for their ability to capture the spatial variability of effective hydraulic conductivity: the constant head (CH) method; the tension infiltration (TI) method; and the mini‐rainfall simulation (RS) method. The three methods yielded different probability distributions of effective hydraulic conductivity and suggested different types of hydrologic response units. Independently from these measurements, the occurrence of infiltration‐excess overland flow was monitored over an area of 6 ha by means of overland flow detectors. The observed pattern of overland flow occurrence did not match any of the patterns suggested by the infiltration measurements. Instead, clusters of spots with overland flow were practically independent from field borders. Geostatistical analysis of the overland flow confirmed the absence of spatial correlation for distances over 40 m. The RS method yielded the pattern closest to the observations, probably because the method simulated better the processes that trigger infiltration‐excess overland flow, i.e. soil sealing and infiltration through macroporosity. The RS hydrologic response unit correlated significantly with observed overland flow frequency. However, the location of clusters and ‘hot spots’ of overland flow remained largely unexplained by land use, soil and topographic variables. It is concluded that using such landscape variables to define hydrologic units may create artificial boundaries that do no correspond to physical realities, especially if the stochastic component within hydrologic units is neglected. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

10.
The resolution of a digital elevation model (DEM) is a crucial factor in watershed hydrologic and environmental modelling. DEM resolution can cause significant variability in the representation of surface topography, which further affects quantification of hydrologic connectivity and simulation of hydrologic processes. The objective of this study is to examine the effects of DEM resolution on (1) surface microtopographic characteristics, (2) hydrologic connectivity, and (3) the spatial and temporal variations of hydrologic processes. A puddle‐to‐puddle modelling system was utilized for surface delineation and modelling of the puddle‐to‐puddle overland flow dynamics, surface runoff, infiltration, and unsaturated flow for nine DEM resolution scenarios of a field plot surface. Comparisons of the nine modelling scenarios demonstrated that coarser DEM resolutions tended to eliminate topographic features, reduce surface depression storage, and strengthen hydrologic connectivity and surface runoff. We found that reduction in maximum depression storage and maximum ponding area was as high as 97.56% and 76.36%, respectively, as the DEM grid size increased from 2 to 80 cm. The paired t‐test and fractal analysis demonstrated the existence of a threshold DEM resolution (10 cm for the field plot), within which the DEM‐based hydrologic modelling was effective and acceptable. The effects of DEM resolution were further evaluated for a larger surface in the Prairie Pothole Region subjected to observed rainfall events. It was found that simulations based on coarser resolution DEMs (>10 m) tended to overestimate ponded areas and underestimate runoff discharge peaks. The simulated peak discharge from the Prairie Pothole Region surface reduced by approximately 50% as the DEM resolution changed from 2 to 90 m. Fractal analysis results elucidated scale dependency of hydrologic and topographic processes. In particular, scale analysis highlighted a unique constant–threshold–power relationship between DEM scale and topographic and hydrologic parameters/variables. Not only does this finding allow one to identify threshold DEM but also further develop functional relationships for scaling to achieve valid topographic characterization as well as effective and efficient hydrologic modelling. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

11.
Amount and composition of dissolved organic matter (DOM) were evaluated for multiple, nested stream locations in a forested watershed to investigate the role of hydrologic flow paths, wetlands and drainage scale. Sampling was performed over a 4‐year period (2008–2011) for five locations with drainage areas of 0.62, 3.5, 4.5, 12 and 79 ha. Hydrologic flow paths were characterized using an end‐member mixing model. DOM composition was determined using a suite of spectrofluorometric indices and a site‐specific parallel factor analysis model. Dissolved organic carbon (DOC), humic‐like DOM and fluorescence index were most sensitive to changes with drainage scale, whereas dissolved organic nitrogen, specific UV absorbance, Sr and protein‐like DOM were least sensitive. DOM concentrations and humic‐like DOM constituents were highest during both baseflow and stormflow for a 3.5‐ha catchment with a wetland near the catchment outlet. Whereas storm‐event concentrations of DOC and humic DOM constituents declined, the mass exports of DOC increased with increasing catchment scale. A pronounced dilution in storm‐event DOC concentration was observed at peak stream discharge for the 12‐ha drainage location, which was not as apparent at the 79‐ha scale, suggesting key differences in supply and transport of DOM. Our observations indicate that hydrologic flow paths, especially during storms, and the location and extent of wetlands in the catchment are key determinants of DOM concentration and composition. This study furthers our understanding of changes in DOM with drainage scale and the controls on DOM in headwater, forested catchments. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

12.
Quantifying the linkages between vegetation disturbance by fire and the changes in hydrologic processes leading to post‐fire erosional response remains a challenge. We measured the influence of fire severity, defined as vegetation disturbance (using a satellite‐derived vegetation disturbance index, VDI), landscape features that organize hydrologic flow pathways (relief and elongation ratios), and pre‐fire vegetation type on the probability of the occurrence of post‐fire gully rejuvenation (GR). We combined field surveys across 270 burned low‐order catchments (112 occurrences of GR) and geospatial analysis to generate a probabilistic model through logistic regression. VDI alone discriminated well between catchments where GR did and did not occur (area under the curve = 0.78, model accuracy = 0.72). The strong effect of vegetation disturbance on GR suggests that vegetation exerts a primary influence on the occurrence of infiltration excess run‐off and post‐fire erosion and that major gully erosion will not occur until fire consumes aboveground biomass. Other topographic and local factors also influenced GR response, including catchment elongation, per cent pre‐fire shrub, mid‐slope riparian vegetation, armoured headwaters, firehose effects, and concentration of severe burn in source areas. These factors highlight the need to consider vegetation effects in concert with local topography and site conditions to understand the propensity for flow accumulation leading to GR. We present a process‐based conceptual hydrologic model where vegetation loss from fire decreases rainfall attenuation and surface roughness, leading to accelerated flow accumulation and erosion; these effects are also influenced by interactions between fire severity and landscape structure. The VDI metric provides a continuous measure of vegetation disturbance and, when placed in a hydrologic context, may improve quantitative analysis of burned‐area susceptibility to erosive rainfall, hazard prediction, ecological effects of fire, landform evolution, and sensitivity to climate change. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

13.
The use of precipitation estimates from weather radar reflectivity has become widespread in hydrologic predictions. However, uncertainty remains in the use of the nonlinear reflectivity–rainfall (Z‐R) relation, in particular for mountainous regions where ground validation stations are often lacking, land surface data sets are inaccurate and the spatial variability in many features is high. In this study, we assess the propagation of rainfall errors introduced by different Z‐R relations on distributed hydrologic model performance for four mountain basins in the Colorado Front Range. To do so, we compare spatially integrated and distributed rainfall and runoff metrics at seasonal and event time scales during the warm season when convective storms dominate. Results reveal that the basin simulations are quite sensitive to the uncertainties introduced by the Z‐R relation in terms of streamflow, runoff mechanisms and the water balance components. The propagation of rainfall errors into basin responses follows power law relationships that link streamflow uncertainty to the precipitation errors and streamflow magnitude. Overall, different Z‐R relations preserve the spatial distribution of rainfall relative to a reference case, but not the precipitation magnitude, thus leading to large changes in streamflow amounts and runoff spatial patterns at seasonal and event scales. Furthermore, streamflow errors from the Z‐R relation follow a typical pattern that varies with catchment scale where higher uncertainties exist for intermediate‐sized basins. The relatively high error values introduced by two operational Z‐R relations (WSR‐57 and NEXRAD) in terms of the streamflow response indicate that site‐specific Z‐R relations are desirable in the complex terrain region, particularly in light of other uncertainties in the modelling process, such as model parameter values and initial conditions. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

14.
15.
We report the complex spatial and temporal dynamics of hyporheic exchange flows (HEFs) and nitrogen exchange in an upwelling reach of a 200 m groundwater-fed river. We show how research combining hydrological measurement, geophysics and isotopes, together with nutrient speciation techniques provides insight on nitrogen pathways and transformations that could not have been captured otherwise, including a zone of vertical preferential discharge of nitrate from deeper groundwater, and a zone of rapid denitrification linking the floodplain with the riverbed. Nitrate attenuation in the reach is dominated by denitrification but is spatially highly variable. This variability is driven by groundwater flow pathways and landscape setting, which influences hyporheic flow, residence time and nitrate removal. We observed the spatial connectivity of the river to the riparian zone is important because zones of horizontal preferential discharge supply organic matter from the floodplain and create anoxic riverbed conditions with overlapping zones of nitrification potential and denitrification activity that peaked 10–20 cm below the riverbed. Our data also show that temporal variability in water pathways in the reach is driven by changes in stage of the order of tens of centimetres and by strength of water flux, which may influence the depth of delivery of dissolved organic carbon. The temporal variability is sensitive to changes to river flows under UK climate projections that anticipate a 14%–15% increase in regional median winter rainfall and a 14%–19% reduction in summer rainfall. Superimposed on seasonal projections is more intensive storm activity that will likely lead to a more dynamic and inherently complex (hydrologically and biogeochemically) hyporheic zone. We recorded direct evidence of suppression of upwelling groundwater (flow reversal) during rainfall events. Such flow reversal may fuel riverbed sediments whereby delivery of organic carbon to depth, and higher denitrification rates in HEFs might act in concert to make nitrate removal in the riverbed more efficient.  相似文献   

16.
17.
The suitability of the physically based model SHETRAN for simulating sediment generation and delivery with a high degree of spatial (20 m) and temporal (sub‐hourly) resolution was assessed through application of the model to a 167‐km2 catchment leading to an estuary in New Zealand. By subdividing the catchment and conducting calculations on a computer cluster for a 6‐month hydrology initialisation period, it was possible to simulate a large rainfall event and its antecedent conditions in 24 h of computation time. The model was calibrated satisfactorily to catchment outlet flow and sediment flux for a large rainfall event in two subcatchments (~2 km2). Validation for a separate subcatchment was successful for flow (Nash–Sutcliff efficiency of 0.84) with a factor 2.1 over‐prediction for sediment load. Validation for sediment at full catchment scale using parameters from the subcatchment scale was good for flow but poor for sediment, with gross under‐estimation of the dominant stream sources of sediment. After recalibration at catchment scale, validation for a separate event gave good results for flow (Nash–Sutcliff efficiency of 0.93) and sediment load within a factor of two of measurements. An exploratory spatially explicit landslide model was added to SHETRAN, but it was not possible to test this fully because no landslides were observed in the study period. Application to climate change highlighted the non‐linear response to extreme rainfall. However, full exploration of land use and climate change and the evaluation of uncertainty were severely constrained by computational limitations. Subdivision of the catchment with separate stream routing is suggested as a way forward to overcome these limitations. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

18.
Lei Yao  Liding Chen  Wei Wei 《水文研究》2016,30(12):1836-1848
Imperviousness, considered as a critical indicator of the hydrologic impacts of urbanization, has gained increasing attention both in the research field and in practice. However, the effectiveness of imperviousness on rainfall–runoff dynamics has not been fully determined in a fine spatiotemporal scale. In this study, 69 drainage subareas <1 ha of a typical residential catchment in Beijing were selected to evaluate the hydrologic impacts of imperviousness, under a typical storm event with a 3‐year return period. Two metrics, total impervious area (TIA) and effective impervious area (EIA), were identified to represent the impervious characteristics of the selected subareas. Three runoff variables, total runoff depth (TR), peak runoff depth (PR), and lag time (LT), were simulated by using a validated hydrologic model. Regression analyses were developed to explore the quantitative associations between imperviousness and runoff variables. Then, three scenarios were established to test the applicability of the results in considering the different infiltration conditions. Our results showed that runoff variables are significantly related to imperviousness. However, the hydrologic performances of TIA and EIA were scale dependent. Specifically, with finer spatial scale and the condition heavy rainfall, TIA rather than EIA was found to contribute more to TR and PR. EIA tended to have a greater impact on LT and showed a negative relationship. Moreover, the relative significance of TIA and EIA was maintained under the different infiltration conditions. These findings may provide potential implications for landscape and drainage design in urban areas, which help to mitigate the runoff risk. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

19.
Heejun Chang 《水文研究》2007,21(2):211-222
This study investigates changes in streamflow characteristics for urbanizing watersheds in the Portland Metropolitan Area of Oregon for the period from 1951 to 2000. The objective of this study was to assess how mean annual runoff ratio, mean seasonal runoff ratio, annual peak runoff ratio, changes in streamflow in response to storm amount, the fraction of time that the daily mean flow exceeds the annual mean flow, 3‐day recession constants, and dry/wet flow ratio vary among watersheds with different degrees of urban development. There were no statistically significant changes in annual runoff ratio and annual peak runoff ratio for the mixed land‐use watershed (Tualatin River watershed) and the urban watershed (Johnson Creek watershed) during the entire study period. The Tualatin River watershed, where most of the urban development occurred in a lower part of the watershed, showed a statistically significant increase in annual peak runoff ratio during the 1976 and 2000 period. The Upper Tualatin River watershed illustrated a significant decrease in annual peak runoff ratio for the entire study period. With significant differences in seasonal runoff ratio, only Johnson Creek exhibited a significant increase in both wet and dry season runoff ratios. Streamflow during storm events declined rapidly in the urban watershed, with a high 3‐day recession constant. At an event storm scale, streamflow in Fanno Creek, which is the most urbanized watershed, responded quickly to precipitation input. The fraction of time that the daily mean flow exceeded the annual mean flow and dry/wet flow ratio are all lower in Johnson Creek. This suggests a shorter duration of storm runoff and lower baseflow in the urbanized watershed when compared to the mixed land use watershed. The findings of this study demonstrate the importance of spatial and temporal scale, climate variability, and basin physiographic characteristics in detecting the hydrologic effects of urbanization in the Pacific Northwest of the USA. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

20.
水资源恶化、水体富营养化严重威胁生态环境健康,农业活动所产生的氮是造成水污染的主要原因之一.本研究以句容水库农业流域为研究对象,基于实地监测数据验证了SWAT模型模拟当地氮污染的适用性,并分析了氮素负荷的时空分布特征及其关键源.结果 表明:硝态氮(NO3-N)和总氮(TN)的年均入库量分别为9.98和27.22 t.时...  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号