Next Article in Journal
Leaves and Twigs Image Recognition Based on Deep Learning and Combined Classifier Algorithms
Previous Article in Journal
Automatic Extraction of Forest Inventory Variables at the Tree Level by Using Smartphone Images to Construct a Three-Dimensional Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multiproxy Approach to Reconstruct the Fire History of Araucaria araucana Forests in the Nahuelbuta Coastal Range, Chile

by
Ariel A. Muñoz
1,2,3,4,*,
Mauro E. González
3,5,6,
Isadora Schneider-Valenzuela
1,2,
Karin Klock-Barría
1,7,
Marcelo Madariaga-Burgos
1,2,
Carmen Gloria Rodríguez
5,
Ana M. Abarzúa
8,
María Eugenia Solari
9,
Alejandra Martel-Cea
8,10,
Bárbara Velásquez
8,
Beatriz Paredes
8,
Fabián Guerrero
11,
Mauricio Montiel
7,
Valeria Tapia-Marzán
1,
Tomás Riquelme
7 and
Paul R. Sheppard
12
1
Laboratorio de Dendrocronología y Estudios Ambientales, Instituto de Geografía, Pontificia Universidad Católica de Valparaíso, Valparaíso 2362807, Chile
2
Centro de Acción Climática, Pontificia Universidad Católica de Valparaíso, Valparaíso 2362807, Chile
3
Center for Climate and Resilience Research, CR2, Santiago 8370449, Chile
4
Fundación San Ignacio del Huinay, Hualaihué 5860000, Chile
5
Laboratorio de Ecología de Bosques, Laboratorio de Dendrocronología y Cambio Global, Instituto de Conservación, Biodiversidad y Territorio, Universidad Austral de Chile, Valdivia 5090000, Chile
6
Center of Fire and Socioecosystem Resilience (FireSES), Valdivia 5090000, Chile
7
Escuela de Graduados, Facultad de Ciencias Forestales y Recursos Naturales, Universidad Austral de Chile, Valdivia 5090000, Chile
8
Laboratorio de Palinología y Reconstrucciones Ambientales, Instituto Ciencias de la Tierra, Universidad Austral de Chile, Valdivia 5090000, Chile
9
Laboratorio de Arqueobotánica e Historia Ambiental, I. de Estudios Antropológicos, Universidad Austral de Chile, Valdivia 5090000, Chile
10
Centro de Estudios Avanzados en Zonas Áridas (CEAZA), La Serena 1780000, Chile
11
Department of Mechanical Engineering, Universidad Técnica Federico Santa María, Avenida España 1680, Valparaíso 2340000, Chile
12
Laboratory of Tree-Ring Research, 407 Bryant Bannister Tree-Ring Building, University of Arizona, Tucson, AZ 85721, USA
*
Author to whom correspondence should be addressed.
Forests 2023, 14(6), 1082; https://doi.org/10.3390/f14061082
Submission received: 1 April 2023 / Revised: 13 May 2023 / Accepted: 14 May 2023 / Published: 24 May 2023
(This article belongs to the Section Forest Ecology and Management)

Abstract

:
Multiproxy reconstructions of fire regimes in forest ecosystems can provide a clearer understanding of past fire activity and circumvent some limitations of single proxy reconstructions. While inferring fire history from scars in trees is the most precise method to reconstruct temporal fire patterns, this method is limited in Araucaria araucana forests by rot after fire injuries, successive fires that destroy the evidence and the prohibition of sample extraction from living Araucaria trees. In this context, dendrochemical studies in Araucaria trees and charcoal analysis from sediment cores can complement and extend the time perspective of the fire history in the relictual Araucaria-Nothofagus forests of the coastal range. We used dendrochemical, fire scar and charcoal records from the Nahuelbuta Coastal Range (37.8° S; 73° W) spanning the last 1000 years to reconstruct the fire history. The results indicate that periods with higher fire activity occurred between 1400 and 1650 AD. Long-term changes in the fire regime are related to increased climate variability over the last 1000 years, and especially with the arrival of settlers to the area after 1860 CE. The most severe fire events in the Nothofagus and Araucaria forests occurred when suitable fire-prone conditions were superimposed with high human densities.

1. Introduction

Fire history is critical for understanding the potential effects of changing climate and human land use on the forest landscape [1,2,3,4]. Under the current climate change scenario, south-central Chile (37°–42° S), characterized by sharp declines in rainfall and increases in the recurrence of drought and heat waves [5,6,7,8], is expected to see an increase in the recurrence and magnitude of forest fires [9,10,11,12,13,14,15,16]. Fire occurrence in this area has also been linked to the climatic effects of Southern Annular Mode (SAM) [17,18] as well as to El Niño Southern Oscillation (ENSO) variability [17,19,20].
Fire is a major disturbance process shaping the Araucaria forest landscape in the Andean and Coastal ranges of north-west Patagonia [18,21,22,23]. Araucaria-Nothofagus forest ecosystems have evolved with fire and some plant species, such as Araucaria araucana (Molina) K. Koch (Araucaria), have developed distinct strategies to cope with this process [24,25,26,27]. The fire regime of Araucaria forests is characterized as mixed severity, where fires may leave behind unburned patches along with patches burned with a low, moderate and high severity (>75% overstory tree kill) [28]. Low-severity surface fires can occur at relatively short intervals and are generally associated with more open sites dominated by Araucaria mixed with Nothofagus antarctica (G. Forster) Oerst. woodlands. Stand-replacing crown fires are common in more productive mesic mixed species stands of Araucaria–Nothofagus forests where emergent, thick-barked Araucaria trees frequently survive [29].
Over recent decades, Araucaria-Nothofagus forests have been affected by large high-severity fires, both in protected areas (national parks and reserves) and on private land [18,21,27,30,31]. Although an increase in the recurrence of fires in Andean Araucaria forests has been documented since the Euro-Chilean period (1883–1970) [18,21], fire events occurring in the last 50 years have been considered exceptional in terms of their frequency, area burned and fire severity [28,30].
Araucaria araucana is a tree species of recognized cultural, scientific and socio-economic importance. This long-lived conifer is endemic to the temperate forests of Chile and Argentina [24,25,32,33], and classified as Vulnerable (VU) in its Andean distribution and Endangered (EN) in its relic distribution in the Nahuelbuta Coastal Range [34,35,36]. It is estimated that the original area of Araucaria forests at the arrival of European colonizers was about 500,000 hectares, but this number was reduced by half by the end of the 20th century, with only 254,000 hectares remaining [37]. The effects of the European-Chilean colonization of southern Chile, and subsequent logging in the Araucanía region in the 20th century, triggered a drastic reduction in forests, leaving those remaining in a critical state [33]. Despite being protected by Chilean law (Decreto N° 43 del Ministerio de Agricultura 1990, which prohibits the cutting and logging of these forests), they are currently at a high risk of degradation, either by illegal logging [19], cattle grazing [26,27,38], the invasion of exotic species [39,40], seed predation by non-native species and the intensive collection of edible seeds [25], or fire-induced disturbances [30,41].
Several studies regarding the fire regime of A. araucana forests have been conducted in the Andean part of its distribution [18,21,28,29], but little is known about the historical range of fire variability in the populations growing in the Coastal Mountain Range, especially in the Nahuelbuta Range. This is of special interest, since this area has a different history of occupation than the Andes, in addition to not having volcanic activity, which can be linked in some degree to the fire activity in the Andes. Knowing the fire regime in Araucaria forests in different areas of its distribution, especially in the Nahuelbuta Coastal Range, would allow us to assess the different degrees of disturbance and fire dynamics in these forests to inform management and restoration plans for these unique ecosystems in protected and private areas.
Fire history has been reconstructed through various environmental proxies, such as fire scars on fire-affected trees, stand age classes, increases or decreases in tree growth ring width, macroscopic charcoal particles and pollen and spores in peatbog sediments, among others [42,43,44,45].
The study of fire scars is based on the analysis of the growth rings of trees in segments of the trunk where the fire burned through the bark causing the death of the cambium and a subsequent scarring of the damaged area [46,47]. These scars make it possible to date fires with an annual or in some cases seasonal resolution throughout the life of the tree [46,47]. Often after fire, an abrupt radial growth response (release) occurs associated with the healing callus [47]. However, due to the low presence of old trees (>500 years), few old fire records can be dated. In addition, fire sensitivity differs amongst trees and tree species, so not all store records of fire in their wood anatomy, and a loss of records after successive large fires over time has been documented, precluding a complete fire reconstruction [21,46,48].
Fossil charcoal, on the other hand, is the most frequently used technique to reconstruct fire history from lake deposits, as it derives from the incomplete combustion of plant tissue [49]. Charcoal dating in lake sediments has been used to date fires in Chile up to 26,000 years old [50]. It can cover larger burned areas than fire scar records, especially because the macroscopic charcoal particles from fires inside a catchment area are deposited just as sediment from that area is deposited in fluvial and lacustrine environments [51]. Thus, this type of record catalogues fire events on larger spatial and temporal scales that are not necessarily close to a given study site (patch scale). The temporal resolution depends on the sedimentation rates of each lake or swamp; however, compared to the ring chronology yearly resolution, it is only possible to date fire events at a decadal to centennial to millennial time resolution [51].
These paleoenvironmental proxies are by definition an approximation of reality and are affected by multiple processes [44]; hence, it is important to obtain multiproxy information to enrich and complement these records and to develop a better understanding of the fire regimes over time.
Trees have been shown to chemically respond to variations in the environment by retaining chemical elements in their growth rings [52], the concentrations of which can be indicators of environmental change over time [53,54]. Dendrochemical studies have been successfully utilized to study trace metal deposition in tree rings and have been widely used to monitor environmental changes in soil and atmosphere [55]. Alternative dendrochemical techniques have been used to study the effects of atmospheric pollution [56,57,58] and volcanic eruptions [59,60,61,62]. Wildfires generate large numbers of pollutants, including trace elements [63,64]. During fire events, several elements (macro and trace nutrients) are released, deposited in the soil [63,64] and then absorbed and stored by trees in their woody tissues [65]. In effect, some studies show a variation in the chemical composition of wood surrounding fire scar samples [60,66,67]. Moreover, changes in the concentration of N and P in tree rings before and after fire events have been described [68], and an increase in the concentrations of various elements has been observed in tree rings for up to nine years after a fire [60]. Increases of B, Al, Mn and Fe have also been found to be useful as indicators of fire occurrence, further corroborating the potential of dendrochemistry for studying fire history [69]. However, the mechanisms by which the elements are fixed inside tree rings are not yet fully understood and show mismatching results [54], making dendrochemistry one of the most challenging lines of dendrochronology at the present time. One way to confirm dendrochemical results is to calibrate them with other proxy records, which is especially needed to assess the potential of dendrochemical results to detect past fires and other disturbance regimes.
The dendrochemical potential of A. araucana has been recently explored to date volcanic eruptions in the Andes Mountains, showing promising results that indicate the species’ ability to record pulses of chemical elements released in eruptions of varying explosivity [70,71]. Ref. [72] point out that evergreen coniferous species may be more sensitive to chemical changes in the environment than deciduous species. Araucaria, being a long-lived evergreen conifer with individuals exceeding 1000 years of age [32], can potentially provide long dendrochemical records of past environmental conditions. Our hypothesis in this work is that the smoke and particles emitted in forest fires could alter A. araucana’s uptake of chemical elements through its leaves or roots and be reflected in the concentrations of these elements in its growth rings through time. This approach could complement the information provided by other environmental proxies used to study fire dynamics, such as fire scars and sedimentary micro-carbon. Furthermore, unlike fire scar sampling, dendrochemistry uses a less invasive technique for sample acquisition; therefore, it can be a very useful tool in those places where scar samples cannot be taken (e.g., national parks).
In this context, this work is the first to evaluate the use of dendrochemical records to study fires in Araucaria araucana forests by comparing them with records from fire scars on trees and charcoal records in sediment cores. This information may be key to achieving high-resolution fire records in Araucaria forests for the last millennium, helping to clarify and better interpret both the impact of the arrival of Spaniard conquistadors on the fire regime in the Araucanía region and to identify the periods of high climatic variability and the influence of high temperatures on the Araucaria forest landscape.

2. Materials and Methods

2.1. Study Area

