Urban Growth Dynamics in Perth, Western Australia: Using Applied Remote Sensing for Sustainable Future Planning

Earth observation data can provide valuable assessments for monitoring the spatial extent of (un)sustainable urban growth of the world's cities to better inform planning policy in reducing associated economic, social and environmental costs. Western Australia has witnessed rapid economic expansion since the turn of the century founded upon extensive natural resource extraction. Thus, Perth, the state capital of Western Australia, has encountered significant population and urban growth in response to the booming state economy. However, the recent economic slowdown resulted in the largest decrease in natural resource values that Western Australia has ever experienced. Here, we present multi-temporal urban expansion statistics from 1990 to 2015 for Perth, derived from Landsat imagery. Current urban estimates used for future development plans and progress monitoring of infill and density targets are based upon aggregated census data and metrics unrepresentative of actual land cover change, underestimating overall urban area. Earth observation provides a temporally consistent methodology, identifying areal urban area at higher spatial and temporal resolution than current estimates. Our results indicate that the spatial extent of the Perth Metropolitan Region has increased 45% between 1990 and 2015, over 320 km 2. We highlight the applicability of earth observation data in accurately quantifying urban area for sustainable targeted planning practices.


Introduction
Over the last 15 years, Perth has experienced exponential economic growth with Gross State Product (GSP) increasing 218% [1].Originally labelled as the 'Cinderella State' due to its remote location and perceived neglect from the rest of Australia, Western Australia (WA) has experienced sustained discovery and extraction of natural resources since the beginning of the 21st century [2].In response to a growing resource sector, the city of Perth has undergone extensive urban expansion at what Dhakal (2014) identified as an unsustainable rate [3].To this end, the Western Australian Planning Commission (WAPC) identified that Perth's urban footprint has increased from 631 km 2 to 870 km 2 in the 10 years between 2002 and 2012 [4,5].However, these figures should be considered with caution as data used in early estimates represent land parcel (Cadastral) valuations only (provided by the Western Australian Value General's Office), with later estimates (from 2009) based on multiple urban zoning classifications, and more recently (from 2010) spatial modelling taking into account land valuation and zoning [6,7].The use of varied data and methods impacts confidence in the ability of the Commission's estimates to represent actual change in urban extent, especially when urban zoning information includes land identified for growth but not necessarily developed.Such inconsistencies could have potential to misinform future development decisions.Consequently, here we present a spatiotemporal assessment of change in areal urban growth based upon medium resolution remote sensing through a single classification model.This provides the first accurate depiction of urban expansion for one of the world's fastest growing cities-Perth, WA.We present our findings and discuss the implications of more accurately classified urban extents in facilitating scientifically evidence-based adaptive and targeted planning policies to help reduce environmental and socio-economic consequences of poorly planned development.

Earth Observation for Monitoring Urban Change
Mapping the spatial extent and temporal profile of urban growth from medium resolution satellite imagery facilitates a consistent, detailed characterisation of the actual urban footprint of a city [8,9].Other conventional spatial datasets such as Cadastral data provide information on freehold and Crown land parcel boundaries including attributes such as ownership and value for a singular temporal period [10].However, attributed data for a singular year provides an ineffective portrayal of actual parcel land cover and temporal change.Thus, the methods and results presented in this study provide foundational information for the development of planning regulations that ensure sustainable growth of our cities, particularly in the reduction of environmental risks from ever-increasing expansion along the wildland-urban interface [11].Specifically, Earth Observation (EO) data allows spatially detailed identification of locations where (un)sustainable urban growth is occurring which enables expansion limits to be imposed through targeted policies [12].In this theme, Schneider et al. (2005) determined the spatial distribution of development zones from 1978 to 2002 in Chengdu, Sichuan province, China in response to the Go West policy of the 1990s, aimed at economically boosting the West of the country [13].Whilst the policy was successful in raising Gross Domestic Product (GDP) levels, urbanisation concurrently increased, generating issues of urban management, including service, infrastructure and resource deficiency.Their results indicated spatial clustering, specialisation of land use and peri urban development (not considered by the original policy) which were subsequently used to tailor policy in remediating issues, facilitating sustainable future urban development [13,14].Similarly, Hepinstall-Cymerman et al (2013) used classified Landsat data to monitor urban growth in regards to imposed growth boundaries in the Central Puget Sound, Washington, USA [15].Surprisingly, more new development occurred outside the growth boundaries than inside within their last time period, illustrating the ineffectiveness of the imposed policy leading to economic and ecological consequences, including a loss of avian diversity in native forest species [15,16].These studies highlight the potential effectiveness of EO data in consistently monitoring the spatiotemporal dynamics of urban development for applied policy outcomes and ensuring sustainable future planning decisions, for which such outputs are unachievable from traditional datasets.

