首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 39 毫秒
1.
Dynamic Hydrologic Simulation of the Bear Brook Watershed in Maine (BBWM)   总被引:2,自引:0,他引:2  
Bear Brook Watershed in Maine (BBWM) consists of a pair of research watersheds, East Bear Brook (EBB) and West Bear Brook (WBB). Years of research and observations have shown both watersheds have high similarity in geographic and hydrologic characteristics; a simple comparison of hydrographs from these two watersheds further substantiates this similarity. The Object Watershed Link Simulation (OWLS) model was developed and used to simulate the hydrological processes within the BBWM. The OWLS model is a 3-dimensional, vector-based, visualized, physically-based, distributed watershed hydrologic model. Simulation results not only provide a close examination of hydrologic processes within a watershed, but also dynamically visualize the processes of flow separations and Variable Source Areas (VSA). Results from flow separations suggest that surface flow from riparian area is the predominate component for the flood rising limb and that macropore flow from riparian area dominates during the falling limb. Soil matrix flow has little effect flood period but is a persistent contributor to base flow. Results from VSA visualization demonstrate 3-D dynamic changes in surface flow distribution and suggest that downstream riparian areas are the major contributing area for peak flow. As water chemistry is highly relevant to the flow paths within a watershed, simulations have provided valuable information about source of stream flow and the water migration dynamics to support the study of watershed chemistry in the BBWM. More specific linkages between the chemistry behavior and the dynamic hydrologic processes should become the next simulation effort in the watershed study. There are many questions that are critical to watershed chemistry studies like: which flow component (surface flow, macropore flow, soil matrix flow) predominates during peak flows? How do the flow components distribute during a flood event? How do flow contributions differ between these two watersheds? Which portion of the watershed contributes the most to the peak flows? These questions remain unknown from previous observations and only can be addressed with a physically-based distributed model.  相似文献   

2.
Chronic elevated nitrogen (N) deposition has altered the N status of temperate forests, with significant implications for ecosystem function. The Bear Brook Watershed in Maine (BBWM) is a whole paired watershed manipulation experiment established to study the effects of N and sulfur (S) deposition on ecosystem function. N was added bimonthly as (NH4)2SO4 to one watershed from 1989 to 2016, and research at the site has studied the evolution of ecosystem response to the treatment through time. Here, we synthesize results from 27 years of research at the site and describe the temporal trend of N availability and N mineralization at BBWM in response to chronic N deposition. Our findings suggest that there was a delayed response in soil N dynamics, since labile soil N concentrations did not show increases in the treated watershed (West Bear, WB) compared to the reference watershed (East Bear, EB) until after the first 4 years of treatment. Labile N became increasingly available in WB through time, and after 25 years of manipulations, treated soils had 10× more extractable ammonium than EB soils. The WB soils had 200× more extractable nitrate than EB soils, driven by both, high nitrate concentrations in WB and low nitrate concentrations in EB. Nitrification rates increased in WB soils and accounted for ~?50% of net N mineralization, compared to ~?5% in EB soils. The study provides evidence of the decadal evolution in soil function at BBWM and illustrates the importance of long-term data to capture ecosystem response to chronic disturbance.  相似文献   

3.
The paired watershed experiment at the Bear Brook Watershed in Maine (BBWM) provided an opportunity to study changes in forest soil O horizon properties as a result of experimental, chronic N additions. The West Bear brook watershed received elevated N and S inputs since November 1989 as bimonthly applications of (NH4)2SO4. Forest floor samples (O horizon) were collected in July of 1992 from three dominant stand and five soil types at BBWM. The (NH4)2SO4 amendments in the treated watershed (West Bear) stimulated potential net nitrification, but significant increases were found only in hardwood O horizons after three years of treatment. Hardwood stand forest floor soil materials had the lowest C:N ratios (mean=23), compared with mixedwood (mean=27) and softwood stands (mean=33). NH4-N accounted for over 95% of the inorganic N in the forest floor. The lack of a strong relationship between soil type and potential net N mineralization at BBWM, coupled with conflicting results in the literature, suggested that stand characteristics were more important than conventional soil nomenclature based on pedogenetic features, or 2.5 years of treatments, in defining differences in soil N dynamics and responses to increased N inputs.  相似文献   

