Bioclimatic Extremes Drive Forest Mortality in Southwest, Western Australia

Extreme and persistent reductions in annual precipitation and an increase in the mean diurnal temperature range have resulted in patch scale forest mortality following the summer of 2010–2011 within the Forest study area near Perth, Western Australia. The impacts of 20 bioclimatic indicators derived from temperature, precipitation and of actual and potential evapotranspiration are quantified. We found that spatially aggregated seasonal climatologies across the study area show 2011 with an annual mean of 17.7 °C (± 5.3 °C) was 1.1 °C warmer than the mean over recent decades (1981–2011,- 16.6 °C ± 4.6 °C) and the mean has been increasing over the last decade. Compared to the same period, 2010–2011 summer maximum temperatures were 1.4 °C (31.6 °C ± 2.0 °C) higher and the annual mean diurnal temperature range (Tmax−Tmin) was 1.6 °C higher (14.7 °C ± 0.5 °C). In 2009, the year before the forest mortality began, annual precipitation across the study area was 69% less (301 mm ± 38 mm) than the mean of 1981–2010 (907 mm ± 69 mm). Using Systeme Pour l'Observation de la Terre mission 5 (SPOT-5) satellite imagery captured after the summer of 2010–2011 we map a broad scale forest mortality event across the Forested study area. This satellite-climatology based methodology provides a means of monitoring and mapping similar forest mortality events- a critical contribution to our understanding the dynamical bioclimatic drivers of forest mortality events.


