Hydrology Author and Earth System Sciences, 8(4), 858859 (2004)
© EGU
Errata The authors regret that errors in the equations in the original manuscripts of two papers published in previous issues of HESS escaped detection at the proof stage so the following Errata should be consulted in using these papers. Wade, A.J., Durand, P., Beaujouan, V., Wessel,W.W., Raat, K.J., Whitehead, P.G., Butterfield, D., Rankinen,K. and Lepisto, A., 2002a. A nitrogen model for European catchments:INCA. new model structure and equations. Hydrol. Earth Syst. Sci. 6, 559582. 1.
Table 3. The symbol for the Base Flow Index is b.
2. Table 4. The unit of soil temperature is degreescentigrade, (oC). 3. C10
Table 5. The in-stream nitrification rate is parameter,
4. Table 5. The in-stream denitrification rate is parameter, C11. 5.
x x .86400 x5 U 3 .100 1 5 10 6 C 7 S1 S 2 Vr x11 Vr x11
7.
x5 10 6 Vr x11
x3 6 Denitrification = C1S1 V x 10 r 11
Change in NO3 mass input into reach, x25 (kg N)
858
x 22 x 23 .86400 x 29
(B.3)
Change in NH4 mass output from reach, x28 (kg N)
dx28 dt
x 22 x 24 .86400 x 29
(B.4)
Change in reach volume, x29 (m3)
dx 29 dt
( S 4 x 22 ).86400
(B.5)
Change in water flow input to reach, x30 (m3)
dx30 dt
S 4 .86400
(B.6)
Change in water flow output from reach, x31 (m3)
dx33 dt
8. Appendix B should read Instream component: equations to track the water volume and N mass balance
S5
dx 27 dt
x 22 .86400
(B.7)
C11 a5,t 1 x 29 1000
(B.8)
Accumulated N mass associated with nitrification, x33 (kg N), is calculated by integrating the change in nitrification, dx33/dt
Eqn. (17) should read
dx 25 dt
Change in NO3 mass output from reach, x27 (kg N)
dx32 dt
Eqn. (16) should read
Immobilisation = C 6 S1
(B.2)
S6
Accumulated N mass associated with denitrification, x32 (kg N), is calculated by integrating the change in denitrification, dx32/dt
x5 x5 10 6 C5 S1100 C 6 S1 10 6 C 4 S1 Vr x11 Vr x11
6.
dx 26 dt
dx 31 dt
Eqn. (8) should read dx5 dt
Change in NH4 mass input into reach, x26 (kg N)
1000
(B.9)
The N mass-balance within the in-stream component is calculated and the results are displayed on the load charts in INCAv1.6. The input to reach, i is calculated as: Inputi
(B.1)
C10 a 6,t 1 x 29
x25,i x26,i
The output from reach, i is calculated as: Outputi x27 ,i x28,i
(B.10)
(B.11)
Errata The N stored in reach, i is calculated as:
Storagei
x23,i x24,i
(B.12)
where the terms are as defined previously. The initial N mass stored in reach, i is calculated as: Initiali x23, 0,i x24,0,i (B.13) where x23,0,i = NO3 stored in the reach at time, t = 0 (kg N); x24,0,i = NH4 stored in the reach at time, t = 0 (kg N); The user supplies all the initial values as input.
The N mass-balance for reach, i is calculated as:
Bali
Initiali Inputi Outputi Storagei (B.14)
Thus, if mass-balance is achieved then the balance will equal zero. The water balance for reach, i is expressed as
Bali
x29, 0,i x30,i x31,i x29,i
(B.15)
where x29,0,i is the initial volume of reach, i, and all the other terms are as previously defined.
Wade, A.J., Whitehead, P.G. and Butterfield, D., 2002d. The Integrated Catchments Model of Phosphorus Dynamics (INCA-P), a new approach for multiple source assessment in heteorogenous river systems: model structure and equations. Hydrol. Earth Syst. Sci. 6, 583 606. 1.
Eqn (28) should read
dx28 dt
1 dx27 Lw dt
(28)
2.Table 8 and Eqn (33). rs is the bed sediment particle density (kg m-3). The symbol is actually given as rs in the description of Eqn (33), even though rs is correctly used in the equation itself. 3.Eqn (37) should read dx34 dt
gain ½ C19 C (a12 a11 ) ® ¾ loss 1000 ¿ ¯ (U 20 ) C22C11T M x31U11C12 a12 Lw C20 a11 (C13 a12 )(C12 x31 )1000 1000 4
(37)
Fortunately, the errors are in these publications only and do not change any of the main results or conclusions in these papers or related work. A. J. Wade Aquatic Environments Research Centre University of Reading Reading RG6 6AB UK Email:
[email protected]
859
