Futures Analysis of Urban Land Use and Wetland Change in Saskatoon, Canada: an Application in Strategic Environmental Assessment

This paper presents a scenario-based approach to strategic environmental assessment (SEA) for wetland trend analysis and land use and land cover (LUC) modeling in an urban environment. The application is focused on the Saskatoon urban environment, a rapidly growing urban municipality in Canada's prairie pothole region. Alternative future LUC was simulated using remote sensing data and city spatial planning documentation using a Markov Chain technique. Two alternatives were developed and compared for LUC change and threats to urban wetland sustainability: a zero alternative that simulated trends in urban development and wetland conservation under a business as usual scenario, in the absence of prescribed planning and zoning actions; and an alternative focused on implementation of current urban development plans, which simulated future LUC to account for prescribed wetland conservation strategies. Results show no improvement in future wetland conditions under the city's planned growth and wetland conservation scenario versus the business as usual scenario. Results also indicate that a blanket wetland conservation strategy for the city may not be sufficient to overcome the historic trend of urban wetland loss; and that spatially distributed conservation rates, based on individual wetland water catchment LUC peculiarities, may be more effective in terms of wetland conservation. The paper also demonstrates the 812 challenges to applied SEA in a rapidly changing urban planning context, where data are often sparse and inconsistent across the urban region, and provides potential solutions through LUC classification and prediction tools to help overcome data limitations to support land use planning decisions for wetland conservation.


