Distant neighbors: recent wildfire patterns of the Madrean Sky Islands of southwestern United States and northwestern Mexico

Background: Information about contemporary fire regimes across the Sky Island mountain ranges of the Madrean Archipelago Ecoregion in the southwestern United States and northern Mexico can provide insight into how historical fire management and land use have influenced fire regimes, and can be used to guide fuels management, ecological restoration, and habitat conservation. To contribute to a better understanding of spatial and temporal patterns of fires in the region relative to environmental and anthropogenic influences, we augmented existing fire perimeter data for the US by mapping wildfires that occurred in the Mexican Sky Islands from 1985 to 2011. Results: A total of 254 fires were identified across the region: 99 fires in Mexico (μ = 3901 ha, σ = 5066 ha) and 155 in the US (μ = 3808 ha, σ = 8368 ha). The Animas, Chiricahua, Huachuca-Patagonia, and Santa Catalina mountains in the US, and El Pinito in Mexico had the highest proportion of total area burned (>50%) relative to Sky Island size. Sky Islands adjacent to the border had the greatest number of fires, and many of these fires were large with complex shapes. Wildfire occurred more often in remote biomes, characterized by evergreen woodlands and conifer forests with cooler, wetter conditions. The five largest fires (>25 000 ha) all occurred during twenty-first century droughts (2002 to 2003 and 2011); four of these were in the US and one in Mexico. Overall, high variation in fire shape and size were observed in both wetter and drier years, contributing to landscape heterogeneity across the region. Conclusions: Future research on regional fire patterns, including fire severity, will enhance opportunities for collaborative efforts between countries, improve knowledge about ecological patterns and processes in the borderlands, and support long-term planning and restoration efforts.


Background
In today's world, understanding between neighboring countries is more vital than ever (Riding 1989). Although the importance of history, culture, and socioeconomic conditions is often emphasized, ecological processes that shape the distribution of species and communities, and thus support resources essential to human livelihoods, are equally critical components of integrated systems that transcend political borders (Liu et al. 2007). Wildland fire is an important ecological process that shapes the distribution of species and communities and associated ecosystem services. Wildfire can enhance or reduce species habitat and connectivity, alter water cycle processes and induce erosion, and influence forest carbon budgets (Cannon et al. 2010, Hurteau and Brooks 2011, Williams et al. 2014, Hutto et al. 2016. Information on the timing, location, and patterns of fires is key to understanding the physical, biological, and ecological processes that generate those patterns, including evolutionary processes that drive biological and ecological diversity upon which people depend (Bond and Keeley 2005, Burton et al. 2009, Whitman et al. 2015. Maps of fire occurrence are a key resource in management and ecological restoration efforts, and can provide insight into the effects of land use, climate, and topography on contemporary fire regimes (Morgan et al. 2001;Stephens and Fulé 2005;Miller et al. 2012a, b). Where available, consistent and accurate data on fire regimes, including frequency, severity, and spatial pattern, facilitate important decisions regarding sustainable forest management and reduction of fuels to control fire risk (Bergeron et al. 2004, Petrakis et al. 2018. Furthermore, fires increase the spatial heterogeneity of vegetation at multiple scales; knowledge of fire regime attributes can be used to understand how fires create, maintain, or otherwise alter wildlife and plant habitats, as well as how fires support or disrupt the connectivity of those habitats (Zozaya et al. 2011). Importantly, fire regimes are expected to change with climate; alterations have already been observed at local to global scales (Hurteau et al. 2014). In places where resources necessary to develop and support fire mapping programs are limited, it is impossible to address these priority issues in a way that is comparable with adjoining countries in which mapping programs are well established.
In the US, federal agencies developed the Monitoring Trends in Burn Severity (MTBS) program to provide a comprehensive, national database of wildfire information about "large" fires (>404 ha in the West and >202 ha in the East) from 1984 to the present (Eidenshink et al. 2007). MTBS data have been used to examine regional and national wildfire trends related to vegetation type, topography, forest management, and climate, providing information to guide policy and management decisions (Wimberly et al. 2009, Dillon et al. 2011, Miller et al. 2012a, Baker 2013. Prior to our efforts, a fire database for northern Mexico that is comparable to MTBS did not exist, and institutional information on wildfires that burned in remote areas was lacking due primarily to complex land ownership (i.e., more communal and private land in forests of Mexico compared to lands administered by federal agencies in the US) and limited forestry and firefighting budgets (Rodríguez-Trejo 2008, Villarreal and. Our work is an important first step in characterizing contemporary fire regimes in a cross-border region with outstanding biodiversity (Whittaker andNearing 1965, Coblentz andRiiters 2004) that is expected to experience the ecological impacts of climate change sooner than other areas of western North America (Notaro et al. 2012).
The objectives of this research were to quantify the spatial and temporal patterns of recent (1985 to 2011) fire occurrence across the Madrean Archipelago Ecoregion of North America, and to analyze these patterns relative to significant climate, topography, ignition sources, and land cover variables. The Madrean Ecoregion provides a unique setting to study spatial patterns of wildfire: mountains in the ecoregion are biogeographically similar, but the region is divided in half by the United States−Mexico international border, where contrasting land uses and forest management approaches often occur. We posed the following questions: 1) What are the basic characteristics of contemporary fire regimes across the study region? 2) How do spatial patterns of fires vary regionally, and within each Sky Island mountain range? 3) How does inter-annual climate variability influence fire patterns? 4) How does underlying spatial heterogeneity in climate, topography, ignition sources, and land cover define burned and unburned environments?