4.
The Bear Brook Watershed in Maine (BBWM), USA, and the Fernow Experimental Forest in West Virginia, USA, represent unique, long-term, paired, whole watershed, experimental manipulations focusing on the effects of nitrogen (N) and sulfur (S) deposition on temperate forests. Both watersheds began whole-ecosystem additions of N and S as (NH4)2SO4 in the fall of 1989, and both are entering their third decade of chronic enrichment of the treated watersheds, while the reference watersheds offer unique opportunities to evaluate forest watershed responses to recovery. Differences between BBWM and Fernow in the history of atmospheric deposition, soil properties, and forest composition all contribute to different response trajectories in stream chemical exports over time. The four watersheds represent a spectrum of N enrichment and retention, ranging from ≈98% N retention in the reference watershed in Maine, to ≈20% N retention in the treated watershed in West Virginia. Despite these differences, there is evidence that mechanisms of response in base cation leaching and other processes are similar among all four watersheds. In both cases, the history to date of two decades of research and monitoring has provided new insights into ecosystem response not evident in more traditional short-term research.  相似文献   

5.
Buried mineral soil-bags and natural solutions were studied as indicators of forest ecosystem response to elevated N and S inputs at the Bear Brook Watershed in Maine (BBWM). The BBWM is the site of a paired watershed manipulation experiment in a northern New England forested ecosystem. The study includes two small (10 ha each) catchments dominated by northern hardwood forests with red spruce in the upper elevations. Treatments consist of (NH4)2SO4 applied to the West Bear watershed six times per year, increasing N and S deposition 3× and 2× above ambient values, respectively. Buried mineral soil-bag changes over time reflected both the native soil environment and the treatments. Most of the treatment effects on mineral soils were evident as higher inorganic S found in the treated watershed soils. Adsorbed SO4 in the buried mineral soil-bags increased by approximately 40% under softwood stands and 50% under hardwood stands over the study period. Hardwood soil solutions responded with significant increases in NO3 and SO4 concentrations that resulted in accelerated cation leaching, primarily Ca and Al. Few differences that could be attributed to treatments were evident in soil solutions under softwoods. No treatment effects were evident in throughfall and stemflow chemistry.  相似文献   

6.
The Bear Brook Watershed in Maine (BBWM) is a long-term, paired watershed experiment that addresses the effects of acid and nitrogen (N) deposition on whole watersheds. To examine stream response at BBWM, we synthesized data on organic matter dynamics, including leaf breakdown rates, organic matter inputs and standing stocks, macroinvertebrate secondary production, and nutrient uptake in treated and reference streams at the BBWM. While N concentrations in stream water and leaves have increased, the input, standing stocks, and breakdown rates of leaves, as well as macroinvertebrate production, were not responsive to acid and N deposition. Both chronic and acute increases of N availability have saturated uptake of nitrate in the streams. Recent experimental increases in phosphorus (P) availability enhanced stream capacity to take up nitrate and altered the character of N saturation. These results show how the interactive effects of multiple factors, including environmental flow regime, acidification, and P availability, may constrain stream response to chronic N deposition.  相似文献   

7.
Dry (NH4)2SO4 (1,800 eq-1 ha-1 yr-1) has been applied to the western of two contiguous 10 ha catchments at the Bear Brook Watershed in Maine (BBWM) since November, 1989. The initial rapid and significant response in both S and N in West Bear, compared to the reference East Bear, slowed after three years. Annual S retention of the total experimental treatment decreased from 86 to 34%, with a seven year cumulative retention of 59%. Hydrology influences the export flux of S; S is retained more in dry seasons and dry years. The annual retention of N has decreased from 96 to 81%, with a cumulative retention of 82%. The export of N from the reference watershed has declined from 178 to 23 eq-1 ha-1 yr-1 during the treatment period. The treatment N (as NH4) initially stimulated nitrification, and caused pre-existing N to be lost in runoff, rather than the treatment N. Retention of the treatment N has decreased to approximately 80%. The majority of the retained N is stored in the soil, but the reasons for the decreased flux from the reference watershed are not known.  相似文献   

