ELECTRONIC SUPPLEMENTAL Materialfor Water Resources Management

ELECTRONIC SUPPLEMENTAL MATERIALfor Water Resources Management

Evaluating the importance of non-unique parameter sets on surface water quality variables under climate change conditions in a mesoscale agricultural watershed

Mehdi B.1,*, Schulz, K.2, Ludwig R.3, Ferber F.3, Lehner B.1

1McGill University, Department of Geography, 805 Sherbrooke Str. West, Montreal, H3A 0B9, Canada.

*Present address: 2University of Natural Resources and Life Sciences -Vienna, Institute for Water Management, Hydrology and Hydraulic Engineering, Muthgasse 18, 1190 Vienna, Austria.

3Ludwig Maximilians University, Department of Geography, Luisenstr. 37, 80333 Munich, Germany.

Corresponding author: B. Mehdi, Email:

1

Figure 1S. Observed and simulated daily streamflow, at Treutchlingen, forthe calibration and validation periods.

1

Figure 2S.Mean monthly streamflow at Treuchtlingen for 30 years, for SWAT simulated with observed climate data from 1970-2000 (blue line), and with the reference climate (1970-2000) and the best run (grey lines), and with the future climate (2041-2070) and the best run (red lines).

Figure 3S.Monthly means for the future period 2041-2070 with SWAT simulated using the best run with the suite of climate simulations (red lines) and using the non-unique parameters sets (black lines).

1

Table 1S. Sources of data used for setting up the Altmühl watershed in ArcSWAT

Data / Source / Database / type of data / Time period / Scale / resolution / Reference
Digital Elevation Model / BavarianGeodetic Survey / Raster / 50 m raster from 1:25000 / DGM25, BayerischeVermessungsverwaltung
Land cover for Europe / European Environmental Agency / Raster / 2006 / 50 m / European Environmental Agency
CORINE land cover v.15,
Protected areas / European Environmental Agency / Shape file / 2000 / Shape file / NATURA 2000, EuropeanEnvironmentalAgency

General land use / Bavarian State Office for the Environment / Shape file / 2010 / Shape file / VEKTOR 500, BayerischeVermessungsverwaltung
Forest cover / Bayerische Landesanstalt für Land und Forstwirtschaft / Shape file / 2009 / Shape file / FFH Lebensraumtypenkarte, Bayerische Forstverwaltung
Agricultural land use / Bayerische Landesanstalt für Landwirtschaft / Shape file / 2008 / Field level / FNN Feldstückgeometrien
Soildata / Bavarian Geodetic Survey;
and
GLOWA-Danube / Shape file;
Raster / Shape file;
1:1000000 / Bodenschätzung Bayern v1.4, Bayerische Vermessungsverwaltung;

River network / Bavarian State Office for the Environment / Shape file / 2005 / 1:25000 / WGN25, Digitales Gewässerverzeichnis,
Hydrological and water quality data / Bavarian State Office for the Environment / Excel / 1948-2011 / Gauges Thann, Treuchtlingen / Wasserwirtscahftsamt Ansbach;
Referat 85, Bayerisches Landesamt für Umwelt
Farm type and size / Bavarian State Office for Statistics and Dataprocessing / Excel / Farm level / Landesamt für Statistik und Datenverarbeitung

Crop management / Bavarian State Office for Agriculture,
Association for Technology and Structures in Agriculture / Report
Reports / 2008
1995, 2009 / Crop level / Bayerische Landesanstalt für Landwirtschaft
KTBL
Fruit management / Bavarian State Office for Viticulture and Horticulture / Report / 2007-2011 / Crop level / Bayerische Landesanstalt für Weinbau und Gartenbau

1

Table 2S. Climate model simulations considered in thisstudy. All simulations are bias corrected. Regional climate models: Rossby Centre Regional Atmospheric Climate Model (RCA), Regional Atmospheric Climate Model (RACMO), Canadian Regional Climate Model (CRCM). Driving Global climate models: Bergen Climate Model (BCM), ECHAM version 5 (members 1, 2 and 3), and Hadley Centre Coupled Model (HadCM3),Canadian General Circulation Model (CGCM).