Introduction
Regional drought and heat stress over the past 30 years has resulted in climate-driven forest and woodland mortality events globally [1].Forests and woodlands that cannot relocate to suitable climatic conditions face a reduction of range, fragmentation, mortality and in extreme cases extinction [2].The warming and drying of Mediterranean biomes increases their vulnerability [3][4][5][6] by reducing water availability and increasing vegetation stress and the frequency of mortality [7].Growing evidence suggests that as climates dry and warm the risks of biodiversity loss [8] and severe wildfires increase [9].Such dramatic landscape scale transformations are not likely to occur rapidly, rather evolve as a response to abiotic [10] and biotic changes and the signs of change are already beginning to emerge [1,11].As with other Mediterranean climates, the vegetation of southwest Western Australia (SWWA) is threatened by these changes.Declines of this kind are not new to the SWWA, for example the declines of the tuart forest (Eucalyptus gomphocephala) on the Swan Coastal Plain (SCP) (refer to Figure 1) in the late 1990s and 2000s [12][13][14].Over the last decade, severe drought has impacted on the climate of southern [15] and southwest Australia resulting in marked declines in forest and woodland condition [16,17].Since the 1970's the southwest has received 15-20% less than average rainfall [17] and a 0.4 °C degree increase in temperature [18] leading to the marked depletion of stream flow through the forested region surrounding Perth [10,17] known as the Northern Jarrah (Eucalyptus marginata) Forest (NJF).For example, Petrone [17] shows that a reduction in the average streamflow in the NJF region (calculated as pre-1989 and post-2008 changes) range between −39% and −78% across the catchments monitoring sites.Petrone [17] also shows that, in the region surrounding the NJF, potential evaporation often exceeds precipitation and only a small portion (3-20%) of runoff is generated directly from rainfall and groundwater levels have been consistently falling [19].Experimental results [10] suggest that reductions in forest cover would result in increased streamflow and silvicultural forest canopy thinning has been used as a solution for increasing streamflow into drinking water catchments in the NJF.
In the context of forest distribution and health, bioclimatic limits (or thresholds) define how the forest is affected by its climate.A bioclimatic limit of the distribution of the Jarrah forest was first proposed by Gentilli [20].At the time, the SCP received around 300 mm less annual rainfall (850 mm) than the Jarrah forest (1100 mm) which includes the NJF.Gentilli [20] estimated that this restricted the Jarrah forests westward expansion towards the SCP on the basis of precipitation and evaporative demand: he defined the forest-sustaining bioclimatic threshold to be 800 mm of rain per annum and showed that a difference from the mean of 165 mm rainfall per annum would be significant to drive changes in the forests distribution across SWWA.
A severe inter-annual change in precipitation of the order of that described by Gentilli [20] came to the region between 2009 and 2011.Specifically, nil rainfall was recorded for 122 days at Perth Airport between 21 November 2009 and 22 March 2010 [21].Additionally, the state area average temperature for summer 2009-2010 was equal highest on record, 1.5 °C above average at 36.7 °C [21].Autumn (March, April and May 2010) temperatures were the sixth warmest on record at 1.4 °C above the long term of state area average of 31.1 °C [22].Furthermore, winter 2010 (June, July and August) was the driest and third coldest on record in the southwest [23].The record winter cold and dry spell was followed by the fifth warmest average summer on record, spatially estimated to be +1.09°C above the states long term average [23].
Recently, the Bureau of Meteorology (2012) noted that the prolonged dry spell through summer 2010-2011 contributed to a series of out-of-control forest fires, in the presence of an abundance of dryfuel, resulted in over 100 homes being destroyed or damaged.Extreme fire events such as these are not unprecedented yet are uncommon and formative in that they illustrate the threat posed by dry vegetation and hot conditions.It is critical to record and understand the spatial and temporal dynamics of these events [24].It is expected that Mediterranean ecosystems will experience more frequent and prolonged climatic variability [25] and the need to understand the spatial extent and consequences of these impacts is central to the success of our adaptation strategies in the long term.
In regions where clear evidence of bioclimatic change is taking place, it is critical to develop suitable methodologies to record the extent and severity of forest and woodland mortality events and to quantify the bioclimatic-dynamics in order to understand the difference between natural variability and change.The only tool for achieving near real time, contiguous landscape and continental scale monitoring is remote sensing.Detection of eucalypt canopy dieback using high spatial resolution remote sensing imagery has been shown by a large number of authors and typically involves using spectral bands focused on chlorophyll pigments of the foliage, leaf structural attributes and defoliation [26][27][28][29][30][31][32].Moderate resolution satellite remote sensing (e.g., 250 m-5 km) has been shown to be useful in detecting abrupt changes in greenness and cover [33][34][35] by extracting a noiseless trend from the raw data once seasonal and high frequency changes are removed [36].Using remote sensing to map in combination with gridded and modeled climatologies [8,37] can provide insight into the temporal and spatial dynamics of forested ecosystems.
The objectives of this study were, firstly to use remote sensing technology to map the extent of "patches of mortality" (POM) resulting from the extreme weather of 2010-2011.Secondly, once spatially mapped, we analyze the bioclimatic conditions preceding the event and identify the underlying drivers of the forest mortality.We conclude with a number of potential future research directions on employing this methodology for understanding the responses of the forest stands to changing climatic conditions.

Overview
Following a severe mortality event (July 2011) the methodology for this manuscript was designed to map and investigate the drivers behind the declines.This involved establishing twenty field sites located across the forest, which a team of forest pathologists established and began monitoring immediately [38].Simultaneously, Système Pour l'Observation de la Terre mission 5 (SPOT-5) Satellite imagery was acquired from PlanetAction.org and a feature extraction model, based on known biophysical relationships of remotely sensed imagery, was developed using the field sites with the aim of identifying POM across the imagery.An independent model evaluation dataset was used on the POM identified by the model.Bioclimatic space and time series analysis of Australian Water Availability Data (AWAP) [39] and Moderate Resolution Imaging Spectroradiometer (MODIS) Fractional Cover were used to evaluate the environmental envelope of the study area.Spatially and temporally averaged statistics are evaluated over the period 1981-2011, hereafter the 30 year study period.

