ISRN Soil Science

Volume 2013 (2013), Article ID 362895, 7 pages

http://dx.doi.org/10.1155/2013/362895

Research Article

Temporal and Spatial Variability of Water Surplus in Ontario, Canada

D. Murray Brown, Humaira Dadfar, David J. Fallow, Robert J. Gordon, John D. Lauzon, and Gary W. Parkin

School of Environmental Sciences, University of Guelph, Guelph, ON, Canada, N1G 2W1

Received 17 September 2013; Accepted 20 October 2013

Academic Editors: L. A. Dawson and R. N. Lerch

Copyright © 2013 D. Murray Brown et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The temporal variability in estimated water surplus in 12 climatic regions of the province of Ontario, Canada, and its spatial distribution throughout most of the province are discussed in this paper. Surplus water is that which results from precipitation that runs off the land surface and that which drains through the soil profile to the water table and through subsurface drainage. A one-dimensional, deterministic model (DRAINMOD) that simulates soil water flow, including plant uptake, evapotranspiration, and freeze/thaw conditions, was used to estimate the water surplus. Simulations were performed using daily climatic data from January 1954 to December 2001 for each region. A reference corn crop and the predominant local soil conditions in each region, with the hydraulic properties for each layer in the soil profile, were used as model inputs. There was considerable year-to-year variability in annual water surplus in all regions caused by both precipitation and soil conditions. It was the least (~150 mm) in three regions and it exceeded 350 mm in another three regions, where winter snowfall is the greatest as a result of these regions being in the lea of one of the Great Lakes. The variability in water surplus generally increased as average water surplus increased.

1. Introduction

At the end of most growing seasons in Ontario, soil water has been partially depleted from the plant root zone and deeper through evaporation from the soil surface and transpiration from plants. The growing season generally ends by mid-November in southwestern Ontario, by the end of October in the rest of southern Ontario, and by mid-October in northern Ontario as discussed in [1, 2]. As a result, from November to March in southwestern Ontario, from November to mid-April in the rest of southern Ontario, and from mid-October to April in northern Ontario, the soil water is replenished by precipitation and a water surplus occurs, as evaporation is less than precipitation. This surplus water either drains through the soil profile to the water table or through tile drainage lines (referred to collectively as deep drainage (DD)) or runs off directly into surface water (referred to as runoff (RO)), especially in the late winter and spring months as discussed in [3–5]. Large annual variability in soil water necessitates the use of long-term studies to adequately describe seasonal changes in water surplus as discussed in [6]. Since long-term field data are generally not available, modeling studies are necessary as discussed elsewhere [7]. The objective of this paper is to determine the variability in water surplus both temporally and spatially across Ontario from long-term modeled data using the soil hydrology model DRAINMOD. In this study, there was no attempt to verify the DRAINMOD estimates of water surplus; nevertheless, verification of DRAINMOD’s satisfactory performance to simulate soil hydrology under the cold conditions of southern Canada has been performed elsewhere in [8, 9]. As well, the emphasis of the study is on comparing estimates of water surplus between different regions of Ontario under different soil profile conditions, not on using model output to truly represent soil hydrology under a specific set of climate, crop, and soil conditions at any given location and time within each region. Nevertheless, the general hypothesis of this research is that water surplus will vary substantially between regions and from year to year.

2. Materials and Methods

Daily climatic data (precipitation, maximum, minimum, and dew point temperatures, sunshine hours or solar radiation, and wind speed) were assembled for ten sites in southern Ontario and two sites in northern Ontario for 48 years from 1954 to 2001 (Table 1 and Figure 1). The sites were selected based on available climatic data, representing one of the climatic regions as defined in [1, 2]. The DRAINMOD 6.0 model as discussed elsewhere [8, 10–12] was applied to the daily climatic data for the 48 years at all 12 sites. Table 1 gives average annual precipitation () and estimated evapotranspiration (ET) for each site.

Table 1: Average annual precipitation (), evapotranspiration (ET), soil order, series, and textural class for local soils in each climatic region for the 12 sites.

Site
Latitude
Longitude
P (mm)
ET (mm)
Soil order, series and texture

Emo
48°38′
93°48′
724
564
Gleysol, Emo clay

Guelph
41°31′
80°14′
858
559
Luvisol, Guelph loam