RCM / Driving GCM / SRES / Grid size
(km) / Name of simulation
RCA / BCM / A1B / 50 / RCA-BCM-50K
RCA / ECHAM5-r3 / A1B / 50 / RCA-ECM-50K
RCA / HadCM3Q3 / A1B / 50 / RCA-HCM-50K
RACMO2 / ECHAM5-r1 / A1B / 50 / RAC-ECM-MB1-50K
RACMO2 / ECHAM5-r2 / A1B / 50 / RAC-ECM-MB2-50K
RACMO2 / ECHAM5-r3 / A1B / 50 / RAC-ECM-MB3-50K
CRCM 4.2.3 / CGCM3 / A2 / 45 / CRC-CGC-45K

Table 3S. Parameter set solutions that i) best meet the objective criterion (“best run”)and ii) that also meet the minimum objective criterion (“non-unique parameter sets”).

Parameter / Best run / Non-unique parameter sets
r__CN2.mgt / 0.05 / 0.03 / 0.02 / 0.05 / 0.06 / 0.06
v__GWQMN.gw / 1.80 / 1.36 / 1.28 / 1.48 / 1.95 / 2.35
v__ESCO.hru / 0.85 / 0.90 / 0.92 / 0.83 / 0.89 / 0.88
v__CH_N2.rte / 0.16 / 0.11 / 0.10 / 0.08 / 0.11 / 0.13
v__CH_K2.rte / 106.38 / 137.41 / 139.06 / 108.10 / 134.33 / 133.39
v__ALPHA_BNK.rte / 0.62 / 0.51 / 0.61 / 0.56 / 0.31 / 0.64
r__SOL_AWC(1).sol / 0.24 / 0.05 / 0.25 / -0.01 / -0.03 / 0.00
v__SURLAG.bsn / 10.13 / 11.51 / 9.54 / 10.46 / 10.49 / 13.13
v__CH_K1.sub / 4.23 / 2.97 / 4.45 / 2.57 / 3.08 / 4.04
v__GW_REVAP.gw / 0.17 / 0.18 / 0.20 / 0.17 / 0.16 / 0.19
v__CANMX.hru / 1.02 / 13.02 / 5.60 / 8.69 / 0.77 / 5.31
v__TIMP.bsn / 0.69 / 0.72 / 0.76 / 0.66 / 0.65 / 0.66
v__SNOCOVMX.bsn / 314.00 / 41.77 / 240.85 / 247.65 / 292.82 / 11.79
v__SNO50COV.bsn / 0.32 / 0.35 / 0.32 / 0.24 / 0.14 / 0.01
v__SFTMP.bsn / 3.99 / 3.68 / 4.25 / -0.44 / -0.96 / 0.77
v__SMTMP.bsn / -0.83 / -2.15 / -3.45 / -2.34 / -3.04 / -1.82
v__SMFMX.bsn / 9.98 / 4.49 / 4.08 / 8.46 / 3.95 / 3.88
v__SMFMN.bsn / 6.87 / 7.63 / 3.97 / 2.64 / 7.47 / 7.41
v__CH_N1.sub / 0.24 / 0.24 / 0.29 / 0.25 / 0.22 / 0.21
v__SOL_ORGP(1).chm / 262.14 / 247.95 / 185.05 / 149.43 / 175.87 / 187.28
v__SOL_SOLP(1).chm / 39.27 / 48.32 / 46.22 / 55.82 / 99.81 / 90.48
v__ERORGP.hru / 1.49 / 1.47 / 1.02 / 1.00 / 1.29 / 1.32
v__PHOSKD.bsn / 174.42 / 178.35 / 140.60 / 158.63 / 159.10 / 192.27
v__K_P.wwq / 0.02 / 0.02 / 0.02 / 0.01 / 0.02 / 0.02
v__RS5.swq / 0.047 / 0.05 / 0.04 / 0.02 / 0.02 / 0.06
v__BC4.swq / 0.35 / 0.39 / 0.40 / 0.51 / 0.38 / 0.48
v__RSDIN.hru / 3657.87 / 2551.62 / 1893.13 / 1517.80 / 5264.56 / 385.21

1

Table 4S.Simulated streamflow (m3/s) at Treuchtlingen. Independent t-tests compare the mean of the future climate simulations to the mean of the reference climate simulations (1970-2000). An asterisk indicates a rejection of the test, the value after the asterisk indicates the difference in m3/s.