Study Area
The NJF directly east of Perth, Western Australia (Figure 1) has a Mediterranean climate, hot, dry summers and wet winters and has a history of tree declines [40].The NJF is classified as a medium open forest of average height between 10-30 m with canopy cover between 30-70% that lies atop the Darling Plateau on a substrate of archaen granite and metamorphic rocks with an average elevation of 300 m above sea level [41,42].Rainfall varies from around 1000 mm on the western side of the scarp to 400 mm to the northern and eastern extremes [42].The region serves as the main surface water supply for Perth as the Swan, Helena, Canning, Serpentine and Murray rivers flow off the Darling Plateau onto the SCP below [17].
The soils vary from lateric gravels to clay and the dominant plant species, Eucalyptus marginata (Jarrah) are distributed in a mosaic across the western side and a mix of species dominate the eastern area with Jarrah co-occurring with Corymbia calophylla.The mid-storey consists of a mix of Banksia grandis, Allocasuarina fraseriana, Persoonia longifolia, Xanthorrhoea preissii and the shrub layer is comprised of a variety of perennial woody species.Across the region this topographical variability, granite substrate and gravel soils serve for rapid changes in floristic distribution from the winter water logged valley floors (drains) to the rocky granite outcrops (uplands).
The NJF is extensively mined for bauxite used in aluminum production, coal and other resources.Over the last three decades considerable emphasis has been placed on restoration of dormant mining sites throughout the forest.A post-mining re-vegetation prescription of 2,500-10,000 stems/ ha is used and designed to produce 350 stems/ha after 12 years of rehabilitation under the caveat that the sites are managed and thinned to this prescription [43].Land use across the forest is thus a mosaic of mining, nature reserves, commercial forestry with areas set aside for Perth's water catchment.

Field Data for Model Development
In July 2011, patches containing mortality were located within the NJF by a team of forest pathologists on a number of field visits.The patches used in this study varied in size from 3,720 to 168,152 m 2 and a complete inventory of the work can be found in Matusick et al. [38].The objective of the field work was not to identify all patches in a region, rather identify a number of sites for detailed investigation of the causes [38] and ultimately to use these for detecting and modeling the landscape scale of the POM in this study.Generally, the sites identified for use in this study were easily identifiable from high-ground and from the road whilst travelling through the forest.Twenty of these plots were located using differential Global Positioning System (dGPS) representing a range of dieback classes.Patches of healthy and dead trees were then delineated using dGPS by trained forest health professionals and a range of indices (not used in this study) were used to determine both the healthy and dead class.

Airborne Imagery for Validation of the Model
To obtain additional information on the extent of the tree mortality across the forest Digital Multi-Spectral Imagery (DMSI) had been acquired by SpecTerra Services (Perth, Western Australia) over small sections of the NJF as part of ongoing research into the hydrology of the forest (Batini personal communications 2011) on the 18 th of November 2010.This imagery provided a clear indication of foliage cover and a means to compare geographical features in the imagery to field made records of the presence and absence of decline since the imagery was taken when the forest did not contain any POM.A 1 km × 1 km grid (refer to Figure 2(b).)was overlaid the imagery and a forest health assessment expert classified where POM were known to have occurred intersecting each grid point as either healthy or containing mortality following the summer of 2010-2011.[44] was used for this exercise.This resulted in a validation area covering 135 km 2 of the NJF.These pixels were then aggregated to produce a confusion matrix of pixels containing greater than 60% mortality and those containing less than 40% mortality, therefore considered healthy.

SPOT-5 Satellite Imagery
SPOT-5 satellite imagery was acquired across a large section of the forest on 29 th June, 2011 and 14 th July 2011 before the first winter rains in order to determine the extent of vegetation declines across the region.The SPOT-5 imagery used was 10 m spatial resolution with four spectral bands green ( 500 to 590 nm), red ( 610 to 680 nm), near-infrared ( 790 to 890 nm) and shortwave-infrared (1,580 to 1,750 nm).The four scenes used to cover the region came with radiometric correction of viewing angle and sensor related distortions in preprocessing level 1A [45].Imagery was orthorectified and converted to surface reflectance.