Harrow
42°02′
82°54′
881
730
Gleysol, Brookston clay loam

Kapuskasing
49°24′
82°26′
844
458
Gleysol, Ryland silty clay loam

London
43°02′
81°09′
942
690
Luvisol, Bryanston silt loam

Mount Forest
43°59′
80°45′
991
586
Luvisol, Listowel silt loam

Ottawa
45°23′
75°43′
865
564
Gleysol, north Gower clay loam

Peterborough
44°14′
78°21′
780
599
Luvisol, Otonabee loam

Renfrew
45°29′
76°42′
803
572
Luvisol, Renfrew silty clay

Smithfield
44°05′
77°40′
931
574
Brunisol, Pontypool sand

Vineland
43°11′
79°24′
836
711
Luvisol, Oneida clay

Wiarton
44°45′
81°06′
989
701
Luvisol, Harkaway silt loam

Figure 1: Location of the 12 stations in the climatic regions of southern and northern Ontario.

Precipitation was obtained from daily climatic data and ET was estimated based on the Thornthwaite (1948) approach for estimating potential evapotranspiration (PET) calibrated to long-term average estimates of actual ET based on the Penman and Monteith method as discussed in [3, 4, 13]. This approach was required because monthly calibration factors required in DRAINMOD for determining Thornthwaite’s PET were not available for each site.

DRAINMOD as discussed in [11, 12] is a field scale deterministic hydrologic model based on the water balance technique for the soil surface and for a section of soil located midway between the adjacent drains extending from the soil surface to an impermeable layer. It estimates infiltration, ET, RO, subsurface drainage, and seepage on an hourly or daily basis for long periods of climatic records (e.g., 50 years). Recently, it was modified to include freezing and thawing, and snowmelt components to enhance its application in cold climates as discussed elsewhere [10, 14]. The model has been successfully tested and used under a wide variety of conditions [15–19] including climates similar to Ontario [8, 9, 20].

In terms of the lower boundary condition chosen for running DRAINMOD 6.0, no vertical or lateral seepage was allowed since the required model inputs were not available. Tile drainage was chosen to represent the lower boundary condition because it is commonly installed in Ontario and hence its required DRAINMOD inputs are well known. Therefore, DD was represented by the model estimated tile drainage. Depth and spacing of the 10 cm diameter tile drains were set as 100 cm and 15 m, respectively, since these are average values for most subsurface drainage systems in Ontario. DRAINMOD simulations were conducted with a representative corn crop at each site over 48 years (1954–2001). Corn was chosen since it is grown throughout most of the study area and its growth parameters required to run DRAINMOD were available. Different crops were not selected according to those typically grown in each region because, as discussed above, the emphasis of this study is on comparing water surplus between regions and not crops. Corn planting dates varied slightly each year depending on springtime temperature as discussed in [21].

A local soil type, including soil horizons and their depth, that occurs in each region (Table 1, Figure 1), where the climatic site data were recorded, was used to determine water surplus (sum of RO and DD) for that site and represent the region. No attempt was made to represent soil spatial variability within a region, only between the 12 different modeled regions. In addition, DRAINMOD 6.0 was applied to clay (Oneida clay), loam (Otonabee loam), and sandy (Pontypool sand) soil profiles at all 12 sites to determine their effect on water surplus. Details on soil layers, hydraulic properties, and climate and crop inputs discussed above are given in [3, 4].

3. Results and Discussion

3.1. Temporal Variations

The low DD with clay or silty clay soils for each site resulted in the standard deviation () of the year-to-year variability being nearly equal or greater than the average annual DD (Table 2). Whereas, when sandy to loamy soils were assumed, the was near equal to or exceeded the average amount of RO. This is best illustrated at Renfrew (a region where silty clay soils are predominant) where DRAINMOD estimated DD and RO were  mm and  mm, respectively (Table 2). Results for the clay soil at Renfrew produced DD and RO estimates of  mm and  mm, respectively, and, for the sandy soil, DD and RO were  mm and  mm. At Harrow (a region where clay loam soils (Brookston) are predominant), DD and RO were  mm and  mm, respectively. At Smithfield (a region where sandy soils are predominant), DD and RO were  mm and  mm, respectively, for the Pontypool sand. The total water surplus was over twice as much at Smithfield than at Harrow because the annual and ET were  mm and  mm at Smithfield and  mm and  mm at Harrow (Table 1).

Table 2: Annual average ± standard deviation of DRAINMOD estimated deep drainage (DD) and runoff (RO) for a local soil in each region where site is located and for a clay (Oneida clay), loam (Otonabee loam), and sandy (Pontypool sand) soil at 12 sites in Ontario.

Site
Local soil
Clay soil
Loam soil
Sandy soil

DD
RO
DD
RO
DD
RO
DD
RO

mm

Emo
127 ± 98
46 ± 44
40 ± 30
139 ± 82
104 ± 83
63 ± 52
171 ± 121
27 ± 38

Guelph
275 ± 106
25 ± 39
81 ± 36
218 ± 80
200 ± 89
89 ± 53
303 ± 116
23 ± 35

Harrow
31 ± 14
114 ± 59
10 ± 14
68 ± 52
20 ± 34
49 ± 48
119 ± 72
34 ± 39

Kapuskasing
154 ± 52
244 ± 77
100 ± 37
296 ± 88
246 ± 98
128 ± 67
404 ± 120
26 ± 41

London
107 ± 68
150 ± 77
71 ± 43
198 ± 95
185 ± 106
73 ± 68
302 ± 134
19 ± 38

Mount Forest
322 ± 99
89 ± 69
103 ± 39
297 ± 105
266 ± 94
123 ± 77
421 ± 137
12 ± 22

Ottawa
90 ± 40
211 ± 71
64 ± 38
233 ± 80
150 ± 91
137 ± 74
269 ± 115
54 ± 59

Peterborough
101 ± 86
90 ± 70
43 ± 34
153 ± 82
101 ± 86
90 ± 70
176 ± 107
52 ± 53

Renfrew
18 ± 15
214 ± 62
35 ± 28
186 ± 63
113 ± 69
78 ± 66
179 ± 100
68 ± 69

Smithfield
318 ± 111
50 ± 55
75 ± 38
262 ± 74
217 ± 84
109 ± 60
318 ± 111
50 ± 55

Vineland
32 ± 26
101 ± 58
32 ± 26
101 ± 58
65 ± 67
54 ± 49
154 ± 101
32 ± 36

Wiarton
130 ± 47
162 ± 78
73 ± 27
215 ± 86
186 ± 72
84 ± 63
354 ± 107
5 ± 12

3.2. Variability in Water Surplus

The year-to-year variability for each of the water balance components is shown in Figures 4 to  10 in [4] for each year at 7 sites and for another 5 sites in Figures 5 to  9 in [3] and tabulated for each year in Appendix B in both of these open access publications. In this paper the number of years in which DRAINMOD estimated DD and RO exceeded specified levels at each site will be summarized. The specified levels were somewhat arbitrary, but the risk of surface and groundwater contamination, from crop production systems, would be most likely at these levels. That is, the sites with a higher frequency and greater amount of DD are seen as having a higher risk of groundwater contamination in the absence of tile drainage (i.e., nutrients, pesticides, and pathogens) as discussed in [22]. Therefore, in addition to average  ± of the water balance components derived from DRAINMOD for the local soil and a sandy soil at each of the 12 sites, the number of years in which DD exceeded 300 and 400 mm and the number of years in which RO exceeded 150 and 200 mm are summarized in Table 3. Additionally, the number of years in which annual precipitation exceeded 1000 mm and the years in which it was less than 800 mm are also provided.

Table 3: Number of years that measured precipitation () and DRAINMOD estimated deep drainage (DD) and runoff (RO) exceeded the specified limits at 12 Ontario sites for the local and sandy (Pontypool sand) soil types for the 1954 to 2001 (48 years) period.

Site
Deep drainage
Runoff        
Precipitation

Local soil
Pontypool

sand
Local soil
Pontypool sand

300 mm
400 mm
300 mm
400 mm
150 mm
200 mm
150 mm
200 mm
<800 mm
>1000 mm

Emo
2
1
4
2
0
0
0
0
37
1

Guelph
15
8
17
11
1
0
1
0
14
8

Harrow
0
0
1
0
11
4
1
0
14
11

Kapuskasing
0
0
37
25
44
35
0
0
17
5

London
0
0
21
11
23
14
2
0
5
16

Mount Forest
25
8
38
23
5
3
0
0
4
19

Ottawa
0
0
19
6
39
24
2
2
15
7

Peterborough
0
0
8
2
10
4
3
1
29
3

Renfrew
0
0
5
1
41
29
6
2
21
2

Smithfield
26
13
26
13
0
0
0
0
4
14

Vineland
0
0
2
1
7
3
1
0
21
6

Wiarton
0
0
33
12
25
15
0
0
2
23

Estimated DD was >300 mm in numerous years for the local soil profile and frequently exceeded 400 mm at Guelph, Mount Forest, and Smithfield. However, when the sandy soil profile was assumed in DRAINMOD at all 12 sites, DD exceeded these limits in numerous years at 7 of the 12 sites. At four of these seven sites there were more years when exceeded 1000 mm and fewer years when it was less than 800 mm. The relation between DD and is illustrated in Figure 2, demonstrating the effect of texture in the soil profile on DD and how it increases with . The values indicate that the positive linear relationship between DD and precipitation weakens as the soil texture becomes finer. As well, the value for the slope of the regression line is only significant for sand (<0.05). This trend is probably related to the greater ability of clay soils to retain water than loam or sandy soils and the greater chance of runoff occurring on the clay soils. Hence, when a event occurs, drainage is more likely to happen in response to it as soil texture becomes coarser and hydraulic conductivity increases.

Figure 2: Comparison between average annual deep drainage (mm) and precipitation (mm) for the clay, loam, and sandy soil profiles used in this study at the 12 sites. indicates significance level of slope <0.05.

Analyzing the model estimates of large annual RO elucidates the high risk of surface water contamination at each site. Estimated RO exceeded 150 mm in more than 22 of the 48 years at five of the 12 sites for the local soil profile and for 14 or more years exceeded 200 mm. Mean annual DD exceeded RO at sites where the soil profile was loamy or sandy and RO exceeded DD where the soil in the region was predominately clay or silty clay. Both climatic and soil conditions had considerable influence in determining the surplus water and the relative amounts of DD and RO. More DD than RO occurred in the summer and fall seasons, perhaps due to the relatively flat surface and daily time step used in the model runs as given in [3, 4]. For the sandy soil profile there were very few years that RO exceeded 150 mm, so the risk of surface water contamination is less likely in regions where sandy soils occur, but there is a greater likelihood of groundwater contamination.

3.3. Timing of Water Surplus

As expected, the most water surplus occurred in late winter and spring seasons. For example, at Mount Forest the estimated mean water surplus in March, April, and May was 103, 169, and 46 mm, respectively, whereas it did not exceed 30 mm in any other month as given in [4]. The water surplus (DD + RO) in each of the four seasons are provided in Appendix E for 7 of the 12 sites in [4] and for five sites in [3].

3.4. Spatial Differences in Water Surplus in Ontario

A map of average annual water surplus was created for southern Ontario and part of northern Ontario (Figure 3), using the local soil type data and previously published maps (Figure 42 in [1]; Figure 39 in [2]; and Figures 9 and  10 in [23]). Their maps were derived using the Thornthwaite approach [24], and they assumed a perennial forage crop growing season for annual water balance and the 1931 to 1960 average monthly temperature and precipitation data for over 160 climate stations.

Figure 3: Isolines of average annual water surplus (mm), with ± standard deviation, (mm) of water surplus shown for Mount Forest only ( for all regions given in Table 2), in southern and part of northern (inset) Ontario.

The average annual water surplus (Figure 3) shows considerable differences across Ontario as discussed here using the climatic regions shown in Figure 1. It was the least (~150 mm) in the Leamington region; in the Niagara Fruit Belt; in an area south of Barrie in the Simcoe-Kawartha Lakes and South Slopes regions; and in the Rainy River-Thunder Bay region of northwestern Ontario. Water surplus was the greatest (≥350 mm) in the Dundalk Uplands, east of Georgian Bay in the Muskoka region, and north and east of Lake Superior in the Height of Land region of northern Ontario. Mount Forest near the centre of the Dundalk Uplands region in southern Ontario averaged  mm and Kapuskasing in the Northern Clay Belt region averaged  mm over the 48-year period as given in [25] and in Appendix in [3]. In the eastern part of southern Ontario water surplus increased gradually from ≤200 mm in the Renfrew region toward the Algonquin Park region, southeastern Haliburton Slopes region, and eastern part of the Eastern Counties region, where it exceeded 350 mm. The variability as depicted by the generally increased as the water surplus increased.

4. Summary and Conclusions

The hydrologic model (DRAINMOD) that simulates soil water flow, including plant uptake, ET, and freeze/thaw conditions, was used to estimate water surplus. Simulations were performed using daily climatic data from 1954 to 2001 for 12 sites in Ontario, each one in a different climatic region. A local soil profile for each region was evaluated along with similar clay, loam, and sandy soil types including the hydraulic properties for each soil profile. The seasonal growth characteristics for a corn crop in each region were also used as input for the model. Although the corn crop was assumed in this study, the average annual water surplus for all 12 sites matched the average water surplus, estimated using the Thornthwaite method, that uses monthly climate values and assumed a perennial crop growing season, in previous studies [1, 2, 23].

The major findings of this study were as follows.

(1)There was considerable annual variability in water surplus at all sites. The varied from 14 to 111 mm for DD and from 39 to 78 mm for RO for the local soil.(2)Mean annual DD exceeded RO at sites where the local soil was loamy or sandy and RO exceeded DD where the soils were predominately clay or silty clay.(3)Both precipitation and soil texture had a strong influence on the relative amounts of DD and RO.(4)Estimated DD was >300 mm in over half of the years at Smithfield and Mount Forest for the local soil profile and less frequently was >400 mm at certain sites (i.e., Guelph, Mount Forest, and Smithfield).(5)When a sandy soil was assumed for all 12 sites, DD was >300 and 400 mm in numerous years at 7 of the 12 sites. At four of the 7 sites there were more years when exceeded 1000 mm and fewer years when it was less than 800 mm.(6)Estimated RO was >150 mm in >22 of the 48 years at 5 of the 12 sites for the local soil profile and for 14 years exceeded 200 mm for all sites. However, for the sandy soil profile, there were very few years that RO exceeded these limits, so the risk of surface water contamination is less likely in regions where sandy soils occur and a greater likelihood of groundwater contamination.(7)Average annual water surplus ranged from 150 mm in the extreme southwest and in the Rainy River district of northwestern Ontario to over 400 mm at Mount Forest, Parry Sound/Muskoka district, and east of Lake Superior and was nearly 400 mm in the Kapuskasing/Timmins area.

The overall conclusion is that surplus water varies substantially between regions and from year to year in Ontario, which must be considered when developing policy and guidelines for source water protection.

Acknowledgments

The financial support provided by the Ontario Ministry of Agriculture and Food and Rural Affairs, Pioneer Hi-Bred Limited, Syngenta Crop Protection, and the Canadian Water Network is gratefully acknowledged. The authors’ thanks are extended to Marie Puddister, Geography Department, University of Guelph for drafting the maps.

References

D. M. Brown, G. A. McKay, and L. J. Chapman, “The climate of Southern Ontario,” Climatological Studies 5, Department of Transport, Meteorological Branch, Toronto, Canada, 1968.

L. J. Chapman and M. K. Thomas, “The climate of Northern Ontario,” Climatological Studies 6, Department of Transport, Meteorological Branch, Toronto, Canada, 1968.

H. Dadfar, D. J. Fallow, D. M. Brown, G. W. Parkin, J. D. Lauzon, and R. J. Gordon, “Estimation of annual and seasonal water surplus for five new regions in Ontario to identify critical regions for water quality monitoring,” Technical Memo 2010-1, University of Guelph, School of Environmental Sciences, 2010, http://hdl.handle.net/10214/2231.

D. J. Fallow, D. M. Brown, G. W. Parkin, J. D. Lauzon, and C. Wagner-Riddle, “Identification of critical regions for water quality monitoring with respect to seasonal and annual water surplus,” Technical Memo 2003-1, University of Guelph, Department of Land Resource Science, 2003, http://hdl.handle.net/10214/2204.

G. W. Parkin, C. Wagner-Riddle, D. J. Fallow, and D. M. Brown, “Estimated seasonal and annual water surplus in Ontario,” Canadian Water Resources Journal, vol. 24, no. 4, pp. 277–292, 1999.

R. De Jong and A. Bootsma, “Estimated long-term soil moisture variability on the Canadian prairies,” Canadian Journal of Soil Science, vol. 68, no. 2, pp. 307–321, 1988.

S. E. Allaire, E. van Bochove, J.-T. Denault et al., “Preferential pathways of phosphorus movement from agricultural land to water bodies in the Canadian Great Lakes basin: a predictive tool,” Canadian Journal of Soil Science, vol. 91, no. 3, pp. 361–374, 2011.

S. Dayyani, C. A. Madramootoo, P. Enright et al., “Field evaluation of DRAINMOD 5.1 under a cold climate: simulation of daily midspan water table depths and drain outflows,” Journal of the American Water Resources Association, vol. 45, no. 3, pp. 779–792, 2009.

C.-C. Yang, S. O. Prasher, S. Wang et al., “Simulation of nitrate-N movement in southern Ontario, Canada with DRAINMOD-N,” Agricultural Water Management, vol. 87, no. 3, pp. 299–306, 2007.

W. Luo, R. W. Skaggs, A. Madani, S. Cizikci, and A. Mavi, “Predicting field hydrology in cold condition with DRAINMOD,” Transactions of the American Society of Agricultural Engineers, vol. 44, no. 4, pp. 825–834, 2001.

R. W. Skaggs, “A water management model for shallow water table soils,” Tech. Rep. 134, Water Resources Research Institute of the University of North Carolina, North Carolina State University, Raleigh, NC, USA, 1978.

Soil and Water Management Group, DRAINMOD. North Carolina State University, Department of Bio and Ag Engineering, 2011, http://www.bae.ncsu.edu/soil_water/index.html

R. G. Allen, L. S. Pereira, D. Raes, and M. Smith, “Crop evapotranspiration: guidelines for computing crop requirements,” FAO Irrigation and Drainage Paper 56., FAO, Rome, Italy, 1998.

W. Luo, R. W. Skaggs, and G. M. Chescheir, “DRAINMOD modifications for cold conditions,” Transactions of the American Society of Agricultural Engineers, vol. 43, no. 6, pp. 1569–1582, 2000.

S. Ale, L. C. Bowling, J. R. Frankenberger, S. M. Brouder, and E. J. Kladivko, “Climate variability and drain spacing influence on drainage water management system operation,” Vadose Zone Journal, vol. 9, no. 1, pp. 43–52, 2010.

M. A. Brevé, R. W. Skaggs, J. E. Parsons, and J. W. Gilliam, “DRAINMOD-N, a nitrogen model for artificially drained soils,” Transactions of the American Society of Agricultural Engineers, vol. 40, no. 4, pp. 1067–1075, 1997.

M. A. Brevé, R. W. Skaggs, J. W. Gilliam et al., “Field testing of DRAINMOD-N,” Transactions of the American Society of Agricultural Engineers, vol. 40, no. 4, pp. 1077–1085, 1997.

G. P. Gupta, S. O. Prasher, S. T. Chieng, and I. N. Mathur, “Application of DRAINMOD under semi-arid conditions,” Agricultural Water Management, vol. 24, no. 1, pp. 63–80, 1993.

M. B. Shukla, S. O. Prasher, A. Madani, and G. P. Gupta, “Field validation of DRAINMOD in Atlantic Canada,” Canadian Agricultural Engineering, vol. 36, no. 4, pp. 205–213, 1994.

S. Wang, S. O. Prasher, R. M. Patel et al., “Fate and transport of nitrogen compounds in a cold region soil using DRAINMOD,” Computers and Electronics in Agriculture, vol. 53, no. 2, pp. 113–121, 2006.

D. M. Brown and A. Bootsma, “Crop heat units for corn and other warm-season Crops in Ontario,” Ministry of Agriculture and Food Factsheet Agdex 111/31, 1993.

J. L. Hutson, “Applying one-dimensional deterministic chemical fate models on a regional scale,” Geoderma, vol. 60, no. 1–4, pp. 201–212, 1993.

M. E. Sanderson and D. W. Philips, “Average annual water surplus in Canada,” Climatological Studies 9, Department of Transport, Meteorological Branch, Ontario, Canada, 1967.

C. W. Thornthwaite, “An approach toward a rational classification of climate,” Geographical Review, vol. 38, pp. 55–94, 1948.

D. M. Brown, “Average Annual Water Surplus in Ontario,” 2010, https://atrium.lib.uoguelph.ca/xmlui/handle/10214/2230