The Case of Perth
Perth's dramatic urban expansion can be attributed to Australia's minerals and energy boom commencing at the turn of the century.Queensland (QLD) and WA were at the forefront of the boom contributing the largest proportion of the nation's resources output, valued at 3.3% of GDP [1].In WA, mining and petroleum extraction dominate exports, peaking at 95% of the state's export earnings between 2010 and 2011 [17].The increase in extraction was predominantly attributable to greater demand for raw materials from China, resulting in steady growth of the WA mineral and petroleum industry from AUD 4.7 billion in 1996 to a peak of AUD 121.6 billion mid-2013.However, in 2009, a 10.3% reduction in the overall value of mineral and petroleum resources resulted from falling commodity prices and the 2007-2009 global financial crisis [17].Again in 2012, a further 9% reduction in resource value was observed as uncertainty in global economic conditions increased [17].The largest decline to date occurred between 2014 and 2015, with an additional 22% reduction in the value of mineral and petroleum resources as a result of surplus capacity, decreased demand, and decline in the value of the Australian dollar [17].The temporal trend in resource value indicates a stagnation and decline since late 2013 (Figure 1).and decline in the value of the Australian dollar [17].The temporal trend in resource value indicates a stagnation and decline since late 2013 (Figure 1).Perth is described as one of the most isolated cities in the world (pop.> 1 million) and was Australia's fastest growing metropolis between 2007 and 2014; however, subsequent to a decline in natural resource value, a slowdown in population expansion soon followed (Figure 1) [2].As a result, 2015 population statistics highlight the lowest population increase since records began with a 0.5% increase from the previous year [1,2].In comparison to other Australian state capitals, based on the Australian Bureau of Statistics (ABS) 2011 population grid, Perth exhibits a relatively sparse spatial distribution of population with a maximum population density of only 3662 people per square kilometre (Melbourne 10,827; Sydney 14,747).Such low density population has generated high demand for dispersed housing, amenities and services, and has influenced changes to Perth's land use patterns in a non-strategic, "lot-by-lot fashion" based on a car-dependent lifestyle [3].Anthropogenic modifications of the landscape from vegetation cover to human-made impervious surfaces represent a critical driving force in both local and global environmental change [18,19].For example, abrupt, poorly planned and uncontrolled urban expansion can lead to environmental impacts which degrade ecological systems including habitat fragmentation and socio-economic issues that deteriorate efficiency of amenity provisioning, both of which can exacerbate localised climate change [11,20].Identifying impacts of Land Use and Land Cover (LULC) change on socioecological systems is vital for future sustainable urban development; as reflected in the "sustainable cities and communities" 2030 sustainable development goal and the effective land use planning criteria of the City Resilience Framework (CRF) [19,21].It is essential for Perth to adapt current practices of outward suburban expansion to achieve more sustainable urban growth and become citysmart for accommodating the predicted additional half a million new residents by 2031, which will result in an overall population exceeding 2.2 million [5].Perth is described as one of the most isolated cities in the world (pop.> 1 million) and was Australia's fastest growing metropolis between 2007 and 2014; however, subsequent to a decline in natural resource value, a slowdown in population expansion soon followed (Figure 1) [2].As a result, 2015 population statistics highlight the lowest population increase since records began with a 0.5% increase from the previous year [1,2].In comparison to other Australian state capitals, based on the Australian Bureau of Statistics (ABS) 2011 population grid, Perth exhibits a relatively sparse spatial distribution of population with a maximum population density of only 3662 people per square kilometre (Melbourne 10,827; Sydney 14,747).Such low density population has generated high demand for dispersed housing, amenities and services, and has influenced changes to Perth's land use patterns in a non-strategic, "lot-by-lot fashion" based on a car-dependent lifestyle [3].Anthropogenic modifications of the landscape from vegetation cover to human-made impervious surfaces represent a critical driving force in both local and global environmental change [18,19].For example, abrupt, poorly planned and uncontrolled urban expansion can lead to environmental impacts which degrade ecological systems including habitat fragmentation and socio-economic issues that deteriorate efficiency of amenity provisioning, both of which can exacerbate localised climate change [11,20].Identifying impacts of Land Use and Land Cover (LULC) change on socio-ecological systems is vital for future sustainable urban development; as reflected in the "sustainable cities and communities" 2030 sustainable development goal and the effective land use planning criteria of the City Resilience Framework (CRF) [19,21].It is essential for Perth to adapt current practices of outward suburban expansion to achieve more sustainable urban growth and become city-smart for accommodating the predicted additional half a million new residents by 2031, which will result in an overall population exceeding 2.2 million [5].