8.
At the Bear Brook Watershed in Maine (BBWM), the forest tree composition was characterized and the effects of the chronic ammonium sulfate ((NH4)2SO4) treatment on basal area growth, foliar chemistry, and gas exchange were investigated on forest species. The BBWM is a paired watershed forest ecosystem study with one watershed, West Bear (WB), treated since 1989 with 26.6 kg N ha???1 year???1 and 30 kg S ha???1 year???1applied bimonthly as (NH4)2SO4, while the other watershed, East Bear (EB), serves as a reference. Tree species richness, density, and mortality were found to be similar between watersheds. Basal area increment was estimated from red spruce and sugar maple, showing that, for the first 7 years of treatment, it was significantly higher for sugar maple growing in WB compared to EB, but no differences were observed for red spruce between watersheds. However, the initial higher sugar maple basal area growth in WB subsequently decreased after 8 years of treatment. Foliar chemical analysis performed in trees, saplings, and ground flora showed higher N concentrations in the treated WB compared to the reference EB. But, foliar cation concentrations, especially Ca and Mg, were significantly lower for most of the species growing in WB compared with those growing in EB. For sugar maple, foliar N was higher on WB, but there were no differences in foliar Ca and Mg concentrations between treated and reference watersheds. In addition, only sugar maple trees in the treated WB showed significantly higher photosynthetic rates compared to reference EB trees.  相似文献   

9.
Despite decades of research about episodic acidification in many regions of the world, the understanding of what controls the transient changes in stream water chemistry occurring during rain and snow melt events is still limited. Here, we use 20 years of hydrological and stream chemical data from the paired watershed study at Bear Brook Watershed in Maine (BBWM), USA to improve the understanding of the effects of acid deposition on the causes, drivers, and evolution of episodic acidification. The long-term experimental study at BBWM includes 18 years of chemical treatment of the West Bear Brook (WB) watershed with (NH4)2SO4. East Bear Brook (EB) serves as reference. The treatment started in 1989 following a 2-year pretreatment period. We analyzed 212 hydrological episodes using an episode model that can separate and quantify individual drivers of the transient change in acid-neutralizing capacity (ANC) during hydrological events. The results suggest that 18 years of N and S addition have not affected the natural drivers of episodic acidification of base-cation dilution, marine sea salt episodes, or organic acidity during rain and snow melt events. The contribution of SO $_{4}^{2-}$ to the ANC decline in WB has been increasing linearly since the beginning of watershed treatment, while the role of NO $_{3}^{-}$ has remained relatively constant after an initial increase. This is contradictory to many previous shorter-term studies and illustrates the need for a more mechanistic understanding of the causes and drivers of episodic acidification during rain- and snow melt-driven hydrological events.  相似文献   

10.
Atmospheric deposition of nitrogen (N) and sulfur (S) containing compounds affects soil chemistry in forested ecosystems through (1) acidification and the depletion of base cations, (2) metal mobilization, particularly aluminum (Al), and iron (Fe), (3) phosphorus (P) mobilization, and (4) N accumulation. The Bear Brook Watershed in Maine (BBWM) is a long-term paired whole-watershed experimental acidification study demonstrating evidence of each of these acidification characteristics in a northeastern U.S. forested ecosystem. In 2003, BBWM soils were studied using the Hedley fractionation procedure to better understand mechanisms of response in soil Al, Fe, and P chemistry. Soil P fractionation showed that recalcitrant P was the dominant fraction in these watersheds (49%), followed by Al and Fe associated P (24%), indicating that a majority of the soil P was biologically unavailable. Acidification induced mobilization of Al and Fe in these soils holds the potential for significant P mobilization. Forest type appears to exert important influences on metal and P dynamics. Soils supporting softwoods showed evidence of lower Al and Fe in the treated watershed, accompanied by lower soil P. Hardwood soils had higher P concentrations in surface soils as a result of increased biocycling in response to N additions in treatments. Accelerated P uptake and return in litterfall overshadowed acidification induced P mobilization and depletion mechanisms in hardwoods.  相似文献   