Best run / Non-unique parameter sets
Referenceclimate
(1970-2000) / Future climate
(2041-2070) / Future climate
(2041-2070)
Mean
μ0 / Std. dev. / Percentile / Mean
μ1 / Std. dev. / Min / Max / Percentile / t-test μo=μ1
(p<0.05) / Mean
μ2 / Std. dev. / Min / Max / Percentile / t-test μo=μ2
(p<0.05)
10th / 50th / 90th / 10th / 50th / 90th / 10th / 50th / 90th
J / 9.5 / 5.4 / 2.8 / 8.9 / 16.4 / 9.5 / 4.9 / 0.22 / 27.7 / 2.9 / 9.3 / 16.2 / 9.5 / 5.0 / 0.17 / 32.9 / 2.9 / 9.3 / 16.1
F / 10.7 / 5.2 / 4.6 / 9.6 / 17.9 / 11.0 / 5.2 / 0.27 / 25.0 / 4.5 / 10.7 / 18.0 / 10.9 / 5.3 / 0.11 / 31.0 / 4.3 / 10.4 / 17.9
M / 9.9 / 5.0 / 4.6 / 9.3 / 15.6 / 9.8 / 4.5 / 0.12 / 25.5 / 3.7 / 9.8 / 14.8 / 9.3 / 4.5 / 0.09 / 25.5 / 3.3 / 9.2 / 14.8
A / 6.5 / 3.5 / 2.4 / 6.4 / 11.4 / 6.9 / 4.3 / 0.05 / 21.8 / 1.9 / 6.5 / 12.6 / 6.3 / 4.1 / 0.03 / 21.8 / 1.5 / 5.8 / 11.9
M / 3.6 / 3.7 / 0.4 / 2.2 / 8.5 / 4.7 / 3.8 / 0.13 / 19.4 / 0.8 / 3.6 / 10.2 / *1.1 / 4.3 / 3.7 / 0.06 / 19.4 / 0.7 / 3.3 / 9.8 / *0.73
J / 2.7 / 2.7 / 0.4 / 1.8 / 6.7 / 3.3 / 3.4 / 0.07 / 15.8 / 0.5 / 2.1 / 8.4 / *0.6 / 3.1 / 3.2 / 0.04 / 17.7 / 0.4 / 1.8 / 7.5
J / 2.6 / 2.9 / 0.3 / 1.5 / 5.7 / 2.9 / 3.7 / 0.12 / 29.5 / 0.3 / 1.7 / 7.6 / 2.6 / 3.5 / 0.08 / 30.0 / 0.3 / 1.3 / 6.7
A / 2.0 / 2.4 / 0.2 / 1.0 / 5.6 / 2.0 / 2.8 / 0.03 / 21.5 / 0.2 / 1.1 / 5.3 / 1.9 / 2.6 / 0.01 / 21.5 / 0.2 / 0.9 / 4.9
S / 1.7 / 2.4 / 0.2 / 0.8 / 5.0 / 1.4 / 2.1 / 0.05 / 14.8 / 0.2 / 0.7 / 3.7 / 1.3 / 2.0 / 0.02 / 16.3 / 0.1 / 0.6 / 3.7 / *-0.39
O / 2.1 / 2.6 / 0.2 / 1.0 / 5.7 / 2.0 / 2.3 / 0.06 / 11.0 / 0.3 / 1.1 / 5.6 / 1.9 / 2.1 / 0.03 / 11.7 / 0.2 / 1.0 / 5.2
N / 3.3 / 3.2 / 0.4 / 2.2 / 8.3 / 3.3 / 3.3 / 0.06 / 19.6 / 0.4 / 2.2 / 7.2 / 3.2 / 3.2 / 0.03 / 21.5 / 0.4 / 2.1 / 7.0
D / 6.8 / 4.7 / 1.8 / 6.2 / 13.1 / 6.4 / 4.6 / 0.07 / 25.7 / 1.2 / 5.5 / 13.2 / 6.3 / 4.7 / 0.05 / 28.5 / 1.2 / 5.4 / 13.0

1

Table5S.Quantiles for simulated streamflow (m3/s) at Treuchtlingen using the mean of 7 climate simulationscomparing the best run approach for the reference climate simulations, the best run approach for future climate simulation, and the non-unique approach for the future climate simulations.