Our study area was the Nahuelbuta Coastal Range (37°–38°50′ S, Figure 1), located between the Bío-Bío and Imperial rivers, with an extension of 175 km in length [73] and a maximum elevation of 1525 m a.s.l. (Figure 1A). The Nahuelbuta Coastal Range is a biogeographically relevant place since it remained free of periglacial processes during the glacial events of the Quaternary, developing a more stable biogeographic history compared to the Andes [73]. This allowed a rich diversity of narrow-range endemic flora and fauna to concentrate in this mountain range, with some vascular plants whose only distribution is in this mountain range [73]. Moreover, it represents the predominant distribution of the Araucaria forest on the coast. This mountain range is characterized by a maritime climate with a Mediterranean influence, reflected in a distribution of maximum rainfall in winter and relatively dry and hot conditions during the summer months (December–March) [35,74]. The amounts of annual precipitation vary between 1500 and 3000 mm, influenced by snow at high altitudes, while the average annual temperature values have a wide variability, with minimum temperatures in July and maximum temperatures in February [38]. In general, soils derive from metamorphic material and in some sectors from granite [38,75]. Unlike the Andes, this mountain range has no volcanic influence, which is why the fire signals can only be due to forest fires, allowing us to remove some of the noise in the records.
During prehistoric times, the valleys of the Nahuelbuta Coastal Range presented different socioeconomic strategies of land use, from the nomadic hunter-gatherers of the early and middle Holocene to the complex agro-pottery systems of the late Holocene (1000–1500 CE), with the presence of ceremonial mounds (kuel) in valleys such as Purén and Lumaco [76]. This gradual sociocultural complexification of lowland areas parallels the continuous use of high-altitude ecosystems and their Araucaria forests, which for the Mapuche people ancestrally correspond to a culturally and ritually significant space in which the gathering of species such as Araucaria fruit (piñón) and hunting endure over time [77,78].
Environmental history shows that since 1550 CE, with the arrival of the Spaniards in the territory occupied by the Mapuche people, a colonization process began in which the high-altitude ecosystems were used by the Mapuches as shelter during wars and skirmishes against the Europeans (16th-early–19th centuries) and, in turn, as grazing areas for introduced Euro-Mediterranean livestock [79,80,81]. Since the end of the 19th century, the timber exploitation of alpine forests has enabled the construction of settlements in the valleys and the construction of a large-scale economy based on wheat cultivation [82,83,84].

2.2. Dendrochronology Field Sampling