11.
The Bear Brook Watershed in Maine (BBWM), USA is a paired watershed study with chemical manipulation of one watershed (West Bear = WB) while the other watershed (East Bear = EB) serves as a reference. Characterization of hydrology and chemical fluxes occurred in 1987–1989 and demonstrated the similarity of the ca. 10 ha adjacent forested watersheds. From 1989–2010, we have added 1,800 eq (NH4)2SO4 ha???1 y???1 to WB. EB runoff has slowly acidified even as atmospheric deposition of SO $_{4}^{2-}$ has declined. EB acidification included decreasing pH, base cation concentrations, and alkalinity, and increasing inorganic Al concentration, as SO $_{4}^{2-}$ declined. Organic Al increased. WB has acidified more rapidly, including a 6-year period of increasing leaching of base cations, followed by a long-term decline of base cations, although still elevated over pretreatment values, as base saturation declined in the soils. Sulfate in WB has not increased to a new steady state because of increased anion adsorption accompanying soil acidification. Dissolved Al has increased dramatically in WB; increased export of particulate Al and P has accompanied the acidification in both watersheds, WB more than EB. Nitrogen retention in EB increased after 3 years of study, as did many watersheds in the northeastern USA. Nitrogen retention in WB still remains at over 80%, in spite of 20+ years of N addition. The 20-year chemical treatment with continuous measurements of critical variables in both watersheds has enabled the identification of decadal-scale processes, including ecosystem response to declining SO $_{4}^{-2}$ in ambient precipitation in EB and evolving mechanisms of treatment response in WB. The study has demonstrated soil mechanisms buffering pH, declines in soil base saturation, altered P biogeochemistry, unexpected mechanisms of storage of S, and continuous high retention of treatment N.  相似文献   

12.

Regime shifts of major salinity constituents (Ca, Mg, Na, K, SO4, Cl, HCO3, and NO3) in the lower Salinas River, an agricultural ecosystem, can have major impacts on ecosystem services central to continued agricultural production in the region. Regime shifts are large, persistent, and often abrupt changes in the structure and dynamics of social-ecological systems that occur when there is a reorganization of the dominant feedbacks in the system. Monitoring information on changes in the system state, controlling variables, and feedbacks is a crucial contributor to applying sustainability and ecosystem resilience at an operational level. To better understand the factors driving salinization of the lower Salinas River on the central coast of California, we examined a 27-year record of concentrations of major salinity constituents in the river. Although limited in providing an understanding of solute flux behavior during storm events, long-term “grab sampling” datasets with accompanying stream discharges can be used to estimate the actual history of concentrations and fluxes. We developed new concentration–discharge relationships to evaluate the dynamics of chemical weathering, hydrological processes, and agricultural practices in the watershed. Examinations of long-term records of surface water and groundwater salinity are required to provide both understanding and perspective towards managing salinity in arid and semi-arid regions while also enabling determination of the influence of external climatic variability and internal drivers in the system. We found that rock weathering is the main source of Ca, Mg, Na, HCO3, and SO4 in the river that further enables ion exchange between Ca, Mg, and Na. River concentrations of K, NO3, and Cl were associated with human activities while agricultural practices were the major source of K and NO3. A more direct anthropogenic positive trend in NO3 that has persisted since the mid-1990s is associated with the lag or memory effects of field cropping and use of flood irrigation. Event to inter-year scale patterns in the lower Salinas River salinity are further controlled by antecedent hydrologic conditions. This study underscores the importance of obtaining long-term monitoring records towards understanding watershed changes-of-state and time constants on the range of driving processes.

  相似文献   

13.
This paper presents the spatial and temporal characteristics of petroleum hydrocarbons (PH) in surficial sediments from the Songhuajiang River using the method of petroleum chemical fingerprinting. Twenty-four surface sediment samples were collected at 17 sites along the river from upstream to downstream in flood season (August 2005) and icebound season (December 2005) and were analyzed for PH including n-alkanes (C(16)-C(33)), isoprenoid alkanes (pristane and phytane), and unresolved complex mixture (UCM). The concentration of PH varied from 22.64 to 91.45 μg g(?-1) dry sediment. n-alkanes with a carbon number from 16 to 33 were detected in all samples, and the UCM was the dominant composition for PH. The variability of the PH concentration was mostly influenced by external conditions, such as seasonal change and industrial area position, as well as internal sediment physicochemical properties, such as organic carbon and grain size. The concentration of hydrocarbons is higher in flood season than in icebound season and is higher upstream than downstream. The diagnostic ratios of specific hydrocarbons showed that the PH pollution comes from a combination of biogenic and petrogenic sources, and petrogenic input is dominant in icebound season relative to flood season. It also indicates that there is a clear terrigenous input of n-alkanes in flood season. Principal components analysis was used to study the composition and characteristics of PH in Songhuajiang River sediments and to assess the spatial and temporal distribution of their natural and anthropogenic sources.  相似文献   