Season / Std. dev., Quantiles / Reference climate / Future climate
Best run
(1970-2000) / Best run
(2041-2070) / Non-unique
(2041-2070)
DJF / Std. dev. / 5.33 / 5.28 / 5.37
25 / 5.11 / 4.75 / 4.67
50 / 8.10 / 8.57 / 8.40
75 / 11.77 / 12.30 / 12.20
Interquartile range / 6.66 / 7.55 / 7.53
MAM / Std. dev. / 4.84 / 4.68 / 4.60
25 / 2.92 / 3.22 / 2.79
50 / 6.27 / 6.51 / 5.90
75 / 9.56 / 10.39 / 9.70
Interquartile range / 6.64 / 7.17 / 6.91
JJA / Std. dev. / 2.67 / 3.36 / 3.15
25 / 0.57 / 0.60 / 0.54
50 / 1.43 / 1.58 / 1.32
75 / 3.42 / 3.58 / 3.21
Interquartile range / 2.85 / 2.98 / 2.67
SON / Std. dev. / 2.84 / 2.72 / 2.61
25 / 0.44 / 0.44 / 0.43
50 / 1.27 / 1.15 / 1.14
75 / 3.13 / 3.07 / 2.92
Interquartile range / 2.69 / 2.63 / 2.49

1

Table 6S.Simulated nitrate nitrogen (kg/ha). Independent t-tests compare the mean of the future climate simulations to the mean of the reference climate simulations (1970-2000). An asterisk indicates a rejection of the test, the value after the asterisk indicates the difference in kg/ha.

Best run / Non-unique parameter
Referenceclimate
(1970-2000) / Future climate
(2041-2070) / Future climate
(2041-2070)
Mean
μ0 / Std. dev. / Percentile / Mean
μ1 / Std. dev. / Min / Max / Percentile / t-test μo=μ1
(p<0.05) / Mean
μ2 / Std. dev. / Min / Max / Percentile / t-test μo=μ2
(p<0.05)
10th / 50th / 90th / 10th / 50th / 90th / 10th / 50th / 90th
J / 1.73 / 0.90 / 0.60 / 1.68 / 3.00 / 2.34 / 1.21 / 0.02 / 6.9 / 1.03 / 2.22 / 3.89 / *0.62 / 2.28 / 1.19 / 0.02 / 7.4 / 0.98 / 2.13 / 3.80 / *0.55
F / 1.64 / 0.77 / 0.67 / 1.62 / 2.71 / 2.13 / 0.96 / 0.04 / 5.3 / 0.99 / 2.05 / 3.33 / *0.49 / 2.05 / 0.93 / 0.01 / 6.3 / 0.95 / 1.98 / 3.26 / *0.41
M / 1.97 / 0.90 / 0.86 / 1.91 / 3.13 / 2.10 / 1.02 / 0.02 / 6.1 / 0.92 / 2.13 / 3.42 / 1.93 / 1.01 / 0.01 / 6.5 / 0.76 / 1.88 / 3.17
A / 1.49 / 1.04 / 0.32 / 1.34 / 3.05 / 1.42 / 1.15 / 0.00 / 5.8 / 0.24 / 1.20 / 3.06 / 1.25 / 1.07 / 0.0 / 5.8 / 0.18 / 0.99 / 2.71 / *-0.24
M / 0.56 / 0.79 / 0.04 / 0.24 / 1.56 / 0.79 / 1.05 / 0.01 / 6.6 / 0.06 / 0.38 / 2.03 / *0.23 / 0.73 / 1.02 / 0.0 / 6.6 / 0.05 / 0.35 / 1.75 / *0.18
J / 0.24 / 0.33 / 0.02 / 0.13 / 0.57 / 0.40 / 0.69 / 0.0 / 4.3 / 0.03 / 0.16 / 1.13 / *0.17 / 0.37 / 0.67 / 0.0 / 4.5 / 0.02 / 0.13 / 1.05 / *0.14
J / 0.21 / 0.29 / 0.02 / 0.10 / 0.52 / 0.34 / 0.67 / 0.01 / 4.9 / 0.01 / 0.11 / 0.93 / *0.14 / 0.32 / 0.66 / 0.0 / 5.4 / 0.01 / 0.08 / 0.83 / *0.11
A / 0.21 / 0.36 / 0.01 / 0.07 / 0.64 / 0.29 / 0.51 / 0.0 / 3.0 / 0.01 / 0.07 / 0.92 / 0.30 / 0.53 / 0.0 / 3.9 / 0.01 / 0.06 / 0.92 / *0.08
S / 0.26 / 0.45 / 0.01 / 0.07 / 0.77 / 0.28 / 0.48 / 0.0 / 3.1 / 0.01 / 0.06 / 0.87 / 0.30 / 0.49 / 0.0 / 3.2 / 0.01 / 0.07 / 1.00
O / 0.42 / 0.54 / 0.01 / 0.21 / 1.17 / 0.66 / 0.85 / 0.0 / 5.5 / 0.02 / 0.35 / 1.76 / *0.24 / 0.67 / 0.84 / 0.0 / 6.1 / 0.02 / 0.40 / 1.74 / *0.25
N / 0.72 / 0.69 / 0.05 / 0.49 / 1.74 / 1.13 / 1.02 / 0.0 / 5.1 / 0.06 / 0.86 / 2.54 / *0.41 / 1.11 / 0.97 / 0.0 / 5.5 / 0.08 / 0.86 / 2.49 / *0.39
D / 1.41 / 0.95 / 0.36 / 1.33 / 2.58 / 1.82 / 1.23 / 0.01 / 7.3 / 0.45 / 1.53 / 3.41 / *0.41 / 1.76 / 1.21 / 0.0 / 8.1 / 0.46 / 1.46 / 3.31 / *0.35