Feature Extraction Model Development
Since, the 10 m SPOT-5 imagery was insufficient to see individual crowns an object based 'feature extraction' approach was developed to identify POM across the study area located within NJF [46].Feature extraction was selected because the sites varied in size, structure and landscape position.Feature extraction involves first delineating areas of similar spectral and textural likeness by manipulating only spatial parameters, thus creating a mosaic of landscape features of similar spectral and textural likeness.The second step was to extract the features of interest (i.e., the POM) by classifying them by their spatial, textural and radiometric properties, however, given the distinct differences between healthy, green forest and the POM, only radiometric and spatial properties were needed.
Three simple rules were identified by the object based image analysis systems to distinguish POM across the otherwise healthy forest; (1) the mean of ratio of the near infra-red and red bands (i.e., mean patch greenness) of the imagery being less than 0.83 (unit less) was identified as a suitable threshold in a histogram analysis of agreement with the field data.(2) The POM size being between 5,000 m 2 , the Food and Agriculture Organization (FAO) [47] definition of a minimum forest patch size and the upper size limit of 325,000 m 2 being twice the size of the limited area field data and representative of POM identified across the entire study area.(3) The mean of the green band was zero, indicative of non-photosynthetic vegetation and the absence of greenness, thus complementing rule 1.
Using these three rules, the model was visually calibrated on the 20 healthy and 20 mortality field sites identified and resulted in a total of n = 1,853 POM across the study area.The accuracy of the object based classification was then assessed by comparing a 135 km 2 forest health assessment grid of 1 km 2 cells with the classified objects in a standard confusion matrix allowing overall, user and producers' accuracy to be calculated.

MODIS Fractional Cover Dataset
We use the MODIS-based vegetation fractional (MODFC) cover dataset developed by Guerschman [33] to determine changes in fractional cover between 2000 and 2012 in order to establish the natural variability experienced by the forest against which we compare the impact of the 2010-2011 event.MODFC was selected because it has been used for national monitoring of ground cover [48] and is currently the only estimate of photosynthetic vegetation (f pv ), non-photosyntheic vegetation (f npv ) and bare soil (f bs ) data validated for Australia.We address the issue of spatial and temporal disparity between remotely sensed greenness (SPOT-5 at 10 m), cover observations (MODFC 0.5 km) and climatologies by aggregating the fractional cover from its native resolution down to that of the climatologies used in our bioclimatic analyses and specifically for comparison across scale.In order to determine significant trends versus signal noise we use the Seasonal Trend Decomposition using Loess (STL) package in the R-Project framework after Cleveland et al. [36].
The STL method provides a loess smoothed trend by calculating and extracting a seasonal cycle and high frequency variability and subtracting these components from the raw f pv , f npv and f bs signal.The resulting trend can then be used to reliably determine significant trends (i.e., positive or negative departures) from the time series mean (M) and standard deviation (SD).We define a significantly positive trend as being when the loess trend is greater than the sum of M and SD, for example, a significant growth trend is said to occur when f pv is greater than the sum of the M and SD (MSD) and a significant decline in cover when f pv is less than MSD.Additionally, we calculate a rolling M and SD for the raw data using a moving window of three months.This provides an estimate of the local variance of the SD of the raw data relative to the STL trend, useful for visually assessing the temporal variance in both the raw and trend signals simultaneously.

