小清河流域非点源污染负荷与降水关系研究

流域非点源污染造成了严重的水环境问题,为准确核算典型流域非点源污染负荷及为流域水环境治理提供支持和依据,以山东省小清河流域为研究区,采用SWAT模型及相关分析方法,研究了流域总氮(TN)、总磷(TP)营养物输出负荷时空变化规律以及营养物输出负荷与降水间的关系。结果表明:SWAT模型对小清河流域非点源污染模拟具有较好的适用性。汛期(7—9月)营养物输出负荷最高,占全年比重的50%以上,而TN是其中主要非点源污染物。TN、TP输出负荷空间分布相似,负荷较高的区域都主要集中在流域南部;该区域TN、TP输出负荷分别在69.72～235.30 kg/hm~2和0.93～4.73 kg/hm~2范围内。不同的土地利用氮输出负荷较高的依次为林地、耕地和草地,磷输出负荷较高的依次为林地、草地和耕地。流域氮输出强度与降水相关性强的区域主要集中在流域中上游,而磷输出强度与降水相关性强的区域则主要集中在流域中游。     

Using a GIS transfer model to evaluate pollutant loads in the Lake Kinneret watershed, Israel.

Lake Kinneret (Sea of Galilee) is the only large surface water body in Israel, encompassing an area of 167 km2 and supplying some 30% of the country's fresh water. Pollution from anthropogenic sources and water abstraction for domestic and agricultural uses has long been threatening the water quality of the lake. Point-source pollution in the watershed has decreased drastically with the development of wastewater treatment. However, diffuse pollution from agricultural activities is still an unresolved issue. In this paper we present an application of AVGWLF (a GIS-based watershed load model) to the Lake Kinneret watershed. The model allows one to simulate daily stream flows and monthly sediment, nitrogen, and phosphorus loads discharged to the lake from the surrounding watershed. Results from simulations yield a satisfactory correspondence between simulated and measured daily water volume. Partition by source of total phosphorus delivered to the lake in the period of 2000-04 confirms the reduction in point source nutrient contribution due to improvement of wastewater treatment facilities in the area. Future management should focus on reduction of nutrients originating from septic systems (point sources) and pasture and cropland areas (diffuse sources). Results from simulations will enable watershed managers to prioritize effective management alternatives for protecting the water quality in the lake.     

Spatially Explicit Modeling of Land Use Specific Phosphorus Transport Pathways to Improve TMDL Load Estimates and Implementation Planning

Diffuse pollution from urban stormwater and agricultural runoff are among the leading causes of water pollution in the USA. A process-oriented, stakeholder-driven research approach was implemented in the small heterogeneous watershed of St. Albans Bay, Vermont to model the relative load of phosphorus from all sources, including diffuse transport pathways, and compared to goals and assumptions outlined by a Total Maximum Daily Load (TMDL) developed for phosphorus in Lake Champlain. Mass-balance and dynamic landscape simulation models were used to describe the distribution of the average annual phosphorus load to streams (10.57 t/year) in terms of space, time, and transport process. The majority of the phosphorus load comes from two subwatersheds dominated by clay soils, Stevens and Jewett Brooks. Dissolved phosphorus in surface runoff from the agricultural landscape, driven by high soil phosphorus concentrations, accounts for 41% of the total load to watershed streams. Direct discharge from farmsteads and stormwater loads, primarily from road sand wash-off, are also significant sources. Results reported in this study could help target watershed interventions both in terms of the types and locations of recommended best management practices (BMPs). The study offers an approach to attaining TMDL diffuse pollution targets in a cost-effective and participatory manner and could be replicated for other TMDL processes around the country.     

Planning for Implementation of Riparian Buffers in the Feitsui Reservoir Watershed