14.
This study was conducted to determine the response of stream water DOC and organic acidity to increased inputs of ammonium sulfate to a whole catchment. Precipitation, throughfall, soil solutions (from Spodosols) and stream waters were characterized for DOC concentrations and fractions (hydrophobic acids and neutrals, hydrophilic acids, bases, and neutrals) in both the control (East Bear) and the treatment (West Bear) catchments of Bear Brook Watershed, Maine (BBWM), a northern hardwood forest. In all solutions except precipitation, DOC was composed primarily of organic acids, with hydrophobic acids dominating (> 60% of DOC) in forest floor leachates (5000 mol C L-1), and a balance of hydrophobic and hydrophilic acids in deep B horizons and stream waters ( 150 mol C L-1). Stream waters had higher concentrations of DOC during storm or snowmelt events (high discharge), often reaching 300 to 400 mol C L-1. Forest floor leachate C was rapidly attenuated by the mineral soils under all flow conditions, indicating how important mineral soil sorption of DOC was in reducing the loss of C via surface water from BBWM. No differences occurred between control and treatment streams for concentration or composition of DOC due to treatment from 1989 through 1994. In 1995, West Bear Brook had much lower concentrations of DOC than East Bear for the first time. However, this occurred during a year of record low runoff, suggesting that hydrology may have affected export of C. Average annual export of DOC from the catchments was similar (1000 to 2000 mol C ha-1 yr-1). Organic anions in streamwaters increased slightly during high flow events (e.g., East Bear had means of 15 and 19 eq L-1 organic anions during base flow and high discharge in 1995). Treatment of West Bear caused a decrease in organic anions, both in concentration and contribution to overall anion composition (organic anions during high discharge as a percentage of total anions decreased from about 8 to 4% for 1987-89 and 1993-95 samples, respectively). This was probably due to decreased solution pH (greater protonation of organics) and higher concentrations of inorganic anions. Overall, there were no clear, detectable changes in stream water DOC, with only minor changes in organic anions, as a result of treatment with ammonium sulfate.  相似文献   

15.
A distributed hydrologic modeling and GIS approach is applied for the assessment of land use impact in the Steinsel sub-basin, Alzette, Grand-Duchy of Luxembourg. The assessment focuses on the runoff contributions from different land use classes and the potential impact of land use changes on runoff generation. The results show that the direct runoff from urban areas is dominant for a flood event compared with runoff from other land use areas in this catchment, and tends to increase for small floods and for the dry season floods, whereas the interflow from forested, pasture and agricultural field areas contributes to the recession flow. Significant variations in flood volume, peak discharge, time to the peak, etc., are found from the model simulation based on the three hypothetical land use change scenarios.  相似文献   

16.
The Bear Brook Watershed Manipulation in Maine is a paired watershed experiment. Monitoring of the paired catchments (East Bear Brook — reference; West Bear Brook — experimental) began in early 1987. Chemical manipulation of West Bear Brook catchment began in November 1989. Process studies on the watershed, outflow observation and monitoring, and modeling simulations of predicted versus observed response, are yielding a wealth of information about the behavior of the paired catchments and their responses to the deposition of acidifying substances. Results from the studies are providing important information relevant to national policies on emissions controls.  相似文献   

17.
Hydrologic response is an integrated indicator of watershed condition, and significant changes in land cover may affect the overall health and function of a watershed. This paper describes a procedure for evaluating the effects of land cover change and rainfall spatial variability on watershed response. Two hydrologic models were applied on a small semi-arid watershed; one model is event-based with a one-minute time step (KINEROS), and the second is a continuous model with a daily time step (SWAT). The inputs to the models were derived from Geographic Information System (GIS) theme layers of USGS digital elevation models, the State Soil Geographic Database (STATSGO) and the Landsat-based North American Landscape Characterization classification (NALC) in conjunction with available literature and look up tables. Rainfall data from a network of 10 raingauges and historical stream flow data were used to calibrate runoff depth using the continuous hydrologic model from 1966 to 1974. No calibration was carried out for the event-based model, in which six storms from the same period were used in the calculation of runoff depth and peak runoff. The assumption on which much of this study is based is that land cover change and rainfall spatial variability affect the rainfall-runoff relationships on the watershed. To validate this assumption, simulations were carried out wherein the entire watershed was transformed from the 1972 NALC land cover, which consisted of a mixture of desertscrub and grassland, to a single uniform land cover type such as riparian, forest, oak woodland, mesquite woodland, desertscrub, grassland, urban, agriculture, and barren. This study demonstrates the feasibility of using widely available data sets for parameterizing hydrologic simulation models. The simulation results show that both models were able to characterize the runoff response of the watershed due to changes of land cover.  相似文献   