Study area
Covering a geographic area of approximately 74 788 km 2 , the Madrean Archipelago Ecoregion (Omernik 1987) encompasses portions of northern Sonora and Chihuahua in Mexico, and southeastern Arizona and southwestern New Mexico in the US (Fig. 1). Characterized by a basin and range topography with elevations ranging from 600 m in the valley floors (Tepache, Sonora) to 3267 m (Mount Graham, Arizona), the Madrean Archipelago consists of isolated forested mountains (referred to as Sky Islands) surrounded by desert lowlands, semi-arid grasslands, and shrub steppe vegetation (Warshall 1995). The ecoregion contains 39 Sky Island mountain complexes that make up more than half of the total land area (42 729 km 2 ; Fig. 1; Warshall 1995). The international border, established after the Gadsden Purchase in 1853, cuts through the center of the ecoregion at 31°33' N ( Fig. 1). Twenty-two of the Sky Islands fall primarily within the US and 17 in Mexico. Three ranges-the Atascosa-Cibuta, Huachuca-Patagonia, and Peloncillo-Pan Duro-cross the international border. Characteristics of the 30 Sky Islands that had fires during the study period are presented in Table 1. Climate of the Madrean Archipelago is characterized by a bimodal precipitation pattern. Winter (November to March) frontal storms bring snow to higher elevations and rain to the lower elevations. In mid to late summer (July to August), moisture from the south creates a monsoonal weather pattern with thunderstorms. The arid foresummer period (April to June) preceding the monsoonal moisture is typically hot and dry and generally corresponds to the fire season. Annual precipitation varies greatly by elevation, with orographic uplift contributing to increased precipitation at upper elevations. The average annual precipitation ranges from about 232 mm at 820 m above sea level to 1035 mm at 3267 m (Fort Thomas, Arizona, and Mount Graham, Pinaleño Mountains, respectively; Hamann et al. 2013).

Sky Island Vegetation
Upland vegetation of the ecoregion generally consists of Sonoran Desertscrub, Chihuahuan Desertscrub, and Sinaloan Desertscrub, and Semidesert Grasslands in the low-elevation valleys and flats, and Madrean Evergreen Woodlands, Montane Conifer Forests, and Subalpine Conifer Forests at the upper elevations (Brown et al. 1979). The sharp vertical relief and unique geographic position of the Sky Islands combine to create a diverse mosaic of vegetation assemblages stratified along elevation and moisture gradients (Niering and Lowe 1984). The Rocky Mountain Cordillera vegetation influence from the north is more prevalent at higher elevations (above 2000 m) while the Madrean vegetation influence from the south dominates lower elevations (below 2000 m) (McLaughlin 1994). Coniferous vegetation within the Sky Islands is generally divided into five types: pinyon−juniper woodlands, Madrean evergreen woodlands, pine−oak forests, pure pine forests, and mixed-conifer forests (Iniguez et al. 2005).
Pinyon−juniper woodlands typically occur in areas below 1800 m elevation, although they can also be found at higher elevations on drier south-facing aspects.