The Feitsui reservoir is a major water supply source for more than five million people in Northern Taiwan. The reservoir water quality has been good, but is threatened by eutrophication due to excessive nutrient input and siltation due to sediment loads. Recently, the water authorities in Taiwan have made considerable efforts to devise strategies using watershed conservation practices for the protection of Feitsui reservoir water quality. The control of non-point source pollution (NPS) represents one of the major strategies and the use of best management practices (BMPs) is under careful consideration. The objective of this study was to assess the effectiveness and cost of a full Feitsui watershed implementation of riparian buffer strips and other appropriate conservation practices. Based on the use of watershed simulation models and a statistical relationship between pollution reduction rate and the width and slope of a buffer strip, a methodology for the planning and design of riparian buffer strips was addressed. Data from field experiments were used to calibrate the coefficients of the regression equations. Several planning scenarios were evaluated by means of cost-benefit analysis coupled with net present value method. Data on local construction and maintenance costs for the selected design and location of buffer strips were used in the analysis. Based on several cost-benefit analyses, the scenario for installing buffer strips with 30 m width and 5% slope along both sides of tributary streams in the sub-watersheds with high nutrient (phosphorus) loadings was found to be most cost effective.     

An Integrated Approach for Targeting Critical Source Areas to Control Nonpoint Source Pollution in Watersheds

This study presents an integrated approach for targeting critical source areas (CSAs) to control nonpoint source pollution in watersheds. CSAs are the intersections between hydrologically sensitive areas (HSAs) and high pollution producing areas of watersheds. HSAs are the areas with high hydrological sensitivity and potential for generating runoff. They were based on a soil topographic index in consistence of a saturation excess runoff process. High pollution producing areas are the areas that have a high potential for generating pollutants. Such areas were based on simulated pollution loads to streams by the Soil and Water Assessment Tool. The integrated approach is applied to the Neshanic River watershed, a suburban watershed with mixed land uses in New Jersey in the U.S. Results show that several land uses result in water pollution: agricultural land causes sediment, nitrogen and phosphorus pollution; wetlands cause sediment and phosphorus pollution; and urban lands cause nitrogen and phosphorus pollution. The primary CSAs are agricultural lands for all three pollutants, urban lands for nitrogen and phosphorus, and wetlands for sediment and phosphorus. Some pollution producing areas were not classified into CSAs because they are not located in HSAs and the pollutants generated in those areas are less likely to be transported by runoff into streams. The integrated approach identifies CSAs at a very fine scale, which is useful for targeting the implementation of best management practices for water quality improvement, and can be applied broadly in different watersheds to improve the economic efficiency of controlling nonpoint source pollution.     

Nonlinear regression approach to evaluate nutrient delivery coefficient.

Implementation of the Korean Total Maximum Daily Load Act calls for new tools to quantify nutrient losses from diffuse sources at a river basin district scale. In this study, it was elucidated that the nonlinear regression model (NRM) reduces the uncertainty of the boundary conditions of the water quality model. The NRM was proposed to analyse the delivery coefficients of surface waters and retention coefficients of pollutants. Delivery coefficient of pollution load was considered as a function of two variables: the watershed form ratio, S(f), which is a measurable geomorphologic variable and the retention coefficient, phi, which is an empirical constant representing the basin-wide retarding capacity of pollutant wash-off. This model was applied on the Geum River, one of the major basins in South Korea. The QUAL2E was used to simulate stream water quality using NRM. In this paper, we elucidate the possibility to use a nonlinear regression model for delivery and retention of nutrients in a drainage basin characterized as both data-rich and data-poor, and the magnitude of the nutrient loads and sources has been uncertain for a long time.     

An extrapolation method for estimating loads from unmonitored areas using watershed model load ratios

It is important to routinely estimate loads from an entire watershed to describe current conditions and evaluate how watershed-wide management efforts have affected the nutrient and sediment export that affect downstream water quality. However, monitoring in most areas, including the Great Lakes watershed, consists of sampling at a limited number of sites that are only periodically used to estimate total watershed loading. Here, we describe a technique to extrapolate loads measured at a limited number of reference sites to the total load from a large watershed using load ratios between monitored sites and unmonitored areas obtained from a watershed model (i.e., model load ratio, MLR, approach). In this study, modeled nonpoint-source load ratios between monitored tributaries (reference sites) and nearby unmonitored areas and point-source delivery factors for all areas were obtained from a Spatially Referenced Regression On Watershed attributes (SPARROW) model and used to extrapolate the measured loads from an ongoing monitoring program (Great Lakes Restoration Initiative Tributary monitoring program) to the entire Great Lakes watershed. The MLR approach incorporates spatial variability in nonpoint- and point-source delivery, watershed characteristics, and hydrology that are often not considered when estimating loads from unmonitored areas, such as using the unit area load (UAL) extrapolation approach. The MLR approach provided smaller watershed loads than the UAL approach because yields from monitored sites, in general, were larger than from unmonitored areas. When both approaches were used to estimate loads at adjacent monitored sites, the MLR approach provided more accurate estimates than the UAL approach.     

