Feb 24, 2010 - Hunting and Fishing in Quebec, developed by a private mapping company (5). Data from 1999 were used to identify the location of every deer ...
Supporting Information Raudsepp-Hearne et al. 10.1073/pnas.0907284107 SI Text Quantifying Ecosystem Services. Information for crops, pork, and
maple syrup was taken from the 2001 Canadian Agricultural Census, which provides data at the level of municipalities. In the Agricultural Census, these are called Census Consolidated Subdivisions (CCS). CCS represent single municipalities, or in some cases, the consolidation of several municipalities if the municipalities in question are judged to have too few farms to allow the privacy of individual farmers to be protected. To extract information on individual municipalities from CCS, ecosystem service (ES) values were estimated for individual municipalities by weighting the CCS values by area of the consolidated municipalities. Crops (Percentage of Municipal Land Dedicated to Crop Production).
Crop production was estimated using 2001 Canadian Agricultural Census data (1). The category Land Under Crop gives the amount of hectares under cultivation. Land Under Crop was chosen as an indicator, because crop-yield values were less comparable across the two watersheds because of differences in crop plants. Land Under Cultivation is an indicator of how important agriculture is to each municipality in terms of how much land in total is dedicated to this economic activity. Pork (Number of Pigs Produced per Square Kilometer). Pork production data were taken from the 2001 Canadian Agricultural Census (2). Because of privacy protection laws, the number of pork-producing farms but not the number of pigs produced on these farms was available in some cases. In these cases, the average number of pigs per farm in all of the surrounding municipalities (those that touched the municipality in question) was used to estimate the number of pigs per farm and thus, the total number of pigs produced. Drinking Water (Rank 1–5 Including Nonintegers Where Average Values Were Calculated from More than One Sampling Station).
The provision of drinking water was estimated from a provincial government water-quality database (3). The majority of water for drinking in the region is harvested from surface waters where water-quality testing stations are located. We used a water-quality index called the indice de qualité bactériologique et physicochimique (IQBP) that is used by the provincial government to assess the raw water supply intended for consumption. The index is based on conventional physicochemic and bacteriologic waterquality parameters and combines eight variables: phosphorus, fecal coliforms, turbidity, ammonia nitrogen, nitrates/nitrites, chlorophyll a, dissolved oxygen, and pH. Long-term water-quality data from the government database was available throughout the watershed but only on the main stems of rivers and their major tributaries. Values for 2001 were estimated by calculating the average of all water-quality samples taken in that year. To estimate water-quality values for municipalities that touched rivers where samples were taken, IQBP values from stations falling within each municipality were used. If more than one sampling station fell within a municipality, the average values of water quality across sampling stations was calculated. In cases where municipalities did not touch a river where sampling was done, one of two methods was used to estimate a value for water quality: (i) municipalities drawing their drinking water from a known location and treatment plant were assigned water-quality values for the area of river where the water is drawn from, or (ii) municipalities not drawing their Raudsepp-Hearne et al. www.pnas.org/cgi/content/short/0907284107
drinking water from surface water were assigned a neutral value of 3.5 (not considered to be good or bad water quality). Maple Syrup (Number of Maple-Syrup Taps per Square Kilometer).
Maple-syrup production was estimated using the 2001 Canadian Agricultural Census (4). This dataset gives the number of maplesyrup taps and maple-syrup farms for each municipality. Maplesyrup taps were chosen as an indicator of production, because the size of each maple-syrup farm varied widely. Because of privacy protection laws, the number of maple-syrup–producing farms but not the number of maple-syrup taps was available in some cases. In these cases, the average number of taps per farm in all surrounding municipalities (those that touched the municipality in question) was used to estimate the number of maple-syrup taps per farm and thus, the total number of taps. Deer Hunting (Number of Deer Killed per Square Kilometer). Hunting data were obtained from the database called The Guide to Hunting and Fishing in Quebec, developed by a private mapping company (5). Data from 1999 were used to identify the location of every deer killed that year. Hunters are required to mark on a map the location where they shot the deer when they register the dead deer at their local tagging station. In ArcGIS, the number of points where deer were killed was calculated for each municipal polygon. Tourist Attractions (Number of Tourist Attractions per Square Kilometer). The number of tourist attractions per municipality
was estimated using a provincial government database (6). A document accompanying the database lists the criteria used to determine whether or not a site is included in the tourism database (criteria include factors such as having infrastructure aimed at hosting tourists, being registered with the provincial regulating body, conforming with local and provincial laws and regulations, etc. as well as more specific criteria pertaining to types of tourist sites) (7). Estimations of the number of tourist attractions per municipality included only tourist attractions themselves and not related services such as restaurants and lodging. Nature Appreciation (Number of Reported Sightings of Rare Species per Square Kilometer). These data were collected by the Ministry for
Sustainable Development, Environment and Parks (MDDEP) and the Ministry for Natural Resources and Fauna (MRNF) from 1988 to 2007 and were available as a shape file for the region requested (8, 9). The number of points where rare species of flora and fauna were observed was calculated in ArcGIS for each municipal polygon. These data included observations of endangered or threatened floral and faunal species within the two study watersheds [in accordance with classifications under Convention on International Trade in Endangered Species (CITES)]. Observations were not considered to represent an accurate portrait of rare and endangered species in the region; Rather, they are a reflection of areas that have been frequented by scientists and nature enthusiasts, because these species were observed. Summer Cottages (Tax Value of Cottages per Square Kilometer). Data for estimating values of summer cottages were taken from a provincial database on property taxes by category of lodging (10). The summer-cottage category of housing (category K—chalet) was available by municipality. Additional data on the number of summer cottages per municipality were collected from each municipality by calling all municipal offices to check whether or not the tax-value data reflected the extent of this land use. Tax 1 of 3
value of cottages and number of cottages were highly correlated; the tax value was used, because it was considered to better reflect the value or importance of this ES across municipalities. Forest Recreation (Percent of Municipal Land Covered by Forest).
Forest recreation was assessed using forest cover as an indicator. Forest cover was estimated from land-cover maps (11) joined to a layer of municipal boundaries. The areas of all forest fragments were calculated for each municipality in ArcGIS, and resulting values were outputted into a spreadsheet and added to obtain the total forest area per municipality. This value was divided by the area of the municipality. Above-Ground Carbon Sequestration (Kilograms of Carbon per Square Kilometer). Above-ground carbon-sequestration values
were calculated from Moderate Resolution Imaging Spectroradiometer (MODIS) satellite data (NASA) representing net primary productivity values (12). MOD17 produces gross primary production of vegetation every day and sums to net primary production at 8-day intervals. The product is computed with daily MODIS landcover, fraction of photosynthetically active radiation/ leaf area index (FPAR/LAI), and global Global Modeling and Assimilation Office (GMAO) surface meteorology at 1 km for the global vegetated land surface. These variables provide the initial calculation for growing season and carbon-cycle analysis, and they are used for agriculture-, range-, and forest-production estimates. Two tiles of 2001 data were ordered from NASA to cover the area of land within the Richelieu and Yamaska watersheds for every 8-day period in 2001. The files were converted to file formats compatible with IDRISI and then windowed to extract the relevant study area. All pixels classified as water or urban (32761 and 32766) were reclassified as 0. Missing data, presumably caused by cloud cover, were a problem during the months of June and July (classified as 32767), and thus, we calculated summer and nonsummer totals separately. June and July maps where many pixels were assigned the value of 32767 (n = 11) were reclassified to correct for missing values. This was done by taking the average value of net primary productivity (NPP) for each pixel over the weeks in June and July, only counting the values that were not 32767. This average value was then multiplied by 11 to find the summer total NPP for those missing 11 8-day periods. Nearest neighbor reclassification was not appropriate in this instance, 1. Statistics Canada (2001) Land in crops (excluding Christmas tree area) by hectare, by province, Census Consolidated Subdivision (CCS), 2001 Quebec (Table 5). In 2001 Census of Agriculture, land use, by province, Census Agricultural Region (CAR), Census Division (CD) and Census Consolidated Subdivision (CCS). Available at http://www. statcan.gc.ca/ca-ra2001/index-eng.htm. Accessed June 5, 2006. 2. Statistics Canada (2001) Pigs by number of pigs, by province, Census Consolidated Subdivision (CCS), 2001 Quebec (Table 20). In 2001 Census of Agriculture, pigs, by province, Census Agricultural Region (CAR), Census Division (CD) and Census Consolidated Subdivision (CCS). Available at http://www.statcan.gc.ca/ca-ra2001/indexeng.htm. Accessed June 5, 2006. 3. Ministry of Sustainable Development, Environment and Parks, Québec (2008) Database on quality of aquatic environment. Available at http://www.mddep.gouv.qc. ca/eau/inter.htm. Accessed January 11, 2008. 4. Statistics Canada (2001) Number of taps on maple trees, by province, Census Consolidated Subdivision (CCS), 2001 Quebec (Table 18). In 2001 Census of Agriculture, maple tree taps, by province, Census Agricultural Region (CAR), Census Division (CD) and Census Consolidated Subdivision (CCS). Available at http://www.statcan.gc.ca/cara2001/index-eng.htm. Accessed June 5, 2006. 5. SoftMap (2001) The Guide to Fishing and Hunting in Quebec (SoftMap Technologies Inc., Quebec City). 6. Ministry of Tourism, Québec (2007) Database on tourism sites of Quebec. Available at http://www.tourisme.gouv.qc.ca/programmes-servidces/services/promation_touristiqu e/donnees_produits_touristiques.html. Accessed July 20, 2007. 7. Ministry of Tourism, Québec (2007) Criteria for admissibility into Quebec tourism database. Available at http://www.tourisme.gouv.qc.ca/programmes-servidces/services/ promation_touristique/donnees_produits_touristiques.html. Accessed July 29, 2007. 8. Database of Natural Patrimony of Quebec (2007) Extraction for the territory of Montérégie, l’Estrie and Centre-du-Québec (Ministry of Sustainable Development, Environment and Parks, Québec).
Raudsepp-Hearne et al. www.pnas.org/cgi/content/short/0907284107
because the missing data pixels were aggregated together in large areas. The total nonsummer NPP was calculated by adding up all of the values for each pixel over the year. The total nonsummer NPP was added to the total summer NPP to get total NPP values for 2001 at a 1 × 1-km resolution. This map was resampled for municipalities using zonal statistics in ArcGIS to obtain an average NPP value per municipality. Crop yields were calculated for each municipality using the 2001 Canadian Agricultural Census (13). The amount of carbon present in the portion of cultures that was removed from fields was calculated using average values of C per crop plant published for Quebec crops (14) multiplied by the percentage of each crop that is typically removed from the field (15). Final amounts of C that were removed from fields per municipality were subtracted from the NPP totals per municipality to get the final value of carbon sequestration. Phosphorus Retention in Soil (Percentage). Soil samples were taken by the provincial government on every farm in the province between 1995 and 2001 (16). A total of 22,984 soil samples were taken within the two study watersheds. Average values for soil retention of phosphorus for each municipality were available from a government online database (16). These values only represent the phosphorus saturation values on agricultural land and not across all land covers. The retention of phosphorus in the soil is indicated by the phosphorus saturation index, which measures the degree to which soil phosphorus sorption sites have been filled (17). For Quebec soils of the type found in these two watersheds, values above 12% phosphorus saturation are considered to be at high risk for runoff into waterways (18). Organic Matter in Soil (Percentage). Soil samples were taken by the provincial government on every farm in the province between 1995 and 2001 (16). A total of 22,984 soil samples were taken within the two study watersheds. Average values for organic matter in soil for each municipality were available from a government online database (16). These values only represent the organic matter in soil on agricultural land and not across all land covers. Organic matter in soil is widely thought to have a critical level of 3.4%, below which the productive capacity of agriculture is compromised by a deterioration in soil physical properties and the impairment of soil nutrient cycling mechanisms (19). 9. Ministry of Natural Resources and Fauna (2007) Database of Natural Patrimony of Quebec. Extraction for the Territory of Montérégie, l’Estrie and Centre-du-Québec (Ministry of Natural Resources and Fauna, Quebec City). 10. Ministry of Municipal Affairs and Regions (2007) Database of Municipal Tax Evaluations (Ministry of Municipal Affairs and Regions, Quebec City). 11. Ministry of Natural Resources and Fauna, Québec (2005) Quebec topographic database 1:20,000. Available at http://www.geogratis.ca/geogratis/en/product/search.do? id-8147. Accessed June 15, 2006. 12. The Land Processes Distributed Active Archive Center (LP DAAC) at the U.S. Geological Survey (USGS) Earth Resources Observation and Science Center (EROS) (2007) MODIS/Terra Net photosynthesis 8-day L4 global 1 km SIN gridMOD17A2. Version V004. January 1, 2001 to December 31, 2001. Available at http://LPDAAC.usgs.gov. Accessed June 15, 2007. 13. Statistics Canada (2001) Total corn, soya, hay, wheat, beans, oats, barley, mixed grains by hectare, by province, Census Consolidated Subdivision (CCS), 2001 Quebec (Table 13). In Census of Agriculture Hay and field crops, by province, Census Agricultural Region (CAR), Census Division (CD) and Census Consolidated Subdivision (CCS). Available at http://www.statcan.gc.ca/ca-ra2001/index-eng.htm. Accessed June 5, 2006. 14. Institut de la Statistique du Québec (2007) Biofood industry—crop production 2001. Available at http://www.stat.gouv.qc.ca/donstat/econm_finnc/filr_bioal/culture/ index_an.htm. Accessed August 1, 2007. 15. Whalen JK, Quancai H, Aiguo L (2003) Compost applications increase water-stable aggregates in conventional and no-tillage systems. Soil Sci Soc Am J 67:1842–1847. 16. Ministry of Agriculture, Fisheries and Food, Québec (2001) Statistical description of the chemical properties of mineral soils of Quebec: Based on soil analyses completed from 1995 to 2001. Available at http://www.agrireseau.qc.ca/references/6/sols/index. htm. Accessed November 3, 2007. 17. Kleinman PJA, Sharpley AN (2002) Estimating soil phosphorus sorption saturation from Mehlich-3 data. Commun Soil Sci Plant Anal 33:1825–1839.
2 of 3
18. Beauchemin S, Simard RR (2000) Phosphorus status of intensively cropped soils of the St. Lawrence Lowlands. Soil Sci Soc Am J 64:659–670.
19. Loveland P, Webb J (2003) Is there a critical level of organic matter in the agricultural soils of temperate regions: A review. Soil Tillage Res 70:1–18.
Fig. S1. The study site consists of the adjacent Richelieu and Yamaska watersheds (dark lines) directly east of Montreal, QC. The territory intersected by these watersheds is divided into 144 municipalities (gray lines).
Fig. S2. Significant correlations between pairs of ecosystem services and data points (*P < 0.05; **P < 0.01; ***P < 0.001). Red signifies negative correlations, and blue signifies positive correlations.
Raudsepp-Hearne et al. www.pnas.org/cgi/content/short/0907284107
3 of 3
