<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">AS</journal-id><journal-title-group><journal-title>Agricultural Sciences</journal-title></journal-title-group><issn pub-type="epub">2156-8553</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/as.2017.88062</article-id><article-id pub-id-type="publisher-id">AS-78737</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Biomedical&amp;Life Sciences</subject><subject> Earth&amp;Environmental Sciences</subject></subj-group></article-categories><title-group><article-title>
 
 
  Analysis of Impervious Surface Area, and the Impacts on Soil-Based Agriculture and the Hydrologic Cycle: A Case Study in the Agricultural Land Reserve in Metro Vancouver, British Columbia, Canada
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Ashley</surname><given-names>Rose</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Julie</surname><given-names>E. Wilson</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Les</surname><given-names>M. Lavkulich</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib></contrib-group><aff id="aff1"><addr-line>Master of Land and Water Systems, Faculty of Land and Food Systems, The University of British Columbia, Vancouver, Canada</addr-line></aff><pub-date pub-type="epub"><day>08</day><month>08</month><year>2017</year></pub-date><volume>08</volume><issue>08</issue><fpage>837</fpage><lpage>856</lpage><history><date date-type="received"><day>July</day>	<month>14,</month>	<year>2017</year></date><date date-type="rev-recd"><day>Accepted:</day>	<month>August</month>	<year>25,</year>	</date><date date-type="accepted"><day>August</day>	<month>28,</month>	<year>2017</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  The province of British Columbia, Canada, has established an Agricultural Land Reserve (ALR) to protect the most suitable soil landscapes for agriculture. Increases in population and urbanization have resulted in development challenges on ALR lands. The Metro Vancouver area is the most productive agricultural area in British Columbia as well as the most rapidly growing urban region. The increase in impervious areas has decreased the amount of arable land for soil-based agriculture and altered the hydrological cycle. Analysis using a combination of aerial imagery and GIS found that impervious areas comprise about 10 percent of the ALR within Metro Vancouver. Farm residences and greenhouses have the largest effect on reducing the soil surface for water infiltration. This decrease in area has negatively influenced the ecosystem heath of the region, as well as, decreasing the amount of agricultural land for soil based agriculture and both surface and groundwater dynamics.
 