Cultural history
The oldest known inhabitants of the general area were the Hohokam, who, as farmers, mainly occupied the desert valley bottoms near reliable water sources (Waters and Ravesloot 2001). Apaches moved into southwestern Arizona in the 1680s and 1690s and, unlike prior inhabitants of the Southwest, they utilized both valley bottoms and upland forests (Bahre 1991). Spanish settlers began moving into the region in the sixteenth and seventeenth centuries, and by the 1600s, they had established settlements around missions near major rivers including the San Pedro, Santa Cruz, Magdalena, and the Rio Sonora (Clemensen 1987). Demand for timber for construction and fuel in the new settlements had a considerable impact on the surrounding forests (Bahre and Hutchinson 1985). However, the presence of Apache groups restricted Euro-American settlement in and around high-elevation forests of the Sky Islands until the end of the nineteenth century (Bigelow Jr 1968).
After the Gadsden Purchase in 1853, the region was divided and the current political borders were established; around this time, the arrival of the Southern Pacific Railroad and subjugation of the Apaches led to an increase in Euro-American settlement that brought about significant land use changes (Bahre 1991). Between 1902 and 1907, several federal forest reserves were established in the Sky Islands of Arizona and New Mexico (Allen 1989). Following the historic fires of 1910 in the western US, federal land management agencies adopted aggressive fire suppression policies (Pyne 1982). In 1935, the US Forest Service established the "10 a.m. policy," with the goal of suppressing every fire by 10 a.m. the day following its initial report. Fire suppression became even more effective after the Second World War with advancement in firefighting technologies, particularly more abundant aircraft and fire suppressing chemicals. In northern Mexico, many of the Sky Island forests were part of privately owned haciendas until the Mexican Revolution (1910Revolution ( to 1920 and the passing of the Agrarian Code in 1934, which reallocated private lands to local communities, called "ejidos" (Thoms and Betters 1998). Unlike in the United States where much of the forest land was federally managed, under the ejido system, local communities were responsible for managing forests and wildfires; therefore, wildfires in northern Mexico were not suppressed to the same extent that they were in the western US. During the 1980s, the forest service of Mexico, aided by the US Forest Service and US Agency for International Development (USAID), developed a more active fire management program (Rodriguez-Trejo et al. 2011); however, the scale of firefighting in Mexico during that period was much smaller, with Mexico spending around one tenth of the dollar amount per acre on wildland firefighting activities compared to the US (Villarreal and Yool 2008).

Fire history
Prior to the twentieth century, wildfires within the Madrean Sky Islands were relatively common, particularly at higher elevations where greater productivity created abundant fuels and a dry pre-monsoon period allowed these fuels to cure and then burn (Baisan and Swetnam 1990). Tree-ring based fire history studies have found that, prior to the 1900s, pine-dominated sites experienced, on average, large fires every 6 to 15 years , Meunier et al. 2014. Mixed-conifer sites, on the other hand, experienced large fires less frequently at intervals ranging from 8 to 25 years (Swetnam and Baisan 1996, Iniguez et al. 2008, Margolis et al. 2011. While fire severity varied by site, fires in places historically dominated by pines were predominantly low severity , Iniguez et al. 2009). Reduced wildfire activity in the mountains surrounding major human settlements has been documented for the US Sky Islands at the end of the nineteeth century Swetnam 1990, Swetnam et al. 2001). These changes have been associated with a sequence of historical events: the end of the Apache wars, the arrival of the railroad, and the introduction of large-scale livestock grazing. As a result, fine fuels required for fire spread were effectively removed over large areas (Savage andSwetnam 1990, Belsky andBlumenthal 1997).
North of the border, the last large wildfires prior to fire suppression occurred between 1900 and 1910; to the south, this process was delayed until the 1930s, following the Mexican Revolution (1910Revolution ( to 1920. While some authors suggest that the allocation of ejidos marks the onset of fire suppression (Heyerdahl and Alvarado 2003), this is not the case for some sites in the Madrean Archipelago or even farther south (Kaib 1998, Cortés Montaño et al. 2012. However, in Mexico, the exclusion of wildfire was less extensive; in some more remote locations, frequent fires continued up to the present day (Fulé et al. 2011, Meunier et al. 2014. Moreover, some of these sites where fires still persist contain forest conditions similar to those believed to dominate the region prior to fire suppression (Cortés Montaño et al. 2012, Yocum Kent et al. 2017.

Fire identification with satellite imagery
Since information on dates and locations of many past fires was not available for the Sky Islands in Mexico, we used a semi-automated remote sensing fire-mapping approach to identify and delineate fires that occurred within the Mexican Sky Island mountain ranges. Our approach adapted satellite-based fire mapping and burn severity techniques used by MTBS in the US that are typically applied on an ad hoc basis, with knowledge of specific fire date and location (Monitoring Trends in Burn Severity: https://mtbs.gov/mapping-methods). The techniques take advantage of the changes in spectral information observed in remotely sensed imagery.
Post-fire mapping of burned area and severity are typically accomplished either through visual interpretation of aerial photography, or by using band-ratio techniques on satellite spectral information (White et al. 1996, Miller and Yool 2002, Lentile et al. 2006. Landsat multispectral imagery has proven suitable for fire mapping, as the data are relatively high resolution (30 m) with a 16-day return, and contain shortwave infrared (SWIR) and near infrared (NIR) information that can be used to discriminate a range of burn severity classes (White et al. 1996). The MTBS method uses Landsat data to identify fire perimeters and estimate burn severity classes using the differenced Normalized Burn Ratio (dNBR) method (Eidenshink et al. 2007). The Normalized Burn Ratio (NBR) is derived from the near infrared reflectance in Landsat band 4 (0.63 to 0.69 μm) and shortwave infrared reflectance in band 7 (2.08 to 2.35 μm) from a single Landsat image (Cocke et al. 2005, Key 2006, and the dNBR is a simple differencing of pre-and post-fire NBRs. The NBR equation is as follows: An important aspect of the dNBR is careful selection of pre-and post-fire images to reduce error introduced by changes in scene conditions, including vegetation phenology and climate-induced conditions (Key 2006). Land surface and atmospheric conditions such as snow, vegetation change in wet versus dry seasons, shadows, smoke, and clouds all affect the outcome of the dNBR product (Key 2006).
The two steps of our semi-automated process were as follows: the first step involved calculating the NBR for a spring or early summer (assumed to be pre-fire) and a fall or early winter (assumed to be post-fire) Landsat scene of each year (1985 to 2011), calculating a dNBR with the two NBRs, then applying a dNBR threshold to identify locations where vegetation change was likely due to fire and to eliminate many non-fire related pixel changes (Prickett 2013). To calculate the dNBRs, we obtained 1517 cloud free Landsat 5 TM images (1985 to 2011) from the USGS Global Visualization Viewer (GLOVIS) website for seven Landsat path or row designations that cover the Madrean Archipelago Ecoregion. All images were corrected for atmospheric effects and haze, and converted to reflectance using the ATmospheric CORrection (ATCOR) module for ERDAS Imagine. The dNBR-derived detections were then filtered by size (<100 ha) and location (intersecting Sky Island polygons) to remove small non-fire vegetation changes (e.g., changes from productive agricultural fields to fallow), and positive identifications were converted to polygons. An initial accuracy assessment of fire detections tested against MODIS Fire Information for Resource Management System (FIRMS) daily detections (from 2001 to 2011) indicated that the automated process correctly identified about 42% of the fires in Mexico during this period.
The second step involved manually identifying additional fires in Mexico through visual interpretation using a combination of dNBR images and false color composites (displaying bands 6, 4, and 3 in RGB) from the full Landsat archive. Once identified, fire perimeters were digitized on screen using the dNBR and post-fire color composites. To assess omission errors that may have occurred due to data limitations (e.g., too many cloudy scenes during post-fire period to detect fire scar), we validated the burn perimeter data with MODIS FIRMS daily detections from 2001 to 2011. We also validated MTBS fires in the US with the FIRMS dataset to identify any commission errors in that dataset and to better ensure comparability of fire regimes across the region from different data sources and methods.

Derived and ancillary datasets
Using the fire perimeter polygons (rasterized using 100 m [1 ha] cell size), we developed a database of patch statistics for each fire, including those that overlap (Table 2). We selected a set of metrics that were representative of several attributes including size, shape, and orientation. The metrics allow characterization of fire patterns that may reflect underlying aspects of the environment or important ecological processes (Table 2). To represent regional climate conditions through the study time period, we obtained the self-calibrating Palmer Drought Severity Index (scPDSI; van der Schrier et al. 2013). The scPDSI provides a standardized, time-integrated drought index that takes into account precipitation, temperature, and available moisture capacity of soils (van der Schrier et al. 2013). A variation of the original Palmer Drought Severity Index (Palmer 1965, Osborn et al. 2017, scPDSI was developed to provide more robust comparisons across Table 2 We quantified fire spatial patterns using metrics representing fire size, shape, and orientation. The table includes a reference for the program used to calculate the metrics, a general description of each, and a summary of their potential ecological significance in terms of pattern−process relationships of interest. Madrean Ecoregion Sky Islands, US and Mexico, 1985  Greater interior, or core area, reflects a higher proportion of area burned that is isolated from the surrounding unburned matrix and its intact resources. Interior areas are likely to experience very different physical conditions than disturbed areas near intact vegetation (Turner et al. 2001).
The interior area of the fire, greater than 100 m from the fire perimeter as a percentage of patch area Shape Eccentricity (Bui et al. 2012) Ratio of the distance between the two farthest extremities of each fire and the distance of the maximum width of its perpendicular axis (i.e., length-to-breadth) Fire shape can reflect underlying factors such as physiography, hydrography, and surficial deposits (Mansuy et al. 2014), time period, vegetation, climate, and management (Parisien et al. 2006). Shape interacts with fire size and severity to influence rate and direction of successional trajectories (Turner et al. 1998

Orientation
Bearing (Hijmans 2016) Orientation (degrees) of the axis of the two farthest extremities of each fire The orientation of a burn can reflect predominant winds and/or topographic (e.g., watershed) orientation (Barros et al. 2013). Given this potential interaction with topography, bearing may also influence the traversability of the burned area by wildlife and plants.
diverse climates (van der Schrier et al. 2013). Approximately 80% of Sky Island fires occurred between the end of April and the beginning of August (Julian dates 117 to 221), so we selected April to represent pre-fire season conditions. At each fire polygon centroid, we extracted the April scPDSI value for the year that the fire burned. We adopted a standard classification of the index to interpret the results (for drought: extreme < −4; severe < −3; moderate < −2); wetter classes are assigned the same threshold values, but with positive sign (Palmer 1965, Osborn et al. 2017).
In addition, we used several datasets to represent spatial heterogeneity of environmental conditions, including those representing physical and biological factors that influence the occurrence of fire, as well as human presence and activities that affect ignition sources and use and management of fire (Table 3). We obtained a 30 m Shuttle Radar Topography Mission (SRTM; Farr et al. 2007) digital elevation model from Google Earth Engine (Gorelick et al. 2017) and resampled the data to 1 km resolution for computation of terrain ruggedness. The spatial climate data included 26 bioclim layers developed by the ClimateWNA project (Hamann et al. 2013) and were based on 1981 to 2010 climate normals (Table 3). Ignition source data were obtained from National Aeronautics and Space Administration (NASA) (lightning; Albrecht et al. 2016); anthropogenic biomes (Ellis and Ramkutty 2008) and roads were downloaded from Commission for Environmental Cooperation (see CEC 2009). Population density data were obtained from Socioeconomic Data and Applications Center (SEDAC; http://sedac.ciesin.columbia.edu/). The lightning, roads, and population density data were smoothed by applying a uniformly weighted kernel (55 km radius window) before sampling. We used Brown and Lowe's Biotic Communities of the Southwest (CBI 2004; available at http://databasin.org) to characterize variability in vegetation composition and structure (Brown et al. 1979).

Characteristics of fire regimes
We used the fire perimeter database to characterize the regional and temporal distribution of fire frequency and area burned for each island and for each year during the study period. We mapped the variability in the landscape structure of fire patterns, including composition (area burned or fire size) and configuration (shape metrics) for each fire across the ecoregion.

Influence of climate on fire patterns
To understand how regional climate influences fire patterns, we used the R package (R Core Team 2017) ggplot2 (Wickham 2009) to generate a hexagon heatmap of two-dimensional bin counts (pattern metric by year) and graphed the results using a color gradient of scPDSI values. We specified 20 bins in horizontal and vertical dimensions for the plots. The position of each bin represents the center of mass (average of x, y values) for that bin and the mean scPDSI for the x, y values in each bin determines its color. The hexagon plots provided a visual tool to examine how fire size and shape varied in relation to wetter and drier years and whether patterns across years indicated trends during the study period.

Variation in fire patterns in relation to underlying spatial heterogeneity
We randomly sampled the study region within islands that experienced fire to obtain values for all spatial variables (Table 3). We applied a mask of fire perimeters to sample within areas that burned (total n = 6670; US n = 3699; Mexico n = 2971) and areas that did not burn (total n = 29 731; US n = 8759; Mexico n = 20 972). We used the sample data to conduct a graphical analysis of burned and unburned environments.
For categorical variables (anthropogenic biomes and biotic communities; see Table 3), we examined patterns in the data using stacked bar charts, for which height of the bar in each color represents its relative occurrence in burned and unburned samples. For continuous variables, we displayed the distribution of burned samples and unburned samples using smoothed density plots to aid in identifying the environmental space occupied by fires in the US and Mexico. We tested the null hypothesis that samples of burned and unburned environments came from a common population using Anderson-Darling statistics (AD.T; k-sample criterion AD; Scholz and Zhu 2016). Because several of the climate variables exhibited similar patterns, we included five of the 26 climate variables in the results, as well as population density, road density, terrain ruggedness, and lightning flash rate (Table 3).

Validation of fire detections in US and Mexico
A total of 254 fires (>404 ha) were identified across the region from 1985 to 2011, 99 fires in Mexico and 155 in the US. (The fire perimeter database is available from US Geological Survey, see Villarreal and Poitras 2018). Of the 254 fires mapped, 142 occurred between 2001 and 2011 and were validated using MODIS FIRMS detections. All 142 fires were confirmed with FIRMS data (no false positives); however, the FIRMS data set contained six fires that were not mapped by MTBS or captured with our mapping approach for Mexico. Four of the missed fires were located in the US and two in Mexico. In all cases, the fires were not mapped because multiple consecutive Table 3 Spatial data used to describe fire environments, including topographic, climatic, ignition sources (including human influence), and vegetation in the Madrean Ecoregion Sky Islands, US and Mexico. The spatial density of roads, lightning, and population was calculated using a uniformly weighted kernel (~55 km radius moving window) for use in analysis. The bioclim variables graphed in figures are listed in bold. For biotic communities and anthropogenic biomes, only categories present in the study region are given in the  (1989, 1994, 2002, and 2011), each of which included more than 20 large fires and burned more than 75 000 ha yr −1 , with a maximum of 231 249 ha burned in 2011 (Fig. 2). Fire patterns during some large fire years varied across the ecoregion; much of the fire activity and area burned in 1994 occurred in the US (Fig. 2) and, in 2011, the area burned in US Sky Islands was more than twice that of Mexico, with the Horseshoe II Fire (91 476 ha) in the Chiricahua Mountains contributing to that difference (Fig. 2). Over the entire study period, the Pan Duro, Huachuca-Patagonia, and Los Ajos-La Madera experienced the greatest number of fires (37, 25, and 21, respectively; Fig. 2). The Pan-Duro had the most hectares burned (136 530 ha), followed by the Chiricahua (116 004 ha) and Atascosa-Cibuta mountains (89 048 ha).

Variation in spatial and temporal patterns of fire patch metrics
Maps of fire size and core area index (i.e., interior area of the fire, greater than 100 m from the fire perimeter as a percentage of fire size) indicate that many of largest fires and fires with considerable core area (>0.90) occurred along the border in the Atascosa-Cibuta, El Pinito, and Pan Duro mountains, and in the Chiricahua and Santa Catalina mountains in the US (Fig. 3). Eccentricity (i.e., ratio of the distance between the two farthest extremities of each fire and the distance of the maximum width of its perpendicular axis, or length-to-breadth ratio) was generally high for most fires (i.e., >0.6), indicating more elongated shapes. Shape index (i.e., fire shape complexity; the degree of departure from a circular shape; ranges from 1 [round] to >1 [complex]) varied widely by fire and range, but complexity tended to be

Influence of regional climate
The scPDSI patterns indicated that drier and wetter pre-fire seasons generally alternated, but multiple, consecutive dry seasons (scPDSI < −2) became more frequent after year 2000 compared to earlier years (Fig. 4). Across years and in relation to climate, spatial patterns were heterogeneous, with a fairly wide range of metric values. Maximum fire size generally increased through time, and this trend was driven by several large fire events that occurred during the droughts of 2000 to 2003 and 2011 (Fig. 4a). Interior areas of fires reflected a lower range of values in early, wetter years (pre-fire season scPDSI > 2) and a clustering of higher values in recent, drier seasons with scPDSI < −2 (i.e., core area index; Fig. 4b).
Shapes, overall, varied in complexity, but the upper range of values for Perimeter Area Ratio exhibited a decreasing trend across years with variable climate (Fig. 4c). More circular shapes, (i.e., shape index < 2) were common, but fires with highly complex shapes were less frequent across years and across variability in scPDSI (Fig. 4d). Recent, drier seasons (scPDSI < −2) produced some of the most complex shapes, but 2001, a year with wetter values (−2 < scPDSI > 2), particularly for some fires, also included complex shapes. Eccentricity, or the length-to-breadth ratio, indicated that many of the fires tended toward elongated, rather than circular, shapes (index value > 0.6; Fig. 4e). Bearing ranged from 23 to 73 degrees (north-northeast to east-northeast) for all fires (Fig.  4f ); these values indicate the orientation of the two farthest extremities of each fire (Table 2).

Fire patterns and underlying spatial heterogeneity
Biotic communities that experienced fire occurred in different proportions to those in the unburned landscape (Fig. 5, top). A majority of the burned area occurred in Madrean Evergreen Woodlands and Petran Montane Conifer Forests; these communities occurred in greater proportion on burned versus unburned landscapes. Semidesert Grassland occupied a large proportion of unburned area, and a somewhat lesser proportion of burned landscapes. Unburned areas were predominantly Madrean Evergreen Woodlands and Semidesert Grassland with higher proportions of Arizona Upland Subdivision−Sonoran Desertscrub, Chihuahuan Desertscrub, and Sinaloan Thornscrub than burned landscapes. Burned and unburned landscapes also differed in anthropogenic biome composition (Fig. 5, bottom). The most noted contrast was in the higher proportion  1989, 1994, 2002, and 2011 were active fire years. Most of the islands with large amounts of burned area occurred at middle latitudes of burned Remote Forests, relative to its abundance in unburned landscapes. Populated Rangelands and Remote Rangelands were less commonly burned than unburned; Villages and Residential anthromes composed a small proportion of the study region, with no evident difference in their occurrence on burned and unburned areas.
Burned areas occupied distinct environments within their unburned context based on visual interpretation of smoothed density plots (Fig. 6) and results of statistical tests (AD.T P-value < 0.001 for all burned-unburned comparisons; see Appendix).
Results for the climate variables indicated that fire environments occurred at the lower to mid-range of heat moisture in both countries, corresponding to a somewhat higher range of mean annual precipitation and lower range in mean annual temperature (Fig. 6a−c). Fires occurred more often in places with later onset of the frost-free period; this trend was observed in both countries, but the density plot for burned samples was particularly skewed toward later Julian dates for Mexico (Fig. 6d). Burned areas were coincident with a higher range in mean summer precipitation in the US, but the peak in burned sample density occurred at Fig. 3 Distribution of fire patch patterns across the Sky Islands of the Madrean Ecoregion, US and Mexico, between 1985 and 2011: fire size and core area index (top row) and eccentricity and shape index (bottom row). Larger fires tended to have greater core area, but more circular (versus more elongated) patterns occurred across the region, as did fires with more or less core area relative to their size a slightly lower value for mean summer precipitation (MSP) than the unburned maximum density in Mexico (Fig. 6e).
Natural ignitions did not limit the occurrence of fire, based on the smoothed density plots for US and Mexico (Fig. 6f ). Specifically, burned environments were located in the lower range of values for lightning flash rate, as well as in places where rates of lightning strikes were higher. Areas in Mexico with the highest flash rates had lower density of area burned or remained unburned. Population density smoothed plots for burned−unburned samples indicated that fires occurred in very remote areas, as well as more populated areas (Fig. 6g). In Mexico, greater density of unburned samples occurred in sparsely populated areas, compared to burned samples. Burned environments included places with a wide range of road density, especially in the moderate to high range, but in both countries, fires were limited where roads were most prevalent (Fig. 6h). Burned environments were located in flatter, as well as more complex terrain; burned sample distributions tended to encompass places with higher values of terrain ruggedness index for both countries, compared to unburned (Fig. 6i).

Characteristics of contemporary fire regimes
Recent fires in the Madrean Ecoregion exhibited several distinctive spatial and temporal patterns. Area burned varied by year and across islands (Fig. 2), but generally increased over the study period, culminating with the 2011 fire year, when the Horseshoe II Fire burned much of the Chiricahua Mountain range in Arizona (91 476 ha), and the El Pinito (28 519 ha) and Murphy Complex fires (26 044 ha) burned along the international border. Fire size in the US was more variable than in Mexico, and much of that variability was related to several large twenty-first century fires. In general, Sky Islands near the border had the greatest number of fires during the period (Pan Duro,37;25;21), and many of the large border fires exhibited complex shapes. As proportion of Sky Island size, the fires in the US tended to burn more area, especially in the Animas, Chiricahua, Huachuca, and Santa Catalina mountains. Many of the southern islands in Mexico had only a few small fires during the study period. The pattern of larger fires north of the international border may be related to observed increases in forest  (Table 2) by year for fires that burned between 1985 and 2011 in the Madrean Ecoregion Sky Island, US and Mexico. Metrics include: a) Fire size, b) Core area index, c) Perimeter area ratio, d) Shape index, e) Eccentricity and f) Bearing. The position of each bin represents the average of the x, y values assigned to the bin; bin color corresponds to the mean scPDSI for fires in the bin (orange to green scale represents drier to wetter conditions). Due to outliers, fire size is presented in log scale on the y-axis of Fig. 4a. Fires exhibited heterogeneity in spatial patterns across years and climates. In early time periods, fires occurred in a mixture of wetter and drier years, but after 2005, a majority of fires burned in years with moderate to severe drought conditions (scPDSI < −2) fuel loads at the stand level and related fuel homogeneity at landscape levels due to a longer period of fire exclusion (Moore et al. 2004).

Climate and topographic influence on fire patterns
We found that, during our study period, the number of fires and area burned were greatest during drought years identified in other studies (e.g., Dillon et al. 2011;Fig. 2). In particular, 1989, 1994, 2002, and 2011 stand out as large fire years, with at least 20 fires (>404 ha) and more than 75 000 ha burned each year. However, we found that fire size and complexity were highly variable in relation to scPDSI, suggesting that fires increased landscape heterogeneity through a range of variation in fire shape and size in years with a mixture of drier and wetter pre-fire seasons (Fig. 4). But we also observed that larger fire sizes followed trends of drought and increased temperature across the region, a trend that is already affecting the distribution of plants in the Sky Islands along an elevational gradient (Brusca et al. 2013). The strong relationship between drought and fire was notable for the Sky Islands in the southwestern part of the ecoregion in Mexico: 12 of the 13 total fires in four Sky Islands (Aconchi, El Carmen-Verde, La Huerta, and La Madera-Cucurpe) occurred during 1989, 1994, 2002; however, these fires were generally of moderate size (μ = 2911 ha).
A number of studies in the American Southwest have documented the strong historical relationship between wildfire and antecedent climate Betancourt 1998, Crimmins andComrie 2004). Although similar studies that include the Mexican Sky Islands are not as common, a recent study by Meunier et al. (2014) in the Sierra San Luis found that, prior to 1886, drought years during the fire year were most conducive for fires. After 1887, in remote areas that continued to burn, the pattern appears to shift and wet conditions the year prior to the fire become a more important factor (Meunier et al. 2014). This suggests that wet years result in increased biomass and fuel production that then cures a year later and supports the ignition and spread of fires. Because patterns of scPDSI varied considerably across our study region, the relationship between scPDSI and fires was complex. There is a need to  , US andMexico, between 1985 and. Comparison of color bar height for each class allows identification of cases where fire occurred disproportionately with regard to biotic and human communities, relative to unburned landscapes more closely investigate these relationships to better understand how antecedent climate and drought influence the geographic patterns of fire across the Madrean Archipelago region.
Vegetation type and fuel moisture content are driven by elevation and topography that create gradients of temperature, precipitation, and evaporation, thus defining unique fire environments and associated fire regimes (Mansuy et al. 2014, Meunier et al. 2014, Perez-Verdin et al. 2014, Whitman et al. 2015. Our results suggest that many of the fire patterns observed in the Madrean Sky Islands occurred at higher elevations with cooler temperatures, received slightly more precipitation, and had fewer frost-free days than sites that remained unburned (Fig. 6a−d). These areas are dominated by Madrean Evergreen Woodlands and Petran Montane Conifer Forests biotic communities that burned in greater proportion on burned versus unburned landscapes (Fig. 5). Historically, these vegetation types produced enough fuels to support surface fires every 8 to 15 years in these middle-elevation forests; higher-elevation forests characterized by longer fire intervals occurred in proportions that were comparable to that of unburned landscapes (Fig. 5). The largest fires we observed were on the US side of the border, including the Horseshoe II Fire and the 2003 Aspen Fire in the Santa Catalinas (32 231 ha). These large fires may represent the interactive effects of regional climate and, in some places, a longer history of fire suppression on fuel accumulation and connectivity (Miller et al. 2009).
Human influence on the regional variability of spatial patterns The importance of understanding the various influences of human activity on fire regimes is increasingly recognized (Mann et al. 2016, Parisien et al. 2016) and takes on a unique relevance in areas where livelihoods are closely dependent on the use of fire (Perez-Verdin et al. 2014). Our geographic analysis revealed patterns of fire Fig. 6 The smoothed density plots illustrate the distribution of area burned (orange) relative to unburned environments (gray) in the Madrean Ecoregion Sky Islands that experienced fire in the time period of the study (1985 to 2011). Variables considered were: a) Annual heat moisture index, b) Mean annual precipitation, c) Mean annual temperature, d) Beginning of frost-free period, e) Mean summer precipitation, f) Lightning flash rate, g) Population density, h) Road density, and i) Terrain ruggedness index. The top graph for each variable was constructed from random points in US burned−unburned locations; the bottom graph used a set of random points for Mexico. Test comparisons indicated rejection of the null hypothesis that burned−unburned samples came from a common distribution (AD.T P-value < 0.001) activity, size, and shape that are likely related to social and cultural dynamics. First, the data indicate greater fire activity surrounding the international border. Specifically, number of fires and area burned was greatest on some of the Sky Islands that traverse the border or are adjacent to it, particularly the Peloncillo-Pan Duro, Huachuca-Patagonia, and Atascosa-Cibuta ranges. This pattern is likely associated with greater human-caused ignitions along major travel corridors. Two of the largest fires that burned during dry 2011, the Horseshoe II Fire and Murphy Complex, were determined to be human-caused by the US Forest Service. Some unique fire patches and patterns occurred near the international border where lack of fuel continuity along the border access road or direct fire suppression resulted in long, straight fire perimeter edges adjacent to the boundary line. These uniquely shaped border fires were characterized by high eccentricity index values (i.e., long and thin) and had strong east to west bearings. Alternatively, fires can be oriented in a particular direction due to prevailing winds during burning, or corresponding to topographic characteristics (Barros et al. 2013).
Secondly, population density was generally similar between burned and unburned sites, with a greater distribution of fires in the low population density range, but with a somewhat bimodal distribution showing a slight increase in fires occurring in the high population density range, indicating the influence of land use and anthropogenic activities on fires for some Sky Islands (Fig. 6g). Remote forests experienced proportionately more fire, suggesting that wildfires dominate the fire regime, rather than human use of fire in agricultural practices (i.e., in Populated Forests; Fig. 5; Perez-Verdin et al. 2014). Furthermore, natural ignition source (i.e., lightning strikes; Fig. 6f) did not limit the occurrence of fire, but areas in the southern Sky Islands with highest strike rate experienced less fire, possibly due to more mesic vegetation types and higher fuel moisture in those islands that receive more precipitation from the south.