Gridded Climatologies and Bioclimate Analysis
The AWAP dataset [39] was used because it is recognised as an accurate resource of the state and trends of the terrestrial water balance of the Australian continent at a spatial resolution of 0.05 degrees (nominally 5 km) [49].Monthly climatologies of minimum and maximum temperature, precipitation, actual (total) evaporation and potential evaporation were spatially extracted from AWAP for the study period.Spatially explicit climatologies of the study area and the NJF were then aggregated for the analysis across seasonal and inter-annual timescales necessary to understand the climatic envelope of the forested study area.
In the NJF, subsurface through-flows are coupled to soil moisture recharge [10] and there is an annual moisture deficit [17] consistent with a Mediterranean climate [50].The net result is that wet season (winter) precipitation that is not intercepted or stored is lost to runoff and stream flow.In addressing this issue for similar climates, Prentice et al. [50] proposed a scalar, α renamed to the Mean Annual Aridity Index (MAAI) for this study, designed to capture the seasonal variability of soil moisture independently of vegetation interception and considered it indicative of the growth-limiting drought stress on vegetation.
We calculate the MAAI as the ratio of actual (total) evaporation (E a ) to potential evaporation (E p ) derived from the AWAP climatologies [39] which models evapotranspiration approximate to the methods of Priestley and Taylor [51].We assume that actual evapotranspiration (E a ) can never be greater than potential evapotranspiration (E p ) and that E a constrains both soil evaporation and transpiration from the vegetation.We calculate MAAI for each NJF pixel and compare mean annual precipitation (annual totals for each year of the study period) and summer diurnal temperature ranges (i.e., for December, January, February) for each summer over the study period.We classify our MAAI in the same way as Middleton's [52] aridity ratio of precipitation to potential evaporation for the purposes of comparison (Table 1).Detailed or gridded meteorological observations of surface energy fluxes do not currently exist across the forest.As such, the best available alternative are to surrogate bioclimatic indices with known biological significance [53] derived using the best available data, in this case-AWAP [39].Using the AWAP climatologies [39] 19 bioclimatic indices (Table 2) were derived using the R-Project species distribution modeling package 'dismo' [54] and combined with MAAI to make a total of 20 indices.The bioclimatic indices summarise the total energy and water inputs across the NJF at seasonal and annual timescales.They characterise the climatic envelope and identify any thresholds contributing to the POM occurrences following the summer of 2010-2011.The aggregate of monthly inputs into the seasonal variables enable the detection of inter-annual shifts in seasonality; BIO15 relates to shifts in seasonal precipitation occurrence, known to be directly related to vegetation drought stress [55] and BIO14 focuses specifically to precipitation in the dry season, in the NJF this is in the summer months (December, January and February).
The mean summer diurnal temperature range, and annual precipitation, used in the climate space analysis, provides insight into the range and magnitude of the energy inputs of temperature and specifically on precipitation preceding the mortality event.The approach assisted in identifying the key variables in the bioclimatic space, preceding and thus contributing to the mortality event.This deterministic methodology is consistent with others which use the response of plant functional traits to model their distribution [56,57] and persistence under variable climate scenarios [58,59].