Applicability of modelling tools in watershed management for controlling diffuse pollution.

Diffuse pollution is hard to analyze, control and manage by its nature. Watershed models and Geographical Information Systems (GIS) are recently developed tools that aid analysis of diffuse sources of pollution. However, their applications are not always easy and straightforward. Turkey is a typical example of a mountainous country rich in rivers and streams. Due to the complex geomorphology, land-use and agricultural practices in most of the watersheds in Turkey, modelling, analyzing and managing diffuse pollution has been a challenge. The complex watershed structure forces the modellers to work with spatially high resolution data. Apart from the data, the models themselves may also cause operational problems. These issues and their probable solutions form the basis of the discussions in this paper. It acts as a guideline for modelling and analyzing diffuse pollution by emphasizing the referred problems and difficulties. Design of an Information Technology-based system tool for watershed and/or water quality modelling, which would be suitable for countries having watersheds with similar structure and problems to those of Turkey, is also outlined.     

Evaluating diffuse and point source phosphorus inputs to streams in a cold climate region using a load apportionment model

Variation in the timing and quantity of diffuse versus point-source inputs of phosphorus (P) to streams can be evaluated by examining P concentration-flow relationships. Diffuse load inputs usually increase with stream flow (due to increased delivery caused by precipitation); whereas, point-source concentrations decrease with rising river flow (due to increased dilution). This study tested the suitability of a load apportionment model (LAM), a power-law function of flow, to estimate contributions of diffuse and point inputs to P loads for eight sub-watersheds in the Red River Valley, a cold-climate rural region of Manitoba, Canada. For all but two sub-watersheds, annual and seasonal (snowmelt and summer) models of P concentration versus flow best fit a strictly diffuse source contribution. The models identified significant point-source inputs (in addition to diffuse sources) in two sub-watersheds, during summer in both watersheds (consistent with the fact that wastewater from sewage lagoons is discharged to upstream reaches between June and September) and during snowmelt for one watershed. Application of a LAM proved to be a simple and rapid method for nutrient source apportionment as well as detection of unknown sources for cold-climate, rural sub-watersheds. Such information is critical for developing the most effective mitigation strategies to reduce P concentrations and eutrophication risk.     

Spatial Optimization of Best Management Practices to Attain Water Quality Targets

Diffuse nutrient loads are a common problem in developed and agricultural watersheds. While there has been substantial investment in best management practices (BMPs) to reduce diffuse pollution, there remains a need to better prioritize controls at the watershed scale as reflected in recent US-EPA guidance for watershed planning and Total Maximum Daily Load development. We implemented spatial optimization techniques among four diffuse source pathways in a mixed-use watershed in Northern Vermont to maximize total reduction of phosphorus loading to streams while minimizing associated costs. We found that within a capital cost range of 138 to 321 USD ha^-1 a phosphorus reduction of 0.29 to 0.38 kg ha^?1 year^?1, is attainable. Optimization results are substantially more cost-effective than most scenarios identified by stakeholders. The maximum diffuse phosphorus load reduction equates to 1.25 t year^?1using the most cost-effective technologies for each diffuse source at a cost of $3,464,260. However, 1.13 t year^?1 could be reduced at a much lower cost of $976,417. This is the practical upper limit of achievable diffuse phosphorus reduction, above which additional spending would not result in substantially more phosphorus reduction. Watershed managers could use solutions along the resulting Pareto optimal curve to select optimal combinations of BMPs based on a water quality target or available funds. The results demonstrate the power of using spatial optimization methods to arrive at a cost-effective selection of BMPs and their distribution across a landscape.     

Simulation of Agricultural Management Alternatives for Watershed Protection