1

Table 7S. Quantiles for simulated nitrate nitrogen (kg/ha) at Treuchtlingen, comparing the best run approach for the reference climate simulations, the best run approach for future climate simulations, and the non-unique approach for the future climate simulations.

Season / Std. dev., Quantiles / Reference climate / Future climate
Best run
(1970-2000) / Best run
(2041-2070) / Non-unique
(2041-2070)
DJF / Std. dev. / 0.89 / 1.16 / 1.14
25 / 0.94 / 1.12 / 1.19
50 / 1.51 / 2.00 / 1.91
75 / 2.16 / 2.76 / 2.66
Interquartile range / 1.22 / 1.64 / 1.47
MAM / Std. dev. / 1.10 / 1.20 / 1.15
25 / 0.33 / 0.42 / 0.35
50 / 1.11 / 1.21 / 1.03
75 / 2.05 / 2.22 / 2.00
Interquartile range / 1.72 / 1.8 / 1.65
JJA / Std. dev. / 0.33 / 0.63 / 0.62
25 / 0.04 / 0.03 / 0.03
50 / 0.10 / 0.11 / 0.09
75 / 0.25 / 0.35 / 0.32
Interquartile range / 0.21 / 0.32 / 0.29
SON / Std. dev. / 0.60 / 0.88 / 0.86
25 / 0.05 / 0.05 / 0.06
50 / 0.23 / 0.33 / 0.36
75 / 0.66 / 1.04 / 1.05
Interquartile range / 0.61 / 0.99 / 0.99

1

Table 8S.Simulated total phosphorus (kg/ha). Independent t-tests compare the mean of the future climate simulations to the mean of the reference climate simulations (1970-2000). An asterisk indicates a rejection of the test, the value after the asterisk indicates the difference in kg/ha.