Research implications and research needs
Dendrochronological fire scar evidence from many sites in northwestern Mexico shows that some forests have maintained historical fire regimes (Kaib 1998, Stephens et al. 2003, Yocum Kent et al. 2017) and can become reference sites for Western forests in the US that experienced long-term suppression (Leopold 1937, Meunier et al. 2014, potentially guiding ecological restoration of altered sites (Stephens and Fulé 2005). Future research on modern fires in the Madrean Sky Islands could more directly consider the role of historical land use, forest management, and climate on fire attributes in this region. For example, the fire regional severity data (i.e., dNBR) that we developed to map fire perimeters, as measured using pre-and post-fire differencing of satellite imagery, could be used to examine the role of climate, fire weather, management history, and active fire management (e.g., back burns) in resulting spatial heterogeneity of severity patterns. Research into land use history of individual Sky Islands would help to elucidate these possible effects, as would more stand-level field data describing historical fire regimes and current forest structure in sites with different fire histories in northern Mexico.

Conclusions
The purpose of this research was to fill a significant information gap on recent wildfire activity across the bi-national Madrean Archipelago Ecoregion of the United States and Mexico. Past fire research in this region, in particular tree-ring research, suggests that different fire regimes observed in mountains of the US and Mexico were linked to historical land use and fire suppression activities. Understanding the links between these environmental histories and contemporary fire regimes of the Sky Islands can help to guide regional fire planning, fuels management, and habitat conservation efforts on both sides of the US−Mexico border.
Our results show increasing fire sizes across the region, particularly in Sky Islands of the US, with large fires occurring during recent twenty-first century droughts. Area burned generally increased over the study period, culminating with the 2011 fire year, when the Horseshoe II Fire burned much of the Chiricahua Mountain range in Arizona. Number of fires and area burned were greatest on some of the mountains that traverse the border or are adjacent to it, suggesting that recent fire patterns in this area are strongly influenced by patterns of human activity. Two of the largest fires that burned during our study period, the Horseshoe II Fire and Murphy Complex, were determined to be human-caused by the US Forest Service and occurred during the dry 2011 fire season. Given projections of increasing temperature and decreasing precipitation (Wilder et al. 2013), along with projections of continued human population growth  in the borderlands region, large and more frequent fires may become more common in the future. Further analyses of fire occurrence and fire severity in the context of interactions among past and current land uses, regional and local climate variables, and fuel conditions can provide additional important information when considering scenarios of possible future fire regimes in this biologically and culturally diverse region.
Our research provides insight into the modern fire regimes of mountains in the Madrean Archipelago Ecoregion, an area that previously lacked cross-border fire information. A major obstacle to our understanding of the connectivity of ecological, hydrologic, and socio-cultural systems in this region is geospatial data that adhere to political boundaries rather natural ecological or physical boundaries. Recent research on ecosystem service flows in the US−Mexico borderlands highlights the challenges and importance of a transboundary approach (López-Hoffman et al. 2010, Norman et al. 2012. Transboundary data are particularly important to consider when managing habitat for migratory species and large predators (e.g., Panthera onca, Linnaeus; McCain and Childs 2008, Atwood et al. 2011), and for many other ecological processes that are not confined to political boundaries. An increased focus on developing and sharing transboundary datasets can contribute to collaborative research efforts between countries that improve knowledge about ecological patterns and processes in the borderlands, and support long-term regional planning and restoration efforts. Table 4 Anderson-Darling statistical output for comparisons of burned and unburned samples (see Fig. 6 and Table 3), for the Madrean Ecoregion Sky Islands, US and Mexico. The standard deviation for the Anderson-Darling test statistic (AD.T) was 0.76, and the null hypothesis that all samples came from a common population was rejected in all comparisons (P < 0.001). Variables with density plots in Fig. 6 are shown in bold. Number of ties refers to the number of times identical values occur in the data. If too many ties occur, there will be a rejection of the null hypothesis and the result will suggest non-normality Publisher's Note