The Bosque River Watershed in Texas is facing a suite of water quality issues including excess sediment, nutrient, and bacteria. The sources of the pollutants are improperly managed cropland and grazing land, dairy manure application, and effluent discharge from wastewater treatment facilities. Several best management practices (BMPs) have been proposed for pollution reduction and watershed protection. The overall objectives of this study were to demonstrate a modeling approach using Soil and Water Assessment Tool (SWAT) model to simulate various BMPs and assess their long-term impacts on sediment and nutrient loads at different spatial levels. The SWAT model was calibrated and validated for long-term annual and monthly flows at Valley Mills and for monthly sediment, total nitrogen (TN) and total phosphorus (TP) at Hico and Valley Mills monitoring locations. The BMPs including streambank stabilization, gully plugs, recharge structures, conservation tillage, terraces, contour farming, manure incorporation, filter strips, and PL-566 reservoirs were simulated in the watershed areas that met the respective practice’s specific criteria for implementation. These BMPs were represented in the pre- and post-conditions by modifying one or more channel parameters (channel cover, erodibility, Manning’s n), curve number (CN), support practice factor (P-factor), filter strip width, and tillage parameters (mixing efficiency, mixing depth). The BMPs were simulated individually and the resulting Hydrologic Response Units (HRUs), subwatershed, and watershed level impacts were quantified for each BMP. Sensitivity of model output values to input parameters used to represent the BMPs was also evaluated. Implementing individual BMPs reduced sediment loads from 3% to 37% and TN loads from 1% to 24% at the watershed outlet; however, the changes in TP loads ranged from 3% increase to 30% decrease. Higher reductions were simulated at the subwatershed and HRU levels. Among the parameters analyzed for sensitivity, P-factor and CN were most sensitive followed by Manning’s n. The TN and TP outputs were not sensitive to channel cover. This study showed that the SWAT modeling approach could be used to simulate and assess the effectiveness of agricultural best management practices.     

Application of the Soil and Water Assessment Tool for six watersheds of Lake Erie: Model parameterization and calibration

The Soil and Water Assessment Tool (SWAT), a physically-based watershed-scale model, holds promise as a means to predict tributary sediment and nutrient loads to the Laurentian Great Lakes. In the present study, model performance is compared across six watersheds draining into Lake Erie to determine the applicability of SWAT to watersheds of differing characteristics. After initial model parameterization, the Huron, Raisin, Maumee, Sandusky, Cuyahoga, and Grand SWAT models were calibrated (1998-2001) and confirmed, or validated (2002-2005), individually for stream water discharge, sediment loads, and nutrient loads (total P, soluble reactive P, total N, and nitrate) based on available datasets. SWAT effectively predicted hydrology and sediments across a range of watershed characteristics. SWAT estimation of nutrient loads was weaker although still satisfactory at least two-thirds of the time across all nutrient parameters and watersheds. SWAT model performance was most satisfactory in agricultural and forested watersheds, and was less so in urbanized settings. Model performance was influenced by the availability of observational data with high sampling frequency and long duration for calibration and confirmation evaluation. In some instances, it appeared that parameter adjustments that improved calibration of hydrology negatively affected subsequent sediment and nutrient calibration, suggesting trade-offs in calibrating for hydrologic vs. water quality model performance. Despite these considerations, SWAT accurately predicted average stream discharge, sediment loads, and nutrient loads for the Raisin, Maumee, Sandusky, and Grand watersheds such that future use of these SWAT models for various scenario testing is reasonable and warranted.     

未来气候变化影响下的流域面源污染负荷特征响应评估

针对流域面源污染负荷在未来气候变化影响下的变化特征,以我国新安江上游率水流域为例,使用通用流域污染负荷模型(GWLF),对其2000-2013年的水量及总氮、总磷面源污染负荷通量进行了模拟,并解析了其负荷来源分配。在此基础上,基于政府间气候变化专业委员会(IPCC)的气候变化评估报告结果,利用GWLF模型分析了到21世纪20年代、50年代、80年代在A1FI(最高排放)和B1(最低排放)情景下,率水流域的水文及总氮、总磷面源污染负荷特征变化。结果表明:未来气候变化对流域水文及面源污染负荷特征均有一定影响。年水资源量先减少后增加,地表径流量和蒸发量逐渐上升而地下水量逐渐下降。到2080s,A1FI情景比B1情景有更多的水资源量。年总氮通量先增加后减少,在2050s最高,而年总磷通量则逐渐增加,且两种污染物均在A1FI排放情景下有更高的污染负荷量,表明人类温室气体的排放会潜在地增加流域水体面源污染负荷。     

Landscape and flow path-based nutrient loading metrics for evaluation of in-stream water quality in Saginaw Bay,Michigan