</p></abstract><kwd-group><kwd>Agricultural Land Reserve</kwd><kwd> Impervious Surfaces</kwd><kwd> Watershed Health</kwd><kwd>  Stormwater Runoff</kwd><kwd> Hydrology</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>Prime agricultural land is commonly located adjacent to urban centers, as historic human settlement and urbanization was associated with local land and water resources that could provide requisites of food and shelter. Urbanization poses an ongoing threat to farmland; every year an estimated 65,000 km<sup>2</sup> of cropland is lost worldwide to urban expansion for uses such as housing, industry and infrastructure [<xref ref-type="bibr" rid="scirp.78737-ref1">1</xref>] . Between 2001 and 2011 nearly one million hectares (ha) of dependable agricultural land was taken out of cultivation in Canada and subsumed, mostly by urban development [<xref ref-type="bibr" rid="scirp.78737-ref2">2</xref>] . In addition to the physical loss of agricultural land, large-scale agriculture has resulted in an increase in farm building infrastructure and impervious areas [<xref ref-type="bibr" rid="scirp.78737-ref3">3</xref>] . At a time when there is growing pressure for local food security, prime agricultural lands in the peri-urban environment are becoming “sealed” from food production by the impervious surfaces. Global estimations indicate 0.43% of the world’s land area is covered with impervious surfaces [<xref ref-type="bibr" rid="scirp.78737-ref4">4</xref>] .</p><p>Arable agricultural land is a scarce resource in Canada, as only 5% of it is suitable for crop production and free from limitations such as climate, topography, water availability and soil quality [<xref ref-type="bibr" rid="scirp.78737-ref5">5</xref>] . According to the 1972 Canada Land Inventory (CLI), which classifies soils based on their capability for agriculture from Class 1 to Class 7, arable agricultural land is defined as Class 1 - 3 [<xref ref-type="bibr" rid="scirp.78737-ref6">6</xref>] . Class1 has the highest capability agricultural soil with no significant limitations for a range of climatically suitable crops, while Class 7 is the lowest and non- arable [<xref ref-type="bibr" rid="scirp.78737-ref7">7</xref>] . In British Columbia (BC) soils are assessed for capability and limitations through the “Land Capability Classification for Agriculture in British Columbia” and Class 1 - 4 is considered suitable for a range of crops [<xref ref-type="bibr" rid="scirp.78737-ref8">8</xref>] .</p><p>In the late 1960s and early 1970s it was recognized that BC was losing an estimated 6000 ha of prime farmland to urban development and other uses each year [<xref ref-type="bibr" rid="scirp.78737-ref9">9</xref>] . In order to preserve this valuable farmland the provincial government introduced BC’s Agricultural Land Commission Act (the ALC Act) on April 18, 1973. The ALC Act identified 4.7 million ha or 5% of BC’s land base for agriculture based upon the CLI and placed these lands into an Agricultural Land Reserve (ALR) [<xref ref-type="bibr" rid="scirp.78737-ref9">9</xref>] . Of the province’s ~945,000 km<sup>2</sup> land base, approximately 1.1% is considered prime agricultural land (CLI Class 1 - 3) and is primarily located in the Fraser River Valley near Metro Vancouver [<xref ref-type="bibr" rid="scirp.78737-ref10">10</xref>] .</p><p>Metro Vancouver’s ALR contains 1.5% of the agricultural land in BC and is some of the best agricultural land in Canada with its combination of soil, climate and available water [<xref ref-type="bibr" rid="scirp.78737-ref11">11</xref>] . However, there is increasing concern about covering fertile soils with impervious surfaces such as greenhouses, farm structures, poultry and dairy barns, and residential developments in a process also known as “soil sealing”. Permanently covering fertile soils with impervious surfaces diminishes the capacity of the region to produce food and feed the growing population, forecasted to grow by 30% from 2001 to 2025 [<xref ref-type="bibr" rid="scirp.78737-ref12">12</xref>] . Changes in impervious surface cover may also alter the hydrology, stream morphology, water and sediment quality, stream temperature and aquatic habitat, and biodiversity of the region [<xref ref-type="bibr" rid="scirp.78737-ref13">13</xref>] .</p><p>The purpose of this study is to investigate the extent of “soil sealing” on prime agricultural land located in Metro Vancouver’s ALR with impervious surfaces. The results are discussed in the context of two primary concerns: the loss of soil infiltration capacity and the loss of land available for food production. The permitted land uses in the ALR that contribute to imperviousness are discussed in relation to food production and farming opportunities and their impact on the hydrologic cycle. Recommendations to minimize the impact of impervious surfaces and to protect the long-term sustainability of soil-based agriculture are suggested.</p><sec id="s1_1"><title>1.1. Impervious Surfaces</title><p>Impervious area is a useful indicator of watershed conditions, as it relates to the health of aquatic ecosystems [<xref ref-type="bibr" rid="scirp.78737-ref13">13</xref>] . Increased impervious area has been linked with changes in stream biotic assemblages, including increases in algae biomass and decreases in benthic macroinvertebrate and fish abundance [<xref ref-type="bibr" rid="scirp.78737-ref14">14</xref>] [<xref ref-type="bibr" rid="scirp.78737-ref15">15</xref>] . Numerous studies have shown that stream quality starts to degrade when a watershed is more than 10% impervious [<xref ref-type="bibr" rid="scirp.78737-ref16">16</xref>] . For example, Stepenuck et al. [<xref ref-type="bibr" rid="scirp.78737-ref17">17</xref>] observed a threshold of 8 to 12 percent imperviousness across 43 southeastern Wisconsin streams, where stream quality is degraded and macroinvertebrate richness declines. In BC, Coho salmon are rarely found in watersheds exceeding 10% - 15% imperviousness, as their headwaters and spawning grounds become affected [<xref ref-type="bibr" rid="scirp.78737-ref13">13</xref>] . However damage thresholds may vary across different ecosystems and the ten-percent rule is a generalization [<xref ref-type="bibr" rid="scirp.78737-ref18">18</xref>] .</p><p>Impervious surfaces impact the quality and quantity of runoff discharged into receiving streams, often directly through sewer systems [<xref ref-type="bibr" rid="scirp.78737-ref19">19</xref>] . This disrupts habitat, water flow, food systems and biotic and chemical interactions [<xref ref-type="bibr" rid="scirp.78737-ref3">3</xref>] [<xref ref-type="bibr" rid="scirp.78737-ref13">13</xref>] . Stormwater carries with it pollutants that accumulate on hard surfaces, including pathogens, nutrients, pesticides, oil and grease, metals, and sediments that degrade water quality [<xref ref-type="bibr" rid="scirp.78737-ref14">14</xref>] [<xref ref-type="bibr" rid="scirp.78737-ref20">20</xref>] . This can lead to what is known as the “first flush phenomenon”, where pollutants that accumulate during dry periods, are transported after a major rainfall event into the drainage system, typically containing higher concentrations of toxins such as heavy metals, polychlorinated biphenyls (PCBs) and polycyclic aromatic hydrocarbons [<xref ref-type="bibr" rid="scirp.78737-ref20">20</xref>] [<xref ref-type="bibr" rid="scirp.78737-ref21">21</xref>] .</p><p>Impervious surfaces reduce the capacity of soils to store water and peak surface runoff rates increase with increasing precipitation [<xref ref-type="bibr" rid="scirp.78737-ref22">22</xref>] . This intensifies channel and bank erosion and increases vulnerability to flooding and sediment transport, which may decrease water clarity, benthic diversity, and spawning survival, as well as clogging the gills of fish [<xref ref-type="bibr" rid="scirp.78737-ref13">13</xref>] . Groundwater recharge and base flow is also reduced, as less precipitation infiltrates the soil and runs off hard surfaces [<xref ref-type="bibr" rid="scirp.78737-ref16">16</xref>] . Thus impervious surfaces indicate a disruption in the hydrological cycle and increased contamination through runoff.</p><p>Climate change projections for the Metro Vancouver region may exacerbate these effects, as they predict more frequent and extreme weather events [<xref ref-type="bibr" rid="scirp.78737-ref23">23</xref>] . The projections also suggest that peak stream flows will occur earlier in the spring, as more precipitation falls as rain and less as snow in the winter [<xref ref-type="bibr" rid="scirp.78737-ref23">23</xref>] . Flows would also decrease in the summer, as snowpack is reduced in the winter [<xref ref-type="bibr" rid="scirp.78737-ref23">23</xref>] . Increased impervious surface area compounds these effects, as stormwater drainage systems collect and convey runoff, giving it less of an opportunity to infiltrate into the soil and contribute to groundwater recharge and base flow of streams [<xref ref-type="bibr" rid="scirp.78737-ref22">22</xref>] .</p><sec id="s1_1_1"><title>1.1.1. Total versus Effective Impervious Area</title><p>Two categories are commonly used to quantify impervious cover: 1) Total Impervious Area (TIA) and 2) Effective Impervious Area (EIA). TIA is the measurement of all the impervious surfaces in a catchment area that do not allow water to infiltrate the soil including streets, sidewalks, driveways, and rooftops [<xref ref-type="bibr" rid="scirp.78737-ref14">14</xref>] . From a hydrological perspective TIA is limited because it ignores areas that may be vegetated but are too highly compacted to allow water to infiltrate [<xref ref-type="bibr" rid="scirp.78737-ref24">24</xref>] . Secondly, it includes paved areas that may not be directly connected to a waterbody via the storm drain network. For example, if runoff is diverted from a rooftop onto a green space such as a rain garden or a grassed swale then it is not contributing to the surface runoff of the residential area, and is helping to reduce peak runoff downstream. EIA measures the impervious surfaces that are directly connected to downstream drainage and is the portion of runoff that will impact both hydrology and water quality [<xref ref-type="bibr" rid="scirp.78737-ref24">24</xref>] [<xref ref-type="bibr" rid="scirp.78737-ref25">25</xref>] .</p></sec><sec id="s1_1_2"><title>1.1.2. Watershed Health and Agricultural Land Use</title><p>A watershed is considered healthy when it can functionally support terrestrial, riparian, aquatic and wetland habitats and sustain diverse native populations of aquatic and riparian species. The Watershed Health Tracking System (WHTS) (2005) is a proposed means of measuring current and predicted watershed health in Metro Vancouver, wherein total impervious area (TIA) and riparian integrity are used as key indicators of watershed health [<xref ref-type="bibr" rid="scirp.78737-ref26">26</xref>] . TIA is ranked relative to stream health as follows: 1) stressed 1% - 10%, 2) impacted 11% - 25% and 3) degraded 26% - 100% [<xref ref-type="bibr" rid="scirp.78737-ref26">26</xref>] . Riparian integrity is measured using a 30- meter (m) setback or buffer zone on both sides of the watercourse over its entire length. If the 30 m riparian area is in pristine condition it is ranked 100% and if half of it is reduced to 15 m, the riparian integrity is closer to 25% [<xref ref-type="bibr" rid="scirp.78737-ref26">26</xref>] .</p><p>The WHTS does not factor in agricultural land use, which may compound the environmental impacts of impervious areas. For example, clearing of forested lands in the Metro Vancouver region for agriculture has increased peak flows from less than 5 to as much as 50% [<xref ref-type="bibr" rid="scirp.78737-ref27">27</xref>] . Agriculture also extracts large amounts of surface and groundwater, particularly in the hotter and drier months of July and August, which can account for 60% of irrigation water and result in minimum low flows [<xref ref-type="bibr" rid="scirp.78737-ref27">27</xref>] . Water quality may be impacted by nutrient inputs from fertilizers and manure, agricultural chemicals and increases in sediment loading from erosion. Wang et al. [<xref ref-type="bibr" rid="scirp.78737-ref19">19</xref>] measured watershed health using the fish community as an indicator and defined a threshold of 50% agricultural land use when agricultural watersheds move from “excellent” to “good” health [<xref ref-type="bibr" rid="scirp.78737-ref27">27</xref>] .</p><p>Integrating agricultural land use into the WHTS was not deemed to be scientifically defensible, based on the assertion that the effects of agriculture relative to urban land use on watershed health have not been adequately examined or studied [<xref ref-type="bibr" rid="scirp.78737-ref27">27</xref>] . This is in part because the term agriculture encompasses a wide range of activities, which have varying degrees of impact that make establishing an empirical relationship between the extent of agricultural land use and watershed health difficult [<xref ref-type="bibr" rid="scirp.78737-ref27">27</xref>] . It is also challenging to define the extent of stream networks in Metro Vancouver’s agricultural areas, as there is a large number of drainage ditches. Measuring the riparian areas adjacent to these ditches would also be somewhat unpractical, as it is a time consuming process [<xref ref-type="bibr" rid="scirp.78737-ref27">27</xref>] .</p><p>The objectives of this study are to estimate the extent of impervious surfaces within the Metro Vancouver ALR and to provide an analysis of the potential effects on the hydrological cycle and loss of land for soil-based agriculture. Residential developments and greenhouses will be discussed in depth due to the extent of arable agricultural soil they cover in the ALR study area. Two methods are used to calculate the extent of TIA: first is a preliminary assessment conducted in the absence of detailed land cover data using satellite imagery, and the second is using detailed land cover data from Agricultural Land Use Inventories (ALUIs) prepared by the BC Ministry of Agriculture [<xref ref-type="bibr" rid="scirp.78737-ref28">28</xref>] .</p></sec></sec></sec><sec id="s2"><title>2. Methods</title><sec id="s2_1"><title>2.1. Study Area: Metro Vancouver Agricultural Land Reserve</title><p>Metro Vancouver is located in the southwest corner of British Columbia in the Fraser River watershed (49˚04'59.7&quot;N, 122˚58'19.5&quot;W). It covers a total area of 278,700 ha and is comprised of 22 municipalities, one electoral area, and one First Nation [<xref ref-type="bibr" rid="scirp.78737-ref28">28</xref>] . The Agricultural Land Reserve within Metro Vancouver is 60,890 ha or 22% of the land base [<xref ref-type="bibr" rid="scirp.78737-ref28">28</xref>] . The study area for this project includes six municipalities: Langley (Township), Surrey, Maple Ridge, Richmond, Delta and Pitt Meadows (<xref ref-type="table" rid="table1">Table 1</xref>), which support 94% of Metro Vancouver’s 2821 farms [<xref ref-type="bibr" rid="scirp.78737-ref29">29</xref>] . The study focused on 90% of the land in Metro Vancouver’s ALR (approximately 54,570 ha) [<xref ref-type="bibr" rid="scirp.78737-ref28">28</xref>] . Ten percent of the ALR was excluded from the study area due to sparse distribution and site limitations that prevented it from being surveyed in the Metro Vancouver ALUI.</p><sec id="s2_1_1"><title>2.1.1. Physiography</title><p>The study area is located within the Fraser River Lowland physiographic region (<xref ref-type="fig" rid="fig1">Figure 1</xref>). The Fraser Lowland was created 10,000 years ago from glacial and alluvial deposits along the Fraser River [<xref ref-type="bibr" rid="scirp.78737-ref30">30</xref>] . The Fraser River drains into the Strait of Georgia, separating the municipalities of Delta, Surrey, and Langley on the south side, Richmond in the middle, and Pitt Meadows and Maple Ridge on the north side. The south side of the Fraser Lowlands (including Richmond) is comprised of wide valley bottoms with gently rolling or flat-topped uplands and elevations typically lower than 15 m above sea level [<xref ref-type="bibr" rid="scirp.78737-ref27">27</xref>] . The north side has a</p><table-wrap id="table1" ><label><xref ref-type="table" rid="table1">Table 1</xref></label><caption><title> Percentage and area of the agricultural land reserve in six municipalities</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Local Government</th><th align="center" valign="middle" >Geographical Coordinates</th><th align="center" valign="middle" >% in ALR</th><th align="center" valign="middle" >ALR Area (ha)</th></tr></thead><tr><td align="center" valign="middle" >Langley (Township)</td><td align="center" valign="middle" >49˚5'33.60&quot;N 122˚30'31.47&quot;W</td><td align="center" valign="middle" >37%</td><td align="center" valign="middle" >22,280</td></tr><tr><td align="center" valign="middle" >Delta</td><td align="center" valign="middle" >49˚4'47.61&quot;N 123˚4'2.10&quot;W</td><td align="center" valign="middle" >15%</td><td align="center" valign="middle" >8,8430</td></tr><tr><td align="center" valign="middle" >Surrey</td><td align="center" valign="middle" >49˚7'48.71&quot;N 122˚46'52.99&quot;W</td><td align="center" valign="middle" >14%</td><td align="center" valign="middle" >8670</td></tr><tr><td align="center" valign="middle" >Pitt Meadows</td><td align="center" valign="middle" >49˚14'25.69&quot;N 122˚38'58.40&quot;W</td><td align="center" valign="middle" >10%</td><td align="center" valign="middle" >6380</td></tr><tr><td align="center" valign="middle" >Richmond</td><td align="center" valign="middle" >49˚10'13.99&quot;N 123˚1'50.98&quot;W</td><td align="center" valign="middle" >8%</td><td align="center" valign="middle" >4760</td></tr><tr><td align="center" valign="middle" >Maple Ridge</td><td align="center" valign="middle" >49˚11'39.14&quot;N 122˚27'41.75&quot;W</td><td align="center" valign="middle" >6%</td><td align="center" valign="middle" >3630</td></tr><tr><td align="center" valign="middle" >Study Area Total</td><td align="center" valign="middle"  colspan="2"  >90%</td><td align="center" valign="middle" >54,570</td></tr><tr><td align="center" valign="middle" >Other</td><td align="center" valign="middle"  colspan="2"  >4%</td><td align="center" valign="middle" >2810</td></tr><tr><td align="center" valign="middle" >Not Surveyed (limitations)</td><td align="center" valign="middle"  colspan="2"  >6%</td><td align="center" valign="middle" >3510</td></tr><tr><td align="center" valign="middle" >ALR Total</td><td align="center" valign="middle"  colspan="2"  >100%</td><td align="center" valign="middle" >60,890</td></tr></tbody></table></table-wrap><fig id="fig1"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref></label><caption><title> Location of the Metro Vancouver agricultural land reserve study area</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/9-3001734x2.png"/></fig><p>similar topography but is up to 200 m elevation with upland surfaces ranging from flat-topped outwash or deltas to hummocky till and glaciomarine sediments [<xref ref-type="bibr" rid="scirp.78737-ref27">27</xref>] .</p><p>The parent material of the Fraser River Lowland ranges from deposits of alluvium near the mouth of the Fraser River to glacial till at higher elevations [<xref ref-type="bibr" rid="scirp.78737-ref27">27</xref>] . The ALR study area is largely composed of fine silts and clay, which have resulted in soils characterized by impeded drainage, but high agricultural productivity [<xref ref-type="bibr" rid="scirp.78737-ref31">31</xref>] .</p></sec><sec id="s2_1_2"><title>2.1.2. Climate</title><p>Metro Vancouver has a modified maritime climate, characterized by wet, cool winters (October to March) and a pronounced dry period in the summer (July to August). The temperature rarely drops below zero degrees Celsius and the region has one of the longest frost-free periods in Canada [<xref ref-type="bibr" rid="scirp.78737-ref27">27</xref>] .</p><p>The average rainfall varies and is divided into three distinct hydrological zones that increase in precipitation from the south to north towards the Cascade Mountains [<xref ref-type="bibr" rid="scirp.78737-ref32">32</xref>] :</p><p>1) Low rainfall zone (&lt;1350 mm annual): Lowlands of Richmond, Delta and south Surrey;</p><p>2) Medium rainfall zone (1350 mm to 1850 mm annual): Township of Langley and north Surrey; and</p><p>3) High rainfall zone (&gt;1850 mm annual): Maple Ridge and Pitt Meadows.</p></sec></sec><sec id="s2_2"><title>2.2. Analysis</title><sec id="s2_2_1"><title>2.2.1. Imagery Analysis of Impervious Surfaces</title><p>An initial assessment of impervious surface area was conducted in the study area using Google Earth Pro [<xref ref-type="bibr" rid="scirp.78737-ref33">33</xref>] , as a means of estimating Total Impervious Area (%) in the absence of detailed land cover data. The boundaries of the ALR were derived from Keyhole Markup Language (KML) files from the Provincial Agriculture Land Commission [<xref ref-type="bibr" rid="scirp.78737-ref34">34</xref>] .</p><p>Four 65 ha plots were analyzed in each municipality in the study area. To choose the plots, four place marks were randomly pinned within the ALR boundary at each location in Google Earth at a 1:20,000 scale. The boundaries of the plots were digitized following the gridlines of the agricultural fields visible through satellite imagery [<xref ref-type="bibr" rid="scirp.78737-ref33">33</xref>] . The land covers identified as impervious were digitized into Google Earth polygons for each plot including:</p><p>・ Built objects such as residential, commercial, industrial, recreational and institutional structures;</p><p>・ Human created bare sites such as aggregate extraction sites, compacted soil or gravel parking lots and storage areas; and</p><p>・ Paved areas including roads, driveways and parking lots.</p><p>Earth Point [<xref ref-type="bibr" rid="scirp.78737-ref35">35</xref>] was used to calculate the area of the impervious land cover in hectares for each plot. The impervious land covers for each individual plot were summed to determine the percentage of the plot that is impervious, or the TIA.</p></sec><sec id="s2_2_2"><title>2.2.2. Analysis of Agricultural Land Use Inventory Data</title><p>1) Land Cover</p><p>To determine the relationship between impervious surfaces and watershed health, watershed boundaries are usually delineated based on topography and the storm sewer drainage system [<xref ref-type="bibr" rid="scirp.78737-ref24">24</xref>] . However, due to the complex nature of the agricultural drainage systems including pumping and retention systems in Metro Vancouver’s ALR, specific watersheds were not delineated in this study. Imperviousness was calculated as a percentage of land cover within the ALR boundaries of the study area.</p><p>A detailed analysis of land cover data from the Agriculture Land Use Inventories (ALUIs) prepared by the BC Ministry of Agriculture was conducted to provide a more accurate assessment of TIA in the Metro Vancouver ALR study area, as compared to the imagery analysis. Data were collected for the ALUIs with AgFocusin the summer of 2010 and 2011, which employs a “windshield” survey method and Geographic Information Systems (GIS) to identify land use and land cover on legal parcels [<xref ref-type="bibr" rid="scirp.78737-ref28">28</xref>] .</p><p>2) Imperviousness Factors</p><p>Land cover types identified as impervious were assigned an imperviousness factors (IF) (<xref ref-type="table" rid="table2">Table 2</xref>), derived from the Sensitive Habitat Inventory Mapping</p><table-wrap id="table2" ><label><xref ref-type="table" rid="table2">Table 2</xref></label><caption><title> Impervious factors derived from the BC Sensitive Habitat Inventory Mapping Manual</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Impervious Land Covers</th><th align="center" valign="middle" >Impervious Factors (%)</th></tr></thead><tr><td align="center" valign="middle" >Cultivated Field Crops</td><td align="center" valign="middle" >3%</td></tr><tr><td align="center" valign="middle" >Farm Infrastructure</td><td align="center" valign="middle" >80%</td></tr><tr><td align="center" valign="middle" >Greenhouses</td><td align="center" valign="middle" >80%</td></tr><tr><td align="center" valign="middle" >Crop Barn</td><td align="center" valign="middle" >80%</td></tr><tr><td align="center" valign="middle" >Unused forage or Pasture</td><td align="center" valign="middle" >3%</td></tr><tr><td align="center" valign="middle" >Unmaintained field crops</td><td align="center" valign="middle" >3%</td></tr><tr><td align="center" valign="middle" >Unmaintained greenhouses</td><td align="center" valign="middle" >80%</td></tr><tr><td align="center" valign="middle" >Unmaintained crop barn</td><td align="center" valign="middle" >80%</td></tr><tr><td align="center" valign="middle" >Managed vegetation</td><td align="center" valign="middle" >3%</td></tr><tr><td align="center" valign="middle" >Golf fairway/green</td><td align="center" valign="middle" >3%</td></tr><tr><td align="center" valign="middle" >Settlement</td><td align="center" valign="middle" >12%</td></tr><tr><td align="center" valign="middle" >Non Built or Bare</td><td align="center" valign="middle" >3%</td></tr><tr><td align="center" valign="middle" >Built up-Other</td><td align="center" valign="middle" >12%</td></tr><tr><td align="center" valign="middle" >Transportation</td><td align="center" valign="middle" >90%</td></tr><tr><td align="center" valign="middle" >Waterbodies</td><td align="center" valign="middle" >0%</td></tr><tr><td align="center" valign="middle" >Utilities</td><td align="center" valign="middle" >80%</td></tr><tr><td align="center" valign="middle" >Vegetated</td><td align="center" valign="middle" >0%</td></tr><tr><td align="center" valign="middle" >Wetlands</td><td align="center" valign="middle" >0%</td></tr><tr><td align="center" valign="middle" >Waterbodies</td><td align="center" valign="middle" >0%</td></tr></tbody></table></table-wrap><p>(SHIM) manual, a standard for fish and aquatic mapping for watersheds in BC [<xref ref-type="bibr" rid="scirp.78737-ref24">24</xref>] . Total Impervious Area was calculated by multiplying the percentage of land cover in the ALR in each category by the corresponding impervious factor using a spreadsheet.</p><p>The IF for residential footprint in the study area was established independently, as imperviousness can vary depending on housing type and density, design of road networks, and structures on the parcel such as recreational and storage facilities [<xref ref-type="bibr" rid="scirp.78737-ref24">24</xref>] . The IF for residential footprint was determined by developing test sites in three randomly chosen locations in the study area: Richmond, Delta and Pitt Meadows. Ten residential footprints in each of the three locations were mapped in detail using Google Earth Pro satellite imagery. The residential footprint, according to the BC Ministry of Agriculture definition, includes: the main residence, yard, driveway, parking, and any auxiliary buildings or structures. The IF was calculated by dividing the TIA by the total residential footprint to get an overall average for the test areas of 29%, which was extrapolated across the entire study area.</p></sec><sec id="s2_2_3"><title>2.2.3. GIS Analysis</title><p>This study used ArcGIS 10 [<xref ref-type="bibr" rid="scirp.78737-ref36">36</xref>] to calculate the extent of impervious surface area covering Class 1 to Class 4 lands in the ALR. Data sources included:</p><p>・ Spatial land cover polygons and tabular land cover survey data from the Metro Vancouver (ALUI) [<xref ref-type="bibr" rid="scirp.78737-ref28">28</xref>] ;</p><p>・ 1:125,000 maps of the Soil Capability for Agriculture classification (Canada Land Inventory) [<xref ref-type="bibr" rid="scirp.78737-ref6">6</xref>] ;</p><p>・ Boundaries for the ALR based on shape files from the Provincial Agricultural Land Commission [<xref ref-type="bibr" rid="scirp.78737-ref34">34</xref>] .</p><p>The following land cover types were identified as impervious: crop barn, farm infrastructure, greenhouse, non-built or bare, residential footprint, settlement (structures and their associated yards, roads and parking), transportation, and utilities. The extent of impervious surface area for each land cover was calculated for each land capability class including: Classes 1 - 4 representing arable agricultural land; Classes 5 &amp; 6 indicating severe limitations only suitable for few crop types or requires special conservation practices; Class 7, with no capacity for arable culture; and Class 0, organic soil, which is not included in capability classes [<xref ref-type="bibr" rid="scirp.78737-ref7">7</xref>] .</p><p>ArcGIS 10 [<xref ref-type="bibr" rid="scirp.78737-ref36">36</xref>] was used to map the land cover types in the Metro Vancouver ALR, to present the varying levels of imperviousness according to the corresponding IFs.</p></sec></sec></sec><sec id="s3"><title>3. Results</title><sec id="s3_1"><title>3.1. Imagery Analysis</title><p>The initial imagery analysis estimated an average total impervious area of 9.7% in Metro Vancouver’s ALR (<xref ref-type="table" rid="table3">Table 3</xref>). The municipality of Delta had the highest</p><table-wrap id="table3" ><label><xref ref-type="table" rid="table3">Table 3</xref></label><caption><title> Averaged total impervious area (%) estimated for six municipalities in the Agricultural Land Reserve using Google earth (2016) imagery analysis</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Local Government</th><th align="center" valign="middle" >Total Impervious Area (%) per 260 hectares</th></tr></thead><tr><td align="center" valign="middle" >Delta</td><td align="center" valign="middle" >20</td></tr><tr><td align="center" valign="middle" >Langley (Township)</td><td align="center" valign="middle" >6</td></tr><tr><td align="center" valign="middle" >Surrey</td><td align="center" valign="middle" >5</td></tr><tr><td align="center" valign="middle" >Pitt Meadows</td><td align="center" valign="middle" >10</td></tr><tr><td align="center" valign="middle" >Richmond</td><td align="center" valign="middle" >9</td></tr><tr><td align="center" valign="middle" >Maple Ridge</td><td align="center" valign="middle" >8</td></tr><tr><td align="center" valign="middle" >ALR Average</td><td align="center" valign="middle" >10</td></tr></tbody></table></table-wrap><p>TIA, followed by Pitt Meadows. Two of the plots surveyed in Delta had large greenhouses, which covered approximately 20% and 54% of the area (<xref ref-type="fig" rid="fig2">Figure 2</xref>). Three out of the four plots in Pitt Meadows also contained large greenhouses. The remaining plots included a variety of land covers, with no single activity covering the majority of the area. Depending on the random selection of the plots, different land covers where captured, while others were not included in this approximation. The average TIA could appear extremely high or low if, for example, the plot included a greenhouse versus a cultivated field.</p></sec><sec id="s3_2"><title>3.2. Agricultural Land Use Inventory Analysis</title><p>Detailed spatial assessment of land cover data in Metro Vancouver ALUIs showed that approximately 6420 hectares, or 11% of the ALR contains a range of impervious farming, residential, commercial and industrial structures, and extraction or disposal sites that are partially impervious (<xref ref-type="table" rid="table4">Table 4</xref>). The highest percentage of impervious land cover types are residential footprint (1650 ha) and farm infrastructure (1580 ha), collectively making up 6% of the ALR study area.</p><p>The calculation for TIA included all land cover types from the Metro Vancouver ALUIs. The average TIA in the study area was calculated to be 7.2% (<xref ref-type="fig" rid="fig3">Figure 3</xref>). The highest TIA was again found in Delta at 8.1%, followed by Surrey at 7.6%. The TIA results are lower than those of the imagery analysis in 3.1, and are a more accurate representation, as they account for all of the impervious land covers in the ALR.</p><p>Farm infrastructure (2.4%) and transportation (0.9%) land cover types contributed the most to the overall percentage of TIA, but are considered necessary for the development of soil-based agriculture. Residential footprint (0.9%) and greenhouses (0.6%) were the next largest contributors to TIA in the ALR, and impede soil-based production.</p></sec><sec id="s3_3"><title>3.3. GIS Analysis of Soil Capability for Agriculture</title><p>Approximately 5260 ha of Class 2 - 4 lands, considered prime agricultural soils are covered by impervious surfaces in the study area. Transportation and residential footprints cover the largest area of prime soils (1686 ha and 961 ha,</p><fig id="fig2"  position="float"><label><xref ref-type="fig" rid="fig2">Figure 2</xref></label><caption><title> Impervious surface area in Delta, BC Agricultural Land Reserve study area (Google Earth Pro, 2014)</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/9-3001734x3.png"/></fig><fig id="fig3"  position="float"><label><xref ref-type="fig" rid="fig3">Figure 3</xref></label><caption><title> Average total impervious area (%) in the Metro Vancouver Agricultural Land Reserve study area, derived from GIS-based Agricultural Land Use Inventories from 2010 &amp; 2011 (Ministry of Agriculture, 2014)</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/9-3001734x4.png"/></fig><p>respectively) (<xref ref-type="table" rid="table5">Table 5</xref>). The majority of impervious surfaces are on Class 4 soil (3394 ha), representing 50% of the ALR in the study area (<xref ref-type="fig" rid="fig4">Figure 4</xref>). Organic soils have the second largest impervious area (1669 ha), however these are not rated for agricultural capability in the Canada Land Inventory (1972). There are no Class 1 soils in the study area for analysis. These reported values are smaller than reported in the ALUIs for Metro Vancouver, as not all the soil was surveyed in areas potentially covered by a range of land covers such as settlements, transportation and utilities.</p></sec></sec><sec id="s4"><title>4. Discussion</title><p>Approximately 11% of the ALR study area was determined to be impervious land covers to some degree. However, after applying the impervious factors, the</p><table-wrap id="table4" ><label><xref ref-type="table" rid="table4">Table 4</xref></label><caption><title> Total impervious area in the Metro Vancouver ALR study area derived from 2010 &amp; 2011 (Ministry of Agriculture, 2014)</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Metro Vancouver ALR Study Area</th><th align="center" valign="middle" >IF (%)</th><th align="center" valign="middle" >ALR (ha)</th><th align="center" valign="middle" >% of ALR</th><th align="center" valign="middle" >TIA (%)</th></tr></thead><tr><td align="center" valign="middle"  colspan="5"  >Impervious Land Covers</td></tr><tr><td align="center" valign="middle" >Residential footprint</td><td align="center" valign="middle" >30</td><td align="center" valign="middle" >1650</td><td align="center" valign="middle" >3</td><td align="center" valign="middle" >10</td></tr><tr><td align="center" valign="middle" >Farm infrastructure</td><td align="center" valign="middle" >80</td><td align="center" valign="middle" >1580</td><td align="center" valign="middle" >3</td><td align="center" valign="middle" >2.5</td></tr><tr><td align="center" valign="middle" >Non Built or Bare</td><td align="center" valign="middle" >3</td><td align="center" valign="middle" >730</td><td align="center" valign="middle" >1</td><td align="center" valign="middle" >0.03</td></tr><tr><td align="center" valign="middle" >Settlement</td><td align="center" valign="middle" >12</td><td align="center" valign="middle" >720</td><td align="center" valign="middle" >1</td><td align="center" valign="middle" >0.1</td></tr><tr><td align="center" valign="middle" >Transportation</td><td align="center" valign="middle" >90</td><td align="center" valign="middle" >630</td><td align="center" valign="middle" >1</td><td align="center" valign="middle" >1.0</td></tr><tr><td align="center" valign="middle" >Utilities</td><td align="center" valign="middle" >80</td><td align="center" valign="middle" >500</td><td align="center" valign="middle" >0.8</td><td align="center" valign="middle" >0.6</td></tr><tr><td align="center" valign="middle" >Greenhouses</td><td align="center" valign="middle" >80</td><td align="center" valign="middle" >440</td><td align="center" valign="middle" >0.5</td><td align="center" valign="middle" >0.6</td></tr><tr><td align="center" valign="middle" >Built up-Other</td><td align="center" valign="middle" >12</td><td align="center" valign="middle" >140</td><td align="center" valign="middle" >0.23</td><td align="center" valign="middle" >0.03</td></tr><tr><td align="center" valign="middle" >Crop Barn</td><td align="center" valign="middle" >80</td><td align="center" valign="middle" >23</td><td align="center" valign="middle" >0.03</td><td align="center" valign="middle" >0.02</td></tr><tr><td align="center" valign="middle" >Unmaintained greenhouses</td><td align="center" valign="middle" >80</td><td align="center" valign="middle" >2</td><td align="center" valign="middle" >0.003</td><td align="center" valign="middle" >0</td></tr><tr><td align="center" valign="middle" >Unmaintained crop barn</td><td align="center" valign="middle" >80</td><td align="center" valign="middle" >2</td><td align="center" valign="middle" >0.003</td><td align="center" valign="middle" >0</td></tr><tr><td align="center" valign="middle"  colspan="2"  >Total</td><td align="center" valign="middle" >6420</td><td align="center" valign="middle" >11</td><td align="center" valign="middle" >6</td></tr><tr><td align="center" valign="middle"  colspan="5"  >Pervious and Semi-pervious Land Covers</td></tr><tr><td align="center" valign="middle" >Cultivated Field Crops</td><td align="center" valign="middle" >3</td><td align="center" valign="middle" >26,970</td><td align="center" valign="middle" >44</td><td align="center" valign="middle" >1.3</td></tr><tr><td align="center" valign="middle" >Vegetated (natural)</td><td align="center" valign="middle" >0</td><td align="center" valign="middle" >15,860</td><td align="center" valign="middle" >26</td><td align="center" valign="middle" >0</td></tr><tr><td align="center" valign="middle" >Managed vegetation</td><td align="center" valign="middle" >3</td><td align="center" valign="middle" >2920</td><td align="center" valign="middle" >5</td><td align="center" valign="middle" >0.2</td></tr><tr><td align="center" valign="middle" >Wetlands</td><td align="center" valign="middle" >0</td><td align="center" valign="middle" >1960</td><td align="center" valign="middle" >3</td><td align="center" valign="middle" >0</td></tr><tr><td align="center" valign="middle" >Golf fairway/green</td><td align="center" valign="middle" >3</td><td align="center" valign="middle" >1500</td><td align="center" valign="middle" >2</td><td align="center" valign="middle" >0.06</td></tr><tr><td align="center" valign="middle" >Unused forage or pasture</td><td align="center" valign="middle" >3</td><td align="center" valign="middle" >640</td><td align="center" valign="middle" >1</td><td align="center" valign="middle" >0.03</td></tr><tr><td align="center" valign="middle" >Water bodies (anthropogenic)</td><td align="center" valign="middle" >0</td><td align="center" valign="middle" >500</td><td align="center" valign="middle" >0.8</td><td align="center" valign="middle" >0</td></tr><tr><td align="center" valign="middle" >Water bodies (natural)</td><td align="center" valign="middle" >0</td><td align="center" valign="middle" >470</td><td align="center" valign="middle" >0.8</td><td align="center" valign="middle" >0</td></tr><tr><td align="center" valign="middle"  colspan="2"  >Total</td><td align="center" valign="middle" >23,850</td><td align="center" valign="middle" >83</td><td align="center" valign="middle" >1.6</td></tr><tr><td align="center" valign="middle"  colspan="4"  >TIA Total</td><td align="center" valign="middle" >7.2</td></tr></tbody></table></table-wrap><table-wrap id="table5" ><label><xref ref-type="table" rid="table5">Table 5</xref></label><caption><title> Areas of impervious land covers on soil Classes 2 to 7 in the Metro Vancouver Agricultural Land Reserve study area (No Class 1 soil in the study area)</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Land Capability Classes</th><th align="center" valign="middle" >Crop Barn (ha)</th><th align="center" valign="middle" >Farm Infrastructure (ha)</th><th align="center" valign="middle" >Greenhouse (ha)</th><th align="center" valign="middle" >Non-built or Bare (ha)</th><th align="center" valign="middle" >Residential Footprint (ha)</th><th align="center" valign="middle" >Settlement (ha)</th><th align="center" valign="middle" >Transportation (ha)</th><th align="center" valign="middle" >Utilities (ha)</th></tr></thead><tr><td align="center" valign="middle" >Class 2</td><td align="center" valign="middle" >2</td><td align="center" valign="middle" >88</td><td align="center" valign="middle" >43</td><td align="center" valign="middle" >32</td><td align="center" valign="middle" >120</td><td align="center" valign="middle" >110</td><td align="center" valign="middle" >220</td><td align="center" valign="middle" >60</td></tr><tr><td align="center" valign="middle" >Class 3</td><td align="center" valign="middle" >5</td><td align="center" valign="middle" >190</td><td align="center" valign="middle" >65</td><td align="center" valign="middle" >112</td><td align="center" valign="middle" >200</td><td align="center" valign="middle" >140</td><td align="center" valign="middle" >360</td><td align="center" valign="middle" >120</td></tr><tr><td align="center" valign="middle" >Class 4</td><td align="center" valign="middle" >5</td><td align="center" valign="middle" >280</td><td align="center" valign="middle" >130</td><td align="center" valign="middle" >306</td><td align="center" valign="middle" >640</td><td align="center" valign="middle" >570</td><td align="center" valign="middle" >1110</td><td align="center" valign="middle" >350</td></tr><tr><td align="center" valign="middle" >Class 2 - 4</td><td align="center" valign="middle" >12</td><td align="center" valign="middle" >560</td><td align="center" valign="middle" >240</td><td align="center" valign="middle" >450</td><td align="center" valign="middle" >960</td><td align="center" valign="middle" >815</td><td align="center" valign="middle" >1690</td><td align="center" valign="middle" >540</td></tr><tr><td align="center" valign="middle" >Class 5</td><td align="center" valign="middle" >0.0</td><td align="center" valign="middle" >24</td><td align="center" valign="middle" >49</td><td align="center" valign="middle" >110</td><td align="center" valign="middle" >86</td><td align="center" valign="middle" >92</td><td align="center" valign="middle" >300</td><td align="center" valign="middle" >23</td></tr><tr><td align="center" valign="middle" >Class 6</td><td align="center" valign="middle" >0.0</td><td align="center" valign="middle" >0.0</td><td align="center" valign="middle" >0.0</td><td align="center" valign="middle" >0.1</td><td align="center" valign="middle" >4</td><td align="center" valign="middle" >4</td><td align="center" valign="middle" >1</td><td align="center" valign="middle" >3</td></tr><tr><td align="center" valign="middle" >Class 7</td><td align="center" valign="middle" >0.0</td><td align="center" valign="middle" >8</td><td align="center" valign="middle" >26</td><td align="center" valign="middle" >10</td><td align="center" valign="middle" >41</td><td align="center" valign="middle" >33</td><td align="center" valign="middle" >40</td><td align="center" valign="middle" >34</td></tr><tr><td align="center" valign="middle" >Organic</td><td align="center" valign="middle" >3</td><td align="center" valign="middle" >71</td><td align="center" valign="middle" >25</td><td align="center" valign="middle" >35</td><td align="center" valign="middle" >400</td><td align="center" valign="middle" >220</td><td align="center" valign="middle" >500</td><td align="center" valign="middle" >97</td></tr><tr><td align="center" valign="middle" >Total</td><td align="center" valign="middle" >14</td><td align="center" valign="middle" >660</td><td align="center" valign="middle" >340</td><td align="center" valign="middle" >925</td><td align="center" valign="middle" >1490</td><td align="center" valign="middle" >1160</td><td align="center" valign="middle" >2520</td><td align="center" valign="middle" >690</td></tr></tbody></table></table-wrap><fig id="fig4"  position="float"><label><xref ref-type="fig" rid="fig4">Figure 4</xref></label><caption><title> Total impervious land cover (ha) per land capability class in the Metro Vancouver Agricultural Land Reserve study area (No Class 1 soils in the study area)</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/9-3001734x5.png"/></fig><p>Metro Vancouver ALR study area was determined to have a total impervious area of 7.2% (<xref ref-type="fig" rid="fig3">Figure 3</xref>). This may provide an assessment of the amount of precipitation that is no longer intercepted by vegetation or infiltrated into the soil and contributes to runoff. In forested catchments, runoff is reported to double with increases of 10% - 20% impervious surface cover (ISC), triple at 35% - 55%, and increase more than five times at 75% - 100% (ISC) [<xref ref-type="bibr" rid="scirp.78737-ref37">37</xref>] . As the study area is 7.2% impervious, the volume of runoff generated from impervious surfaces would also almost double. This is substantial as Metro Vancouver annually receives &lt;1350 mm in low rainfall zones and &gt;1850 mm in high rainfall zones [<xref ref-type="bibr" rid="scirp.78737-ref32">32</xref>] .</p><p>In the Metro Vancouver ALR, drainage is one of the limiting factors for crop production, as the soils are composed of fine silts and clay, characterized by low infiltration rates [<xref ref-type="bibr" rid="scirp.78737-ref38">38</xref>] . To manage high water tables, agricultural producers install surface and subsurface drainage and construct ditches [<xref ref-type="bibr" rid="scirp.78737-ref38">38</xref>] . Excess runoff from impervious surfaces largely drains into these ditches, which are often connected to natural and channelized streams that drain into larger rivers and the ocean [<xref ref-type="bibr" rid="scirp.78737-ref38">38</xref>] . Increased impervious surface area reduces the capacity of soils to store water, and as the soil becomes saturated peak surface runoff rates increase [<xref ref-type="bibr" rid="scirp.78737-ref22">22</xref>] . Runoff from agricultural fields may carry with it several non-point sources of pollution including pesticides, sediments, fertilizers and manure, which degrades water quality and aquatic health [<xref ref-type="bibr" rid="scirp.78737-ref3">3</xref>] .</p><p>Future climate change scenarios project increases in total precipitation and more frequent and extreme weather events in Metro Vancouver, which will further magnify drainage and flooding issues in the ALR [<xref ref-type="bibr" rid="scirp.78737-ref39">39</xref>] . Projections suggest that peak flows will occur earlier in the spring, as more precipitation falls as rain in the winter. They also indicate that flows will decrease in the summer, as snowpack is reduced in the winter. This would result in reduced groundwater recharge rates and base flow at times of low flow in the late summer [<xref ref-type="bibr" rid="scirp.78737-ref23">23</xref>] . In dry years, base flow helps to sustain stream flow and many streams would not exist without it [<xref ref-type="bibr" rid="scirp.78737-ref22">22</xref>] . Increased impervious surface area would exacerbate these effects, as drainage systems spatially concentrate distributed runoff and diminish opportunities for precipitation to infiltrate into the ground [<xref ref-type="bibr" rid="scirp.78737-ref22">22</xref>] . This compounds the impacts of impervious surfaces and reduces the amount of precipitation that contributes to groundwater recharge and replenishes the base flow of streams.</p><p>The TIA values for the Metro Vancouver ALR study area are below the 10% imperviousness threshold, where it is generally accepted that aquatic health becomes degraded, as a guideline for urban watersheds [<xref ref-type="bibr" rid="scirp.78737-ref18">18</xref>] . Applying this threshold to the ALR does not provide a complete picture from a watershed perspective. The boundaries of the ALR are not compatible with watershed boundaries, as they were established based on agricultural capability and suitability, and not on the basis of topography and how water drains from ridges into rivers or creeks at lower elevations [<xref ref-type="bibr" rid="scirp.78737-ref8">8</xref>] . Moreover, the majority of watersheds in Metro Vancouver are associated with urban land uses, which often have impervious factors ranging from 80% to 90% [<xref ref-type="bibr" rid="scirp.78737-ref24">24</xref>] .</p><p>The ALR contains a variety of different land covers such as cultivated field crops, natural and semi-natural vegetation, residential developments, utilities and transportation. There are several sites in the ALR study area that are largely impervious, with land covers that range from 25% to 100% impervious (<xref ref-type="fig" rid="fig5">Figure 5</xref>). Concentrated areas of ISC that receive high amounts of rainfall can become hot spots for the transportation of agricultural non-point sources (NPS) of pollution, especially when located next to water bodies. In agricultural areas water quality has been shown to decline when agricultural land use exceeds 50% of the watershed area [<xref ref-type="bibr" rid="scirp.78737-ref27">27</xref>] . Considering the ALR study area is 50% farmed, the IBI scores would be expected to show a decline in water quality regardless of the TIA. Metro Vancouver has already experienced issues associated with the management of nutrients including: nitrate contamination of groundwater, eutrophication, pathogens and ammonia in surface water and excessive amounts of potassium in the soil [<xref ref-type="bibr" rid="scirp.78737-ref39">39</xref>] . The Serpentine and Nicomekl Rivers for example, which flow through the Surrey ALR into the Strait of Georgia are at risk of contamination from both urban and agricultural land uses and often fail to meet fecal contamination, turbidity and dissolved oxygen objectives [<xref ref-type="bibr" rid="scirp.78737-ref40">40</xref>] . Increasing impervious surface area in the ALR will only exacerbate potential water quality issues associated with agricultural NPS pollution, and push thresholds towards ecosystem degradation.</p><sec id="s4_1"><title>4.1. Soil Sealing</title><p>There is an estimated 60,940 hectares of land in the Metro Vancouver ALR, 80% of which is considered prime soil for agriculture (Class 2 - 4) [<xref ref-type="bibr" rid="scirp.78737-ref11">11</xref>] . In the study area, approximately 5256 ha or 10% of Class 2 - 4 soils (<xref ref-type="fig" rid="fig4">Figure 4</xref>) are covered by impervious surfaces associated with agricultural, residential, commercial and</p><fig id="fig5"  position="float"><label><xref ref-type="fig" rid="fig5">Figure 5</xref></label><caption><title> Impervious factors of land cover types in six municipalities in the Metro Vancouver Agricultural Land Reserve study area</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/9-3001734x6.png"/></fig><p>industrial development. These soils are now permanently sealed and lost to soil-based agriculture.</p><p>While farm infrastructure and transportation comprise a large TIA, these are considered essential for agricultural enterprise. Residential development and greenhouses are another significant area of TIA and are arguably less essential and of concern from the perspective of soil-based agriculture. Therefore, these are the focus of further discussion.</p><sec id="s4_1_1"><title>4.1.1. Residential Footprint</title><p>One of the concerns in the ALR with regards to residential footprint (RF) is the increased number of estate-sized homes (&gt;465 m<sup>2</sup> or 5000 ft<sup>2</sup>) encroaching into farmland. Test areas determined that residential imperviousness was 29%. This exceeds estimates for rural residential land uses, which range from 8% - 12% impervious [<xref ref-type="bibr" rid="scirp.78737-ref13">13</xref>] [<xref ref-type="bibr" rid="scirp.78737-ref24">24</xref>] . Test-plots in this study included estates with tennis courts, swimming pools, auxiliary buildings and long driveways, not typically expected in rural agricultural areas.</p><p>Residential footprint covers approximately 1650 hectares of land, or 3% of the Metro Vancouver ALR study area, of which 58% (960 ha) is on Class 2 - 4 soils. While residential housing is a necessary accessory to agriculture, it can also hinder the viability of agriculture. This occurs when the primary motivation for purchasing farmland is for residential use and not for agriculture [<xref ref-type="bibr" rid="scirp.78737-ref28">28</xref>] . Currently, only 4780 or 41% of the parcels in the Metro Vancouver ALR are used for farming [<xref ref-type="bibr" rid="scirp.78737-ref28">28</xref>] . When residences are built without the purpose of farming there is often little consideration for agriculture, and homes are built to both maximize views and privacy and are setback from the road [<xref ref-type="bibr" rid="scirp.78737-ref41">41</xref>] . This can increasingly fragment the land, compromising future agricultural use [<xref ref-type="bibr" rid="scirp.78737-ref41">41</xref>] .</p></sec><sec id="s4_1_2"><title>4.1.2. Greenhouses</title><p>Imagery analysis of impervious surfaces in Metro Vancouver’s ALR found an average TIA of 9.7% (<xref ref-type="table" rid="table2">Table 2</xref>). Delta had the highest TIA at 20%, as it was sampled to have a higher occurrence of greenhouses. Two of the plots (160 acres each) in Delta had greenhouses on them, including one that was 53.6% impervious. Pitt Meadows had the second highest TIA (10.3%), which was also largely attributed to greenhouses. This preliminary analysis shows how drastically areas with concentrated greenhouses can increase imperviousness in localized agricultural areas.</p><p>Glass greenhouses occupy 441 ha (not including poly or unmaintained green- houses) of land or less than 1% of Metro Vancouver’s ALR [<xref ref-type="bibr" rid="scirp.78737-ref28">28</xref>] . This represents a large portion of the greenhouses in BC, as 90% of the greenhouses in BC are located in the Lower Mainland [<xref ref-type="bibr" rid="scirp.78737-ref42">42</xref>] . The highest concentration is in Delta, which accounts for 52% of the production in the province [<xref ref-type="bibr" rid="scirp.78737-ref42">42</xref>] . Greenhouses in Delta cover approximately 2% of the ALR (162 ha) and are the most impervious land cover type.</p><p>The GIS analysis showed that approximately 240 ha or 70% of all greenhouses in the study area are on Class 2 - 4 soil (<xref ref-type="table" rid="table5">Table 5</xref>). In areas such as Delta this can account for a large portion of land, as almost all the soils when drained are highly productive and climatically suited for the crops in the area [<xref ref-type="bibr" rid="scirp.78737-ref42">42</xref>] . Covering these soils with greenhouses arguably defeats the purpose of protecting productive agricultural soil in the ALR for food production.</p></sec></sec></sec><sec id="s5"><title>5. Conclusions</title><p>Large areas of fertile soil are lost each year to impervious surfaces, as agricultural land is converted to urban, commercial and industrial uses, decreasing the viability of soil based-agriculture. Impervious surfaces also disrupt the hydrologic cycle, reducing the infiltration of water into the ground and accelerating runoff and the transport of contaminants into waterways. The Agricultural Land Reserve in British Columbia was originally created to protect arable soil from being lost to urban and other development. However, there is a wide array of permitted activities in the ALR that “seal” the soil with impervious surfaces, such as the construction of greenhouses, crop and dairy barns, estate-sized houses, roads, and commercial facilities.</p><p>After applying impervious factors to account for variances in permeability the total impervious area of the Metro Vancouver ALR is approximately 7.2%, which is below the 10% threshold that correlates with changes in ecosystem health. However, this calculation was based on ALR boundaries, not watershed boundaries, so they are not comparable. If watersheds were to be considered, they would contain a higher TIA value by incorporating neighboring urban developments. Moreover, the impacts of agricultural land use in terms of water quality and quantity would be expected to compound the effects of TIA on watershed health. Runoff generated from impervious surfaces in the ALR largely drains into ditches that are often connected to fish bearing streams or waterways that eventually drain into the ocean.</p><p>Approximately 5260 hectares or 10% of the arable agricultural soil (Class 1 - 4) has been covered with impervious surfaces in the Metro Vancouver ALR, with the highest percentage of that being for transportation (1690 ha) and residential development (960 ha). Residential development and greenhouses are considered to be of the most concern, as they inhibit the viability of soil-based agriculture. Approximately 70% of greenhouses (240 ha) and 58% of residential developments are on Class 2 - 4 soil in the ALR, permanently eliminating soil-based production on these prime lands.</p><p>The fertile soils found in the ALR are part of the region’s natural capital, providing nutrients and water for food to grow, and regulating floods and droughts and the movement of contaminants into waterways. Long-term action and planning is suggested between governments, communities, landowners, the Agricultural Land Commission, developers, environmental stewardship and other stake- holder groups for the protection of arable agricultural land from being covered with impervious surfaces and to achieve the goal of the long-term sustainability and viability of British Columbia’s agricultural land base. This could include limiting the construction of impervious structures to areas that are less productive (Class 5 - 7 soil), and restricting the size of residential footprints and regulating the siting of residential uses. Specific targets to limit total impervious area in watersheds to less than 10% would further help protect stream health and conserve valuable soil resources.</p></sec><sec id="s6"><title>Acknowledgements</title><p>The authors would like to thank the British Columbia Ministry of Agriculture for providing the land cover data results from the summer 2010-2011 Agriculture Land Use Inventory and Agricultural Water Demand Projects in BC. We would also like to thank HaojingXie for her assistance in producing the GIS maps.</p></sec><sec id="s7"><title>Cite this paper</title><p>Rose, A., Wilson, J.E. and Lavkulich, L.M. (2017) Analysis of Impervious Surface Area, and the Impacts on Soil-Based Agriculture and the Hydrologic Cycle: A Case Study in the Agricultural Land Reserve in Metro Vancouver, British Columbia, Canada. Agricultural Sciences, 8, 837-856. https://doi.org/10.4236/as.2017.88062</p></sec></body><back><ref-list><title>References</title><ref id="scirp.78737-ref1"><label>1</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>D&amp;#214;&amp;#214;s</surname><given-names> B.R. </given-names></name>,<etal>et al</etal>. (<year>2002</year>)<article-title>Population Growth and Loss of Arable Land</article-title><source> Global Environmental Change</source><volume> 12</volume>,<fpage> 303</fpage>-<lpage>311</lpage>.<pub-id pub-id-type="doi"></pub-id></mixed-citation></ref><ref id="scirp.78737-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Lin, J. (2014) Human Activity and the Environment: Agriculture in Canada. Catalogue No. 16-201-X, Statistics Canada.  
http://www.statcan.gc.ca/pub/16-201-x/16-201-x2014000-eng.pdf</mixed-citation></ref><ref id="scirp.78737-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Allan, J.D. (2004) Landscapes and Riverscapes: The Influence of Land Use on Stream Ecosystems. Annual Review of Ecology, Evolution, and Systematics, 35, 257-284. https://doi.org/10.1146/annurev.ecolsys.35.120202.110122</mixed-citation></ref><ref id="scirp.78737-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">European Commission (2012) In-Depth Report: Soil Sealing. Science for Environmental Policy.  
http://ec.europa.eu/environment/archives/soil/pdf/sealing/Soil%20Sealing%20In-depth%20Report%20March%20version_final.pdf</mixed-citation></ref><ref id="scirp.78737-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Hofmann, N., Filoso, G. and Schofield, M. (2005) The Loss of Dependable Agricultural Land in Canada. Rural and Small Town Canada Analysis Bulletin, 6, Catalogue No. 21-006-XIE.  
http://publications.gc.ca/collections/Collection/Statcan/21-006-X/21-006-XIE2005001.pdf</mixed-citation></ref><ref id="scirp.78737-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Canada Land Inventory (1972) Soil Capability Classification. Department of the Environment, Ottawa.</mixed-citation></ref><ref id="scirp.78737-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Agriculture and Agri-Food Canada (2013) Overview of Classification Methodology for Determining Land Capability for Agriculture.  
http://sis.agr.gc.ca/cansis/nsdb/cli/classdesc.html</mixed-citation></ref><ref id="scirp.78737-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Agricultural Land Commission (ALC) (2013) Agricultural Capability Classification in BC. Province of British Columbia. 
http://www.alc.gov.bc.ca/assets/alc/assets/library/agricultural-capability/agriculture_capability_classification_in_bc_2013.pdf</mixed-citation></ref><ref id="scirp.78737-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Agricultural Land Commission (ALC) (2014a) Agricultural Land Reserve History. Province of British Columbia.  
http://www.alc.gov.bc.ca/alc/content/alr-maps/alr-history</mixed-citation></ref><ref id="scirp.78737-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">Smith, B.E. (2012) A Work in Progress: The British Columbia Farmland Preservation Program. BC Government and Agricultural Land Commission eBook Collection.  
http://www.alc.gov.bc.ca/assets/alc/assets/library/archived-publications/alr-history/a_work_in_progress_-_farmland_preservation_b_smith_2012.pdf</mixed-citation></ref><ref id="scirp.78737-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">Metro Vancouver (2011) Metro Vancouver’s Regional Food System Strategy. Sustainable Region Initiative, Metro Vancouver, Feb. 2011.  
https://www.metrovancouver.org/services/regional-planning/PlanningPublications/RegionalFoodSystemStrategy.pdf</mixed-citation></ref><ref id="scirp.78737-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">BC Ministry of Agriculture and Lands (BCMARR) (2006) B.C.’s Food Self-Reliance. Province of British Columbia, 820105-1. http://www.plancanada.com/BCfood.pdf</mixed-citation></ref><ref id="scirp.78737-ref13"><label>13</label><mixed-citation publication-type="other" xlink:type="simple">Guthrie, R. and Deniseger, J. (2001) Impervious Surfaces in French Creek. BC Ministry of Water, Land and Air Protection, Victoria. 
http://www.env.gov.bc.ca/van-island/es/pdf/Impervious%20Surfaces%20technical%20document.pdf</mixed-citation></ref><ref id="scirp.78737-ref14"><label>14</label><mixed-citation publication-type="other" xlink:type="simple">US Environmental Protection Agency (2016) CADDIS Volume 2: Sources, Stressors and Responses. https://www3.epa.gov/caddis/ssr_urb_is3.html</mixed-citation></ref><ref id="scirp.78737-ref15"><label>15</label><mixed-citation publication-type="other" xlink:type="simple">Walsh, C., Roy, A., Feminella, J., Cottingham, P., Groffman, P. and Morgan, R. (2005) The Urban Stream Syndrome: Current Knowledge and the Search for a Cure. The North American Benthological Society, 24, 706-723.</mixed-citation></ref><ref id="scirp.78737-ref16"><label>16</label><mixed-citation publication-type="other" xlink:type="simple">Chithra, S.V., Nair, M.V.H., Amarnath, A. and Anjana, N.S. (2015) Impact of Impervious Surfaces on the Environment. International Journal of Engineering Science Invention, 4, 27-31.</mixed-citation></ref><ref id="scirp.78737-ref17"><label>17</label><mixed-citation publication-type="other" xlink:type="simple">Stepenuck, K.F., Crunkilton, R.L. and Wang, L. (2002) Impacts of Urban Land Use on Macroinvertebrate Communities in Southeastern Wisconsin Streams. Journal of the American Water Resources Association, 38, 1041-1051.  
https://doi.org/10.1111/j.1752-1688.2002.tb05544.x</mixed-citation></ref><ref id="scirp.78737-ref18"><label>18</label><mixed-citation publication-type="other" xlink:type="simple">Beach, D. (2002) Coastal Sprawl the Effects of Urban Design on Aquatic Ecosystems in the United States. South Carolina Coastal Conservation League: Pew Oceans Commission. 
http://www.pewtrusts.org/en/research-and-analysis/reports/2002/04/08/coastal-sprawl-the-effects-of-urban-design-on-aquatic-ecosystems-in-the-united-states</mixed-citation></ref><ref id="scirp.78737-ref19"><label>19</label><mixed-citation publication-type="other" xlink:type="simple">Wang, L., Lyons, J. and Kanehl, P. (2001) Impacts of Urbanization on Stream Habitat and Fish across Multiple Spatial Scales. Environmental Management, 28, 255-266.</mixed-citation></ref><ref id="scirp.78737-ref20"><label>20</label><mixed-citation publication-type="other" xlink:type="simple">Jartun, M., Ottesen, R.T., Steinnes, E. and Volden, T. (2008) Runoff of Particle Bound Pollutants from Urban Impervious Surfaces Studied by Analysis of Sediments from Stormwater Traps. Science of the Total Environment, 396, 147-163.</mixed-citation></ref><ref id="scirp.78737-ref21"><label>21</label><mixed-citation publication-type="other" xlink:type="simple">Soller, J., Stephenson, J., Olivieri, K., Downing, J. and Olivieri, A.W. (2005) Evaluation of Seasonal Scale First Flush Pollutant Loading and Implications for Urban Runoff Management. Journal of Environmental Management, 76, 309-318.</mixed-citation></ref><ref id="scirp.78737-ref22"><label>22</label><mixed-citation publication-type="other" xlink:type="simple">Shuster, W., Bonta, J., Thurston, H., Warnemuende, E. and Smith, D. (2005) Impacts of Impervious Surface on Watershed Hydrology: A Review. Urban Water Journal, 2, 263-275. https://doi.org/10.1080/15730620500386529</mixed-citation></ref><ref id="scirp.78737-ref23"><label>23</label><mixed-citation publication-type="other" xlink:type="simple">Crawford, E. and MacNair, E. (2012) BC Agriculture Climate Change Adaptation Risk + Opportunity Assessment: Provincial Report. BC Agriculture and Food Climate Action Initiative. 
http://pics.uvic.ca/sites/default/files/uploads/BC%20Agriculture_full%20report.pdf</mixed-citation></ref><ref id="scirp.78737-ref24"><label>24</label><mixed-citation publication-type="other" xlink:type="simple">Mason, B. and Knight, R. (2001) Sensitive Habitat Inventory and Mapping. Community Mapping Network, Vancouver, 315 p.  
http://cmnmaps.ca/cmn/files/methods/shim_manual/shim_10jun02.pdf</mixed-citation></ref><ref id="scirp.78737-ref25"><label>25</label><mixed-citation publication-type="other" xlink:type="simple">Beyerlein, D. (1996) Effective Impervious Area: The Real Enemy. Proceedings of the Impervious Surface Reduction Research Symposium, Olympia.</mixed-citation></ref><ref id="scirp.78737-ref26"><label>26</label><mixed-citation publication-type="other" xlink:type="simple">Kerr Wood Leidal Consulting Engineers (KWL) (2005) Template for Integrated Stormwater Management Planning 2005. Report for Greater Vancouver Regional District, Kerr Wood Leidal Associates Ltd., Burnaby. 
http://www.waterbucket.ca/rm/sites/wbcrm/documents/media/18.pdf</mixed-citation></ref><ref id="scirp.78737-ref27"><label>27</label><mixed-citation publication-type="other" xlink:type="simple">Page, N., Rood, K., Holz, T., Zandbergen, P., Horner, R. and McPhee, M. (1999) Proposed Watershed Classification System for Stormwater Management in the Greater Vancouver Sewerage and Drainage District Area. Greater Vancouver Sewerage and Drainage District, Burnaby. 
http://www.raincoastappliedecology.ca/wp-content/uploads/2012/05/Proposed-Watershed-Classification-System-for-the-GVSDD-Area.pdf</mixed-citation></ref><ref id="scirp.78737-ref28"><label>28</label><mixed-citation publication-type="other" xlink:type="simple">BC Ministry of Agriculture (2014) The Metro Vancouver Regional Report (Summer 2010 &amp; 2011) Agriculture Land Use Inventory. Ref No. 800.510-29.2014.  
http://www2.gov.bc.ca/assets/gov/farming-natural-resources-and-industry/agriculture-and-seafood/agricultural-land-and-environment/strengthening-farming/land-use-inventories/metrovanregional2010_11_aluireport.pdf</mixed-citation></ref><ref id="scirp.78737-ref29"><label>29</label><mixed-citation publication-type="other" xlink:type="simple">Metro Vancouver (2017) About Agriculture in the Region.  
http://www.metrovancouver.org/services/regional-planning/agriculture/about/Pages/default.aspx</mixed-citation></ref><ref id="scirp.78737-ref30"><label>30</label><mixed-citation publication-type="other" xlink:type="simple">Robinson, J.L. (2015) Fraser River Lowland. The Canadian Encyclopedia.  
http://www.thecanadianencyclopedia.ca/en/article/fraser-river-lowland</mixed-citation></ref><ref id="scirp.78737-ref31"><label>31</label><mixed-citation publication-type="other" xlink:type="simple">Luttmerding, H.A. (1981) Soils of the Langley—Vancouver Map Area. RAB Bulletin 18, BC Ministry of Environment, Map 6.</mixed-citation></ref><ref id="scirp.78737-ref32"><label>32</label><mixed-citation publication-type="other" xlink:type="simple">Kerr Wood Leidal Consulting Engineers (KWL) (2015) Options for a Region-Wide Baseline for On-Site Rainwater Management. Prepared for Metro Vancouver, Liquid Waste Services Department, Kerr Wood Leidal Associates Ltd., Burnaby.</mixed-citation></ref><ref id="scirp.78737-ref33"><label>33</label><mixed-citation publication-type="other" xlink:type="simple">Google Earth Pro (2014) Delta Agricultural Land Reserve, BC, Canada (49°04’03.45"N, 123°06’06.68"W, Eye alt. 1.2 km).</mixed-citation></ref><ref id="scirp.78737-ref34"><label>34</label><mixed-citation publication-type="other" xlink:type="simple">Agricultural Land Commission (ALC) (2014b) Agricultural Land Reserve Maps and GIS. http://www.alc.gov.bc.ca/alc/content/alr-maps/maps-and-gis</mixed-citation></ref><ref id="scirp.78737-ref35"><label>35</label><mixed-citation publication-type="other" xlink:type="simple">Earth Point (2016) KML Shapes—Polygon Area, Linestring Length, Placemark Point. http://www.earthpoint.us/Shapes.aspx</mixed-citation></ref><ref id="scirp.78737-ref36"><label>36</label><mixed-citation publication-type="other" xlink:type="simple">Environmental Systems Research Institute (ESRI) (2014) ArcGIS Desktop Version 10.2.</mixed-citation></ref><ref id="scirp.78737-ref37"><label>37</label><mixed-citation publication-type="other" xlink:type="simple">Paul, M.J. and Meyer, J.L. (2001) Streams in the Urban Landscape. Annual Review of Ecology and Systematics, 32, 333-365.</mixed-citation></ref><ref id="scirp.78737-ref38"><label>38</label><mixed-citation publication-type="other" xlink:type="simple">Dobb, A. (2013) BC Farm Practices &amp; Climate Change Adaptation: Drainage. British Columbia Agriculture &amp; Food Climate Action Initiative. 
http://www.bcagclimateaction.ca/wp/wp-content/media/FarmPractices-Drainage.pdf</mixed-citation></ref><ref id="scirp.78737-ref39"><label>39</label><mixed-citation publication-type="other" xlink:type="simple">Fraser Basin Council (2001) Nutrient Management Planning Strategies for the Fraser Valley. BC Ministry of Environment Lands &amp; Parks. 
http://www.fraserbasin.bc.ca/_Library/FVR/report_fvr_nutrient_mgt_2001.pdf</mixed-citation></ref><ref id="scirp.78737-ref40"><label>40</label><mixed-citation publication-type="other" xlink:type="simple">Bull, J. and Freyman, L. (2013) Status of Water Quality Objectives Attainment in the Little Campbell River, Serpentine River and Nicomekl River, 1971-2009. BC Ministry of Environment.</mixed-citation></ref><ref id="scirp.78737-ref41"><label>41</label><mixed-citation publication-type="other" xlink:type="simple">BC Ministry of Agriculture (2011) Bylaw Standard for Residential Uses in the Agricultural Land Reserve. 
http://www2.gov.bc.ca/gov/content/industry/agriculture-seafood/agricultural-land-and-environment/strengthening-farming/local-government-bylaw-standards-and-farm-bylaws/residential-uses-in-the-alr</mixed-citation></ref><ref id="scirp.78737-ref42"><label>42</label><mixed-citation publication-type="other" xlink:type="simple">Zbeetnoff Agro-Environmental Consulting and Quadra Planning Consultants (2011) Delta Agricultural Plan, Phase 1: Delta Agriculture Profile. The Corporation of Delta.  
http://www.delta.ca/docs/default-source/community-planning-and-development/agricultural-plan/delta-agricultural-profile.pdf?sfvrsn=0</mixed-citation></ref></ref-list></back></article>