Data Preprocesing
EO data have been extensively used to monitor the sustainability of urban areas [22,23].However, accurate identification and temporal monitoring of urban land is frequently precluded due to the coarse resolution (300 m-1 km) of a number of commonly used remotely sensed datasets including night time lights (1 km) and the Moderate Resolution Imaging Spectroradiometer (MODIS) land cover product (0.083 • ) [22,24].Whilst 30 m resolution data (e.g., Landsat) are more suitable to detect nuances of urban development the majority of studies and classified products which have used these finer resolution products implement large temporal windows, negating the possibility of detailed temporal urban characterisation e.g., GloeLand30 [25][26][27][28][29].This research provides the first comprehensive temporal evolution analysis quantifying land cover change and associated urban expansion for the Perth Metropolitan Region (PMR) using 30 m Landsat imagery, the longest temporal record of medium spatial resolution imagery, for seven sequential time snapshots between 1990 and 2015.
Cloud free imagery was acquired in or close to the month of July for 1990July for , 2000July for , 2003July for , 2005July for , 2007July for , 2013July for and 2015.Analysis of imagery acquired from WA winter season coincided with peak green-up which provided the greatest contrast between spectrally similar surfaces (e.g., bare earth and urban) [30][31][32].Imagery date selection was founded upon the strong positive relationship between Australian soil moisture (related to rainfall) and the Normalised Difference Vegetation Index (NDVI) [33], which exhibits an approximate one month lag between peak soil moisture and peak NDVI [33].
Productive photosynthesising plants use energy in the visible red (VIS) portion of the electromagnetic spectrum whilst reflecting in the near-infrared (NIR) region.NDVI ((NIR -VIS)/ (NIR + VIS)) is a representative measure of growth allowing for the identification of green, healthy vegetation [33][34][35], as illustrative of Southwest WA's winter months.A total of 14 images from Landsat 5 Thematic Mapper (TM) (eight images), Landsat 7 Enhanced Thematic Mapper Plus (ETM+) (two images) and Landsat 8 Operational Land Imager (OLI) (four images) were acquired for the specified years.Seamless images were produced based on Voroni diagrams that locate the bisector between images; adjacent edges were identified as seamlines constraining effective mosaic polygons that specify inclusion pixels for the final mosaicked product, permitting less visible boundaries through blending overlapping pixels [36].Mosaicked images were subsequently clipped to the original PMR study area boundary.
The atmospherically corrected Landsat data used in this study were obtained from the Landsat Ecosystem Disturbance Adaptive Processing System (LEDAPS) and the Landsat 8 Surface Reflectance (L8SR) algorithm [37,38].Some inherent residual noise remained, for example, due to the differences in modelled atmospheric correction parameters [39].To correct for this, surface reflectance values were standardised as 1 : where p i,b is the standardised pixel value i, from band b based on the original surface reflectance x, standardised through division of a priori specific upper reflectance limit for each band (maxb):  and between (intra) each image (see Figure S1).