18.
Ion-exchange resins (IER) offer alternative approaches to measuring ionic movement in soils that may have advantages over traditional approaches in some settings, but more information is needed to understand how IER compare with traditional methods of measurement in forested ecosystems. At the Bear Brook Watershed in Maine (BBWM), one of two paired, forested watersheds is treated bi-monthly with S and N (28.8 and 25.2kgha−1yr−1 of S and N, respectively). Both IER and ceramic cup tension lysimeters were used to study soil solution responses after ∼11 years of treatment. Results from both methods showed treatments resulted in the mobilization of base cations and Al, and higher SO4—S and inorganic N in the treated watershed. Both methods indicated similar differences in results associated with forest type (hardwoods versus softwoods), a result of differences in litter quality and atmospheric aerosol interception capacity. The correlation between lysimeter and IER data for individual analytes varied greatly. Significant correlations were evident for Na (r=0.75), Al (r=0.65), Mn (r=0.61), Fe (r=0.57), Ca (r=0.49), K (r=0.41) and NO3—N (r=0.59). No correlation was evident between IER and soil solution data for NH4—N and Pb. Both IER and soil solution techniques suggested similar interpretations of biogeochemical behavior in the watershed.  相似文献   

19.
Spring snowmelt is an important period of mercury (Hg) export from watersheds. Limited research has investigated the potential effects of climate variability on hydrologic and Hg fluxes during spring snowmelt. The purpose of this research was to assess the potential impacts of inter-annual climate variability on Hg mobility in forested uplands, as well as spatial variability in hillslope hydrology and Hg fluxes. We compared hydrological flows, Hg and solute mobility from three adjacent hillslopes in the S7 watershed of the Marcell Experimental Forest, Minnesota during two very different spring snowmelt periods: one following a winter (2009-2010) with severely diminished snow accumulation (snow water equivalent (SWE) = 48 mm) with an early melt, and a second (2010-2011) with significantly greater winter snow accumulation (SWE = 98 mm) with average to late melt timing. Observed inter-annual differences in total Hg (THg) and dissolved organic carbon (DOC) yields were predominantly flow-driven, as the proportion by which solute yields increased was the same as the increase in runoff. Accounting for inter-annual differences in flow, there was no significant difference in THg and DOC export between the two snowmelt periods. The spring 2010 snowmelt highlighted the important contribution of melting soil frost in the timing of a considerable portion of THg exported from the hillslope, accounting for nearly 30% of the THg mobilized. Differences in slope morphology and soil depths to the confining till layer were important in controlling the large observed spatial variability in hydrological flowpaths, transmissivity feedback responses, and Hg flux trends across the adjacent hillslopes.  相似文献   

20.
The concentrations of polychlorinated biphenyls (PCBs) in sediments from the Fenhe reservoir and watershed were detected at 28 sites in wet and dry seasons. The ∑(123)PCBs ranged from n.d. to 126.49 ng g(-1) dw. The dominated congeners were tri-PCBs (34.29%) and tetra-PCBs (24.05%). In the Fenhe reservoir, ∑(123)PCBs presented a decreasing trend, while percentages of low chlorinated congeners showed an increasing trend. For the temporal variations, PCBs homologues profiles of sediment samples and spatial distribution of ∑(123)PCBs for the two periods were similar (with CD = 0.021 and r(2) = 0.999 respectively), although PCBs concentrations in the wet season were significantly higher than in the dry season. PCA was applied to analyze the possible sources for PCBs, suggesting that PCBs might be mainly influenced by Aroclor 1016 and 1242. Compared with 3 established sediment quality guidelines, levels of PCBs in sediments of the investigated watershed might have a potential biological impact, especially in the wet season.  相似文献   

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

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