Landscape metrics are often used to model nonpoint source pollution from agricultural and urban surface runoff. By considering topography and the spatial arrangement of land cover, landscape metrics can better account for hydrologic connectivity, loading quantity, and vegetated buffer filtering between nutrient loading sources and streams. For this study we develop a surface runoff nutrient loading metric that considers source (i.e. cropland or developed) loading and buffer filtering along hydrologic transport vectors, or flow paths. We use General Additive Modeling to evaluate the relationship between this metric and in-stream nitrogen (N) and phosphorus (P) concentrations in the Saginaw Bay watershed in Michigan, US and compare the relative predictive power between this metric and other landscape metrics that do not consider hydrologic connectivity. The flow path-based cropland loading metric was a stronger predictor of in-stream NO₃ concentrations than alternative metrics. In-stream P concentrations were best predicted by models that included 48-h antecedent precipitation and catchment-wide proportion of developed landcover. Metric maps reveal high nutrient loading areas where only a small proportion of loading reaches streams via surface runoff, highlighting the need to consider nutrient loading via drainage tiles and other subsurface pathways in efforts to quantify nonpoint source loading to surface waters. The flow path-based loading metric is represented spatially as a gridded dataset showing estimates of nutrient loading adjacent to streams, and with higher resolution stream delineation data it could be used by land managers to target locations for precision buffer placement to intercept surface runoff and reduce nutrient loading.     

Estimation of non-point source pollution loads with flux method in Danjiangkou Reservoir area,China

The estimation of non-point source pollution loads into the Danjiangkou Reservoir is highly significant to environmental protection in the watershed. In order to overcome the drawbacks of traditional watershed numerical models, a base flow separation method was established coupled with a digital filtering method and a flux method. The digital filtering method has been used to separate the base flows of the Hanjiang, Tianhe, Duhe, Danjiang, Laoguan, and Qihe rivers. Based on daily discharge, base flow, and pollutant concentration data, the flux method was used to calculate the point source pollution load and non-point source pollution load. The results show that: (1) In the year 2013, the total inflow of the six rivers mentioned above accounted for 95.9% of the total inflow to the Danjiangkou Reservoir. The total pollution loads of chemical oxygen demand (COD_Mn) and total phosphorus (TP) from the six rivers were 58.20 × 10³ t and 1.863 × 10³ t, respectively, and the non-point source pollution loads were 39.82 × 10³ t and 1.544 × 10³ t, respectively, indicating that the non-point source pollution is a major factor (with a contribution rate of 68.4% for COD_Mn and 82.9% for TP). (2) The Hanjiang River is the most significant contributor of pollution loads to the Danjiangkou Reservoir, and its COD_Mn and TP contribution rates reached 79.3% and 83.2%, respectively. The Duhe River took the second place. (3) Non-point source pollution mainly occurred in the wet season in 2013, accounting for 80.8% and 90.9% of the total pollution loads of COD_Mn and TP, respectively. It is concluded that the emphasis of pollution control should be placed on non-point source pollution.     

Development of a new indicator of pollutant loads and its application to the Chesapeake Bay watershed

Pollutant load reductions are often required to restore aquatic ecosystems experiencing eutrophication. Loads can be estimated using watershed models or data from monitoring stations, however data availability can limit the timeliness or comprehensiveness of the load estimates. We developed an approach to address this challenge that used watershed model results to estimate the proportion of annual nonpoint source nitrogen (N), phosphorus (P) and sediment (Sed) loads derived from unmonitored catchments. This proportion was multiplied by the nonpoint portion of United States Geological Survey (USGS) estimated annual river loads to account for annual variation in hydrologic conditions. Total loads were calculated as the sum of measured river loads, reported point sources from unmonitored areas and the estimated nonpoint source loads from unmonitored catchments. We applied this approach to the Chesapeake Bay because of its socio‐economic and ecological importance. Median watershed loads for N, P and Sed were 140, 6.4 and 3030 Mg year^?1, respectively (1990–2004). Nonpoint source loads from the monitored areas constituted the greatest source of N, P and Sed (55, 47 and 74% respectively) to the Bay. The high N, P and Sed yield rates (7.3, 0.38 and 99 kg ha^?1 year^?1, respectively) from nonpoint loads originating from unmonitored areas near the Bay resulted in 25, 32 and 26% (N, P and Sed, respectively) of the Bay's total loads (excluding direct atmospheric deposition, shoreline erosion and oceanic inputs). Disproportionately high loads of P and Seds were associated with years that experienced elevated discharge whereas N loads were directly related to discharge. Error estimates indicated that our methods were most reliable for N (±6%) but reasonable for P (±22%) and provide an effective technique for the timely estimation of pollutant loads from watersheds with unmonitored catchments. Management strategies that decrease N deposition and reduce runoff to control P and Sed transport will effectively reduce pollutant loads. Published in 2010 by John Wiley & Sons, Ltd.     