Data Classification
The normalised Landsat imagery was classified using the Import Vector Machine (IVM) which builds upon the popular Support Vector Machine (SVM) methodology 2 [41].In order to obtain the optimum classification, the IVM algorithm explores all possible subsets of training data for optimal selection (termed import vectors) which are derived through successively adding training data samples until a given convergence criterion is met [41].Data samples are selected according to their contribution to the classification solution.However, a pure forward system is unable to remove import vectors that become obsolete after addition of other vectors.Therefore the implemented version of IVM utilised here is a hybrid forward/backward strategy that adds import vectors whilst concurrently testing if they can be removed in each step, thus leading to a sparse and more accurate solution [41].Furthermore, the IVM selects data points from the entire distribution resulting in a smoother decision boundary which is based on the optimal separating hyperplane in multidimensional space compared to that of SVM algorithms [42].The benefits of the IVM algorithm have resulted in this approach being successfully applied in a number of studies (e.g., [42][43][44][45]) due to its accuracy and performance advantages over alternative methodologies including SVM and the traditional Maximum Likelihood (ML) classifiers [44,45].
Model training samples were selected using the July 2005 Landsat 5 TM image coinciding with the month post maximum rainfall of all considered Landsat 5 TM and 7 ETM+ to facilitate optimum spectral seperability 3 [33].Land cover was defined as high albedo urban (e.g., concrete), low albedo urban (e.g., asphalt) or other.Two urban classes were initially identified in order to reduce confusion between spectrally similarly classes (e.g., urban and bare earth) being merged post-classification to represent complete urban coverage [46].For each class, 250 pixels were randomly selected as training data, which is consistent with Foody and Mathur (2006) and Pal and Mather (2003) (see Supplementary S2).Training data parametrised the IVM algorithm, creating a classification model of spectral profiles that are compared to Landsat spectral profiles for classification.The classification model was then applied to all Landsat 5 TM and Landsat 7 ETM+ images obtaining similar spectral wavebands, considered to be equivalent [47].However, due to Landsat 8 OLI sampling different spectral regions, a new classification model was developed using the same training areas, as these were deemed to remain representative of the land cover, but with Landsat 8 OLI spectral wavebands [47,48].Validation was performed through an accuracy assessment based on an independent dataset (Google Earth high resolution imagery) consistent with Landsat acquisition months following previously published methods (e.g., [22,[49][50][51][52][53]). For each land use category, 50 pixels per class per year were visually identified and classified based on the majority land cover within the coincident Landsat pixel from Google Earth imagery for the available years: 2000, 2003, 2005, 2007, 2013 and 2015 consistent with recommended land cover accuracy sample size of Congalton (2001) [54].

Results
The spatial footprint of PMR development has increased 45% between 1990 and 2015, over 320 km 2 (Figures 2 and 3), with a 37% increase occurring since 2000.The classification accuracy assessment indicates an average overall accuracy of 84.1% and Kappa Coefficient of 0.73 being comparable to other studies (e.g., [52,[55][56][57]) (see Tables S1 and S2).Urban expansion mirrors population increase and as population growth has slowed, urban development has concurrently exhibited a levelling trend compared to expansion previously observed (Figure 3).LGAs, where land might be developed but less valuable than central LGAs.
For urban estimates post 2005, two urban land zones, urban and urban deferred, are used within the Perth Metropolitan Region Scheme (MRS), the division of the State Planning Policy Framework applicable to the PMR, pursuant to the Planning and Development Act (2005) that inform recent WAPC land parcel based estimates [58,59].Urban land refers to locations where activities in line with 4 Analysed in ArcGIS version 10.2.2 5 Outlines of LGAs are displayed in Figure S2  [58,59].Urban land refers to locations where activities in line with urban development are permitted, but not necessarily constructed (e.g., housing and commercial use) whilst urban deferred represents land suitable for future development with remaining planning, servicing or environmental issues [59,60].For land to be assigned urban deferred, it must obtain characteristics of the urban zone including being able to provide essential services, a logical progression of development, and able to satisfy regional requirements (e.g., roads and open spaces).The 2012 WAPC estimates were derived from stock of land zoned urban or urban deferred, cadastral land plot and value information, conditional subdivision approvals, and ongoing regional rezoning and subdivisions [61].Similarly to 2000, valuation data may misrepresent suburban urban land cover resulting in overestimation.Inclusion of additional variables that are unrepresentative of actual land cover change (e.g., rezoning and conditional approvals) could exacerbate differences between WAPC and EO derived urban estimates (Figure 4b), leading to the potential confounding of errors in WAPC estimates.urban development are permitted, but not necessarily constructed (e.g., housing and commercial use) whilst urban deferred represents land suitable for future development with remaining planning, servicing or environmental issues [59,60].For land to be assigned urban deferred, it must obtain characteristics of the urban zone including being able to provide essential services, a logical progression of development, and able to satisfy regional requirements (e.g., roads and open spaces).The 2012 WAPC estimates were derived from stock of land zoned urban or urban deferred, cadastral land plot and value information, conditional subdivision approvals, and ongoing regional rezoning and subdivisions [61].Similarly to 2000, valuation data may misrepresent suburban urban land cover resulting in overestimation.Inclusion of additional variables that are unrepresentative of actual land cover change (e.g., rezoning and conditional approvals) could exacerbate differences between WAPC and EO derived urban estimates (Figure 4b), leading to the potential confounding of errors in WAPC estimates.

Discussion
WA state government planning documentation states that the majority of new development within the PMR has occurred as low-density suburban growth, responding to consumer preferences and market forces [4].Additionally, sustainable policy objectives suggest that new development should be managed and focused on current communities, making the most efficient use of existing urban areas [4,5].Planning policy research has highlighted issues of outward urban expansion as being costly in economic, environmental and social terms based on dispersed service requirements, habitat fragmentation and neighbourhood segregation [20].Thus, urban expansion in the PMR may result in further economic, social and environmental costs associated with servicing and maintaining low-density lifestyles, owing to the rapid outward urban growth estimates between 2000 and 2007 [11,20].
In contrast, the witnessed slowdown of urban growth, population and natural resource value since 2013 indicates the possibility that the 'boom' of previous years has reached a turning point.Stagnation of urban growth implies that issues associated with spatially distributed urban areas might be contained to the current urban extent.Nevertheless, it is conceivable that prosperous future economic circumstances could initiate growth at a rate previously observed, and that the economic slowdown might be a temporary hiatus responding to current economic conditions [62].For example, in 2014-2015, WA continued to attract the largest proportion of state mineral exploration expenditure at 58%, with QLD (the second ranked stated) obtaining only 20% [17].Furthermore, as of September 2015, WA had an estimated AUD 171 billion in mineral and petroleum projects under construction, with a further AUD 110 billion allocated for future expansion [17].Comparatively, during the peak (mid-2013) in terms of total sales, WA only had an estimated AUD 160 billion worth of projects under construction and a further AUD 108 billion for future development [17].Whilst 2014-2015 observed the greatest decline in total sales of resources, sustained investment and improved global economics could reinvigorate the industry and reinitiate urban expansion within the PMR.
Future development (urban and urban deferred) is guided by Directions 2031 amending the MRS and local planning schemes [5,[63][64][65][66].WAPC aims to achieve 47% of future development as infill and a 50% increase in average residential density by 2050 of 10 dwellings per urban zoned hectare and 15 per new urban zoned hectare [5].In monitoring progress towards the infill target, zoned development land within the PMR is considered, including residential, industrial and commercial land uses [4].Densities are defined as infill or Greenfield if above or below an undocumented residential threshold from census data [60].Initial results from delivering directions 2031, 2014 report indicate the requirement of a significant increase in infill development if the above targets are to be met [67].Similarly, average residential density monitoring has been achieved with land valuation data (from the Valuer General's Office) for major activity centres, being unrepresentative of actual density change and providing an incomplete metropolitan comparison [67].Inclusion of EO data would permit quantitative evidence of urban expansion, infill and density at a higher spatial and temporal frequency than current census based estimates.This would facilitate credible, evidence-based efficient targeted action founded upon improved representative urban area, insuring infill and density attainment.In this theme, Schneider et al. (2005) and Hepinstall-Cymerman et al. (2013) used spatial metrics (e.g., urban area mean patch size) based on classified Landsat data in either pre-defined census units [15] or development corridors [13] to monitor development type (infill or expansion) over time for adapting inappropriate static urban development policy.Using EO derived land cover data in this manner aids in understanding dynamics of the urban environment through monitoring, planning and mitigating land use changes that impact natural assets and increase vulnerability of city systems [14,15,68].Information of this sort aligns with criteria of the CRF in improving city resilience from effective land use planning, possible at lower expense and higher temporal frequencies than in situ measurements [21].
The universal methodology implemented within this research lends itself to credibly inform policy in a similar manner in other global cities through monitoring urban expansion in order to identify rapid, unsustainable development.For example, Jakarta, obtaining the world's second largest urban area with a population of 28 million, has yet to have any quantitative urban area delineation [69,70].Identification of actual urban growth in developing cities is vital to city planners, environment managers and policy makers due to the difference between planned growth and actual growth [14].Such information could be of critical importance for regulating urban expansion due to extreme poverty and high level of risk to environmental hazards, such as that posed from flooding [71,72].EO data presents many opportunities for added value within urban planning policy, and additional analyses could be pursued into specific human-induced environmental issues, such as detecting thermal changes in the urban environment for planning issues associated with urban heat islands (e.g., cooling provisions) and their impact on human health (e.g., air quality).

Recommendations
Consistent and accurate LULC estimates are a vital aspect of sustainable urban development throughout the world, especially considering the predicted additional 2.5 billion city dwellers by 2050.LULC models that require agents that are representative of land use decisions can often fail in practicality due to the difficulty in quantifying driving forces of change and multi-level relationships.Models of this nature are also temporally independent, with each annual iteration implementing new data or data not representative of actual LULC change.EO data provides a replicable detailed representation of the complete urban extent requiring no additional data.The use and application of EO data reported within this paper highlights several improvements to WAPC policy for consistent urban area estimations with associated accuracy measures.Therefore it is recommended that planning authorities, such as WAPC, integrate EO data to achieve the following: (1) provide scientific urban estimates based on a temporally consistent model within future regional structure plans, metropolitan region and local planning schemes; (2) monitor infill development at a higher temporal frequency than census years for policy targeting to meet key goals; (3) monitor urban density through areal urban expansion compared to current metrics using land valuations; and (4) restrict development based on temporal urban analysis that degrades amenity efficiency and ecological systems whilst promoting development in locations to maximise efficiency and long-term sustainability.Additional EO datasets (e.g., finer resolution Sentinel 2 satellite imagery or aerial imagery) could be used to refine planning decisions based on areas of concern identified from Landsat.For example, finer spatial resolution datasets could facilitate enhanced feature extraction, optimising sustainable planning decisions through the identification of candidate infill sites.EO data of this nature provides an essential tool for timely planning policy that is adaptive to changes in urban landscape to mitigate socio-environmental issues associated with poorly planned urban areas for the future sustainability of our cities.

Supplementary Materials:
The following are available online at www.mdpi.com/2073-445X/6/1/9/s1.Additional methodological detail, full accuracy results and Local Government Area (LGA) outlines are reported in the supplementary documentation.

Figure 1 .
Figure 1.Timeline of natural resource value (based on Department of Mines and Petroleum annual reports) fitted with a fourth order polynomial trend line and population (based on Australian Bureau of Statistics data) also indicating key milestones.

Figure 1 .
Figure 1.Timeline of natural resource value (based on Department of Mines and Petroleum annual reports) fitted with a fourth order polynomial trend line and population (based on Australian Bureau of Statistics data) also indicating key milestones.

1
Using the Interactive Data Language (IDL) version 8.3 where ∑ i p i,b is the sum of each standardised pixel across all bands [40].Normalised Landsat data obtained a statistically significant reduction of spectral variation per land cover class within (inter)

2
Using the open source Environmental Mapping and Analysis Program version 2.1.1 (EnMAP) 3 Achieved in the ENvironment for Visualizing Images software version 5.2 (ENVI)

Figure 2 .
Figure 2. Urban expansion within the Perth Metropolitan Region (PMR) between 1990 and 2015.Vast urban growth has been observed in PMR with graduating colours exhibiting outward expansion (a); (b) and (c) exhibit static snapshots of urban extent from 2000 (b) and 2015 (c); whilst (d) depicts percentage of urban change per subnational administrative boundary (Local Government Area;LGA).

Figure 2 .
Figure 2. Urban expansion within the Perth Metropolitan Region (PMR) between 1990 and 2015.Vast urban growth has been observed in PMR with graduating colours exhibiting outward expansion (a); (b) and (c) exhibit static snapshots of urban extent from 2000 (b) and 2015 (c); whilst (d) depicts percentage of urban change per subnational administrative boundary (Local Government Area; LGA).

Figure 3 .
Figure 3.Time line of urban expansion in kilometers squared derived from Earth observation data with associated classification error derived from validation data (points indicating classified image years).Alongside population data in millions per year since 1988 (based on Australian Bureau of Statistics data, 2015 data is projected) with key natural resource milestones indicated, and average annual urban and average annual population growth rate indicated between classified image years.WAPC's urban estimates of the PMR from Directions 2031 (the strategic plan for the Perth and Peel region) were provided for comparison to those produced within this study 4 [5].WAPC's estimates note an expansion from 637 km 2 to 813 km 2 between 2001 and 2012.Our results indicate an expansion of 747 km 2 to 1050 km 2 from 2000 to 2013 illustrating an overall underestimation by WAPC figures (Figure 4).Within suburban areas surrounding the two major cities in the metropolitan region, Perth and Fremantle, WAPC's estimates underrepresent the amount of urban area derived from EO, being more pronounced in 2013 than 2000.The Local Government Area (LGA) 5 of Stirling South Eastern represented the maximum overestimation in 2013 urban area with 34% (2000: 10%) additional urban area per km 2 of LGA established on a difference of 2.89 km 2 , 40.2% (2000: 0.83 km 2 , 15%) between EO data and WAPC's estimates.Outer Northern and Southern LGA WAPC urban values were consistently underestimated, with the LGA of Belmont representing the maximum underestimation of percent per km 2 of LGA in 2013 with 24% (2000: 13%) due to a difference of 9.37 km 2 , 40.39% (2000: 5 km 2 , 26.46%).Prior to 2009, WAPC's estimates were solely based upon land parcel valuations from the Western Australian Value General's Office, consequently valuation thresholds designating land to urban may have been inappropriately applied to outer suburbanLGAs, where land might be developed but less valuable than central LGAs.For urban estimates post 2005, two urban land zones, urban and urban deferred, are used within the Perth Metropolitan Region Scheme (MRS), the division of the State Planning Policy Framework applicable to the PMR, pursuant to the Planning and Development Act (2005) that inform recent WAPC land parcel based estimates[58,59].Urban land refers to locations where activities in line with

Figure 3 .
Figure 3.Time line of urban expansion in kilometers squared derived from Earth observation data with associated classification error derived from validation data (points indicating classified image years).Alongside population data in millions per year since 1988 (based on Australian Bureau of Statistics data, 2015 data is projected) with key natural resource milestones indicated, and average annual urban and average annual population growth rate indicated between classified image years.

Figure 4 .
Figure 4. Percentage differences relative to local government area size, permitting a change metric standardised by Local Government Area (LGA) area between Earth Observation (EO) and the Western Australian Planning Commission's (WAPC) urban estimates for: (a) 2000 (EO) and 2001 (WAPC) and (b) 2012 (WAPC) and 2013 (EO), whilst (c) depicts the percentage difference in the relative urban rate of change (km 2 per LGA area) between 2000 and 2013 (EO) and 2001 and 2012 (WAPC).Positive values indicate underestimation by WAPC whilst negative values represent overestimation by WAPC.

Figure 4 .
Figure 4. Percentage differences relative to local government area size, permitting a change metric standardised by Local Government Area (LGA) area between Earth Observation (EO) and the Western Australian Planning Commission's (WAPC) urban estimates for: (a) 2000 (EO) and 2001 (WAPC) and (b) 2012 (WAPC) and 2013 (EO), whilst (c) depicts the percentage difference in the relative urban rate of change (km 2 per LGA area) between 2000 and 2013 (EO) and 2001 and 2012 (WAPC).Positive values indicate underestimation by WAPC whilst negative values represent overestimation by WAPC.
Figure S1 Inter year classification reflectance variation categorised by classified output for each spectral band for: pre (a) and post (b) normalisation correction.
Table S1 Classification accuracy and associated Kappa Coefficient per year of classified Landsat.Table S2 Producer's and User's accuracy per year of classified Landsat imagery.Figure S2 Local Government Areas (LGAs) located in Perth Metropolitan Region (a); with (b) exhibiting LGAs South and West of the Swan River and (c) LGAs North and East of the Swan River.The classified data reported in this paper (doi:10.1594/PANGAEA.871017) are archived at https://doi.pangaea.de/10.1594/PANGAEA.871017, the pangea open access data publisher for earth and environmental science.