Introduction
The majority of wetland alteration is driven by, or associated with, changes to land use and land cover (LUC) resulting from human-induced surface disturbance [1].In recent years, urban growth and regional development have become significant drivers of LUC change, due to the direct and indirect effects of surface disturbances to wetland hydrology, geomorphology, and ecology [2,3].Spatial planning, including urban planning and development zoning, is a key policy instrument designed to direct future land use and development actions [4], but such planning and policy instruments do not always provide for adequate consideration of the potential impacts of urban LUC change to wetlands.Strategic Environmental Assessment (SEA), the assessment and integration of environmental and sustainability considerations in policy, plan, or program (PPP) development and decision making [5,6], has gained considerable momentum in recent years as a proactive and spatially appropriate framework for assessing and shaping land use policies, plans, and development actions [7,8], including urban and regional land use plans [9,10], and managing their impacts on wetlands [11].
The basic premise of SEA is that it helps ensure that PPPs are developed and operationalized in an environmentally sensitive way, and that land uses and development actions are implemented within a sustainability framework.SEA is about understanding the context of a PPP or strategy being developed and assessed, identifying and understanding key trends, and assessing future environmental and sustainable outcomes to help achieve strategic objectives [12].Its application typically focuses on the identification and assessment of trends in environmental baselines, and the analysis of future outcomes or scenarios under different land use or disturbance conditions [13,14], thus identifying desired outcomes and what is required to achieve those outcomes [15].
The development and assessment of current and future land uses via scenario analysis and LUC modeling are foundational to SEA in the context of spatial planning and development [16,17].Such scenario-based approaches, which produce visions of future conditions with and without currently planned development actions or initiatives [18][19][20], are widely promoted as good-practice in the impact assessment literature and recognized as a means to help overcome the challenges associated with predicting future outcomes under uncertain planning conditions [13,16,[21][22][23].However, despite the variety of LUC models applicable to land use analysis that currently exist [24,25], as well as a number of available scenario-based approaches [21][22][23], their combined application in SEA is limited [26], particularly within the context of land use planning and assessment for managing impacts to urban wetlands [27,28].Noble et al. [29] report that the majority of assessment methods and techniques applied in SEA are limited to a number of common, qualitative-based approaches and argue the need for more analytical-based, quantitative assessment methods to support SEA application in regional land use contexts.
This paper presents a scenario-based approach to SEA for wetland area change trend analysis and LUC modeling in an urban environment using a Markov Chain technique.The application is focused on the Saskatoon urban environment, a rapidly growing urban municipality in Canada's prairie pothole region.Specifically, this paper: (i) demonstrates a scenario-based methodology for assessing the impact of a business as usual scenario through simulation of future LUC, based on trend analysis of historical wetland and urban conditions; (ii) assesses the impact of the City's proposed urban growth and development plan; and (iii) undertakes a comparison of futures with respect to the implications for urban wetland conservation.In doing so, the objective is to advance SEA design and quantitative approaches for wetland conservation planning in urban environments.
The following section provides a brief review of the wetland assessment and LUC modeling for futures and trend analysis in SEA application for urban planning and development.This is followed by a description of the study area, data, and methods used for wetland baseline and futures assessment; finally, the study results are presented.The paper concludes with a discussion of the results, challenges, and further research to advance strategic assessment methods for land use planning assessment in an urban wetland context.

SEA Context: Wetlands Assessment and LUC Modeling
The environmental significance of wetlands is well recognized: wetlands provide habitat for a large number of species, serve as flood control areas, and support ecosystems in terms of water quality maintenance, nutrient cycling, and carbon sequestration [1].Despite this, wetlands are decreasing more rapidly than any other type of ecosystem [30,31].Currently, only about 5-8 percent of Earth's land surface is covered by wetlands-about 14 percent of Canada is wetlands [32,33].In the context of the Canadian prairie region, a reported 71 percent of wetlands have been lost, due primarily to agricultural land conversion but also due to road development and urban growth [3,11,31].
According to Rubec and Hanson [34], the Canadian Federal Policy on Wetland Conservation (FPWC) is intended to ensure an appropriate level of wetland conservation, the mitigation of environmental impacts to wetlands, and to sustain wetland functions.The FPWC applies to federal land, all activities on it, and all programs, expenditures, and decisions under federal jurisdictions.However, the majority of wetlands in Canada are on privately owned agricultural land [35], or in and surrounding urban municipalities, which makes both federal and provincial wetland protection policy largely ineffective.
Further, despite numerous policy and regulatory initiatives in Canada that relate to wetlands, there are no standardized or formal methods for the assessment and mitigation of the environmental impacts of development and land use decisions on wetlands [3,34].The tool most often applied to wetland impact assessment and mitigation planning in Canada is project based Environmental Assessment (EA).However, Noble et al. [33], Seitz et al. [36], and Nielsen et al. [3] report that a considerable number of proposed activities, which potentially affect wetlands directly and/or indirectly, including urban growth and expansion, are either not classified as significant enough to trigger EA or are not subject to formal EA regulation.When EA is conducted, assessments are often of the "screening-type", designed for formulaic projects with seemingly predictable impacts, and EA is often conducted perfunctorily, which is both spatially and temporally restrictive [11,33].
The development and application of tools for the assessment of LUC change and potential impacts on wetlands can methodologically improve EA for the sustainable management of wetlands in urban environments.LUC change is broadly defined as a spatiotemporal iteration between biophysical and human related drivers [37].A variety of LUC change models are available for representation of LUC dynamics, its causes, and possible consequences.Schrojenstein Lantman et al. [38] propose four core principles of LUC change modeling, where any given LUC change model is based on at least one principle: Historical trends.This principle is based on extrapolation of past LUC change, assuming permanency of the factors that underline the trend.(ii).
Suitability.The suitability principle describes LUC change based on the specific characteristics of a parcel of land (e.g., spatial, biophysical, or socio-economical preference) in terms of quantitative or qualitative based minimization of costs or maximization of profit.(iii).
Neighborhood interaction.The basis of this principle is that the possibility of land change depends on surrounding land uses and use characteristics.LUC change drivers in this case can be biophysical or socio-economic.(iv).
Actor interaction.In this principle, LUC change depends on the decision-making process, which is the result of actors' interaction.The drivers are socio-economic values and development policy.
Regardless of the principles behind a LUC change model, there are several considerations for any LUC model application [24,25], namely: land use history, temporal dynamics of land use, and the trajectory of land use changes; representative driving factors for a particular LUC change model; land use patterns, their spatial interaction and neighborhood effects; level of analysis and its complexity; scale, including spatial and temporal scale, model granularity, and analytical dimensions.The weight of a consideration given to any set of considerations can vary depending on the particular LUC change model, the goals of the application, and practical constraints.
There are several underlying analytical approaches that are used most often in LUC change model applications, including: cellular automata, statistical analysis, Markov Chain, artificial neural networks, economic-based models, and agent-based model [38].Any one of these analytical approaches can be based on one or more of principles described above and a number of them can be combined in any given LUC model.In the context of wetlands assessment, a Markov Chain approach, which is based on a historical trend principle, is a convenient and accurate model for LUC simulation and has been adopted in several simulation-based studies (see for example Zhang et al. [39], Arsanjani et al. [40], and Ma et al. [28]).However, the benefits of such an approach to assessing scenarios of wetland change in urban development have not been explored within the context of SEA application.

Study Area
The city of Saskatoon is located in the province of Saskatchewan, Canada, on the banks of the South Saskatchewan River.The province of Saskatchewan is part of the prairie pothole region, which covers approximately 480,000 km 2 [1] and contains approximately 11 percent of Canada's wetlands [41].
Saskatoon is the largest city in the province of Saskatchewan with an estimated population of 254,000 as of 30 June 30 2014 [42].From the first permanent settlements in 1883, the city of Saskatoon has grown steadily [43]; 1951-2011 was historically the period of fastest growth, with an average annual rate of population increase of 2.6% [44,45], and with a projected population of 387,742 (with annual growth rate of 2.5%) to 2032 [46].
The study area consists of five development sectors situated in the Suburban Development Areas of the city of Saskatoon.The development sectors are identified in the City of Saskatoon's Official Community Plan Bylaw No. 8769 [47] and with planning documentation (Sector Plans) publicly available [48].However, publicly available Sector Plans do not fully reflect the most current planning strategy for the city.As such, the Sector Plan boundaries for each development area were updated based on discussions with the City of Saskatoon Planning and Development Department, taking into consideration the City's most current development strategy.Next, following Sizo et al. [49], the Sector Plans were then adjusted to the next closest (in terms of spatial extent) water catchments (land area that drains to a common waterbody, e.g., river or wetland)-referred hereinafter as assessment areas (Figure 1).These are considered the smallest geographically and ecologically meaningful regions for a regional level environmental assessment [50].A water catchment was chosen for this assessment as it represents the minimum hydrological unit in terms of the ability of a wetland's system to maintain its functions and stability over the long-term [2,51].The assessment areas contain approximately 1870 ha of wetlands, including about 506 ha of wetlands within the updated Sector Plans' area [49].The future growth of the city will have both direct and indirect impacts on urban lands, including wetland areas.Similar to other urban municipalities across Canada, the impacts of urban development on wetlands under planned urban development activity is not subject to SEA under current federal or provincial laws or regulations [33].

Assessment Approach
Baseline assessment is a key step in any SEA application that serves to establish the regional context for an assessed policy, plan, or program, and consists of an analysis of baseline conditions and changes across space and/or over time.Baseline assessment is arguably the most data and ancillary information intense stage of the SEA process, and is fundamental to the overall analytical quality of SEA and the reliability of SEA results [6,8,[52][53][54][55][56]. Using the information from a baseline assessment, the scenario assessment and evaluation stage of SEA identifies and evaluates the potential impacts of strategic alternatives, or development scenarios, on current baseline conditions.A zero alternative or business as usual scenario describes a future without planned action and serves as a future baseline condition for comparison with potential impacts from planned activities, and is a required component of good SEA practice [6,8,[18][19][20]55,56].Two alternative scenarios were developed and assessed for the study area: (i) a zero alternative, based on a trend analysis of historical wetland and urban change data; and (ii) a current development plan (CDP) scenario, which was based on considering the current development plans and growth strategy for the City of Saskatoon.A discrete-time Markov Chain technique was applied for simulation of future LUC within the study area under each scenario.

Data
The assessment incorporated historic LUC change data and data from a landscape based temporal analysis of wetland change in the Saskatoon urban development region [49].The assessment of LUC change was limited to three classes: (i) wet area; (ii) built-up area; and (iii) non built-up area.Wet area was used as a proxy for wetlands and consists of open water or saturated areas without taking into account morphological properties of the soil [57].Built-up area captured urban development features and associated structures (e.g., housing units, schools) [58].All other areas, which were not identified as wet area or built-up area, were classified as non built-up areas.
Historic LUC change data covered 1985-2011, and three Landsat 5 TM Climate Data Record Surface Reflectance datasets [59] were used as remote sensing (RS) data sources.The selection of RS data was based on image quality, availability, and seasonality.Spring images were preferred, as suggested by Dahl and Watmough [60] for wetland identification in Canada prairie region.The acquisition dates for RS images were: May 11, 1985, April 28, 2006, and May 19, 2011.The post classification comparison method was used for LUC change data extraction in the form of a change detection table (Table 1).RS unsupervised classification method was used for the image classification with a resulting overall accuracy of 92.4% for 1985, 93.6% for 2006, and 92.0% for 2011, and Khat coefficients of 0.884, 0.903, and 0.878, respectively [49].Change detection data was used as a LUC area transition matrix for subsequent Markov Chain models.

Simulation of Future LUC: Markov Chain Technique
The Markov Chain technique is a stochastic model and can be defined as a set of states, St = {St0, St1, St2, …, Stn}, where the modeled process moves from one state to another in a series of steps with a denoted transition probability pij.When a transition matrix P = [pij] is defined, each element of the matrix, pij, shows the probability of land use area change from type i in time period n-1 to type j in time period n, with R total number of land use types: Using the probability matrix in the initial state, it is possible to calculate the state transition probabilities from the initial state to the n th state: The probability matrix was calculated based on the LUC area transition matrix (Table 1), using the following formula: where Aij is the area of land use type that has been in state i in period t-1 and is in state j in period t.
The initial state matrix, defined as St(0), identifies the starting situation-the current state of the baseline environment.The future land use distribution was predicted based on the initial state matrix and transition probability, using a Markov Chain simulation model St(n) In Equation ( 4), n is the relation between desired time period and observation length: where Ld is the desired time period (for a land use prediction) and Lo is observation length.
For computation of the probability matrix P (n) the spectral decomposition approach was applied, where matrix P was represented by its eigenvalues (EigVal) and eigenvectors (EigVect): and respectively: where and λm is the m th eigenvalue and associated with m column of EigVect.The negative eigenvalues when appeared were replaced by zero [61].The initial state matrix St(0), based on LUC data for the assessment areas, was defined as: where A is the area for a respective LUC class.

Markov Chain Models and LUC Alternatives
Two Markov Chain LUC simulation models were developed using 1985 LUC data based St( 0

Quality Assessment of Markov Chain Simulation
The assessment was performed for the built-up class, as wetness condition (and respectively wet area change tendency) changed significantly between 1985 and 2011 [62]

City Planning and Neighborhood Urban Development Evaluation
The City of Saskatoon's Official Community Plan Bylaw No. 8769 identifies the temporal scope of CDP application not by year but with a population threshold of 500,000 [47].As such, to calculate an end date for the CDP for modeling purposes, the City of Saskatoon Population Projection data [46] was used.The projection report covered 2012-2032; the report's population estimations for the years 2012, 2017, 2022, 2027, and 2032 (assuming annual growth rate 2.5%) were extrapolated using a 2nd order polynomial regression to the population level of 500,000, which resulted in the year 2043.The expected scope of planned development over the next 30 years was confirmed by discussions with a senior planner at the Planning and Development Branch, the City of Saskatoon.
To assess future LUC for the CDP alternative, the following steps were performed (Figure 3): first, the updated Sector Plan area (the area of proposed urban development) was removed from the assessment areas and reclassified; second, the MC1985-2006 model was applied to the rest of the assessment area to simulate future LUC to the end date of the expected scope of planned development; and, finally, reclassification results of the Sector Plans and the MC1985-2006 simulation for the remaining assessment areas were combined as the CDP alternative.Please see details below.The planning and development documentation for the City (bylaws and Sector Plans) were prepared in different years and with different levels of detail contained in each plan.Some planning areas (e.g., Blairmore, University Heights, and Holmwood) contained only neighborhood outlines; four neighborhoods (Kingston, Aspen Ridge, Brighton, and Rosewood) had very detailed Sector Plans, containing street details and neighborhood housing/land use block boundaries; the North Sector plan for the city contained the lowest level of planning detail, where only the boundary of planned development with the area was identified (Figure 4).
To compensate for the difference in the level of planning details provided across the Sector Plans, and to include the most current planning strategy in the LUC simulation, LUC within the Sector Plans was reclassified into three classes: wet areas, built-up, and non built-up.The residential, business, and industrial areas were classified as built-up.Water and storm water ponds and constructed wetlands (where specified) were classified as wet areas.Wet areas also included: RS identified wet areas within development sectors where no land alteration had been planned for the next 30 years; and wet areas located in west, north, and north-east swales that had been directly identified for conservation.All other areas were classified as non built-up.The reclassification scheme was reviewed and discussed with the City of Saskatoon Planning Department for confirmation of the most current planning strategy and data available.As a result, a number of suggestions were added to the reclassification results (Figure 5): (i) In the neighborhoods with detailed plans, open space areas (e.g., parks) were identified as non built-up; constructed wetlands, conserved wetlands with original boundaries, and water/storm ponds were identified as wet areas; all other areas were identified as built-up.(ii) Neighborhood areas without detailed plans were classified as built-up with a variable rate of non built-up: (1) 0% and 5% for non residential and (2) 10% for residential use.The neighborhoods without wet areas that were directly prescribed for conservation (e.g., west and north swale) were identified for conservation of 30% of existing wet areas per neighborhood.

Results
The comparison of simulated wet areas with historical data (Figure 6) shows a steady decrease in wet area for most of the study area, and simulated values are below the 5% average wetland area threshold for temperate-zone watersheds, suggested by Mitsch and Gosselink [63].Wetland area is a major indicator for wetland habitat evaluation, as it is reasonable to assume that the loss or degradation of wetland area negatively affects the ability of wetlands in a watershed to maintain their functions within a water catchment, including adequate flood control, water quality maintenance, and support efficient wildlife habitat [60,63].
The CDP and zero alternatives identify and estimate future LUC change in the study area.The wet areas simulated for the zero alternative scenario considered the historical trend of wetland/urban change in the study area, while the simulation for the CDP alternative included the most recent development Sector Plans and the City's current development strategy for the next 30 years.The comparison of wet area simulation between the two alternatives shows that five of the nine assessment areas (i.e., adjusted Sector Plan areas) have a less than 0.1% difference in wet area (# 1, 3, 4, 7, and 8); three assessment areas show a difference between scenario conditions of 0.15% and 0.3% (# 2, 6, and 9); and assessment area # 5 resulted in the most difference between scenarios, with a 0.67% difference in wet area.For the entire study area, the 2043 wet area under the zero alternative is 2.59% and for the CDP alternative is 2.56%.Overall, the 2043 simulation results for wet areas are close, which may indicate that the wetland conservation strategy, identified in the CDP, will have little overall effect on wet area when compared to business as usual under the zero alternative.Part of the challenge is that the City's planning initiatives, including wetland conservation strategies, are bounded by administrative units, which is often the case for urban development [64], as opposed to capturing wetland functional scales [17], which were considered in this analysis and may be considered minimal spatial units in terms of ensuring the ability of wetlands to maintain their long-term functions and resiliency in urban environments [2,51].The Markov Chain technique itself is not an explicit spatial analysis method.However, it is possible to examine the spatial variability of wet area conservation trends among assessment areas.It can be done through an investigation of the spatial distribution of differences between the CDP and the zero alternative with regard to wet areas Markov Chain simulation, despite the close simulation result overall (for the study area).This allows identification of areas where the current development strategy either benefits (assessment areas with positive difference, i.e., # 2, 3, 4, and 9) or does not benefit (with negative difference, i.e., # 1, 5, 6, 7, and 8) from the city's proposed wetland conservation strategy, compared to the historical trend (Figure 7).The variability in net benefits indicates that, in terms of ensuring wetland conservation improvements, a unified strategy as identified under the CDP may not be effective for all urban planning units within city boundaries.Rather, what may be required is a diversity of conservation plans within the study area, where the planning strategy respects overall wetland vulnerability by water catchment, e.g.considering existing wetland/urban patterns, overall wetland area, size, and distribution.

Discussion
Appropriately managing the potential impacts of urban development on wetlands requires an understanding of trends in environmental baselines and some forethought and foresight about future urban growth and resulting LUC conditions.Scenario-based assessment is foundational to SEA [16,[21][22][23]; however, its application, quantitatively, and supporting methods and tools have been relatively limited [29].Consideration of a zero alternative (business as usual scenario) can provide a vision of a future without currently planned development actions, and support the comparison of the impacts of future outcomes with those of planned activities [18][19][20].It could be argued that there is a need for consistent and complete baseline data to support sophisticated future trends analysis in SEA applications.Schrojenstein Lantman et al. [38] noted that consistent data is a precondition for good LUC modeling.Unfortunately, as was found during this research, consistent and comprehensive datasets are rarely available in actual urban planning contexts.Publicly available planning data is often not up-to-date, particularly in fast-growing municipalities such as the one in this study and as such may not actually reflect the most current planning and development strategy.González et al. [65] reported that impact assessment with the use of GIS techniques is widely constrained by the level of detail of available datasets.However, Therivel [6] argued that "not all the baseline data must be available for an SEA to proceed;" João [66] agreed, in that SEA can be completed with missing baseline data and argues that there is a need to reach a balance between data collected and data needed.The current study used development plans proposed to establish a development strategy for the city of Saskatoon for the next 30 years.Each of the neighborhood or land use Sector Plans, however, contained varying levels of planning detail.As such, to facilitate trends analysis and scenario development, it was necessary to supplement and standardize the granularity of data and land use classification based on more qualitative data, based on discussion and recommendations from the representatives of the City's Development and Planning Department.
A Markov Chain, a stochastic method, was then used for LUC simulation.A stochastic approach is recommended against others to account for the uncertainty and variability associated with future conditions (see for example MacDonald [67]).However, as with all futuring methods, the Markov Chain has its limitations [68].The probability distribution of the next period's state depends on the current state and does not depend on the states the chain had previously passed through.That means the variation of the probability of change between two states is not a part of the Markov Chain model and is not reflected in the simulation of future.Also, as the transition probability determination plays an important role in Markov Chain modeling, the accuracy of the input data needs careful consideration.In terms of the current study, the probability matrix was calculated from the area transition matrix, based on RS change analysis of the study region, with an overall accuracy of 92.4% for 1985, 93.6% for 2006, and 92.0% for 2011, and Khat coefficients 0.884, 0.903, and 0.878, respectively [49].If feasible, the quality assessment of the Markov Chain model should be performed for the area of application.For the study area, the MC1985-2011 model was evaluated with a good overall average and median accuracy of 66.3% and 69.1%, respectively.Another assumption of the Markov Chain technique is that the probability law relating the next period's state to the current state does not change.That means the change remains stationary over time after the evaluated states.In the current study, this feature of the Markov Chain was used for the zero alternative of LUC change, i.e., for the historic trend description; while the planned change of the LUC condition, i.e., planned urban development, was used for the CDP alternative.
The assessment of LUC change in this research was limited to the three basic classes: a wet areas class was adopted as a proxy for wetlands and built-up and non built-up classes were used for urban development.The LUC simulation was based on the change of RS identified wet areas, though some amount of dried wetlands could have be classified as non built-up class due to the difference in environmental conditions over time.However, the intention was to show the applicability of the Markov Chain approach for scenario development and analysis in SEA practice.For practitioner's needs, the classification scheme can be adjusted depending on desired information resolution and can include, for example, dried wetlands and agriculture field classes separately.
The zero alternative simulated the tendency of urban development before prescribed urban planning and development actions have taken place, while the CDP alternative simulated future LUC in order to account for prescribed development plans and wetland conservation initiatives.Overall, results for the year 2043 under the CDP and zero LUC simulations were similar.In the context of wetlands management, this suggests that the application of the current wetland conservation strategy, only within administrative based planning units, will not be sufficient to compensate for the historic trend of wetland area decrease in the city.However, the variability analysis of difference between the current development strategy and the historical trend wet areas simulations (the CDP and the zero alternatives) identified areas of positive and negative difference.The areas with positive differences are located at the north-east of the city and might experience improvements in wetland coverage due to conservation initiatives; the rest of the Study Area does not have a wetland benefit in terms of wetland conservation.The variability is most likely a reflection of a difference in landscape wetland/urban spatial patterns that was not fully considered in development planning.Overall, for wetland conservation policy improvement, it is suggested that: (i) any wetland conservation strategy should consider wetland functional scale and ecologically meaningful units, e.g., water catchments, as opposed to solely administrative based planning units; and (ii) a variable conservation rate based on consideration of the particular qualities and variability of the LUC (e.g., existing wetland area and wetland/urban spatial patterns) is preferred to a standard conservation prescription across the urban area.
The scenario-based approach presented in this work provides a basis for wetland trends and future LUC analysis, which, in future practice, may be extended by using a range of alternatives; for example, based on modeled planning documentation with respective desired outcomes that may involve different urban development plans, different patterns of urban growth and density, and alternatives wetland conservation policies or targets.Despite the frequent lack of consistency in input data for undertaking such quantitative-based analysis in urban wetland contexts, the approach presented here utilizes a reliable and replicable method for futures analysis that can be contextualized for wetland conservation.This research may be used to support urban planners and wetland policy makers in the development and comparison of completing planning, development, and land use zoning options to help ensure that urban wetland conservation goals are achieved.

Conclusions
This paper presented an approach for the assessment of urban planning futures in the context of urban wetland conservation.The research was designed as a scenario planning exercise, where two alternatives, a zero alternative and current development plan alternative, described potential wetland futures based on past trends of urban development and wetland change.The study addresses the currently recognized need for contextualized methodology in SEA to support scenario development and LUC modeling [26][27][28], particularly within the context of wetland conservation.The scenario exercise did reveal some shortcomings and difficulties in on-the-ground application, specifically related to assessment contextualization, data quality, data availability, and data consistency.The approach presented in this paper may be useful for SEA practitioners, whose work is related to wetland and urban growth analysis, in environmental assessment and leverage data and information that can be helpful in development plan design, in the comparison and revising of wetland conservation strategies, as well as in providing a sound basis for decision and policy making in urban environmental management.

Figure 1 .
Figure 1.Saskatoon urban environment study area and adjusted sector plan areas (assessment areas).
): (i) MC1985-2006, using a 1985-2006 probability matrix, and (ii) MC1985-2011, using a 1985-2011 probability matrix.The MC1985-2011 model was used for a quality assessment of the Markov Chain technique, while the MC1985-2006 model was used for scenario development.Two alternative scenarios were developed for the future analysis: (i) The zero alternative, which simulates business as usual, where future LUC is based on historical trends in wetland/urban change using the MC1985-2006 model for the entire study area and describing future conditions without the introduction of planned urban development actions.(ii) The current development plans (CDP) alternative, which describes future conditions assuming the City's development plans and strategy are implemented, using the MC1985-2006 model for the area which lies outside the spatial extent of the proposed Sector Plans of urban development.
. The simulated MC1985-2011 model LUC area for 2006 was compared with results of 2006 RS images classification.The Relative Error (RE) was calculated for each assessment area as follows: = | − | (10) where Xo are the original values (the result of the 2006 RS classification) and Xp are modeled values (the results of the MC1985-2011 simulation for 2006).The accuracy assessment for the model yielded good overall average and median accuracy of 66.3% and 69.1%, respectively (Figure2), with one outlier in assessment area #1 (more than two standard deviations).

Figure 2 .
Figure 2. Quality assessment of Markov Chain simulation (built-up class).

Figure 4 .
Figure 4. Sector Plans' LUC, publicly available for the study area (overview, please see [48] for original maps).

Figure 6 .
Figure 6.CDP and zero alternatives wet areas comparison.

Figure 7 .
Figure 7. Difference of wet areas Markov Chain simulation between the CDP and the zero alternatives.

Table 1 .
Land use and land cover (LUC) change data (area transition matrix).