Dendrochemistry: A total of 28 Araucaria araucana trees were sampled in April 2021 from the Nahuel site in the Nahuelbuta Coastal Range (NHL, Figure 1B). From each tree, two cores were extracted using 5 mm increment borers (https://haglofsweden.com/project/increment-borers/ (accessed on 30 March 2023)). The chosen trees were healthy, had no fire scars and were growing in shallow soil or rocky substrate. Additionally, data on the diameter at breast height (DBH) and a GPS point at each tree were taken.
Fire scars: Forest patches in the two sampling sites (Figure 1A–C) were intensively searched for fire-scarred A. araucana and Nothofagus sp. trees and, whenever possible, samples were collected in clusters of several trees to improve the chances of obtaining the most complete fire record possible and increase the likelihood of precisely dating the fires [85]. On standing or fallen dead trees, cross-sections of single and multiple fire scars were extracted with a chainsaw (Figure 1C). For live trees, 2 to 4 cores were extracted close to the fire scar with an increment borer to determine fire dates [86]. In addition, the species, DBH, number of visible fire scars, scar height and azimuth of the scar face and GPS point of each tree were recorded. In total, for the two sampling sites, COI and NHL (Figure 1), 10 and 12 fire-scarred trees were collected, respectively.

2.3. Tree Ring Sample Preparation, Dating and Measurement

Following standard dendrochronological procedures, tree cores and cross-section samples were left to dry, then mounted and sanded with sandpaper of different granulometry (100 to 1000 units) [87,88]. Tree ring samples were then visually cross dated using a stereomicroscope, following the Schulman convention for the Southern Hemisphere [89]. Ring width was measured using a Velmex measuring system (±0.001 mm) (https://www.velmex.com/Products/Pre-configured_Systems/Tree_Ring_Measuring_System.html (accessed on 30 March 2023), and ring width series cross dated using COFECHA software to ensure correct dating and measurement [90]. Ring width series were standardized by applying a negative exponential curve to remove age-related growth effects while building a tree ring chronology for the site using ARSTAN software [91]. This tree ring chronology served as a reference series to precisely date the fire scars of dead tree samples. These scars were identified by the characteristic formation of a cambial death lesion and the healing patterns of radial tree ring growth (Figure 1C) [42].

2.4. Climate-Growth Analysis

In order to identify relationships between climate variables (precipitation and temperature) and tree growth, the resulting chronology NHL was compared with precipitation values obtained from the TerraClimate data set [92] for the period 1958–2020 and mean temperature retrieved from KNMI Climate Explorer (https://climexp.knmi.nl/start.cgi (accessed on 6 November 2022), available for the period 1950–2020. Correlation coefficients were calculated for moving periods, from 1 to 12 months, starting from the second previous growing season to the current dendrochronological year in the Southern Hemisphere.

2.5. Chemical Measurements

The growth series from living trees that made up the NHL chronology (Figure 1D,F) with the best crossdating statistics were selected for chemical measurements in order to develop a dendrochemical chronology for the period 1351–2020. Due to the small size of the Araucaria tree rings (Figure 1D) and minimum mass requirements for chemical analysis (50 mg of sample), this was performed by pooling two consecutive annual rings for the entire period and by creating a composite chemical chronology of selected series. Although the individual signal of each tree is lost with this approach, it allows obtaining a common signal of the group of trees, as well as reducing costs and covering longer periods of time [62]. Thus, we selected 32 tree ring series that jointly covered the entire period under study to obtain a 2-year resolution chemical chronology.
Under a binocular microscope, each core was cut every two years using a sterile non-metallic, ceramic knife over an acrylic table. In order to avoid possible cross-contamination, all segments were then polished with a diamond tool to remove all the surface that was in contact with extraction and sanding tools. Isopropyl alcohol was used to clean all materials. All 2-year samples from the 32 series were then placed in labeled tubes to form composite samples for each biennium (Figure 2A), obtaining a total of 335 samples covering the 1351–2020 period. Processed samples were subsequently sent to the Arizona Laboratory for Emerging Contaminants (ALEC) at the University of Arizona for chemical measurements (https://www.alec.arizona.edu/ (accessed on 15 March 2022)) (Figure 2).
The samples were prepared via acid digestion. Wood samples were weighed in precleaned, pre-weighed, trace metal-free polypropylene centrifuge tubes. Samples were then predigested with 3 mL of concentrated Optima grade nitric acid (HNO3) at room temperature for 24 h and then heated at 90 °C in a digestion block for 2 h to complete the digestion (Figure 2B). Sample tubes were then reweighed to calculate exact dilution factors. A 1.5 mL aliquot of digestate was gravimetrically diluted by a factor of approximately 10 with ultrapure 18.2-MΩ/cm water from a Milli-Q system (Merck, Darmstadt, Germany) (Figure 2C,D). The liquid samples (Figure 2E) were analyzed by an inductively coupled plasma mass spectrometry (ICP-MS; 7700) system (Agilent, Santa Clara, CA, USA). Major cations (potassium (K), calcium (Ca), sodium (Na), magnesium (Mg)) and trace metals (beryllium (Be), aluminum (Al), manganese (Mn), nickel (Ni), copper (Cu), zinc (Zn), molybdenum (Mo), silver (Ag), cadmium (Cd), tin (Sn), antimony (Sb), barium (Ba) and lead (Pb)) were measured, along with standards (Pine Needles (SRM 1575a) and one blank sample. The calibration of ICP-MS was performed using Multi-element Standard Solution of Spex CertiPrep (Metuchen, NJ, USA) as well as an independent quality control check from High-Purity Standards (Charleston, SC, USA). Results were expressed as μg g−1 of dry weight of wood.

2.6. Dendrochemical Data Analysis

The biennial time series of chemical element concentrations recorded in A. araucana tree rings (Figure S1A,B) were individually detrended and standardized using a spline curve with a 50% frequency cutoff at 70 years in order to remove the magnitude effect of the chemical concentrations and to be able to compare the time series concentrations between them. Superposed Epoch Analysis (SEA) was used to determine the relationship between fire events and chemical element concentrations in the years before, during and after the event. The mean values of these data were calculated for windows of 4 biennia, i.e., 8 years, including the biennium of the fire event. The mean values of the standardized chemical series during the years of the fire event were compared with variations in the entire record by performing 1000 Monte Carlo simulations that randomly select years, calculate expected averages and provide 95% confidence intervals [93,94]. In each case, the number of randomly selected years is equal to the number of actual fire years. In this way, we compare the chemical signals to a subselection of fire dates recorded on at least two site trees from the sampling sites. All statistical analyses were conducted using R [95], RStudio [96] and the dplR package [97,98].
Finally, in order to propose new fire dates prior to the records obtained from fire scars, we selected the chemical elements that showed a clearer response to the fire from the Superposed Epoch Analysis. The proposed fire dates were selected as those in which at least 3 elements of the subset of previously selected elements had a positive chemical pulse above the 90th percentile, either during the same biennium or after two biennia. In addition, we developed a Principal Component Analysis (PCA) for eight selected elements (Figure S2) according to the SEA results to compare with carbon records from sediment cores. In addition, we used a temperature [99] and PDSI [100] reconstruction to compare the fire record estimate from dendrochemistry and charcoal analyses.

2.7. Charcoal Record

Two sediment cores were collected during January 2021 from Aguas Calientes (ACA) peatbog using a Russian corer (Figure 1E). Sediments were described and dated using 14C. Age–depth model performed by Bayesian analyses in “rPlum” package [101,102] in R platform [95]. Six radiocarbon dates (bulk, plant tissue, charcoal) were analyzed in the University of Arizona LTRR and AMS Laboratory and calibrated with the SHcal20 curve [103] (Figure S3). One date was not included in the age–depth model due to its chronostratigraphic error (Table S1).
To infer fire regime changes at Nahuelbuta National Park we used macroscopic charcoal particles (>125 μm for fire peak detection in a radius of <3 km from coring site; [104]). For macroscopic charcoal counts, 2 cm3 of sediment was taken at 1 cm intervals for the length of the core, sieved through 125 μm screen and tallied under a stereomicroscope at 10–40× following the methodology outlined by [51]. Statistical analyses were performed to detect fire peaks using CHARanalysis software [105]. For that purpose, we interpolated samples at 20 year intervals and calculated the background component using a moving average window with a width of 1000 year and a threshold ratio value of 1.1 to identify fire episodes above background.

3. Results and Discussion

3.1. Fire History of Araucaria-Nothofagus Forests

The fire history of the two sampling sites (NHL and COI) was reconstructed from 19 cross-dated samples, recording 18 fire years in nearly 500 years. Most fire-scarred samples (58%) recorded a single fire scar, while 26% recorded two–three fire scars. The fire chronology developed in COI spanned 468 years, and NHL spanned 478 years. The oldest dated fire was in 1738 and the most recent one occurred in 1981 (Figure 3). The frequency of fire in the last three centuries has increased; in the 18th, 19th and 20th centuries, four, five and nine fire dates, respectively, have been recorded from at least one tree (Figure 3). In this record, seven fires were documented by at least two trees; the majority of those events (five fires) occurred after 1900 (Figure 3).
The patterns found in these two sites are like other studies carried out in Araucaria forests in Chile and Argentina. In the Andean Araucaria distribution, during the period of the Euro-Chilean settlements, fire frequency increased significantly [18,21,22,106]. Perhaps counterintuitively, the establishment of protected areas such as the Nahuelbuta National Park, established in 1939, and the Malleco National Reserve, created in 1907, have not reduced fire frequency in these ecosystems. Only after the creation of the National Forest Service (CONAF) in the 1970s, and the implementation of a fire suppression policy, did fire frequency in these protected areas diminish [107]. A similar result was also observed in Tolhuaca and Villarrica National Parks [18,21]. On the other side of the Andes Range, in Argentina, more effective fire suppression policies were implemented in Araucaria forests after the creation of protected areas [22]. Despite the differing fire management strategies on either side of the border, some fire years during the twentieth century, such as the 1944 and 1959 fires at the Nahuelbuta sites and periods of high fire frequency in 1890–1900 and 1760–1770, have also been found in other Araucaria forests in the Chilean and Argentine Andean distribution, suggesting a climatic link between some of these fires [18,22]. On the other hand, the fires identified in the mid-nineteenth century could be related to the recurrent clashes in this area between the military and Mapuche people, who refused to incorporate their territory into the Chilean State [108]. The Nahuelbuta Mountain Range was very important from a military strategic point of view as it constituted a natural boundary from the settlements, and its thick forests and surrounding areas provided refuge for Mapuche people. In the clashes between indigenous people and the Chilean military, the use of fire was a constant (fires in agricultural fields, pastures and settlements in general) and was at its peak between 1852 and 1883, which may partly explain the occurrence of fires in this area, considering that the dates of the fire scars do not coincide with generalized climatic events [108,109,110].

3.2. Ring Width Chronology of Araucaria araucana

The Ring Width Index chronology of the NHL site used to correctly date fire scars included a total of 28 trees of A. araucana, spanning 921 years for the period 1100–2020 with an interseries correlation of 0.501 (Figure 4A). This correlation value is similar to those seen in other chronologies of this species [22,111,112]. The biggest increases in tree growth were observed in 1414 and between 1956 and 1958, while its lowest growths were observed in 1627 and 1896. A very noticeable feature in this chronology is the positive growth trend between the years 1946 and 1956, which is the most important in the whole chronology (Figure S1 and Figure 1D). This release was recorded in most of the trees and appears to be peculiar to this particular site when compared with other chronologies of the species in Nahuelbuta National Park (Figure S4) and seems to be related to fire occurrence (Figure 4A) and logging activities documented in this area for this period. At the NHL site, we also found an increase in variability and releases at the beginning of the chronology, around 1400–1550, which were not associated with decreases or important variations in sample size (Figure 4A).
The growth patterns of A. araucana in the NHL site show a positive correlation with the precipitation values of the summer prior to ring formation and also with the current growing season. Growth shows a negative correlation with temperature in the spring and summer of the current growing season (Figure 4B). The highest values of correlation between ring width and a climatic variable at NHL were with accumulated precipitation between past November and current January and with the average temperature of August to October in the current year (r = 0.33 (1958–2020) and r = −0.31 (1950–2020)) (Figure 4B). This climate–growth relationship follows the regional and also local tree-growth patterns shown by [111,112] for the Araucaria population across the distribution of Araucaria araucana in Chile and Argentina, including the Piedra del Aguila and Pichinahuel sites, also located in the Nahuelbuta National Park. The Nahuelbuta National Park Araucaria ring width chronologies are some of the longest in the Araucaria network in Chile and Argentina. A SEA using the ring width chronology and selected fire dates revealed negative growth anomalies one year before the fire events (Figure 4C), suggesting that dry conditions, resulting in the formation of small rings, are likely to promote fire events.

3.3. Dendrochemistry and Fire Scar Dates

The concentrations of various chemical elements were measured and were used to determine if there was a relationship with the other fire study proxies, such as fire scars and charcoal in sediments. Figure 5 and Figure 6 show fire scar dates and the detrended and standardized chemical time series since 1700 for all chemical elements analyzed except Be, which was not detected in any sample. The responses of the chemical elements to fire events varied between elements and fire dates recorded in two or more trees. By comparing the dendrochemical data set with these fire dates (1820, 1854, 1900, 1925, 1942, 1944, 1981), we can see that:
1. Some elements increased in concentration during the year of the fire or 2 to 8 years after the fire dates. Positive anomalies (defined as concentration values of >=5 elements above the 90th percentile) were observed in 11 elements (Ba, Ca, Mg, Mn, Cu, Na, Ni, Al, Mo, Pb, Zn and Sb) for the biennial period 1947–1948 (Figure 5 and Figure 6), after the fire events recorded at the NHL site in 1942 and 1944, which could be affecting the chemical responses of elements in tree rings at the site. An increase in tree growth was also observed in the NHL site for this period (Figure 4A). For the fire recorded in 1900 at the COI site, increasing peaks of Cu, Na, Ni, Al and Pb were found in the same biennial period of the fire. On the other hand, the fire that occurred in 1820, also at the COI site, showed positive anomalies two biennial periods after the fire events in Ba, Ca, Mg, Mn, Ni and Zn. Nevertheless, we could not find significant positive anomalies for more than five elements at the same time for the fires occurring in 1854, 1925 and 1981 (Figure 5 and Figure 6).
2. A few elements increased their concentrations in the years prior to fire events. This is true for the year 1900, where Cu, Na, Ni, Al and Pb increased their concentration 1 biennium prior to the fire date. Other fire dates that show this behavior in trace elements are 1912 (same elements as 1900), 1869 (Ba, Ca, Mg, Cu and K) and 1738 (Ba, Ca, Mg, Mn, Cu, Ni and Mo), corresponding to fire dates detected only in one tree.
The SEA carried out to assess the response pattern of significant changes in chemical elements related to the fire events recorded by two or more trees (1820, 1854, 1900, 1925, 1942, 1944, 1981, Figure 3) showed positive significant (p > 0.05) anomalies for Mg, Ca, Sn, Mo and Ba after the fire dates for up to 3 biennia (Figure 7). Sn (p > 0.05), Ag (p > 0.05), Mn (p > 0.05), Ba (p > 0.05), Na (p > 0.1) and Pb (p > 0.1) showed positive significant anomalies in the same biennial period of the fire event (Figure 7). On the other hand, Al and Zn showed positive anomalies a two-year period before the fire (p > 0.05), as well as Ag (p > 0.1) (Figure 7). The remaining elements (i.e., K, Ni, Cu, Cd and Sb) did not show significant anomalies related to fire events (Figure 7).
Positive anomalies in chemical series and SEA somehow demonstrate that the main response of most elements in tree rings after fire is an increase in concentrations following fire events. It is widely known that fire remobilizes trace elements in the environment, increasing their concentration in soils and sediments [63,64]. Through the uptake, discrimination and fixation mechanisms [65] can be different; these trace elements can be incorporated into the wood tissues during the same year of growth or be delayed by a few years because of differential incorporation rates [65]. Furthermore, trace elements have been shown to remobilize and translocate between tree rings [53,113]; this could also be an explanation for why some elements show increases in the years prior to the fire events. Another explanation for prior chemical increases could be fire dates potentially were not captured by the fire scar sampling in Nahuelbuta forest, especially because currently it is not permitted to take fire scar samples from living trees, as has been specified previously for other Araucaria forests. This could explain why we did not find fire records around the end of 19th century that were found in other fire histories of Araucaria forests in Chile and Argentina [18,22]. On the other hand, a study conducted in an Araucaria-Nothofagus forest in the Tolhuaca National Park in the Andes Mountains, affected by a high-severity fire in 2002 [114], found that there were elevated concentrations of Mg, Ba, Mn and Ca three biennia before the fire related to the presence of cells occluded with phenolic compounds, typical of the effect of chemical compartmentalization. This chemical response coincides with those found by [66,67] near the healing callus where these occlusions occur as resin soaks before and after fires [115]. However, in the current study, we tried to avoid erroneous indicators of possible fire events from the Araucaria dendrochemical records by: (a) selecting only trees without fire scars, and (b) selecting fire dates where at least three or more chemical elements showed a clear increase in the chemical concentration, which is best explained in Section 3.4 of this study.
Some elements in tree rings have been proven to record chemical changes associated with natural or anthropogenic activities in the same site and also close to the site where the trees came from [56,62,116]. In the case of our dendrochemical results in the Nahuelbuta Range, we found an interesting match between the increase in concentrations of several elements in the Araucarias of the NHL site and the presence of fire scars there and in a nearby site (COI) as well as with carbon abundance from nearby lake sediments cores (Figure 1A,B). Dendrochemical records integrate chemical changes in the local area and also a wide geographical range surrounding the sample site and its environment. It is possible this is due to the potential of wildfires to mobilize chemical elements in ecosystems that are then spread by smoke far downwind [64]. However, several fires were not detected by concentrations of elements in tree rings. A possible explanation could be a low fire intensity or severity. It is highly probable that more severe fires are more easily detected in dendrochemical records than small ones because of the magnitude of remobilized trace elements in the environment [63]. Another explanation for the variation in element concentrations in trees that survive fire could be explained by variability in which types of elements are mobilized in each individual fire [66,117]. It is known that fires significantly increase nutrient availability, such as Ca, Mg, K and Na, associated with ash deposition and the combustion of organic matter, which increase their concentrations in soils with increasing fire intensity [118]. However, these increases in concentration of exchangeable cations can be short-lived in soils due to the high volatilization and evapotranspiration rates associated with new recruitment and remnant plants after fires [118]. Furthermore, [119] analyzed the concentrations of several trace elements in soils and ash after forest fires in eucalyptus and pine plantations in Portugal, finding high levels of V, Mn, Ni, Cd and Pb in burned soils, while Co and Cu concentrations did not change after fire. Even though increases were reported following fire events, some of them (i.e., Mn and Cd) showed fast reductions after the first rain events, while V, Co and Ni increased through the first eight months after the fire event, and Cu and Pb showed small changes through time. Soil acidity also plays an important role in controlling element availability for plants and other ecosystem processes. Ref. [120] reported Ca, Mg and Mn were significantly correlated to the acidity status of the soil. They suggest that if the elements correlate significantly to a common factor, then they might also correlate with each other. This also can explain the common signal between some elements after fires in the Araucaria Nahuelbuta forest.
Dendrochemical studies have been widely used to reconstruct trace metal deposition from various pollution sources, and its use in dendrochemical monitoring seems promising [55]. There is evidence that the chemical composition of wood is affected by fire, showing changes in element concentrations surrounding fire scar samples [60,66,67], as is the case for the present study. Hence, it may be possible to use some of these elements found in this study to extend the fire history of the Nahuelbuta Coastal Range based on tree ring dendrochemical records. However, the chemical responses to fire dates were not always the same, so the interpretation of these conclusions and those of dendrochemical studies in general should be cautiously applied [121] due to the recurring dendrochemical uncertainties described over the years by several authors (e.g., [53,54,113,121]). These are, but are not limited to: (1) the variety of uptake paths for chemical elements, uses and movement within the xylem and the trunk [53,121], (2) the variability between trees, sites and species [67,122], (3) the climate and environmental conditions potentially related to the uptake of some chemical elements [123,124] and (4) the age of the trees, structure and composition of the forest and canopy [125] and the pH of soil, which play a significant role as well [54].

3.4. Extending the Fire History of Nahuelbuta Range Using Dendrochemistry and Charcoal Records from Sediment Cores

Using a set of chemical elements with pulses of concentration over the 90th percentile, which also match with three other elements in the same condition, we propose a new set of fire periods (Figure 8). The new fire dates proposed include, in most of the cases, chemical elements considered relatively less mobile in the trunk, and avoiding groups only represented by chemical elements described as easily mobile (K, Mn, Cu, Zn) [72]. From this analysis, at least three periods where various elements presented pulses of chemical elements that exceeded their p90 were identified. The first, between 1363 and 1404, presents pulses in all its elements at least once, except Na. The second period extends between 1455 and 1552, and a third period between 1601 and 1700, with concentration pulses of all chemical elements above their 90th percentile (Figure 8). These results suggest the occurrence of periods of greater fire activity in the past than in the calibration period of scars and fires (1700–2020; Figure 8).
These results were also compared with regional climate reconstructions and charcoal records from sediment cores obtained in nearby lakes to the dendrochemical site and the fire scars. This analysis suggests that there are important coincidences between periods of greater presence of carbon in sediments (Figure 9A), periods of high temperatures and dry conditions (Figure 9B), where there are also pulses of chemical elements in Araucaria trees (Figure 9C), demonstrating the consistency of these records and corroborating the influence of the climate on the occurrence of fires. This information constitutes the first time that dendrochemical records calibrated with fire scars have provided information to reconstruct fire regimes, which in turn coincide with charcoal records in lake sediments (Figure 9).
The results from Aguas Calientes charcoal records show that the greatest accumulation occurs between the years 1200 and 1650 AD, distributed with three maximums or peaks recorded in the years 1200, 1380 and 1600 AD (Figure 9A). These charcoal particle maximums coincide with warm and dry periods, from regional climate reconstructions developed with tree rings (Figure 9B), and also with increases in the concentrations of chemical elements in Araucaria trees at the NHL site, which shows its greatest activity (higher chemical pulses) between 1400–1650, confirming greater climatic variability in these periods.
These small differences in the timing of fire activity may be associated with the sedimentation rates and chronological control of the Aguas Calientes records; one radiocarbon date of 525 ± 20 years 14C AP located in the AC2101AT1 core, 33 cm deep, supports the age–depth model for this section of the sediments. This dating is associated with the maximum amount of charcoal registered in the year 1400 AD, which allows us to establish that this moment was the largest fire registered in Nahuelbuta in the last thousand years. It should also be noted that, given the nature of the charcoal data, these fire events are expressed in the sediments as a Gaussian curve, typically caused by the resuspension and reaccumulation of sediments after one fire event. In this sense, wetlands or swamps are excellent sites for the detection of fire events given the little remobilization of sediments that occurs there [126]. More radiocarbon dating in the sediments, and a large number of trees, sites and chemical elements, could improve the statistical relationship between these two proxies (dendrochemistry and charcoal record). Moreover, other paleoenvironmental indicators, such as palynology and geochemical analysis in sediment cores, can also provide most indicators to analyze together.
The beginning of the period with the highest fire activity in Nahuelbuta (around 1200–1450 AD; Figure 9) is also consistent with a period of higher temperatures [99] and lower-than-average precipitation around 1200 CE at 41° S [127]. Thereafter, multi-decadal dry spells were recurrent during the first half of the last millennium (11th, 13th, 16th centuries) in south-central Chile [5], which also supports the case for extended fire activity around 1100–1450 CE. For the Araucaria distribution, warm atmospheric temperatures superimposed with dry years and/or summers are also a keystone in promoting widespread fires mainly driven by the desiccation of dead biomass [18,19,20].
Several authors have shown a correlation between fire years in Patagonia and south-central Chile with droughts derived by large-scale atmospheric anomalies associated with SAM and/or ENSO forcings [17,18,19,20,128], including at multi-centennial scale [126,129,130,131,132]. Our paleoenvironmental reconstructions show that the major biomass burning period around 1200–1450 AD coincided very nearly with El Niño-like and positive SAM-like conditions [133,134]. Today, the co-occurrence of the positive SAM phase and El Niño events during summers has led to extremely warm and dry conditions across Patagonia and promoted widespread fires on both sides of the Andes [17,22,28,128]. This study, however, is the first time this pattern has been proved in the mid-latitudes coastal areas. On the other hand, the cultural history of the study area is characterized, starting around 1300 BP, by the presence of groups that have a diversified adaptive transition strategy between hunter-gatherers and the first pottery manifestations (Pitren), with high mobility, which allowed them in their seasonal itineraries to access resources from the sea, lakes and mountain lagoons, as well as collect the fruit of Araucaria trees [135]. Further investigation into the anthropic signal associated with the incendiary activities for that period may be possible with information from future archaeological studies in the highlands of the Nahuelbuta Range.
The documentary use of archives, complemented by the analysis of archaeological sites, allows us to understand how the cultural practices associated with the use of fire may have left their mark on the study area in the last 1000 years. Fire management through time in these territories is directly related to use intensification by human communities, both in prehistoric and historical times, and the changes in the economic patterns involved (gathering/horticulture-slash-and-graze agriculture/Euro-Mediterranean livestock farming/strategies of warfare/the installation of new European economies), resulting, thereby, in a mosaic-type landscape in the territory [136,137], which is congruent with cleared areas observed since the 16th century by Spaniards [138].
On the other hand, the maximum of charcoal observed around the year 1600 AD is very interesting (Figure 9A), where the regional climatic trend is associated with higher amounts of precipitation and lower temperatures [99,139], and just a few dendrochemical records match with the driest events in the PDSI reconstruction in the years 1638–1639, but most of their pulses occur in the wet and cold previous period (Figure 9B,C). Because of that, we suppose it could be related to increased human activity in the study area. Historical records from central and southern Chile document that, at the time of the arrival of the Spanish, the Nahuelbuta Mountain Range would have been used as an area of refuge for indigenous communities that were located in lowland areas, such as Lanalhue or the Purén-Lumaco Valley. Since the 17th century, the Mapuche have become an itinerant and mainly livestock-focused society, adopting Euro-Mediterranean animals such as cows, horses and sheep [108,140]. At the same time, the war between the Spanish and the Mapuche accelerated the abandonment of the valleys until the end of the 19th century [141,142,143,144].
Most of the observed patterns in this multiproxy study of the Nahuelbuta Range confirm the mixed contribution of both climate and humans to the fire regime through time. Future studies should add more paleoenvironmental records to better understand the fire regime changes in this relictual Araucaria forest.

4. Conclusions

Our data indicate that dendrochemical records can contribute as a proxy of fire activity, complementing fire history studies developed by other techniques.
The main result of this study is the increase in chemical concentration of some elements such as Ca, Mg, Ba among others after fires, which were calibrated with fire scar tree ring data and compared with charcoal records from sediment cores. At the same time, our results also show that dendrochemical approaches require calibration with other proxies and a careful interpretation of chemical patterns. Additional studies of dendrochemical response to fires are also needed to improve understanding of the potential capacities and limitations of these techniques to provide longer records of fire history studies in these forest ecosystems. More information about the chemical mobilization in the trunk from different elements and the chemical characterization of forest composition is needed to best understand sources and concentration changes in this forest ecosystem.
Through the multiproxy approach used in this study, for the first time we presented a millennial high-resolution fire history. These results show that an apparent higher fire activity occurs in the period 1200–1650, especially after the years around 1200, 1380 and 1600. On the other hand, a progressively increasing frequency of low-severity fire was also observed after 1700, concomitant with other Araucaria fire histories developed for the Andean range in Chile and Argentina. This information creates a new perspective about fire temporal patterns before and after European colonization in the coastal mountains of south-central Chile.
Finally, climate variability and human land use influence forest fire dynamics in Araucaria forests. More information about past human activities, climate and fire relationships before and after the American colonization are needed to understand current landscape changes and their consequences for forest ecosystems.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/f14061082/s1, Figure S1A,B: Raw element concentrations derived from dendrochemical analyses in the NHL site; Figure S2: Result of Principal Component Analysis (PCA) to be compared with carbon records, temperature and PDSI reconstructions. The elements that contribute the most to PC1 (fuchsia) and those that contribute the least to PC1 (blue) are observed; Figure S3: Age–depth model performed by Bayesian analyses. Six radiocarbon dates (bulk, plant tissue and charcoal) were analyzed and calibrated with the SHcal20 curve; Figure S4: Ring width index chronology of three sites in Nahuelbuta coastal range: NHL (fuchsia), PAG (green), NAH (blue) (The orange vertical line indicates the year prior to release in NHL (1945)); Table S1: Table with dates used in age–depth model. One date was not included in the age–depth model due to its chronostratigraphic error, marked with an asterisk (*).

Author Contributions

Conceptualization, A.A.M., M.E.G. and A.M.A.; methodology, P.R.S., K.K.-B., I.S.-V., M.M.-B., M.M., C.G.R., V.T.-M., B.P., A.M.-C., B.V., A.A.M., M.E.G., A.M.A., M.E.S. and T.R.; formal analysis, I.S.-V., M.M.-B., K.K.-B., C.G.R., B.P., A.M.-C. and B.P.; investigation, A.A.M., M.E.G., A.M.A. and M.E.S.; resources, A.A.M., M.E.G., A.M.A. and M.E.S.; writing—original draft preparation, I.S.-V., A.A.M., M.M.-B., V.T.-M., M.E.G. and K.K.-B.; writing—review and editing, A.A.M., M.E.G., I.S.-V., A.M.A., M.M.-B., K.K.-B., M.E.S., P.R.S. and F.G.; supervision, A.A.M., M.E.G., I.S.-V. and K.K.-B.; project administration, M.E.G.; funding acquisition, M.E.G., A.A.M., A.M.A. and M.E.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Fondecyt/ANID grant N°1201528 and FONDAP-ANID 1522A0001.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

We thank Mary Kay Amistadi from the Arizona Laboratory for Emerging Contaminants. We thank the Corporación Nacional Forestal (CONAF) for facilitating access to the study sites and research permits. The authors also thank Triana Anderson for the English revision and improving of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Swetnam, T.W.; Farella, J.; Roos, C.I.; Liebmann, M.J.; Falk, D.A.; Allen, C.D. Multiscale perspectives of fire, climate and humans in western North America and the Jemez Mountains, USA. Philos. Trans. R. Soc. B Biol. Sci. 2016, 371, 20150168. [Google Scholar] [CrossRef] [PubMed]
  2. Okin, G.S.; Dong, C.; Willis, K.S.; Gillespie, T.W.; MacDonald, G.M. The Impact of Drought on Native Southern California Vegetation: Remote Sensing Analysis Using MODIS-Derived Time Series. J. Geophys. Res. Biogeosci. 2018, 123, 1927–1939. [Google Scholar] [CrossRef]
  3. Jones, M.W.; Smith, A.; Betts, R.; Canadell, J.G.; Prentice, I.C.; Le Quéré, C. Climate change increases the risk of wildfires. Sci. Brief Rev. 2020, 116, 117. Available online: https://www.preventionweb.net/files/73797_wildfiresbriefingnote.pdf (accessed on 13 December 2022).
  4. Xu, R.; Yu, P.; Abramson, M.J.; Johnston, F.H.; Samet, J.M.; Bell, M.L.; Haines, A.; Ebi, K.L.; Li, S.; Guo, Y. Wildfires, Global Climate Change, and Human Health. N. Engl. J. Med. 2020, 383, 2173–2181. [Google Scholar] [CrossRef] [PubMed]
  5. Garreaud, R.D.; Alvarez-Garreton, C.; Barichivich, J.; Boisier, J.P.; Christie, D.; Galleguillos, M.; LeQuesne, C.; McPhee, J.; Zambrano-Bigiarini, M. The 2010–2015 megadrought in central Chile: Impacts on regional hydroclimate and vegetation. Hydrol. Earth Syst. Sci. 2017, 21, 6307–6327. [Google Scholar] [CrossRef]
  6. Garreaud, R.D.; Boisier, J.P.; Rondanelli, R.; Montecinos, A.; Sepúlveda, H.H.; Veloso-Aguila, D. The Central Chile Mega Drought (2010–2018): A climate dynamics perspective. Int. J. Clim. 2019, 40, 421–439. [Google Scholar] [CrossRef]
  7. Demortier, A.; Bozkurt, D.; Jacques-Coper, M. Identifying key driving mechanisms of heat waves in central Chile. Clim. Dyn. 2021, 57, 2415–2432. [Google Scholar] [CrossRef]
  8. Jacques-Coper, M.; Veloso-Aguila, D.; Segura, C.; Valencia, A. Intraseasonal teleconnections leading to heat waves in central Chile. Int. J. Clim. 2021, 41, 4712–4731. [Google Scholar] [CrossRef]
  9. Liu, Y.; Stanturf, J.; Goodrick, S. Trends in global wildfire potential in a changing climate. For. Ecol. Manag. 2010, 259, 685–697. [Google Scholar] [CrossRef]
  10. González, M.E.; Lara, A.; Urrutia, R.; Bosnich, J. Cambio climático y su impacto potencial en la ocurrencia de incendios forestales en la zona centro-sur de Chile (33°–42° S). Bosque 2011, 32, 215–219. [Google Scholar] [CrossRef]
  11. González, M.E.; Gómez-González, S.; Lara, A.; Garreaud, R.; Díaz-Hormazábal, I. The 2010–2015 Megadrought and its influence on the fire regime in central and south-central Chile. Ecosphere 2018, 9, e02300. [Google Scholar] [CrossRef]
  12. Yoon, J.-H.; Kravitz, B.; Rasch, P.J.; Simon Wang, S.-Y.; Gillies, R.R.; Hipps, L. Extreme Fire Season in California: A Glimpse Into the Future? [in “Explaining Extremes of 2014 from a Climate Perspective”]. Bull. Am. Meteorol. Soc. 2015, 96, S5–S9. [Google Scholar] [CrossRef]
  13. Úbeda, X.; Sarricolea, P. Wildfires in Chile: A review. Glob. Planet. Chang. 2016, 146, 152–161. [Google Scholar] [CrossRef]
  14. Turco, M.; von Hardenberg, J.; AghaKouchak, A.; Llasat, M.C.; Provenzale, A.; Trigo, R.M. On the key role of droughts in the dynamics of summer fires in Mediterranean Europe. Sci. Rep. 2017, 7, 81. [Google Scholar] [CrossRef]
  15. Turco, M.; Rosa-Cánovas, J.J.; Bedia, J.; Jerez, S.; Montávez, J.P.; Llasat, M.C.; Provenzale, A. Exacerbated fires in Mediterranean Europe due to anthropogenic warming projected with non-stationary climate-fire models. Nat. Commun. 2018, 9, 3821. [Google Scholar] [CrossRef] [PubMed]
  16. Sun, Q.; Miao, C.; Hanel, M.; Borthwick, A.G.L.; Duan, Q.; Ji, D.; Li, H. Global heat stress on health, wildfires, and agricultural crops under different levels of climate warming. Environ. Int. 2019, 128, 125–136. [Google Scholar] [CrossRef] [PubMed]
  17. Holz, A.; Paritsis, J.; Mundo, I.A.; Veblen, T.T.; Kitzberger, T.; Williamson, G.J.; Aráoz, E.; Bustos-Schindler, C.; González, M.E.; Grau, H.R.; et al. Southern Annular Mode drives multicentury wildfire activity in southern South America. Proc. Natl. Acad. Sci. USA 2017, 114, 9552–9557. [Google Scholar] [CrossRef]
  18. González, M.E.; Muñoz, A.A.; González-Reyes, Á.; Christie, D.A.; Sibold, J. Fire history in Andean Araucaria–Nothofagus forests: Coupled influences of past human land-use and climate on fire regimes in north-west Patagonia. Int. J. Wildland Fire 2020, 29, 649–660. [Google Scholar] [CrossRef]
  19. González, M.E.; Veblen, T.T. Climatic influences on fire in Araucaria araucana–Nothofagus forests in the Andean cordillera of south-central Chile. Écoscience 2006, 13, 342–350. [Google Scholar] [CrossRef]
  20. Urrutia-Jalabert, R.; González, M.E.; González-Reyes, Á.; Lara, A.; Garreaud, R. Climate variability and forest fires in central and south-central Chile. Ecosphere 2018, 9, e02171. [Google Scholar] [CrossRef]
  21. González, M.E.; Veblen, T.T.; Sibold, J.S. Fire history of Araucaria-Nothofagus forests in Villarrica National Park, Chile. J. Biogeogr. 2005, 32, 1187–1202. [Google Scholar] [CrossRef]
  22. Mundo, I.A.; Kitzberger, T.; Roig Juñent, F.A.; Villalba, R.; Barrera, M.D. Fire history in the Araucaria araucana forests of Argentina: Human and climate influences. Int. J. Wildland Fire 2013, 22, 194. [Google Scholar] [CrossRef]
  23. Dickson, B.; Fletcher, M.; Hall, T.L.; Moreno, P.I. Centennial and millennial-scale dynamics in Araucaria–Nothofagus forests in the southern Andes. J. Biogeogr. 2020, 48, 537–547. [Google Scholar] [CrossRef]
  24. Veblen, T.T.; Burns, B.R.; Kitzberger, T.; Lara, A.; Villalba, R. The ecology of the conifers of southern South America. In Ecology of the Southern Conifers; Enright, N., Hill, R., Eds.; Melbourne University Press: Melbourne, Australia, 1995; pp. 120–155. [Google Scholar]
  25. González, M.E.; Cortés, M.; Gallo, L.; Bekessy, S.; Echeverría, C.; Izquierdo, F.; Montaldo, P. Coníferas chilenas: Araucaria araucana. In Las Especies Arbóreas de Los Bosques Templados de Chile y Argentina: Autoecología; Donoso, C., Ed.; Marisa Cuneo Ediciones: Valdivia, Chile, 2013; pp. 36–53. [Google Scholar]
  26. Arroyo-Vargas, P.; Fuentes-Ramírez, A.; Muys, B.; Pauchard, A. Impacts of fire severity and cattle grazing on early plant dynamics in old-growth Araucaria-Nothofagus forests. For. Ecosyst. 2019, 6, 44. [Google Scholar] [CrossRef]
  27. Fuentes-Ramirez, A.; Salas-Eljatib, C.; González, M.E.; Urrutia-Estrada, J.; Arroyo-Vargas, P.; Santibañez, P. Initial response of understorey vegetation and tree regeneration to a mixed-severity fire in old-growth Araucaria–Nothofagus forests. Appl. Veg. Sci. 2020, 23, 210–222. [Google Scholar] [CrossRef]
  28. González, M.E.; Szejner, M.; Muñoz, A.A.; Silva, J. Incendios catastróficos en bosques andinos de Araucaria-Nothofagus: Efecto de la severidad y respuesta de la vegetación. Bosque Nativ. 2010, 46, 12–17. Available online: https://www.yumpu.com/es/document/view/14189903/incendios-catastroficos-en-bosques-andinos-de-araucaria- (accessed on 27 January 2023).
  29. González, M.E.; Veblen, T.T.; Sibold, J.S. Influence of fire severity on stand development of Araucaria araucana-Nothofagus pumilio stands in the Andean cordillera of south-central Chile. Austral Ecol. 2010, 35, 597–615. [Google Scholar] [CrossRef]
  30. González, M.E.; Lara, A. Large fires in the Andean Araucaria forests: When a natural ecological process becomes a threat. Oryx 2015, 49, 394. [Google Scholar] [CrossRef]
  31. Urrutia-Estrada, J.; Fuentes-Ramírez, A.; Hauenstein, E. Diferencias en la composición florística en bosques de Araucaria-Nothofagus afectados por distintas severidades de fuego. Gayana Botánica 2018, 75, 625–638. [Google Scholar] [CrossRef]
  32. Aguilera-Betti, I.; Muñoz, A.A.; Stahle, D.; Figueroa, G.; Duarte, F.; González-Reyes, Á.; Christie, D.; Lara, A.; González, M.E.; Sheppard, P.R.; et al. The First Millennium-Age Araucaria araucana in Patagonia. Tree-Ring Res. 2017, 73, 53–56. [Google Scholar] [CrossRef]
  33. Premoli, A.; Quiroga, P.; Gardner, M. Araucaria araucana. IUCN Red List. Threat. Species 2013, e.T31355A2805113. [Google Scholar] [CrossRef]
  34. Molina, J.R.; Martín, Á.; Drake, F.; Martín, L.M.; Herrera, M.Á. Fragmentation of Araucaria araucana forests in Chile: Quantification and correlation with structural variables. iForest-Biogeosci. For. 2015, 9, 244–252. [Google Scholar] [CrossRef]
  35. MMA. Decreto 79 Aprueba y Oficializa Clasificación de Especies Según Estado de Conservación, Decimocuarto Proceso [Approves and Makes Official Classification of Species According to Conservation Status, Fourteenth Process]. 2018. Available online: https://bcn.cl/2gaep (accessed on 27 January 2023).
  36. Lara, A.; Solari, M.E.; Rutherford, P.; Thiers, O.; Trecamán, R. Cobertura de la Vegetación Original de la Ecoregión de los Bosques Valdivianos de Chile Hacia 1550. [Original Vegetation Cover of the Valdivian Forest Ecoregion of Chile around 1550]; Technical Report Project FB 49; World Wildlife Found/Universidad Austral de Chile: Valdivia, Chile, 1999. [Google Scholar]
  37. González, M.; Cortés, M.; Izquierdo, F.; Gallo, L.; Echeverría, C.; Bekkesy, S.; Montaldo, P. Araucaria araucana (Molina) K. Koch Araucaria (ó), Pehuen, Piñonero, Pino Araucaria, Pino chileno, Pino del Neuquén, Monky puzzle tree. In Las Especies Arbóreas de los Bosques Templados de Chile y Argentina: Autoecología; Donoso, C., Ed.; María Cuneo Ediciones: Valdivia, Chile, 2006; pp. 36–53. [Google Scholar]
  38. Zamorano-Elgueta, C.; Cayuela, L.; Gonzalez-Espinosa, M.; Lara, A.; Parra-Vazquez, M.R. Impacts of cattle on the South American temperate forests: Challenges for the conservation of the endangered monkey puzzle tree (Araucaria araucana) in Chile. Biol. Conserv. 2012, 152, 110–118. [Google Scholar] [CrossRef]
  39. Urrutia, J.; Pauchard, A.; García, R.A. Diferencias en la composición vegetal de un bosque de Araucaria araucana (Molina) K.Koch y Nothofagus antarctica (G. Forst.) Oerst: Asociadas a un gradiente de invasión de Pinus contorta Douglas ex Loudon. [Differences in the vegetal composition of a forest of Araucaria araucana (Molina) K.Koch and Nothofagus antarctica (G. Forst.) Oerst. associated with an invasion gradient of Pinus contorta Douglas ex Loudon]. Gayana Botánica 2013, 70, 92–100. [Google Scholar] [CrossRef]
  40. Cóbar-Carranza, A.J.; García, R.A.; Pauchard, A.; Peña, E. Effect of Pinus contorta invasion on forest fuel properties and its potential implications on the fire regime of Araucaria araucana and Nothofagus antarctica forests. Biol. Invasions 2014, 16, 2273–2291. [Google Scholar] [CrossRef]
  41. Assal, T.J.; González, M.E.; Sibold, J.S. Burn severity controls on postfire Araucaria-Nothofagus Regeneration in the Andean Cordillera. J. Biogeogr. 2018, 45, 2483–2494. [Google Scholar] [CrossRef]
  42. Agee, J.K. Fire Ecology of Pacific Northwest Forests; Island Press: Washington, DC, USA, 1993; 493p. [Google Scholar]
  43. Whitlock, C.; Anderson, R.S. Fire history reconstructions based on sediment records from lakes and wetlands. In Fire and Climatic Change in Temperate Ecosystems of the Western Americas; Veblen, T.T., Baker, W.L., Montenegro, G., Swetnam, T.W., Eds.; Springer: New York, NY, USA, 2003; pp. 3–31. [Google Scholar]
  44. Conedera, M.; Tinner, W.; Neff, C.; Meurer, M.; Dickens, A.F.; Krebs, P. Reconstructing past fire regimes: Methods, applications, and relevance to fire management and conservation. Quat. Sci. Rev. 2009, 28, 555–576. [Google Scholar] [CrossRef]
  45. Baker, W.L. Fire Ecology in Rocky Mountain Landscapes; Island Press: Washington, DC, USA, 2009; 605p. [Google Scholar]
  46. Kitzberger, T.; Veblen, T.T.; Villalba, R. Métodos dendrocronológicos y sus aplicaciones en estudios de dinámica de bosques templados de Sudamérica [Dendrochronological methods and their applications in studies of dynamics of temperate forests in South America]. In Dendrocronología en América Latina; Roig, F., Ed.; EDIUNC: Mendoza, Argentina, 2000; pp. 17–78. [Google Scholar]
  47. Smith, K.T.; Sutherland, E.K. Terminology and biology of fire scars in selected central hardwoods. Tree-Ring Res. 2001, 57, 141–147. Available online: https://www.fs.usda.gov/research/treesearch/18749 (accessed on 27 January 2023).
  48. Guyette, R.P.; Dey, D.C.; Stambaugh, M.C.; Muzika, R. Fire scars reveal variability and dynamics of eastern fire regimes. In Proceedings of the Fire in Eastern Oak Forests: Delivering Science to Land Managers, Proceedings of a Conference; Columbus, OH, USA, 15–17 November 2005, Dickinson, M.B., Ed.; General Technical Report NRS-P-1; U.S. Department of Agriculture, Forest Service, Northern Research Station: Newtown Square, PA, USA, 2006; pp. 20–39. Available online: https://www.nrs.fs.usda.gov/pubs/gtr/gtr_nrs-p1/guyette_p1_20.pdf (accessed on 27 January 2023).
  49. MacDonald, G.M.; Larsen, C.P.S.; Szeicz, J.M.; Moser, K.A. The reconstruction of boreal forest fire history from lake sediments: A comparison of charcoal, pollen, sedimentological, and geochemical indices. Quat. Sci. Rev. 1991, 10, 53–71. [Google Scholar] [CrossRef]
  50. Abarzúa, A.M.; Pinchicura, A.G.; Jarpa, L.; Martel-Cea, A.; Sterken, M.; Vega, R.; Pino, Q.M. Environmental Responses to Climatic and Cultural Changes. In The Teleoscopic Polity; Dillehay, T., Ed.; Springer International Publishing: Berlin/Heidelberg, Germany, 2014; pp. 123–141. [Google Scholar] [CrossRef]
  51. Whitlock, C.; Larsen, C. Charcoal as a fire proxy. In Tracking Environmental Change Using Lake Sediments: Terrestrial, Algal, and Siliceous Indicators; Smol, J.P., Birks, H.J.B., Last, W.M., Eds.; Springer: Dordrecht, The Netherlands, 2001; pp. 75–97. [Google Scholar] [CrossRef]
  52. Nabais, C.; Freitas, H.; Hagemeyer, J. Tree rings and dendroanalysis. In Metals in the Environment. Analysis by Biodiversity, New York. Prasad, M.N.V. (2001). Metals in the Environment; Prasad, M., Ed.; CRC Press: Boca Raton, FL, USA, 2001; pp. 367–400. [Google Scholar] [CrossRef]
  53. Watmough, S.A. An evaluation of the use of dendrochemical analyses in environmental monitoring. Environ. Rev. 1997, 5, 181–201. [Google Scholar] [CrossRef]
  54. Binda, G.; Di Iorio, A.; Monticelli, D. The what, how, why, and when of dendrochemistry: (Paleo) Environmental information from the chemical analysis of tree rings. Sci. Total Environ. 2021, 758, 143672. [Google Scholar] [CrossRef]
  55. Ballikaya, P.; Marshall, J. Cherubini, PCan tree-ring chemistry be used to monitor atmospheric nanoparticle contamination over time? Atmos. Environ. 2022, 268, 118781. [Google Scholar] [CrossRef]
  56. Muñoz, A.A.; Klock-Barría, K.; Sheppard, P.R.; Aguilera-Betti, I.; Toledo-Guerrero, I.; Christie, D.A.; Gorena, T.; Gallardo, L.; González-Reyes, Á.; Lara, A.; et al. Multidecadal environmental pollution in a mega-industrial area in central Chile registered by tree rings. Sci. Total Environ. 2019, 696, 133915. [Google Scholar] [CrossRef] [PubMed]
  57. Cocozza, C.; Alterio, E.; Bachmann, O.; Guillong, M.; Sitzia, T.; Cherubini, P. Monitoring air pollution close to a cement plant and in a multi-source industrial area through tree-ring analysis. Environ. Sci. Pollut. Res. 2021, 28, 54030–54040. [Google Scholar] [CrossRef]
  58. Ballikaya, P.; Song, W.; Bachmann, O.; Guillong, M.; Wang, X.; Cherubini, P. Chemical Elements Recorded by Quercus mongolica Fisch. ex Ledeb. Tree Rings Reveal Trends of Pollution History in Harbin, China. Forests 2023, 14, 187. [Google Scholar] [CrossRef]
  59. Pearson, C.; Manning, S.W.; Coleman, M.; Jarvis, K. Can tree-ring chemistry reveal absolute dates for past volcanic eruptions? J. Archaeol. Sci. 2005, 32, 1265–1274. [Google Scholar] [CrossRef]
  60. Pearson, C.L.; Manning, S.W.; Coleman, M.; Jarvis, K. A dendrochemical study of Pinus sylvestris from Siljansfors Experimental Forest, central Sweden. Appl. Geochem. 2006, 21, 1681–1691. [Google Scholar] [CrossRef]
  61. Alfaro-Sánchez, R.; Camarero, J.J.; Querejeta, J.I.; Sagra, J.; Moya, D.; Rodríguez-Trejo, D.A. Volcanic activity signals in tree-rings at the treeline of the Popocatépetl, Mexico. Dendrochronologia 2020, 59, 125663. [Google Scholar] [CrossRef]
  62. Bertin, L.J.; Christie, D.A.; Sheppard, P.R.; Muñoz, A.A.; Lara, A.; Alvarez, C. Chemical Signals in Tree Rings from Northern Patagonia as Indicators of Calbuco Volcano Eruptions since the 16th Century. Forests 2021, 12, 1305. [Google Scholar] [CrossRef]
  63. Odigie, K.O.; Khanis, E.; Hibdon, S.A.; Jana, P.; Araneda, A.; Urrutia, R.; Flegal, A.R. Remobilization of trace elements by forest fire in Patagonia, Chile. Reg. Environ. Chang. 2015, 16, 1089–1096. [Google Scholar] [CrossRef]
  64. Fernandez-Marcos, M.L. Potentially Toxic Substances and Associated Risks in Soils Affected by Wildfires: A Review. Toxics 2022, 10, 31. [Google Scholar] [CrossRef] [PubMed]
  65. Balouet, J.-C.; Oudijk, G.; Smith, K.T.; Petrisor, I.; Grudd, H.; Stocklassa, B. Applied Dendroecology and Environmental Forensics. Characterizing and Age Dating Environmental Releases: Fundamentals and Case Studies. Environ. Forensics 2007, 8, 1–17. [Google Scholar] [CrossRef]
  66. Smith, K.T.; Arbellay, E.; Falk, D.A.; Sutherland, E.K. Macroanatomy and compartmentalization of recent fire scars in three North American conifers. Can. J. For. Res. 2016, 46, 535–542. [Google Scholar] [CrossRef]
  67. Abraham, J.; Dowling, K.; Florentine, S. The Unquantified Risk of Post-Fire Metal Concentration in Soil: A Review. Water Air Soil Pollut. 2017, 228, 1–33. [Google Scholar] [CrossRef]
  68. McBride, J.R. Analysis of tree rings and fire scars to establish fire history. Tree-Ring Bull. 1983, 43, 51–67. Available online: https://repository.arizona.edu/handle/10150/261216 (accessed on 27 January 2023).
  69. Wei, M.; Guo, X.; Ma, Y.; Tigabu, M.; Zheng, W.; Liu, M.; Guo, F. Impact of forest fire on radial growth of tree rings and their element concentrations of Pinus sylvestris and Larix gmelinii in northern China. Front. For. Glob. Chang. 2023, 6, 1136039. [Google Scholar] [CrossRef]
  70. Aguilera-Betti, I.; Lucas, C.; Ferrero, M.E.; Muñoz, A.A. A Network for Advancing Dendrochronology, Dendrochemistry and Dendrohydrology in South America. Tree-Ring Res. 2020, 76, 94. [Google Scholar] [CrossRef]
  71. Figueroa, G.; Muñoz, A.A.; Aguilera-Betti, I.; Toledo-Guerrero, I.; Puchi, P.; Sheppard, P. Effects of eruptions Villarrica Volcano observed in the chemical variability of Araucaria araucana tree-rings records. In Proceedings of the Cities on Volcanoes 9. Understanding Volcanoes a Society: The Key for Risk Mitigation, Puerto Varas, Chile, 20–25 November 2016. [Google Scholar]
  72. Cutter, B.E.; Guyette, R.P. Anatomical, Chemical, and Ecological Factors Affecting Tree Species Choice in Dendrochemistry Studies. J. Environ. Qual. 1993, 22, 611–619. [Google Scholar] [CrossRef]
  73. Smith-Ramírez, C. The Chilean coastal range: A vanishing center of biodiversity and endemism in South American temperate rainforests. Biodivers. Conserv. 2004, 13, 373–393. [Google Scholar] [CrossRef]
  74. Di Castri, F.; Hajek, E. Bioclimatología de Chile [Bioclimatology of Chile]; Vicerrectoría Académica de la Universidad Católica de Chile: Santiago, Chile, 1976; 129p. [Google Scholar]
  75. Donoso, C.; Gerding, V.; Olivares, B.; Real, P.; Sandoval, V.; Schlatter, R.; Schlegel, F. Antecedentes Para el Manejo del Bosque Nativo de Forestal Arauco [Background for the Management of the Native Forest of Forestal Arauco]; Informe de Convenio N° 74; Universidad Austral de Chile: Valdivia, Chile, 1984; 183p. [Google Scholar]
  76. Dillehay, T. La Organización Política Temprana de los Mapuche. Materialidad y Patriarcado Andino [The Early Political Organization of the Mapuche. Materiality and Andean Patriarchy]; Pehuén: Santiago, Chile, 2017. [Google Scholar]
  77. Aldunate, C.; Villagrán, C. Recolectores de los bosques templados del cono-sur americano [Collectors of the temperate forests of the American southern cone]. In Botánica Indígena de Chile; de Moesbach, W., Ed.; Andrés Bello: Santiago, Chile, 1992; pp. 23–38. [Google Scholar]
  78. Zavala, J.M. La cordillera de Nahuelbuta, ¿lugar de naturaleza o de cultura?: Concepción del espacio cordillerano en fuentes españolas del siglo XVI [The Nahuelbuta mountain range, a place of nature or culture?: Conception of the mountain range space in Spanish sources of the 16th century]. Rev. Cuhso 2008, 15, 17–24. [Google Scholar]
  79. de Góngora Marmolejo, A. Historia de Chile Desde su Descubrimiento Hasta el año de 1575. Colección de Historiadores de Chile y de Documentos Relativos a la Historia Nacional [History of Chile from Its Discovery to the Year 1575. Collection of Chilean Historians and Documents Related to National History]; Impr. del Ferrocarril: Santiago, Chile, 1862. [Google Scholar]
  80. González de Nájera, A. Desengaño y Reparo de la Guerra del Reino de Chile Donde se Manifiestan las Principales Ventajas que en Ella tienen los Indios a Nuestros Españoles [Disappointment and Repair of the War in the Kingdom of Chile Where the Main Advantages That the Indians Have over Our Spaniards in It Are Manifested]; Impr. Ercilla: Santiago, Chile, 1889. [Google Scholar]
  81. Fitz-Roy, R. Viajes del “Adventure” y el “Beagle”: Diario (Biblioteca Darwiniana); Catarata: Madrid, Spain, 2013. [Google Scholar]
  82. Lara, H. Crónica de la Araucanía: Descubrimiento y Conquista, Pacificación Definitiva i Campaña de Villa-Rica (Leyenda Heroica de Tres Siglos) [Chronicle of the Araucanía: Discovery and Conquest, Final Pacification and Campaign of Villa-Rica (Heroic Legend of Three Centuries)]; Impr. El Progreso: Santiago, Chile, 1888. [Google Scholar]
  83. Mansoulet, J. Guía-Crónica de la Frontera Araucana de Chile: Años 1892–1893. Apuntes Históricos, Topográficos, Geográficos, Descriptivos y Estadísticos de sus Poblaciones, de su Comercio, Industrias y Agricultura [Guide-Chronicle of the Araucanian Border of Chile: Years 1892–1893. Historical, Topographical, Geographical, Descriptive and Statistical Notes of Its Populations, Its Commerce, Industries and Agricultura]; Impr. Barcelona: Santiago, Chile, 1893. [Google Scholar]
  84. Domeyko, I. Araucanía y sus Habitantes: Recuerdos de un Viaje Hecho en las Provincias Meridionales de Chile en los Meses de Enero y Febrero 1845 [Araucanía and Its Inhabitants: Memories of a Trip Made in the Southern Provinces of Chile in the Months of January and February 1845]; Editorial Francisco de Aguirre: Buenos Aires, Argentina, 1971. [Google Scholar]
  85. Arno, S.F.; Sneck, K.M. A Method for Determining Fire history in Coniferous Forests of the Mountain West; General Technical Report INT-42; USDA Forest Service, Intermountain Forest Experiment Station: Tucson, AZ, USA, 1977. [Google Scholar]
  86. Barrett, S.W.; Arno, S.F. Increment-Borer Methods For Determining Fire History in Coniferous Forests; General Technical Report INT-244; US Department of Agriculture, Forest Service, Intermountain Research Station: Ogden, UT, USA, 1988; Volume 244, 15p.
  87. Stokes, M.A.; Smiley, T.L. An Introduction to Tree-Ring Dating; University of Arizona Press: Tucson, AZ, USA, 1996; 73p. [Google Scholar]
  88. Speer, J.H. Fundamentals of Tree Ring Research; The University of Arizona Press: Tucson, AZ, USA, 2010. [Google Scholar]
  89. Schulman, E. Dendroclimatic Change in Semiarid America; University of Arizona Press: Tucson, AZ, USA, 1956; 142p. [Google Scholar]
  90. Holmes, R.L. Computer-assisted quality control in tree-ring dating and measurement. Tree-Ring Bull. 1983, 43, 69–78. [Google Scholar]
  91. Cook, E.R.; Holmes, R.H. ARSTAN Program and Users’ Manual. Laboratory of Tree-Ring Research; University of Arizona: Tucson, AZ, USA, 1984; 15p. [Google Scholar]
  92. Abatzoglou, J.T.; Dobrowski, S.Z.; Parks, S.A.; Hegewisch, K.C. TerraClimate, a high-resolution global dataset of monthly climate and climatic water balance from 1958–2015. Sci. Data 2018, 5, 170191. [Google Scholar] [CrossRef]
  93. Mooney, C.Z.; Duval, R.D. Bootstrapping: A Nonparametric Approach to Statistical Inference; Sage Publications, Inc.: Thousand Oaks, CA, USA, 1993. [Google Scholar] [CrossRef]
  94. Grissino-Mayer, H.D. Tree-Ring Reconstructions of Climate and Fire History at El Malpais National Monument, New Mexico. Ph.D. Thesis, University of Arizona, Tucson, AZ, USA, 1995; 407p. [Google Scholar]
  95. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2022; Available online: https://www.R-project.org/ (accessed on 25 February 2023).
  96. RStudio Team. RStudio: Integrated Development Environment for R. RStudio; PBC: Boston, MA, USA, 2022; Available online: http://www.rstudio.com/ (accessed on 25 February 2023).
  97. Bunn, A.G. A dendrochronology program library in R (dplR). Dendrochronologia 2008, 26, 115–124. [Google Scholar] [CrossRef]
  98. Bunn, A.; Korpela, M.; Biondi, F.; Campelo, F.; Mérian, P.; Qeadan, F.; Zang, C. dplR: Dendrochronology Program Library in R. R Package Version 1.7.4. 2022. Available online: https://CRAN.R-project.org/package=dplR (accessed on 25 February 2023).
  99. Lara, A.; Villalba, R.; Urrutia-Jalabert, R.; González-Reyes, A.; Aravena, J.C.; Luckman, B.H.; Cuq, E.; Rodríguez, C.; Wolodarsky-Franke, A. +A 5680-year tree-ring temperature record for southern South America. Quat. Sci. Rev. 2020, 228, 106087. [Google Scholar] [CrossRef]
  100. Morales, M.S.; Cook, E.R.; Barichivich, J.; Christie, D.A.; Villalba, R.; LeQuesne, C.; Srur, A.M.; Ferrero, M.E.; González-Reyes, Á.; Couvreux, F.; et al. Six hundred years of South American tree rings reveal an increase in severe hydroclimatic events since mid-20th century. Proc. Natl. Acad. Sci. USA 2020, 117, 16816–16823. [Google Scholar] [CrossRef] [PubMed]
  101. Aquino-López, M.A.; Blaauw, M.; Christen, J.A.; Sanderson, N.K. Bayesian Analysis of 210Pb Dating. J. Agric. Biol. Environ. Stat. 2018, 23, 317–333. [Google Scholar] [CrossRef]
  102. Blaauw, M.; Christen, J.A.; Aquino-López, M.A. Rplum: Bayesian Age-Depth Modelling of ‘210Pb’-Dated Cores. R Package Version 0.1.4. 2020. Available online: https://CRAN.R-project.org/package=rplum (accessed on 30 January 2023).
  103. Hogg, A.G.; Heaton, T.J.; Hua, Q.; Palmer, J.G.; Turney, C.S.; Southon, J.; Bayliss, A.; Blackwell, P.G.; Boswijk, G.; Bronk Ramsey, C.; et al. SHCal20 Southern Hemisphere Calibration, 0–55,000 Years cal BP. Radiocarbon 2020, 62, 759–778. [Google Scholar] [CrossRef]
  104. Higuera, P.E.; Whitlock, C.; Gage, J.A. Linking tree-ring and sediment-charcoal records to reconstruct fire occurrence and area burned in subalpine forests of Yellowstone National Park, USA. Holocene 2011, 21, 327–341. [Google Scholar] [CrossRef]
  105. Higuera, P.E.; Brubaker, L.B.; Anderson, P.M.; Hu, F.S.; Brown, T.A. Vegetation mediated the impacts of postglacial climate change on fire regimes in the south-central Brooks Range, Alaska. Ecol. Monogr. 2009, 79, 201–219. [Google Scholar] [CrossRef]
  106. Mundo, I.A.; Wiegand, T.; Kanagaraj, R.; Kitzberger, T. Environmental drivers and spatial dependency in wildfire ignition patterns of northwestern Patagonia. J. Environ. Manag. 2013, 123, 77–87. [Google Scholar] [CrossRef]
  107. Iglesias, W. Reconstrucción de 500 Años de Incendios en Bosques de Araucaria araucana en la Primera Área Protegida por el Estado: Reserva Nacional Malleco, Región de la Araucanía [Reconstruction of 500 Years of Fires in Araucaria araucana Forests in the First Area Protected by the State: Malleco National Reserve, Araucanía Region]. 2011. Available online: https://www.curriculumnacional.cl/estudiante/621/articles-262032_recurso_01.pdf (accessed on 27 January 2023).
  108. Bengoa, J. Historia del Pueblo Mapuche: Siglo XIX y XX [History of the Mapuche People: XIX and XX Centuries]; LOM: Santiago, Chile, 2000; 430p. [Google Scholar]
  109. Saavedra, C. Documentos Relativos a la Ocupación de Arauco [Documents Related to the Occupation of Arauco]; Biblioteca Fundamentos de la Construcción de Chile: Santiago, Chile, 2009; 392p. [Google Scholar]
  110. Villalobos, S. Incorporación de La Araucanía. Relatos Militares [Incorporation of La Araucanía. Military Accounts]; Catalonia: Santiago, Chile, 2018; pp. 1822–1883. 327p. [Google Scholar]
  111. Muñoz, A.A.; Barichivich, J.; Christie, D.A.; Dorigo, W.; Sauchyn, D.; González-Reyes, Á.; Villalba, R.; Lara, A.; Riquelme, N.; González, M.E. Patterns and drivers of Araucaria araucana forest growth along a biophysical gradient in the northern Patagonian Andes: Linking tree rings with satellite observations of soil moisture. Austral Ecol. 2014, 39, 158–169. [Google Scholar] [CrossRef]
  112. Mundo, I.A.; Roig Juñent, F.A.; Villalba, R.; Kitzberger, T.; Barrera, M.D. Araucaria araucana tree-ring chronologies in Argentina: Spatial growth variations and climate influences. Trees 2012, 26, 443–458. [Google Scholar] [CrossRef]
  113. Hagemeyer, J. Trace metals in tree rings: What do they tell us. In Trace Elements—Their Distribution and Effects in the Environment; Elsevier: Amsterdam, The Netherlands, 2000; pp. 375–385. [Google Scholar] [CrossRef]
  114. Montiel, M.; González, M.E.; Muñoz, A.A.; Christie, D.A.; Sheppard, P.R. Dendroecological and Dendrochemical Responses after Two Successive Fires in an Araucaria araucana Forest; Escuela de Graduados, Facultad de Ciencias Forestales y Recursos Naturales, Universidad Austral de Chile: Valdivia, Chile, 2023. [Google Scholar]
  115. Austruy, A.; Yung, L.; Ambrosi, J.P.; Girardclos, O.; Keller, C.; Angeletti, B.; Dron, J.; Chamaret, P.; Chalot, M. Evaluation of historical atmospheric pollution in an industrial area by dendrochemical approaches. Chemosphere 2019, 220, 116–126. [Google Scholar] [CrossRef] [PubMed]
  116. Agbeshie, A.A.; Abugre, S.; Atta-Darkwa, T.; Awuah, R. A review of the effects of forest fire on soil properties. J. For. Res. 2022, 33, 1419–1441. [Google Scholar] [CrossRef]
  117. Campos, I.; Abrantes, N.; Keizer, J.J.; Vale, C.; Pereira, P. Major and trace elements in soils and ashes of eucalypt and pine forest plantations in Portugal following a wildfire. Sci. Total Environ. 2016, 572, 1363–1376. [Google Scholar] [CrossRef] [PubMed]
  118. Houle, D.; Tremblay, S.; Ouimet, R. Foliar and wood chemistry of sugar maple along a gradient of soil acidity and stand health. Plant Soil 2007, 300, 173–183. [Google Scholar] [CrossRef]
  119. Bukata, A.R.; Kyser, T.K.; Al, T.A. Nitrogen Isotopic Composition and Elemental concentration of Tree-rings May Help Map the Extent of Historic Fire Events. Fire Ecol. 2008, 4, 101–107. [Google Scholar] [CrossRef]
  120. Pearson, C.L.; Dale, D.; Lombardo, K. An investigation of fire scars in Pseudotsuga macrocarpa by Scanning X-ray Fluorescence Microscopy. For. Ecol. Manag. 2011, 262, 1258–1264. [Google Scholar] [CrossRef]
  121. Scharnweber, T.; Hevia, A.; Buras, A.; van der Maaten, E.; Wilmking, M. Common trends in elements? Within- and between-tree variations of wood-chemistry measured by X-ray fluorescence—A dendrochemical study. Sci. Total Environ. 2016, 566–567, 1245–1253. [Google Scholar] [CrossRef]
  122. Gavrikov, V.L.; Fertikov, A.I.; Sharafutdinov, R.A.; Vaganov, E.A. Species-specific and Non-species-specific Elemental Trends in Tree Rings. Dokl. Earth Sci. 2021, 496, 96–99. [Google Scholar] [CrossRef]
  123. Sánchez-Salguero, R.; Camarero, J.J.; Hevia, A.; Sangüesa-Barreda, G.; Galván, J.D.; Gutiérrez, E. Testing annual tree-ring chemistry by X-ray fluorescence for dendroclimatic studies in high-elevation forests from the Spanish Pyrenees. Quat. Int. 2019, 514, 130–140. [Google Scholar] [CrossRef]
  124. Gavrikov, V.; Fertikov, A.; Sharafutdinov, R.; Pyzhev, A.; Vaganov, E. Dendrochemical Challenge in Climate Science: Whether Chemical Elements in Wood Reflect the Fluctuations in Weather Parameters. Plants 2022, 11, 3240. [Google Scholar] [CrossRef] [PubMed]
  125. Peckham, M.A.; Gustin, M.S.; Weisberg, P.J. Assessment of the Suitability of Tree Rings as Archives of Global and Regional Atmospheric Mercury Pollution. Environ. Sci. Technol. 2019, 53, 3663–3671. [Google Scholar] [CrossRef] [PubMed]
  126. Holz, A.; Haberle, S.; Veblen, T.T.; De Pol-Holz, R.; Southon, J. Fire history in western Patagonia from paired tree-ring fire-scar and charcoal records. Clim. Past 2012, 8, 451–466. [Google Scholar] [CrossRef]
  127. Lamy, F.; Hebbeln, D.; Röhl, U.; Wefer, G. Holocene rainfall variability in southern Chile: A marine record of latitudinal shifts of the Southern Westerlies. Earth Planet. Sci. Lett. 2001, 185, 369–382. [Google Scholar] [CrossRef]
  128. Garreaud, R.D. A plausible atmospheric trigger for the 2017 coastal El Niño. Int. J. Clim. 2018, 38, e1296–e1302. [Google Scholar] [CrossRef]
  129. Abarzúa, A.M.; Moreno, P.I. Changing fire regimes in the temperate rainforest region of southern Chile over the last 16,000 yr. Quat. Res. 2008, 69, 62–71. [Google Scholar] [CrossRef]
  130. Fletcher, M.S.; Moreno, P.I. Have the Southern Westerlies changed in a zonally symmetric manner over the last 14,000 years? A hemisphere-wide take on a controversial problem. Quat. Int. 2012, 253, 32–46. [Google Scholar] [CrossRef]
  131. Moreno, P.I.; Vilanova, I.; Villa-Martínez, R.P.; Francois, J.P. Modulation of fire regimes by vegetation and site type in southwestern Patagonia since 13 ka. Front. Ecol. Evol. 2018, 6, 34. [Google Scholar] [CrossRef]
  132. Whitlock, C.; Moreno, P.I.; Bartlein, P. Climatic controls of Holocene fire patterns in southern South America. Quat. Res. 2007, 68, 28–36. [Google Scholar] [CrossRef]
  133. Yan, H.; Sun, L.; Wang, Y.; Huang, W.; Qiu, S.; Yang, C. A record of the Southern Oscillation Index for the past 2000 years from precipitation proxies. Nat. Geosci. 2011, 4, 611–614. [Google Scholar] [CrossRef]
  134. Dätwyler, C.; Neukom, R.; Abram, N.J.; Gallant, A.J.E.; Grosjean, M.; Jacques-Coper, M.; Karoly, D.J.; Villalba, R. Teleconnection stationarity, variability and trends of the Southern Annular Mode (SAM) during the last millennium. Clim. Dyn. 2018, 51, 2321–2339. [Google Scholar] [CrossRef]
  135. Adán, L.; Mera, R. Acerca de la distribución espacial y temporal del Complejo Pitrén. Una reevaluación a partir del estudio sistemático de colecciones [About the spatial and temporal distribution of the Pitrén Complex. A reassessment from the systematic study of collections]. Boletín Soc. Chil. Arqueol. 1997, 24, 33–37. [Google Scholar]
  136. Lara, A.; Solari, M.E.; Prieto, M.R.; Peña, M.P. Reconstrucción de la cobertura de la vegetación y uso del suelo hacia 1550 y sus cambios a 2007 en la ecorregión de los bosques valdivianos lluviosos de Chile (35°–43°30′ S) [Reconstruction of vegetation cover and land use around 1550 and its changes to 2007 in the Valdivian rain forests ecoregion of Chile (35°–43°30′ S)]. Bosque 2012, 33, 13–23. [Google Scholar] [CrossRef]
  137. Luebert, F.; Pliscoff, L. Sinopsis Bioclimática y Vegetacional de Chile [Bioclimatic and Vegetational Synopsis of Chile]; Editorial Universitaria: Santiago, Chile, 2017. [Google Scholar]
  138. Camus, P. Bosques y tierras despejadas en el período de la conquista de Chile [Forests and cleared lands in the period of the conquest of Chile]. In Estudios Coloniales II Santiago. Editorial Biblioteca Americana, Universidad Andrés Bello; Retamales Ávila, J., Ed.; RIL Editores: Santiago, Chile, 2002; pp. 159–180. [Google Scholar]
  139. Le Quesne, C.; Acuña, C.; Boninsegna, J.A.; Rivera, A.; Barichivich, J. Long-term glacier variations in the Central Andes of Argentina and Chile, inferred from historical records and tree-ring reconstructed precipitation. Palaeogeogr. Palaeoclimatol. Palaeoecol. 2009, 281, 334–344. [Google Scholar] [CrossRef]
  140. Bengoa, J. Historia de los Antiguos Mapuches del sur [History of the Ancient Mapuches of the South]; Editorial Catalonia: Santiago, Chile, 2003. [Google Scholar]
  141. León, L. Maloqueros y Conchavadores en Araucanía y las Pampas, 1700–1800 [Maloqueros and Conchavadores in Araucanía and the Pampas, 1700–1800]; Serie Quinto Centenario; Ediciones Universidad de La Frontera: Araucanía, Chile, 1991. [Google Scholar]
  142. León, L.; Herrera, P.; Parentini, L.C.; Villalobos, S. Araucanía: La Frontera Mestiza, Siglo XIX [Araucanía: The Mestizo Border, 19th century]; Ed. Lom: Santiago, Chile, 2004. [Google Scholar]
  143. Torrejón, F.; Cisternas, M. Alteraciones del paisaje ecológico araucano por la asimilación mapuche de la agroganadería hispano-mediterránea (siglos XVI y XVII) [Changes in the Araucanian ecological landscape due to the Mapuche assimilation of Spanish-Mediterranean farming (16th and 17th centuries)]. Rev. Chil. Hist. Nat. 2002, 75, 729–736. [Google Scholar] [CrossRef]
  144. Torrejón, F.; Cisternas, M. Impacto ambiental temprano en la Araucanía deducido de crónicas españolas y estudios historiográficos [Early environmental impact in Araucanía deduced from Spanish chronicles and historiographical studies]. Bosque 2003, 24, 45–56. [Google Scholar] [CrossRef]
Figure 1. (A) Study area: the entire extension of the Nahuelbuta Coastal Range is shown, along with the distribution of Araucaria araucana in this area. (B) Sampling locations for the different methodological approximations: dendrochemistry (NHL), fire scars (NHL, COI) and charcoal record (ACA). (C) Scar sample. (D) Segment of processed tree ring core. (E) Sediment core. (F) A. araucana landscape in the Nahuelbuta Coastal Range.
Figure 1. (A) Study area: the entire extension of the Nahuelbuta Coastal Range is shown, along with the distribution of Araucaria araucana in this area. (B) Sampling locations for the different methodological approximations: dendrochemistry (NHL), fire scars (NHL, COI) and charcoal record (ACA). (C) Scar sample. (D) Segment of processed tree ring core. (E) Sediment core. (F) A. araucana landscape in the Nahuelbuta Coastal Range.
Forests 14 01082 g001
Figure 2. Dendrochemical processing: digestion and dilution of wood samples in Arizona Laboratory for Emerging Contaminants (ALEC). (A) Composite sample. (B) Acid predigestion of composite sample. (C) Dilution of 1.5 mL aliquot of digestate with Mili-Q water. (D) Samples in a digestion block to complete digestion. (E) Samples ready to be analyzed by ICP-MS.
Figure 2. Dendrochemical processing: digestion and dilution of wood samples in Arizona Laboratory for Emerging Contaminants (ALEC). (A) Composite sample. (B) Acid predigestion of composite sample. (C) Dilution of 1.5 mL aliquot of digestate with Mili-Q water. (D) Samples in a digestion block to complete digestion. (E) Samples ready to be analyzed by ICP-MS.
Forests 14 01082 g002
Figure 3. Fire-scar chronologies from two sites in Nahuelbuta Coastal Range: COI and NHL. Horizontal lines represent individual trees and fire-scar events are shown by short vertical lines. Fire events and dates (≥1 or ≥2 fire-scarred trees) are indicated below next to the x-axis.
Figure 3. Fire-scar chronologies from two sites in Nahuelbuta Coastal Range: COI and NHL. Horizontal lines represent individual trees and fire-scar events are shown by short vertical lines. Fire events and dates (≥1 or ≥2 fire-scarred trees) are indicated below next to the x-axis.
Forests 14 01082 g003
Figure 4. (A) Ring width index chronology (RWI) of NHL site for the period 1100–2020 (fuchsia line). Vertical dotted line indicates where EPS > 0.85. Right vertical axis shows the number of samples (gray line). Orange triangles indicate years where ≥2 trees presented fire scars. (B) Correlations between accumulated precipitation (blue, 1958–2020) and average temperature (orange, 1950–2020) in periods of three months with the RWI chronology. Horizontal dashed lines indicate critical r values and asterisks indicate significant correlations. (C) Superposed Epoch Analysis between RWI chronology in the period 1350–2020 and fire dates registered in ≥2 trees. The x axis represents a 7-year period, from 1 year before the fire to 5 years after the fire. The dark-green colored box indicates a 95% confidence level for the analysis.
Figure 4. (A) Ring width index chronology (RWI) of NHL site for the period 1100–2020 (fuchsia line). Vertical dotted line indicates where EPS > 0.85. Right vertical axis shows the number of samples (gray line). Orange triangles indicate years where ≥2 trees presented fire scars. (B) Correlations between accumulated precipitation (blue, 1958–2020) and average temperature (orange, 1950–2020) in periods of three months with the RWI chronology. Horizontal dashed lines indicate critical r values and asterisks indicate significant correlations. (C) Superposed Epoch Analysis between RWI chronology in the period 1350–2020 and fire dates registered in ≥2 trees. The x axis represents a 7-year period, from 1 year before the fire to 5 years after the fire. The dark-green colored box indicates a 95% confidence level for the analysis.
Forests 14 01082 g004
Figure 5. Standardized series of element concentrations in tree rings over time. Vertical lines extending from the x-axis indicate fire events recorded by fire scars in 1 tree (gray) and are shown in the upper x-axis, and ≥2 trees (black) are shown in the lower x-axis. Due to the extent of scar-dated events, the chemical record is shown from 1700 to 2020. Asterisks (*) indicate elements above p90 when ≥5 elements are over p90. RWI = Ring Width Index.
Figure 5. Standardized series of element concentrations in tree rings over time. Vertical lines extending from the x-axis indicate fire events recorded by fire scars in 1 tree (gray) and are shown in the upper x-axis, and ≥2 trees (black) are shown in the lower x-axis. Due to the extent of scar-dated events, the chemical record is shown from 1700 to 2020. Asterisks (*) indicate elements above p90 when ≥5 elements are over p90. RWI = Ring Width Index.
Forests 14 01082 g005
Figure 6. Standardized series of element concentrations in tree rings over time. Vertical lines extending from the x-axis indicate fire events recorded by fire scars in 1 tree (gray) and ≥2 trees (black). Due to the extent of scar-dated events, the chemical record is shown from 1700 to 2020. Asterisks (*) indicate elements above p90 when ≥5 elements are over p90. RWI = Ring Width Index.
Figure 6. Standardized series of element concentrations in tree rings over time. Vertical lines extending from the x-axis indicate fire events recorded by fire scars in 1 tree (gray) and ≥2 trees (black). Due to the extent of scar-dated events, the chemical record is shown from 1700 to 2020. Asterisks (*) indicate elements above p90 when ≥5 elements are over p90. RWI = Ring Width Index.
Forests 14 01082 g006
Figure 7. Superposed Epoch Analysis (SEA) comparing standardized “departures” of element concentrations contained in biannual portions of A. araucana rings, with fire dates recorded in two or more trees in the period 1700–2020. The x-axis represents 5 biennia, equivalent to 10 years, from 1 biennium before the fire event to 3 biennia after. The black bars and yellow bars represent significant confidence levels at 95% and 90%, respectively.
Figure 7. Superposed Epoch Analysis (SEA) comparing standardized “departures” of element concentrations contained in biannual portions of A. araucana rings, with fire dates recorded in two or more trees in the period 1700–2020. The x-axis represents 5 biennia, equivalent to 10 years, from 1 biennium before the fire event to 3 biennia after. The black bars and yellow bars represent significant confidence levels at 95% and 90%, respectively.
Forests 14 01082 g007
Figure 8. Standardized series of element concentrations in tree rings over time. Vertical lines extending from the x-axis indicate fire events recorded by fire scars in ≥2 trees (solid black lines), 1 tree (dotted black lines). Red asterisks (*) show all years in which standardized element concentrations exceed their 90 percentile. Possible fire events or periods are shown in yellow rectangles when ≥3 element concentrations are above p90 at the same biennium or after one or two biennia. Black border rectangles show three periods separated for 51 and 49 years, respectively.
Figure 8. Standardized series of element concentrations in tree rings over time. Vertical lines extending from the x-axis indicate fire events recorded by fire scars in ≥2 trees (solid black lines), 1 tree (dotted black lines). Red asterisks (*) show all years in which standardized element concentrations exceed their 90 percentile. Possible fire events or periods are shown in yellow rectangles when ≥3 element concentrations are above p90 at the same biennium or after one or two biennia. Black border rectangles show three periods separated for 51 and 49 years, respectively.
Forests 14 01082 g008
Figure 9. (A) Results of CHAR analysis of the Aguas Calientes sediment record for the last 1000 years. In the top graph, the MCHAR (charcoal accumulation rates) values are shown in gray, the interpolated Background curve in orange. The red bars in the lower graph shows the magnitude of the statistically significant fire events detected by the CHAR analysis marked with the red crosses. (B) Climatic records from the temperature reconstruction (°C) developed by [99] (blue) and PDSI reconstruction developed by [100] (dark red). (C) Dendrochemical records. The upper graph shows the result of the PC1 for eight selected elements according to the results of the SEA (Mg, Na, Ca, Mo, Ag, Sn, Ba and Pb). Ba, Mg and Ca were the elements that contributed the most to PC1 (see Figure S2).
Figure 9. (A) Results of CHAR analysis of the Aguas Calientes sediment record for the last 1000 years. In the top graph, the MCHAR (charcoal accumulation rates) values are shown in gray, the interpolated Background curve in orange. The red bars in the lower graph shows the magnitude of the statistically significant fire events detected by the CHAR analysis marked with the red crosses. (B) Climatic records from the temperature reconstruction (°C) developed by [99] (blue) and PDSI reconstruction developed by [100] (dark red). (C) Dendrochemical records. The upper graph shows the result of the PC1 for eight selected elements according to the results of the SEA (Mg, Na, Ca, Mo, Ag, Sn, Ba and Pb). Ba, Mg and Ca were the elements that contributed the most to PC1 (see Figure S2).
Forests 14 01082 g009
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Muñoz, A.A.; González, M.E.; Schneider-Valenzuela, I.; Klock-Barría, K.; Madariaga-Burgos, M.; Rodríguez, C.G.; Abarzúa, A.M.; Solari, M.E.; Martel-Cea, A.; Velásquez, B.; et al. Multiproxy Approach to Reconstruct the Fire History of Araucaria araucana Forests in the Nahuelbuta Coastal Range, Chile. Forests 2023, 14, 1082. https://doi.org/10.3390/f14061082

AMA Style

Muñoz AA, González ME, Schneider-Valenzuela I, Klock-Barría K, Madariaga-Burgos M, Rodríguez CG, Abarzúa AM, Solari ME, Martel-Cea A, Velásquez B, et al. Multiproxy Approach to Reconstruct the Fire History of Araucaria araucana Forests in the Nahuelbuta Coastal Range, Chile. Forests. 2023; 14(6):1082. https://doi.org/10.3390/f14061082

Chicago/Turabian Style

Muñoz, Ariel A., Mauro E. González, Isadora Schneider-Valenzuela, Karin Klock-Barría, Marcelo Madariaga-Burgos, Carmen Gloria Rodríguez, Ana M. Abarzúa, María Eugenia Solari, Alejandra Martel-Cea, Bárbara Velásquez, and et al. 2023. "Multiproxy Approach to Reconstruct the Fire History of Araucaria araucana Forests in the Nahuelbuta Coastal Range, Chile" Forests 14, no. 6: 1082. https://doi.org/10.3390/f14061082

APA Style

Muñoz, A. A., González, M. E., Schneider-Valenzuela, I., Klock-Barría, K., Madariaga-Burgos, M., Rodríguez, C. G., Abarzúa, A. M., Solari, M. E., Martel-Cea, A., Velásquez, B., Paredes, B., Guerrero, F., Montiel, M., Tapia-Marzán, V., Riquelme, T., & Sheppard, P. R. (2023). Multiproxy Approach to Reconstruct the Fire History of Araucaria araucana Forests in the Nahuelbuta Coastal Range, Chile. Forests, 14(6), 1082. https://doi.org/10.3390/f14061082

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop