Effects of Aggregate Extraction on Water Storage in the Pepin Creek Watershed , British Columbia , Canada

Aggregate extraction in the Lower Fraser Valley of British Columbia, Canada, has affected the soil water storage of the Pepin Creek watershed. Although local government has set regulations for aggregate extraction projects to avoid negative environmental impacts, the gradual loss of soil materials and associated changes in vegetative cover has led to an alteration of the water balance within the watershed, which may affect surface or groundwater levels, and aquatic habitats. The study assessed the effects of aggregate extraction on the water storage of the Pepin Creek watershed and estimated that 25% of the surface area of the Canadian portion of the watershed has been affected by aggregate mining with an estimated loss of water storage of 10%. Evapotranspiration has decreased as a result of the removal of the vegetative cover. Precipitation has remained relatively constant over the study period but the annual discharge measured at Pepin Creek has decreased. Recommendations for enhancing environmental monitoring to better measure and understand ecological functions of the watershed during aggregate extraction are provided.


Introduction
Aggregate materials (e.g., rock, gravel, sand and stone) are important for construction, transportation and urban structure maintenance activities.It has been estimated that each Canadian uses, indirectly, 10 to 16 tonnes of aggregate per year [1].As the demand for aggregates continues to increase, particularly by the construction industry, aggregate mining in the Fraser Valley and Metro Vancouver Regional Districts is also increasing [2].Thus, there is a concern that the process of mining aggregate can have a potentially significant impact on the surrounding environment [3].The Fraser Valley Regional District FVRD immediately identifies the major concerns as: noise, dust, blasting and vibration, traffic and road safety, agricultural and public health, sediment erosion, drainage, and impacts on water supplies [2].Aggregate extraction may disrupt the water balance and impact the hydrological cycle of a watershed.In addition, the watershed may decrease its capacity of assimilating contaminants as a result of the changes of the stream regime, such as flow and temperature [4].From the perspective of water resources, aggregate extraction may also impact water quality on and nearby extraction sites, such as residential wells, as deleterious substances may enter the watershed through groundwater and surface run off [5].These in turn, may negatively affect wetlands and on-site lakes, with reduced water levels and increased summer water temperatures, which may be detrimental to cool water fish.
However, aggregates are in high demand in Metro Vancouver, which suggests that aggregate mining will continue within the study area.The goal of this study was to provide an assessment of the effects of aggregate extraction on the calculated water balance within a watershed with a focus on the effects on water storage.The Pepin Creek watershed, which contains several active mine sites for sand and gravel extraction, was selected as the case study to assess the effects of aggregate extraction on water balance dynamics.
The objectives of the study were to assess the effects of aggregate mining on the water balance, notably the water storage of Pepin Creek watershed, to identify potential environmental impacts of aggregate removal, and to suggest mitigation measures to decrease potential negative ecological effects.

Study Area
Pepin Creek watershed, located in the Lower Fraser Valley, was selected as a case study as it has a well-documented history of aggregate extraction (Figure 1).Pepin Creek is located within the Bertrand Creek watershed, south of the City of Abbotsford, British Columbia, Canada and flows into Whatcom County, Washington, USA (49˚00'13"N, 122˚28'18"W).The Canadian portion of the Pepin Creek watershed is approximately 16.5 square km in area, and 7.5 km in length at its longest point [6].
Aggregate extraction within the district started in the early 1960s, and the scale of extraction has increased relatively rapidly since that time [2].An example of the area affected by aggregate mining is shown in Figure 2. The surficial deposits are described as Sumas Drift [7].The dominant materials are mediumtextured aeolian materials overlying several meters of glaciofluvial sands and gravels [8].The soils are well drained and possess good permeability.Hydrometric data were obtained from the hydrometric station located at the Water storage (S) was calculated by Equation (1) [21]: where: S = the water storage of material removed (m 3 ), V = the volume of materials removed (m 3 ), ρ = the bulk density of materials removed (kg/m 3 ), WHC = the water holding capacity of materials removed (kg/kg), and 1000 mass (kg) of water to volume (m 3 ).
The volume of water lost as a result of the removal of the soil material, or storage (S), was estimated from the water balance Equation (2) presented in [21]: where: Q = the water discharge monitored for Pepin Creek (m 3 ), INPUT = the water input into the watershed (precipitation) (m 3 ), ΔS = the change in water storage within the watershed, ET = the evapotranspiration from vegetation (m 3 ).
Equation (2) assumes there is no surface runoff, nor water entering from adjacent sources.
To calculate actual evaporation, Equation (3) from [22] was used: where: EA = actual evaporation, and ET = the potential evaporation from a grass surface.
The calculation is based on the assumption that the soil is saturated, or not below field capacity (FC) and is equal to soil moisture (SM) at the beginning of the growing season.C 2 and C 4 are soil moisture loss values that change if SM becomes less than FC and are constant at a site.Kristensen and Jensen [22]

Water Balance
The study focused on a 26-year time interval, from 1985 to 2011, which had the most complete data available.Data availability was limited by the frequency of hydrometric data collection on Pepin Creek.The discharge measured on Pepin Creek and total precipitation, from May to October for the time period, are plotted in Figure 3.The total precipitation and potential evapotranspiration is shown in Figure 4.The interval linear trend lines plotted in the figures are from regression calculation and presented to serve as general guides.The variability over the relatively short time frame does not lend itself to statistical trend analysis.Separating the results for potential evapotranspiration and total precipitation within the watershed by growing season shows that the largest change has been in evapotranspiration in the dry season over the 26-year period (Figure 5).
There are three possible causes of the reduction in the apparent discharge, which are discussed in the following section: 1) A change in climate; 2) The loss of surface vegetation and evapotranspiration due to land clearing for aggregate extraction; and 3) The loss of water storage due to aggregate extraction.

Change in Climate
A comparison of total volume of precipitation and discharge between May and October is shown in Figure 6 and Figure 7.There is large variation in both precipitation and discharge measured annually.The total precipitation data suggest a slight negative trend but overall remain relatively stable, while the trend of total discharge exhibits a more pronounced negative trend.Since the reduction of the total discharge appears to be of greater magnitude than total precipitation, the reduction of discharge over the 6 months cannot be attributed exclusively to climate.

Loss of Vegetation (Bare Soil)
As a result of aggregate mining, there has been a progressive increase in nonvegetated surface area, and therefore a potential decrease in potential plant ET and total ET in the watershed.This is partially offset by the increase in evaporation from the soil.The overall effect of vegetation removal can be analyzed using Equation (3) from [22].The calculation is based on the assumption that soil is initially saturated and soil moisture (SM) remains at field capacity (FC).C 2 and C 4 are soil moisture loss constants that change if SM becomes less than FC.If we assume that SM does not fall below FC then EA/ET = 1, which means the evapotranspiration of bare soil (EA), equals the evapotranspiration of soil with grass (ET).Therefore, under this assumption, which is conservative, the loss of vegetation should not cause a large change in discharge.In estimating evapotranspiration from land surfaces in British Columbia, Spittle house [23] suggested that: where: α = a dimensionless constant, and E eq = the equilibrium evaporation rate.
Values for α in the Lower Fraser Valley are 1.26 ± 0.1 for pasture and 1.27 ± 0.1 for bare soil, which suggests that in the region there is little difference in ET as a result of changes in the soil vegetative cover [23].However, in the dry season, the actual soil moisture in unirrigated soil could be below field capacity; that is, the EA will be slightly below the ET.If this is correct, then the extra water (ET -EA) could be retained in the soil and eventually contribute to discharge.

Loss of Water Storage
To ascertain if a vegetative cover is significant we calculated the effects of two hypothetical scenarios.The first was to assume all the area was planted with sweet corn, the dominant cash crop in the region, and the second scenario with grass, which is common practice used in restoration projects.Assigning the land use in Pepin Creek watershed to sweet corn or grass, the calculated water storage (S) was then estimated with Equation (5) [24]: If sweet corn is selected as an example vegetative cover, the calculation for evapotranspiration of sweet corn may be obtained by use of Equation ( 6): where: ET c = evapotranspiration of sweet corn (mm/day), ET 0 = reference evapotranspiration (grass) (mm), and K c = single crop coefficient.K c is an empirical crop coefficient that has been adapted to the Pepin Creek area [21].The growing stage and corresponding coefficient was obtained from [25] and [26] and is presented in Table 1.The growing season was assumed to be 120 days, followed by 14 days of harvest, and no growth (fallow) for the remainder of the season (50 days), when the land is assumed to be bare and thus equal to the reference evapotranspiration rate (ET 0 ).After mapping the active extraction sites in the study area by year, the bare soil area and area covered with vegetation was calculated using ArcGIS [13].Combining the area values with the evapotranspiration data from [26], total evapotranspiration of both sweet corn and grass were calculated and the results are presented in Table 2.
The water storage of Pepin Creek watershed was calculated by Equation ( 7): Figure 8 shows that the change in water storage of the study area is decreasing over time.Combining both the climate and hydrological data indicates that the water balance has been affected by the removal of the aggregates from the site, regardless of whether the cover is sweet corn or grass.In the year of 1997, there was extraordinarily high precipitation, which resulted in the unusual low storage value in Figure 8.However, when the 1997 data is exempted, the trend lines for both scenarios still have negative trends.
The data for 1997 necessitates further discussion.As noted in Table 3, that year reported an exceptionally high volume of precipitation (over 50% higher than the average for the eight years).If one ignores that single year, the overall trend remains slightly negative.
This observation helps to explain the observed trends between precipitation and discharge during August to October.In July, the dry season, the discharge is high and probably influenced by the base flow from the previous winter.As a result of the reduced annual change in water storage, in the remainder of the dry season (August and September), the site did not store as much water, thus the excess contributed to discharge.In conclusion, the major impact has been the Table 1.Growing stage and single crop coefficient of BC sweet corn (adapted from [25] [26]).

Quantifying the Change in Annual Water Storage
A series of maps identifying the areas under aggregate extraction were constructed using ArcGIS, with examples from 1984, 1999 and 2011 shown in Figure 9.The images confirm that the Pepin Creek watershed has experienced a marked increase in the area of extracted aggregates, which is summarized in Table 4.The gravel sites were divided into 6 districts.The depth of each pit was obtained from Abbotsford Council Reports for soil removal and the amount of material removed was estimated by extrapolation from the reported duration of construction outlined in these reports [14]- [20].Bulk density and water holding capacity information of the removed material was derived from [8] to calculate total water stored in Table 4.
Figure 10 presents the actual loss of water storage in the Pepin Creek watershed under the scenario that all vegetation is grass or sweet corn with no irrigation.The positive trends suggest that the aggregate extraction is progressing relatively rapidly.The compensation value of re-establishing a vegetative cover approximately equals the total actual water storage of the watershed and the water storage of material extracted; that is, the value represents the amount of water that could be held.
Interestingly, the loss of water storage under grass cover is relatively unchanged over time (Figure 10), suggesting that change in annual storage over Comparing the balance of the compensation value of vegetation (whether grass or corn) has a positive trend.There is a suggestion that sweet corn is a better crop for retaining water and mitigating the negative impact of aggregate extraction on the water balance.However, sweet corn is harvested annually and the water contained in the corn is removed from the water cycle from the study area.
That is, the positive trend of the red line would decrease and suggest less loss of water storage over time.
The water balance of the Pepin Creek Watershed is altered and the water storage (S) is gradually decreasing.Since water storage of the removed material is a significant part of that loss, this suggests that aggregate extraction has a marked negative effect on the water balance within the watershed.

Discussion
Pepin Creek is a transboundary watershed, with a relatively mild climate that experiences little snowfall during winter, and is dependent on rainfall as the dominant form of precipitation.The watershed is rural and effects of urbanization

Potential Impacts
The A significant remaining question is concerns of a potential negative impact of lower discharge to Pepin Creek and the potential negative effects on cold water fish in the Creek [28], as the increased water temperature can change the metabolic rates and biological activity of fish [29].In addition, as it is well accepted that since warmer waters contain less dissolved oxygen, most aquatic species, including some aquatic plants, experience increased respiration rates as water temperature increases, which denatures cellular enzymes leading to an irreversible reduction of biological activity of fish, if there is an oxygen deficiency [30].

Conclusions
There have been changes in the overall water balance in Pepin Creek during the time frame of this study.The analysis suggests that climate change is not the main influence on the discharge of the Pepin Creek system over the 26 years of analysis.Total precipitation, evapotranspiration and discharge, although they exhibit minor trends, are neither synchronous nor statistically significant.This may, in part, be the result of the large variation in data from year to year and the Current environmental assessments provided by extraction companies usually focus on avoiding negative impact on the environment over a short time frame, such as controlling deleterious materials entering water bodies and the disruption of groundwater systems [31].This study suggests that regulators should conduct assessments of the impact of aggregate extraction over the length of excavation time and identify the accumulated impact on the environment, including water quantity, water quality and habitat of wildlife on a regular schedule.The indicators being monitored by Climate Canada and Water Survey Canada [11] [12], including water discharge, minimum, mean and maximum precipitation (rainfall and snowfall) and temperature, need to be continued and evaluated on a more frequent basis in an attempt to identify trends.To assess more completely the environmental and ecological effects of aggregate mining within the watershed, it is recommended that data on water temperature, water quality and dissolved oxygen, as well as groundwater dynamics be initiated.

Figure 1 .
Figure 1.Boundary of Pepin Creek watershed (within Canada) and location of climate and hydrometric stations.

Figure 2 .
Figure 2. Excavated soil and surficial materials at an aggregate extraction site in Pepin Creek watershed (Image source: L. M. Lavkulich).

Figure 3 .
Figure 3.Total precipitation in Pepin Creek watershed and discharge rate of Pepin Creek during the growing season from 1985 to 2011.(Note: Linear trend line shows general direction of trend; not statistically significant).

Figure 4 .Figure 3
Figure 4. Potential evapotranspiration and total precipitation for the Pepin Creek watershed.(Note: Linear trend line shows general direction of trend; not statistically significant).

Figure 5 .
Figure 5. Potential evapotranspiration and total precipitation by season (wet and dry) and the complete year for Pepin Creek from 1985 to 2011.(Note: Linear trend line shows general direction of trend; not statistically significant).

Figure 6 .
Figure 6.Volume of total annual precipitation in m 3 (May-Oct.)calculated for the Pepin Creek watershed.(Note: Linear trend line shows general direction of trend; not statistically significant).

Figure 7 .
Figure 7. Volume of total annual discharge in m 3 (May-Oct.)in Pepin Creek.(Note: Linear trend line shows general direction of trend; not statistically significant).

Figure 8 .
Figure 8. Loss of water storage (S) of Pepin Creek watershed calculated under grass or sweet corn vegetative cover.(Note: Linear trend line shows general direction of trend; not statistically significant).

Table 4 .Figure 10 .
Figure 10.Change in water storage (S) in Pepin Creek watershed, under different vegetative cover and bare soil scenarios.(Note: Linear trend line shows general direction of trend; not statistically significant).
Abbotsford[27], include detailed information about extraction activity, including types of soil, estimated depth of topsoil, estimated volume of soil removed and detailed groundwater surveys.Topsoil is stockpiled for replacement on excavated areas to be used for rehabilitation after active mining.Rehabilitation is required after the extraction activity ends.In the Pepin Creek watershed, of the 421 ha cumulative extracted area between 1984 and 2011, approximately 110 ha, or 26% of the mine site areas have undergone rehabilitation.This equates to approximately 6.7% of the Canadian portion of the watershed.There is currently no reliable information on the types or amounts of material used in rehabilitation of these aggregate extraction sites, save for the topsoil reserved and returned to the land.This amount of material is likely small in comparison to the total extracted volume of aggregate.Replanting surface vegetation and returning the site to viable agricultural production are the ultimate goals of the rehabilitation of these sites, not restoring the hydrological characteristics of the landscape.Rehabilitating these sites may result in a net increase in surface ET, but not in the water storage component.
relatively short time frame of analysis.However, the yearly trend patterns for precipitation and discharge are synchronous.The removal of aggregate has resulted in a loss of surface vegetative cover of about 20% and a loss of 10% water storage, and although the effects of climate change cannot be excluded; the loss of water storage as a result of aggregate extraction appears to be the major factor in affecting the water balance in the Pepin Creek watershed.Limitations of the analysis that require additional study include an analysis of groundwater information and its contribution to the hydrological dynamics of the watershed; the decreasing volumes of discharge and to what extent these dynamics might be influenced and impacted by groundwater from adjacent watersheds.The assumption that the surficial geographical boundaries reflect the groundwater aquifer Y. N. Wang et al.DOI: 10.4236/nr.2017.87030475 Natural Resources boundary needs further attention.
define C 2 as the evaporation taking place, regardless of vegetation density and soil 4 is the amount of soil water that should be removed from the upper soil layer to reduce the evaporation from bare soil to the equivalent of C 2 .If SM remains above FC, EA/ET equals 1, or evapotranspiration from bare soil, EA, equals the evapotranspiration of the soil with grass, ET.

Table 2 .
Total evapotranspiration of Pepin Creek watershed calculated under grass or sweet corn vegetative cover.

Table 3 .
Water storage change of Pepin Creek watershed calculated under grass or sweet corn vegetative cover.Figure9shows the extent of aggregate extraction areas in the Pepin Creek watershed since 1984.At the beginning of the study period (1984), there were only a few parcels under extraction, with a total area of 116 hectares.In 1999, the area of extracted area had more than doubled to 326 ha, and by 2011, had in- The depth of aggregate extraction is limited to within 15 -18 meters from the surface, depending on specific circumstances, in order to reduce negative influence on groundwater.This study shows that by 2011, about 25% of the area within Pepin Creek watershed has been under aggregate extraction, and the water storage has decreased by 10% -20%.The loss of the water storage will likely affect the contribution of the lost stored water to groundwater and surface water supplies in the watershed.A decrease in discharge during dry the season, may result in higher water temperatures, and lower volumes for other purposes in Pepin Creek.Although there is some evidence that climate change, notably temperature, may have an effect on discharge, its effects and that of changing vegetative cover have a lower impact on the water balance than the removal of the materials by aggregate mining.However, as there was no information available on groundwater dynamics in the study area and it is possible that the groundwater aquifer boundary is not congruous with the surficial boundary of the watershed, the effects and contribution of groundwater to the Pepin Creek discharge needs further evaluation.