Best run / Non-unique parameter
Referenceclimate
(1970-2000) / Future climate
(2041-2070) / Future climate
(2041-2070)
Mean
μ0 / Std. dev. / Percentile / Mean
μ1 / Std. dev. / Min / Max / Percentile / t-test μo=μ1
(p<0.05) / Mean
μ2 / Std. dev. / Min / Max / Percentile / t-test μo=μ2
(p<0.05)
10th / 50th / 90th / 10th / 50th / 90th / 10th / 50th / 90th
J / 0.11 / 0.12 / 0.005 / 0.060 / 0.27 / 0.06 / 0.06 / 0.0 / 0.37 / 0.007 / 0.040 / 0.16 / *-0.045 / 0.06 / 0.07 / 0.0 / 0.54 / 0.008 / 0.040 / 0.15 / *-0.044
F / 0.11 / 0.11 / 0.010 / 0.077 / 0.25 / 0.07 / 0.09 / 0.0 / 0.47 / 0.003 / 0.041 / 0.18 / *-0.038 / 0.07 / 0.09 / 0.0 / 0.63 / 0.005 / 0.038 / 0.19 / *-0.037
M / 0.08 / 0.12 / 0.006 / 0.040 / 0.19 / 0.04 / 0.05 / 0.0 / 0.39 / 0.005 / 0.024 / 0.10 / *-0.034 / 0.04 / 0.05 / 0.0 / 0.45 / 0.005 / 0.024 / 0.09 / *-0.037
A / 0.02 / 0.02 / 0.002 / 0.014 / 0.041 / 0.03 / 0.03 / 0.0 / 0.20 / 0.005 / 0.019 / 0.07 / *0.011 / 0.03 / 0.03 / 0.0 / 0.25 / 0.005 / 0.020 / 0.06 / *0.010
M / 0.02 / 0.03 / 0.003 / 0.014 / 0.049 / 0.03 / 0.04 / 0.0 / 0.22 / 0.005 / 0.022 / 0.08 / *0.010 / 0.03 / 0.03 / 0.0 / 0.23 / 0.005 / 0.022 / 0.07 / *0.010
J / 0.02 / 0.03 / 0.003 / 0.013 / 0.049 / 0.03 / 0.03 / 0.0 / 0.15 / 0.003 / 0.014 / 0.07 / 0.02 / 0.03 / 0.0 / 0.25 / 0.015 / 0.015 / 0.06
J / 0.02 / 0.03 / 0.003 / 0.012 / 0.045 / 0.03 / 0.06 / 0.0 / 0.42 / 0.003 / 0.011 / 0.07 / 0.03 / 0.04 / 0.0 / 0.42 / 0.002 / 0.011 / 0.05 / *0.005
A / 0.01 / 0.02 / 0.002 / 0.007 / 0.031 / 0.01 / 0.02 / 0.0 / 0.23 / 0.001 / 0.006 / 0.03 / 0.01 / 0.02 / 0.0 / 0.23 / 0.001 / 0.006 / 0.03
S / 0.01 / 0.02 / 0.001 / 0.005 / 0.021 / 0.01 / 0.01 / 0.0 / 0.10 / 0.001 / 0.005 / 0.02 / 0.01 / 0.01 / 0.0 / 0.12 / 0.001 / 0.005 / 0.02
O / 0.01 / 0.02 / 0.0009 / 0.007 / 0.029 / 0.02 / 0.03 / 0.0 / 0.18 / 0.001 / 0.008 / 0.04 / *0.005 / 0.02 / 0.03 / 0.0 / 0.23 / 0.001 / 0.008 / 0.04 / *0.004
N / 0.02 / 0.024 / 0.001 / 0.008 / 0.038 / 0.02 / 0.03 / 0.0 / 0.22 / 0.002 / 0.009 / 0.05 / 0.02 / 0.03 / 0.0 / 0.24 / 0.002 / 0.009 / 0.04
D / 0.06 / 0.078 / 0.006 / 0.031 / 0.17 / 0.05 / 0.06 / 0.0 / 0.37 / 0.004 / 0.023 / 0.12 / *-0.015 / 0.05 / 0.07 / 0.0 / 0.53 / 0.004 / 0.025 / 0.13 / *-0.013

1

Table 9S.Quantiles for simulated total phosphorus (kg/ha) at Treuchtlingen, comparing the best run approach for the reference climate simulations, the best run approach for future climate simulations, and the non-unique approach for the future climate simulations.

Season / Std. dev., Quantiles / Reference climate / Future climate
Best run
(1970-2000) / Best run
(2041-2070) / Non-unique
(2041-2070)
DJF / Std. dev. / 0.11 / 0.072 / 0.078
25 / 0.019 / 0.014 / 0.015
50 / 0.055 / 0.035 / 0.034
75 / 0.13 / 0.077 / 0.077
Interquartile range / 0.11 / 0.063 / 0.062
MAM / Std. dev. / 0.074 / 0.042 / 0.038
25 / 0.0069 / 0.0098 / 0.010
50 / 0.018 / 0.022 / 0.022
75 / 0.040 / 0.040 / 0.039
Interquartile range / 0.033 / 0.030 / 0.029
JJA / Std. dev. / 0.024 / 0.034 / 0.032
25 / 0.0047 / 0.0043 / 0.0042
50 / 0.011 / 0.0098 / 0.010
75 / 0.021 / 0.024 / 0.023
Interquartile range / 0.016 / 0.020 / 0.019
SON / Std. dev. / 0.021 / 0.027 / 0.025
25 / 0.0026 / 0.0030 / 0.0032
50 / 0.0059 / 0.0073 / 0.0075
75 / 0.015 / 0.017 / 0.016
Interquartile range / 0.012 / 0.014 / 0.013

1