Intermountain basins use in subtropical regions and their influences on benthic fauna

Conversion of high mountain forests into pastures or croplands and changes in land uses have produced changes in hydrology, vegetation cover and terrestrial aquatic linkages, resulting in increased nutrient loads on water bodies. The developing countries today need to obtain an indicator response signature under multiple stressor scenario, to identify causes of water bodies impairment. Intermountain land use impact on a Subtropical Andean River was studied in Tucumán province, Northwestern Argentina. Benthic macroinvertebrates were used as indicators of diffuse alterations produced by anthropic activities on a small valley at ca. 1000 m.a.s.l. The objectives were: (a) test the effect of land use on macroinvertebrate assemblages; (b) test the response of benthic fauna to diffuse organic pollution and (c) test the accuracy of the current biological water quality indices. Nutrient levels were higher than international standards. First axis in multivariate analysis shows an important upstream–downstream gradient meanwhile the second axis shows seasonal effects. These effects are specially produced by monsoonal regime of rainfall (ca. 80% in summer). No effect of the land use or diffuse pollution was detected on the macroinvertebrate assemblage structure, probably due to the good oxygen availability. Of the several biotic indices (BI) and metrics tested, only one worked on acceptable levels being ‘affected’ by small oxygen variations. Dry winter, dry and warm spring with high sunlight exposure plus high levels of nutrients could be an undesirable combination for standing waters of a reservoir projected in this watershed. A correct monitoring of water quality and alternatives as management of reforested margins is the right beginning. These actions and more studies represent interesting alternatives for poor regions of undeveloped countries. Copyright © 2008 John Wiley & Sons, Ltd.     

妫水河流域农业非点源污染负荷估算与分析

本文利用改进的输出系数模型对妫水河流域农业非点源污染负荷进行了估算,并且分别通过三分法,即分区、分类、分期的方法对结果进行了分析。结果表明:(1)2017年妫水河流域农业非点源污染物TN和TP的负荷量分别为1 402 214.9 kg/a,279 629.1 kg/a;(2)流域内各乡镇TN和TP的负荷量差异较大,旧县镇、永宁镇所贡献的TN、TP负荷量最多且单位负荷强度最大,属于重点治理区域;(3)流域内对TN的贡献最大的污染源为农业种植,对TP的贡献最大的污染源为农村生活污染,这两个污染源属于优先控制的污染源;(4)TN、TP负荷与降雨量呈现正相关关系,降雨冲刷是导致的污染物负荷增加的一个主要因素。农业非点源污染的三分研究结果为妫水河流域非点源污染治理与农业结构调整提供了科学依据。     

Pollution from animal husbandry in China: a case study of the Han River Basin

Animal husbandry is one of the major agricultural pollution sources in China. The Xiangyang Reach of the Han River Basin was used as a case study to identify pollutants from animal rearing. The gross amount of pollutants from livestock and poultry rearing in the Xiangyang Reach was estimated using two empirical models with different data sets. The pig, cattle, sheep, and poultry population in 2009 amounted to 2.6, 0.6, 0.5, and 39.2 million head, respectively. The total annual pollutant loads generated from the feces and urine of livestock and poultry were 270,400 t of chemical oxygen demand; 228,900 t of biochemical oxygen demand; 26,500 t of ammonia nitrogen; 16,500 t of total phosphorus; and 63,900 t of total nitrogen. Approximately 12% of these pollutant loads were estimated to enter the Han River through the watershed outlet. Animal breeding has been one of the main pollution sources in this area, followed by domestic sewage and industrial wastewater. Cattle produced the most pollution, with the heaviest pollution load in downtown Xiangyang City. Several recommendations are presented to control the pollution caused by livestock and poultry breeding.