Hydrochemical modelling of a stream dominated by organic acids and organically bound aluminium

An existing hydrogeochemical model, the Birkenes model (BIM), has been extended to include production of dissolved organic carbon (DOC) in the soil/water system and its chemical interactions with H⁺ and A1³⁺. The model has been calibrated and verified using precipitation and runoff data from the Svartberget catchment in northern Sweden. The catchment is impacted by moderate amounts of sulphur deposition (6 kg S ha^?1 a^?1; the stream is high in DOC (7–35 mg L^?1) and experiences episodically low pH-values (～4.0). The refined BIM is able to simulate the main variations in major chemical components of the stream water (H⁺, Ca²⁺, Mg²⁺, Al_i, Al₀, Cl^?, SO₄ ^?2 and DOC), as well as stream discharge and ¹⁸O, while Na⁺ is not well reproduced. Although very simplified relationships are used for the chemical interactions between DOC and H⁺ and aluminium, the model shows that these interactions are essential to stream acidity and aluminiumspecies content. Some of the model parameters are poorly determined by calibration with runoff data only. Soil and soil water investigations may contribute to the evaluation and development of the model structure and the representation of chemical processes. Further improvements of the model should emphasize DOC-production/absorption, detailed studies on DOC chemical behaviour and hydrological structure.     

Transient acid surges in an upland stream

This paper examines the hydrological and chemical processes controlling transient acid surges in a small, upland catchment in the Cairngorm Mountains of Scotland. The major episodes of increased acidity are produced by preferential elution of pollutants from snow at the onset of melt and by heavy autumn rainstorms when the catchment is saturated or frozen. Records of flow and water quality over a 3 yr period have been used to calibrate a lumped conceptual hydrochemical model. Using this model it is shown that concentration and distribution of input pollution by hydrological processes is the major control on episodic variations of stream water quality in this catchment.     

ORYZA2000模型与水稻群体茎蘖动态模型的耦合

以水稻群体茎蘖动态模拟为例,采用单向耦合方法,将ORYZA2000模型模拟的逐日发育进程和生物量作为水稻群体茎蘖动态模型的输入,驱动茎蘖动态模型模拟,尝试将水稻群体茎蘖动态模型引入ORYZA2000模型中,作为该模型的子模块以扩展ORYZA2000模型对茎蘖动态、籽粒灌浆动态和叶龄动态等水稻生长要素模拟的功能。为检验耦合模型（ORY-TIL）,选用2012年杂交籼稻两优培九和扬稻6号5个播期的大田观测数据,在两模型参数定标的基础上,模拟2个播期水稻群体茎蘖动态,并对模拟结果进行误差分析。结果显示,ORY-TIL模型能较好模拟水稻发育速率、地上部总生物量和群体茎蘖动态,各项模拟值与实测值的相关系数均达0.95以上,且通过0.01水平显著性检验;扬稻6号和两优培九群体茎蘖密度模拟值与实测值的均方根误差分别为24.3条?m-2和34.9条?m-2,误差主要出现在茎蘖消亡阶段。总体来看,ORY-TIL模型具有较好模拟性能,提出的耦合方法对扩展水稻生长模型具有一定参考价值。     

A critique of models for freshwater and soil acidification

Four types of models quantifying effects of acid deposition on freshwaters are reviewed. These include Henriksen's empirical model, an adsorption isotherm model, soil-oriented charge balance models (Reuss-Johnson, Birkenes, MAGIC, ILWAS), and the Trickle Down model. Emphasis is on an assessment of critical assumptions; no attempt has been made to run the various models and compare results. The models range from simple to very complex and from empirical to highly process oriented. The various types have all proven useful and there has been a significant convergence concerning key processes. The importance of anion mobility, sulfate adsorption, ion exchange, dissolution of Al bearing minerals and weathering seems to be accepted by most workers. Future model improvement, however, relies to a large extent on further checking against observations.     

不同土壤物理属性空间分布对水文过程模拟的影响

基于不同空间分辨率的土壤物理属性数据(土壤田间持水量、孔隙度和饱和水力传导度),运用WATLAC分布式水文模型模拟了西苕溪流域2005-2010年的水文过程。对比评价了土壤物理属性空间分布对地下水补给、土壤蒸发、径流量及其组分的影响。结果显示,利用不同空间分辨率土壤物理属性数据模拟的流量过程与实测值都匹配的很好,模拟精度相当,更详尽的土壤属性空间分布信息未能明显提高模型模拟的精度;而对于地下水补给量,高分辨率的土壤属性空间分布会使其模拟结果大幅减小,但对土壤蒸发量则影响较小,两种数据模拟的结果及其空间分布都很接近;同时,虽然土壤物理属性空间分布的分辨率对模拟的径流总量影响甚微,但对基流与地表径流的分割却产生了较大影响。     

基于MODIS遥感数据计算无定河流域日蒸散

为研究无定河流域日蒸散分布规律，应用遥感数据、农业气象站测量数据及Nishida模型等对该流域日蒸散进行了模拟。首先用2001～2002年晴天中国科学院禹城生态试验站Lysimeter测量日蒸散验证模型，模拟与测量的日蒸散相关系数达到0.61。随后，用该模型计算了无定河流域日蒸散，发现无定河流域日蒸散存在较为明显的空间分布规律：2001～2003年连续3年的8月份日蒸散都表现为东北部蒸散明显小于西南部，这是因为东北部基本是荒漠而东南部多是农田，且8月份日蒸散基本在2～5 mm之间变化；从2001年8月份第222 d日蒸散空间分布看，无定河主干道两边蒸散显著高于其他位置，这是由于8月份无定河流域为多雨季节，河谷土壤水分较高的缘故；从2002年内变化来看，不同的土地利用/覆被类型日平均蒸散差别不显著。     

Effects of cold and warm years on the water chemistry at the Birkenes catchment,Norway

During the last years, several scientists have been concerned with the consequences of global warming. In the Birkenes catchment, situated in southernmost Norway, several warm years have occurred during the last decade. In this climatic region, warm years primarily mean warmer winters with mean air temperatures close to or higher than 0 °C. Small temperature changes are therefore decisive for the quality of precipitation, i.e. whether the precipitation enters the catchments as snow or rain, which accordingly has large consequences for the residence time of water and its pathways through the catchment. Another important effect of higher winter-temperatures at Birkenes, is that stronger and more frequent seasalt episodes seem to occur. This may enhance the negative effects to aquatic organisms, because the most extreme concentrations of acute toxic Al-forms in streamwater are related to seasalt episodes during warm winters, when the precipitation comes as rain and thus leaves the catchment relatively fast through the uppermost soil layers. When seasalt-enriched water enters into the uppermost soil horizons, a large amount of sodium will be temporarily retained due to cation exchange processes, where sodium may substitute for H⁺ and cationic Al. At Birkenes, but also in similar acidified areas with high seasalt input, seasalt episodes are probably of greater importance for the periodic variations in streamwater H⁺ and Alⁿ⁺ than commonly recognized. Higher winter-temperatures may also change the timing of the highest streamwater fluxes of nitrate. This may have local consequences for the primary-production in the fjords, both quantitatively and qualitatively. Thus, global warming may have negative consequences for the surface water chemistry in near-coastal areas that also receive high amounts of acidifying atmospheric compounds, and may change the eutrophication status of many fjords.     

Calibrating and validating the LISEM model for two data sets from the Netherlands and South Africa

Within the framework of the GCTE Soil Erosion Network the Limburg Soil Erosion Model (LISEM) has been tested and validated in two catchments in South-Limburg (the Netherlands) and Zululand (South Africa). The calibration and validation of the Green–Ampt version of the LISEM model using 10 storms from the Catsop catchment in the Netherlands shows that differences in measured and simulated hydrographs and sediment loads can be large. These differences may be caused by the sensitivity of the model to some of the input variables, such as saturated hydraulic conductivity and the initial soil moisture content. Given the uncertainty in the input maps of these variables and a limited number of point data that is used to create these maps, a large part of the differences between measured and simulated data can be explained by these uncertainties. Thus, it is clear that detailed process-based models such as LISEM require very detailed and high-resolution input data in order to produce quantitative reliable results. The runoff processes in the Zululand catchment appeared to be dominated by slow throughflow and groundwater flow, which are processes that are not incorporated in LISEM.     

An atmosphere-ocean model for integrated assessment of global change

This paper describes the atmosphere-ocean system of the integrated model IMAGE 2.0. The system consists of four linked models, for atmospheric composition, atmospheric climate, ocean climate and for ocean biosphere and chemistry. The first model is globally averaged, the latter are zonally averaged with additional resolution in the vertical. The models reflect a compromise between describing the physical, chemical and biological processes and moderate computational requirements. The system is validated with direct observations for current conditions (climate, chemistry) and is consistent with results from General Circulation Model experiments. The system is used in the integrated setting of the IMAGE 2.0 model to give transient climate projections. Global surface temperature is simulated to increase by 2.5 K over the next century for socio-economic scenarios with continuing economic and population growth. In a scenario study with reduced ocean circulation, the climate system and the global C cycle are found to be appreciably sensitive to such changes.     

Dynamic Modelling of Spatially Variable Catchment Hydrochemistry for Critical Loads Assessment

Concern about acidification in upland areas has brought about the need to model the stream hydrochemical response to deposition and land-use changes and calculate critical loads. Application of dynamic models such as MAGIC are preferable to steady-state methods, since they are able to produce an estimate of the time scale required to meet some water chemistry target given a reduction in acid deposition. These models typically consider annual changes in stream chemistry at one point. However, in order to protect biota from 'acid episodes', quantification of temporal variability needs to encompass event responses; in addition spatial variability across the catchment also needs to be considered. In this paper, modelling of both spatial and temporal variability is combined in a new framework which enables quantification of catchment hydrochemical variability in time and space. Both low and high flow hydro-chemical variability are quantified in terms of statistical distributions of ANC (Acid Neutralisation Capacity). These are then input as stochastic variables to an EMMA (End-Member Mixing Analysis) model which accounts for temporal variability and ANC is hence predicted as a function of time and space across the whole catchment using Monte-Carlo simulation. The method is linked to MAGIC to predict future scenarios and may be used by iteration to calculate critical loads. The model is applied to the headwaters of the River Severn at Plynlimon, Wales, to demonstrate its capabilities.     

The Impact of Climate Change on Metal Transport in a Lowland Catchment

This study investigates the impact of future climate change on heavy metal (i.e., Cd and Zn) transport from soils to surface waters in a contaminated lowland catchment. The WALRUS hydrological model is employed in a semi-distributed manner to simulate current and future hydrological fluxes in the Dommel catchment in the Netherlands. The model is forced with climate change projections and the simulated fluxes are used as input to a metal transport model that simulates heavy metal concentrations and loads in quickflow and baseflow pathways. Metal transport is simulated under baseline climate (“2000–2010”) and future climate (“2090–2099”) conditions including scenarios for no climate change and climate change. The outcomes show an increase in Cd and Zn loads and the mean flux-weighted Cd and Zn concentrations in the discharged runoff, which is attributed to breakthrough of heavy metals from the soil system. Due to climate change, runoff enhances and leaching is accelerated, resulting in enhanced Cd and Zn loads. Mean flux-weighted concentrations in the discharged runoff increase during early summer and decrease during late summer and early autumn under the most extreme scenario of climate change. The results of this study provide improved understanding on the processes responsible for future changes in heavy metal contamination in lowland catchments.     

Investigating the Uncertainty in Predicting Responses to Atmospheric Deposition using the Model of Acidification of Groundwater in Catchments (MAGIC) within a Generalised Likelihood Uncertainty Estimation (GLUE) Framework

This study investigates the uncertainty associated with the modelled response of a catchment to historic and predicted future acidic deposition for the period 1851–2041. The MAGICmodel is applied within a GLUE framework to the 3.88 km² Afon Gwy catchment at Plynlimon, Wales. Nine million Monte Carlo simulations resulted in 5700 being accepted as behaviouralas defined by a fuzzy measure comparing observed to simulated variables. Model output and parameter sensitivity analysis indicate that, for this example where weathering rates are low,model dynamics are limited compared to control exerted by modelinitial conditions and by the specified acidic deposition boundary conditions. The results show that despite the small number of behavioural simulations, they are widely spread acrossthe ranges for most of the parameters varied. The GLUE methodology allows simulated prediction ranges for important variables to be presented as quantitative likelihood weighteduncertainty estimates rather than a single prediction for eachvariable over time.     

Simulating acidification and recovery processes in experimental catchments with the ILWAS model

The Integrated Lake Watershed Acidification Study (ILWAS) model was used to simulate soil discharge chemistry at two neighboring experimental catchments. One catchment underwent deacidification because of the artificial application of deacidified precipitation whereas the other catchment received unaltered acidic precipitation. Simulated results reproduce the observed seasonal dynamics in the concentrations of base cations, NO ₃ ^? , Al, and H₄SiO ₄ ⁰ in soil discharges for both catchments. Simulated results also indicate that the export flux of base cations was decreased by 30% at the deacidification catchment in response to the decrease in acid deposition. However, simulated SO ₄ ^2? concentrations show decreases that are about 40% more rapid than were observed. Simulated organic acid concentrations were also substantially lower than those observed at the deacidification catchment, indicating that organic matter decomposition processes were not correctly simulated. Acid-base budgets for both 5 and 50-yr simulations indicate that acid displacement by base cations through ion exchange is the principal process delaying recovery of runoff alkalinity, whereas SO ₄ ^2? desorption has a minor role. Silicate weathering is the dominant acid-consuming process at both catchments. Criteria proposed here for assessing forecast reliability include reproducing seasonal dynamics in discharge chemistry, providing numerically accurate chemical concentrations when compared to monitoring data, and correctly predicting deacidification rate and extent. The ILWAS model generally meets these criteria, indicating that the model can produce a reliable forecast of the effects of acid deposition on the acid-base chemistry of surface waters given sufficient temporal data for confident optimization of the calibrated variables in the model.     

Modelling of snowmelt erosion and sediment yield in a small low-mountain catchment in Germany

Temporal variability and spatial heterogeneity of surface runoff generation triggers the dynamics of source areas of sediment and sediment-associated nutrient transport. Reliable modelling of hydrological special situations i.e. snowmelt is of high importance for the quality of erosion and sediment yield modelling. Data from the research catchment Schäfertal demonstrate the individuality of snowmelt events in terms of runoff coefficient and delivery ratio. This 1.44 km² low mountain catchment is characterised by a high portion of arable land with a winter grain/winter rape crop rotation. The integrated winter erosion and nutrient load model (IWAN) considers these dynamic aspects by coupling a hydrological model with a sediment load model. Cell size of this raster-based approach is 10 × 10 m². Additionally, snowmelt rill erosion is simulated with a newly developed physically based model that is firstly applied on a catchment scale. A sensitivity analysis of this model system component demonstrates the plausibility of the model approach and the overall robustness of the model system IWAN. The results of the long-term hydrological modelling from 1991 to 2003 are reliable and form the basis for the simulation of six snowmelt events which were observed in the Schäfertal catchment. The estimated total runoff volumes for these events match the observations well. The modelled overland runoff coefficients vary from 0.001 to 0.72. The mean values of cell erosion, which were modelled with one set of parameters for all six events range from 0.0006 to 0.96 t ha^− 1. The total modelled erosion for the events with unfrozen soil and low amount of surface runoff is of a factor 50 below those with partly frozen soil. In addition to these distinctions, the major differences are caused by flow accumulation in shallow depressions in variable parts of the catchment. However, the validation of these results on the single event scale is restricted due to limited spatial data. Total simulated sediment yield at the catchment outlet was as high as 13.84 t which underestimates the observed values, with the exception of one event. Oversimplification of the modelled channel processes may be a reason. The temporal variability and spatial heterogeneity of the surface roughness parameter, which was identified to be sensitive, also causes uncertainty in the parameter estimation. Despite these findings, the model system IWAN was applied successfully on the catchment scale and the simulated results are reliable.     

Atmospheric deposition effects on the chemistry of a stream in Northeastern Georgia

A modeling study of the Allt a Mharcaidh catchment in the Cairngorm region of Scotland has been undertaken to investigate long term trends in acidification and model sensitivity to soil physical and chemical characteristics. The MAGIC model (Model of Acidification of Groundwater In Catchments) is used to demonstrate that the sulphate adsorption ability of the soil and quality and quantity of rainfall inputs have significant effects on model output. Optimal weathering rates and predicted present day ion concentrations in streamwater compare well with measured and observed values. The analysis shows that the catchment has become progressively acidified since pre-industrial times but major changes in stream acidity have yet to occur.     

Modelling of snowmelt erosion and sediment yield in a small low-mountain catchment in Germany

Temporal variability and spatial heterogeneity of surface runoff generation triggers the dynamics of source areas of sediment and sediment-associated nutrient transport. Reliable modelling of hydrological special situations i.e. snowmelt is of high importance for the quality of erosion and sediment yield modelling. Data from the research catchment Schäfertal demonstrate the individuality of snowmelt events in terms of runoff coefficient and delivery ratio. This 1.44 km² low mountain catchment is characterised by a high portion of arable land with a winter grain/winter rape crop rotation. The integrated winter erosion and nutrient load model (IWAN) considers these dynamic aspects by coupling a hydrological model with a sediment load model. Cell size of this raster-based approach is 10 × 10 m². Additionally, snowmelt rill erosion is simulated with a newly developed physically based model that is firstly applied on a catchment scale. A sensitivity analysis of this model system component demonstrates the plausibility of the model approach and the overall robustness of the model system IWAN. The results of the long-term hydrological modelling from 1991 to 2003 are reliable and form the basis for the simulation of six snowmelt events which were observed in the Schäfertal catchment. The estimated total runoff volumes for these events match the observations well. The modelled overland runoff coefficients vary from 0.001 to 0.72. The mean values of cell erosion, which were modelled with one set of parameters for all six events range from 0.0006 to 0.96 t ha^− 1. The total modelled erosion for the events with unfrozen soil and low amount of surface runoff is of a factor 50 below those with partly frozen soil. In addition to these distinctions, the major differences are caused by flow accumulation in shallow depressions in variable parts of the catchment. However, the validation of these results on the single event scale is restricted due to limited spatial data. Total simulated sediment yield at the catchment outlet was as high as 13.84 t which underestimates the observed values, with the exception of one event. Oversimplification of the modelled channel processes may be a reason. The temporal variability and spatial heterogeneity of the surface roughness parameter, which was identified to be sensitive, also causes uncertainty in the parameter estimation. Despite these findings, the model system IWAN was applied successfully on the catchment scale and the simulated results are reliable.     

基于SWAT模型的径流与土壤侵蚀过程模拟

将可表征土地利用变化对水文响应影响的分布式水文模型SWAT应用于密云水库潮河子流域,对模型的敏感性进行了分析,并采用1986～1991年下会水文站数据对其径流量和泥沙负荷进行了率定和验证,模型评估结果表明该模型对径流的模拟效果非常好,对泥沙的模拟效果较好。通过计算34个子流域内的侵蚀模数,对流域内的土壤侵蚀强度进行了分级,并分析了不同土地利用方式对产沙的影响。结果表明,流域内的土壤侵蚀主要发生在农田,为保障密云水库的水质应对其采取有效的水土保持措施。     

Remote-sensing data as an alternative input for the ‘STREAM’ runoff model

Water erosion of cropland constitutes an issue for natural environments along runoff flowpaths due to property damage by soil-laden water and the associated transfer of nutrients and pesticides. In the Pays de Caux region of northwestern France, the silty soils with crusting properties induce a high risk of runoff and erosion. Changes in agricultural practices, land use and landscape patterns appear to have increased the occurrence of erosion and mud flows over the past few decades.A runoff and erosion model called STREAM, applicable to single rainfall events at catchment scale, has been developed to simulate the impacts of land-use modifications. The model takes into account processes that degrade surface states when calculating infiltration rates, as well as agricultural aspects when computing the runoff circulation network. STREAM is based on an expert-system approach that focuses on the dominant processes whilst having only a few input parameters: three of these are used to determine the runoff circulation network, and the other four to calculate infiltration rates. Input nevertheless requires field observations, which restricts application of the model to small catchments.Satellite data covering large areas is considered as an alternative input for such a model, the main objectives being to adapt STREAM accordingly, and to compare the obtained results with field data. In view of previous work involving the extraction and validation of roughness indices using RADARSAT data, this study is based on RADARSAT and LANDSAT TM data collected during the winter of 1998.After adaptation to receive remote-sensing data, the resulting STREAM-TED model requires less input, namely (1) slope and orientation, (2) land-use classification from optical remote-sensing data, (3) roughness indices from radar remote-sensing data, and (4) previous rainfall.Runoff volumes at a gauged catchment outlet (Bourville in Upper Normandy, France) are simulated by four successive versions of the model ranging from the original STREAM to the adapted STREAM-TED. Predictions of the four versions are compared, and performance of the successive simulations is assessed in relation to measured values and according to five statistical indices.Predictions of runoff volume at the catchment outlet using STREAM-TED are similar to those using the original STREAM model, but with a tendency towards overestimation. The final STREAM-TED version is capable of identifying areas sensitive to runoff within a group of catchments and could be used as a planning decision tool in the implementation of conservation practices.     

A pragmatic approach to modelling soil and water conservation measures with a catchment scale erosion model

To reduce soil erosion, soil and water conservation (SWC) methods are often used. However, no method exists to model beforehand how implementing such measures will affect erosion at catchment scale. A method was developed to simulate the effects of SWC measures with catchment scale erosion models. The method was implemented by applying the LISEM model to an agricultural catchment on the slopes of Mt. Kenya. The method consisted of a field scale calibration based on P-factors, followed by application at catchment scale. This calibration included factors such as saturated conductivity, Manning's n, roughness and slope angle. It was found that using data on P-factors, such models can be calibrated to give acceptable predictions at pixel scale. However, P-factors were also found to vary with land use type and storm size. Besides, more data on the physical effectiveness of SWC measures are needed. At catchment scale, the effect of SWC was found to be different from that at pixel scale. Most SWC were simulated to be more effective at catchment scale, indicating additional infiltration during transport through the catchment to the outlet. However, slope corrections in case of terraces were found to be less effective at this scale. Nevertheless, a simulation for current land use with current SWC measures indicated that these SWC measures decrease runoff by 28% and erosion by 60%.     

冬小麦/夏玉米轮作中NO3-N在土壤剖面的累积及移动

通过田间试验研究了冬小麦 /夏玉米轮作中NO- 3 N在土壤剖面的累积及移动 ,结果表明 ,尿素施入旱地土壤后 ,硝化作用一般在 7d之内完成 ,NH 4 N只在施肥后的短期内保持较高浓度 ,其它时期NH 4 N含量基本在 1～ 3mgkg- 1 范围内 ,土壤剖面不同层次NH 4 N一般也低于 4mgkg- 1 ,NH 4 N的含量不能反映土壤有效氮的水平。土壤剖面中的NO- 3 N随施氮量的增加而显著升高。在低施氮量条件下 (N <12 0kghm- 2 ) ,NO- 3 N主要在 0～ 40cm土层内移动 ,但当施氮量高于N 2 40kghm- 2 时 ,冬小麦季即有相当数量的氮移出 0～ 10 0cm土体。NO- 3 N在土体中的移动存在着很大的年际变化 ,在干旱年份 ,即使夏玉米季 ,NO- 3 N向深层移动的可能性也很小。试验年份中 ,除 1999年夏玉米季发生了较严重的气体损失以外 (该季节特别干旱 ) ,其余季节损失的肥料氮主要以NO- 3 N的形式在深层土壤剖面中累积 ,这在两个试验点的结果相当一致。