Detection of POM and the Accuracy Assessment
The spatial extent of the study area was 6074.5 ha (60 km 2 ).The expert assessment of the 135 km 2 grid resulted in greater accuracy in the detection of POM, 85% hit versus 74% missed, as compared to grid cells containing no death at all (Table 3).However, there were more than 2.5 times more cells containing POM than there were healthy ones.No POM were delineated inside the "validation" area (Figure 2(b)) and assessment of the presence and absence was conducted independently by forestry professionals (Batini personal communications.2011).Examples of the comparison of the POM identified in-situ and use for training the model (Figure 2(a(.)reveal the complexity of the patch shapes and their stark visual differences between the surrounding, healthy vegetation.Figure 2 is a false-colour image in which red dominant colours indicate photosynthetically active cover (i.e., brighter red = more leaves), green-grey-black colours indicate non-photosynthetic cover and white colours are generally indicative of the high levels of reflection from bare soil and sands.Figure 3 shows the distribution of the POM across the study area including those areas excluded.

Changes in Fractional Cover
The MODFC data were used to investigate if an event of comparable extent had occurred prior to the summer of 2010-2011 and quantify the decadal natural variability in fractional cover across the study area.At 0.5 km resolution the MODFC data potentially contains information on changes in fractional cover for none, one or many POM (i.e., those detected from 10 m resolution SPOT-5 data).When MODFC is aggregated to 5 km the climate data will contain information relevant to many MODFC pixels.Addressing this spatial disparity the MODFC data were aggregated to 5 km and both the study area and POM location trends are extracted for comparison.This provided spatial parity between the impact of the POM in MODFC and AWAP data at 5 km yet did not yield any significant differences in the occurrence of the MSD exceedance events (data not shown).Considering all pixels across the study area, the MODFC time series (refer Figure 4) shows a simultaneous decrease in f pv (Figure 4(a).) and increase in f npv (Figure 4(b).)occurred following the summer of 2010-2011.Additionally, around 2001, there was a significant drop in f pv , which did not coincide with an increase in f npv.Comparing only pixels containing a modeled estimate POM, the average of all of these pixels (Figure 5) shows an equally clear and significant decline in f pv corresponding with an increase in f npv and importantly there was no significant change in f bs over the summer 2010-2011 period.In contrast to the impacts across the study area (Figure 5) there was no significant decrease in f pv during the summer of 2001-2002.

The Changing Climatic Space
Spatial analysis of the temporal climate data revealed the NJF has been experiencing a sustained net average reduction in precipitation over the last three decades.Annual precipitation in 2009 was 69% less (301 mm ± 38 mm) than the mean of 1981-2010 (907 mm ± 69 mm).In 2010, was 47% less (477 mm ± 38 mm) over the same period (refer Figure 6).Over this two year period, the study area received 1038 mm less precipitation than the long term mean which resulted in the ecosystem shifting to a more arid regime.Figure 6 shows pockets of sub-humid climate retracted to a more semi-arid climate.The study area has also been getting hotter over the study period.Summer 2010-2011 was 1.1 °C (17.7 °C ± 5.3 °C) warmer than the average summer temperatures (16.6 °C ± 4.6 °C) of the last three decades .This rate of warming has been increasing over the last decade.Whilst the annual mean minimum temperature across the study area was 0.2 °C ± 4.7 °C cooler (10.1 °C ± 4.7 °C) than the last three decades (10.3 °C ± 3.6 °C), the summer mean and maximum temperatures were higher.The mean 2010-2011 summer maximum temperature was 1.4 °C (31.6 °C ± 2.0 °C) higher than the previous 30 years (30.2 °C ± 2.0 °C).The 2010 annual mean diurnal range (BIO02) (Figure 7) (T max − T min ) was 1.6 °C higher (14.7 °C ± 0.5 °C) than the mean of the study period (13.1 °C ± 0.5 °C).As with precipitation, this illustrates a shift from a predominately semi-arid envelope towards an increasing arid one-completely departing its sub-humid zones.

Bioclimatic Indices and the Aggregated Extremes
The spatially aggregated indices (Figure 8) show the biologically significant changes in climate across the NJF between 1981-2010.There is a gradually increasing trend (1981-2010) in mean annual temperature (BIO01 refer Figure 8) driven by both an increase in the frequency of mean temperature years (white) over the last decade and a reduction in the frequency of cooler years (blue) between the hotter, more extreme years (red).The temperature indices, BIO05 (maximum temperature of the warmest month) and BIO06 (minimum temperature of the coldest month), show that warmest and coldest months have been less warm during 2000-2008 compared to 1981-1999.BIO05 shows 2009 and 2010 are the warmest of the last decade.BIO06 shows three of the coldest months between 2006-2010 were colder.Simply put, the seasonal temperature extremes are expanding, the coolest months are getting colder and the warmest-hotter.In 2010, the annual range of temperature (BIO07) was the highest of the entire study period.Winter (i.e., June, July and August) in SWWA is usually the wettest part of the year and the mean temperature of the wettest quarter (BIO08) has been declining in intensity since 2005.Summer (i.e., December, January and February) usually falls in the driest and warmest quarter and the mean temperature of the driest quarter (BIO09) has been on a warming trend since 2005 as has the mean of the warmest quarter (BIO10).Although the coldest month (BIO06) is getting colder, the mean of the coldest quarter (BIO011) has been increasing.[53] for the period 1981-2010.Derived from the AWAP [39] using Hijmans et al. [54].Red colours represent quantities greater than the mean (white) and blue less than, intensity represents magnitude.

Discussion
Over the last decade the study area has become drier, and the mean annual diurnal temperature range more extreme, increasing the risk of damage and mortality from a change in the frequency and intensity of heat waves and the potential for frost.As precipitation decreased the forest was exposed to greater mean annual diurnal temperature range which in turn contributed to an increase in evapotranspiration and gradually resulted in a bioclimatic shift towards aridity.These results suggest the bioclimatic factors preconditioning the POM across the study area were: 1.A succession of warm summers with summer maximum temperatures in 2010-2011 1.4 °C hotter than the 30 year average.2. A sustained multi-decadal reduction in wet years over the last 30, with 1982 being the wettest.3. A succession of extremely dry years, 2009 was 69% less than the 30 year average and 2010 was 47% less over the same period.4. A sustained increase in BIO02 over the last decade, peaking in 2010 at 1.6 °C above the 30 year average.
Whilst it is well documented that climatic warming is expected for the southwest of Western Australia, we found that the increased diurnal temperature range (BIO02) was a stronger indicator than an independent change in the maximum temperature-which was expected [55].Minimum temperatures also declined as the temperatures become more extreme (hot and cold) the mean annual diurnal temperature range is increasing, reaching a study period maximum of 14.7 °C ± 0.5 °C in 2010.This greater range increases the chance of frost in winter and heat waves in summer.In combination, the results suggests that in the presence of drought, exposure to both more intense cold and heat compounded vegetation water stress sufficiently to exhaust available energy reserves required for rooting and transpiration defense strategies-an effect found in plantation forests where soil depth and water storage capacity were found to be low [16].
This 2010-2011 event is also similar to the documented genesis of other forest mortality episodes [1] for which it is argued that they will ultimately accelerate species turnover and drive landscape evolution as climates shift under global climate change [60,61].These results also concur with the pioneering bioclimatic analysis of Gentilli [20] who predicted a similar precipitation threshold (+/− 300 mm) as that which sustains the forest in its present geographic location.The spatial analysis shows a reduction of 1038 mm accumulated precipitation (2009 + 2010), over 700 mm in excess of the Gentilli [20] threshold.When compared to the 1981-2011 mean, this equates to 69% less in 2009 and 47% less in 2010.As the climate of SWWA continues to become dryer and warmer we need to improve our understanding of the relationships between key climatic indicators and the spatiotemporal responses of our native forests at the landscape scale.
More recently, numerous studies have identified the bioclimatic envelopes that drive trait based species distribution across landscapes.The underlying objective of these studies is to parameterise models that attempt to predict future climate change vulnerability and the adaptive capacity of specific forest types (e.g., [8,37]).Many of these studies claim there is no clear link between strong-drought induced water stress and plant physiological responses [62][63][64] on the basis that stomatal control limits carbon assimilation during drought [65,66].The precipitation data shown here suggests water limitation (i.e., drought and water stress) has been the primary driver of this mortality event.
Recent studies of Mediterranean climates indicate that severe droughts can result in non-stomatal limitation of photosynthesis relative to leaf water stress [67] when driven by low water availability and high evaporative demand over time [68,69].Such responses confound the assimilatory capacity of vegetation and in severe drought conditions, can exceed physiological water-carbon deficiency thresholds that could in turn trigger reduced recovery hysteresis, cascading effects and mortality [24,[70][71][72].Further confounding this debate is evidence that slightly warmer temperatures within a tree's optimum range could enhance forest growth under enhanced CO 2 [73].
Satellite imagery has become a well-accepted source of information, providing useful information on the spatio-temporal variation of vegetated landscapes.In this study we used 3 different remotely sensed data; (1) DMSI airborne imagery (0.5 m), (2) SPOT-5 high-resolution (10 m) and (3) MODIS moderate resolution (0.5 km) satellite imagery that were further aggregated to 5 km.We combined them to first detect and map a decline event at a single point in time, highlighting the significant spatial extent of the event.Then we go back in time to understand how the past impacted on the changes in photosynthetic cover (i.e., forest canopy) we observed and modeled following the summer of 2010-2011.We found both satellite sensors (SPOT-5 and MODIS) detected changes in photosynthetic and non-photosynthetic cover following the summer of 2010-2011.The MODFC data provided a long-term (2000-2012) baseline for fractional cover suggesting that in the years prior to the 2010-2011 event the POM modeled sites were persisting within their environmental niche-in a "healthy" state with no evidence of widespread stress.In addition to the independent validation, the MODFC analysis provides greater confidence in this multiple sensor methodology.The similar results following spatial down-aggregation of the MODFC data, i.e., to the same resolution as the climatologies, adds yet further confidence to the usefulness of MODFC for mapping such events at lower resolutions when no better data are available.Additionally, this approach demonstrates the potential use of programmable satellites such as SPOT-5 to acquire these data on demand and for broad-scale, near-real time post-mortem POM detection.
The threat of fire in the NJF remains significant as long as the residues of each POM remain scattered across the forest.Should the POM combust, it could potentially create a wildfire in the same way that the 1998 Yellowstone fires evolved from distributed patches into one large conflagration which persisted for months [74] and such an event could have a significant impact on the biodiversity of the flora and fauna of the forests structure and threaten the welfare of the local population.

Conclusion
Expansion of the thermal extremes of minimum and maximum temperature and prolonged drought were major drivers of forest mortality in the NJF following the summer of 2010-2011.This forest scale mortality event was successfully captured using SPOT-5 imagery using a three-rule feature extraction algorithm that produced a map of their distribution across the landscape.Used in this way, the combination of remote sensing and gridded climatologies provides a means of monitoring, mapping and analyzing the drivers of bioclimatic induced change of forested ecosystems for known events and a basis for a new generation models on similar spatial and temporal scales.As such, this approach makes a significant contribution towards a robust methodology for monitoring the impact of climate induced drought and the underlying ecological consequences of extreme climate variability on forest dynamics.

Figure 1 .
Figure 1.Map featuring the study area and the Interim Biogeographical Region of Australia's NJF and SCP surrounding Perth.

Figure 2 .
Figure 2. Système Pour l'Observation de la Terre mission 5 (SPOT-5) infrared imagery (10 m) overlaid with (a) field data and feature extraction modeled patches of mortality (POM) and (b) the model evaluation grid and POM contained within a number of cells.
object based classification of the presence and absence of healthy (n = 38) and patches of mortality (n = 97) across 135 km 2 (n = 135) within in the NJF.

Figure 3 .
Figure 3. SPOT-5 Imagery overlaid with the mining, water and cloud masks used across the study area in modeling POM and the validation grid used.

Figure 4 .
Figure 4. Spatially aggregated 0.5 km study area average Seasonal Trend Decomposition using Loess (STL) trend plots showing the 2000-2012 period with the summer of 2010-2011 marked for each component of MODFC Guerschman [33].Subplot (a) is the fraction (%) of photosynthetic cover (f pv ); (b) the fraction (%) of non-photosynthetic cover (f npv ) and c. the fraction (%) of bare soil (f bs ).

Figure 6 .
Figure 6.Spatially aggregated comparison of the NJF's Mean Annual Aridity Index versus Mean Annual Precipitation (mm/year).Averages over the period 1981-2009 are shown in grey and the 2010-2011 in black.

Figure 7 .
Figure 7. Spatially aggregated comparison of the NJF's Mean Annual Aridity Index versus Summer (December, January and February) Average Diurnal Temperature Range (deg C).Averages over the period 1981-2009 are shown in grey and the 2010-2011 in black.

Figure 9 .
Figure 9. BIO02 is the 'Mean Diurnal Range', the mean of the difference of the monthly maximum and minimum temperature[53], for the period 1981-2010.Derived from the AWAP[39] using Hijmans et al.[54].

Figure 10 .
Figure 10.Spatial aggregated monthly time series of mean daily rainfall (mm) across the study area showing the loess trend, mean and standard deviations and the original monthly data, extracted from AWAP [39] for the period 1981-2010.

Table 1 .
[52]ity-humidity Classification of the value of MAAI, the ratio of actual to potential evaporation (after[52]).