SciELO - Scientific Electronic Library Online

 
vol.44 issue1Planting native forest species in northern Patagonian woodlands: response to canopy opening in contrasting sitesForest Management with Integrated Livestock in Argentine Patagonia: methodological adjustment for determining the baseline in complex ecosystems and heterogeneous landscapes author indexsubject indexarticles search
Home Pagealphabetic serial listing  

Services on Demand

Journal

Article

Indicators

Related links

  • On index processCited by Google
  • Have no similar articlesSimilars in SciELO
  • On index processSimilars in Google

Share


Bosque (Valdivia)

On-line version ISSN 0717-9200

Bosque (Valdivia) vol.44 no.1 Valdivia Apr. 2023

http://dx.doi.org/10.4067/s0717-92002023000100241 

Artículos

Radial growth response to long-term thinning strategies in Nothofagus pumilio forests of Tierra del Fuego

Respuesta del crecimiento radial a estrategias de raleo de largo plazo en bosques de Nothofagus pumilio de Tierra del Fuego

Julián Rodríguez-Souillaa  * 

Juan M Cellinib 

Fidel A Roigc 

María V Lencinasa 

Jimena E Chavesa 

Dardo Paredesd 

Martin Parodid 

Santiago Favorettie 

Pablo L Perif 

Guillermo J Martínez Pastura 

a Centro Austral de Investigaciones Científicas (CADIC- CONICET), Laboratorio de Recursos Agroforestales, Ushuaia, Tierra del Fuego, Argentina, tel.: 54 2901422310 int 165

b Universidad Nacional de la Plata (UNLP), Laboratorio de Investigaciones en Maderas (LIMAD), La Plata, Buenos Aires, Argentina

c Instituto Argentino de Nivología, Glaciología y Ciencias Ambientales (IANIGLA-CONICET), Mendoza, Mendoza, Argentina

d Dirección General de Desarrollo Forestal del Gobierno de la Provincia de Tierra del Fuego, Ushuaia, Tierra del Fuego, Argentina

e Universidad Nacional de Tierra del Fuego (UNTDF), Instituto de Ciencias Polares, Ambiente y Recursos Naturales, Ushuaia, Tierra del Fuego, Argentina

f Universidad Nacional de la Patagonia Austral (UNPA), Instituto Nacional de Tecnología Agropecuaria (INTA), Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), Río Gallegos, Santa Cruz, Argentina

SUMMARY

New silvicultural proposals based on long-term trials are needed to develop forest management strategies. Radial growth response was the preferred proxy for managers to define thinning levels and best years for intervention. We analyzed long-term thinning strategies in secondary Nothofagus pumilio forests (55 years old) in Tierra del Fuego (Argentina) to determine and quantify the effect on radial growth. Samplings were conducted in a trial area established in 1965 where three interventions (17, 33 and 44 years after harvesting) with different thinning strategies were applied. We employed dendrochronological techniques to measure and co-date ring widths, comparing 19 different thinning strategies (n = 100 slices x 2 reading sections). Average ring width (RW), basal area increment (BAI), and periodic annual increment every five years (PAI5) were calculated and compared through univariate comparisons. Radial growth responses were directly related to thinning intensity; however, a differential response was found according to thinning strategy (e.g. selective and systematic cuttings) and always significantly higher than the control (3 vs. 1 mm year-1). The maximum growth occurred between 2 and 4 years after thinning (YAT), and progressively decreased until reaching pre-cutting levels (5 to 8 YAT). The analyses of the thinning strategies showed that initial systematic thinning followed by at least light thinning when growth rates recover their pre-intervention values is the most convenient combination of intermediate treatments for Nothofagus pumilio forests. These results can contribute to better silvicultural decision-making in secondary forests and provide information to analyze the feasibility of their implementation.

Keywords: dendrochronology; secondary forests; selective and systematic cuttings; Patagonia

RESUMEN

Las propuestas silvícolas basadas en ensayos a largo plazo son necesarias para desarrollar estrategias de manejo forestal. La respuesta de crecimiento radial ha sido el proxy preferido por los gestores para definir los niveles de raleo y los años de intervención. Analizamos estrategias de raleo a largo plazo en bosques secundarios de Nothofagus pumilio (55 años) en Tierra del Fuego (Argentina) para determinar y cuantificar el efecto sobre el crecimiento radial. Los muestreos se realizaron en un área de ensayo establecido en 1965 donde se aplicaron tres intervenciones (17, 33 y 44 años después de la cosecha) con diferentes estrategias de raleo. Se emplearon técnicas dendrocronológicas, donde se midieron y co-fecharon anchos de anillos definiendo 19 estrategias de raleos diferentes (n = 100 rodajas x 2 secciones de lectura). Se calculó el ancho medio de los anillos (RW), el incremento del área basal (IAB) y el incremento periódico anual cada cinco años (PAI5), comparando valores mediante pruebas univariadas. Las respuestas de RW estuvieron relacionadas con la intensidad del raleo y se encontró una respuesta diferencial según la estrategia aplicada (por ejemplo, cortas selectivas y sistemáticas), siendo siempre significativamente mayor que el control (3 vs 1 mm año-1). El crecimiento máximo ocurrió entre 2 y 4 años después del raleo (ADR), y disminuyó hasta alcanzar los niveles previos a los raleos (5 a 8 ADR). Los análisis mostraron que una primera intervención sistemática seguida de al menos un raleo suave cuando las tasas de crecimiento recuperan los valores previos a la intervención, es la combinación de tratamientos intermedios más conveniente para los bosques de Nothofagus pumilio. Estos resultados pueden contribuir a una mejor toma de decisiones silvícolas en los bosques secundarios, proporcionando herramientas para analizar la viabilidad de su aplicación.

Palabras clave: dendrocronología; bosques secundarios; raleos selectivos y sistemáticos; Patagonia

INTRODUCTION

Silviculture of secondary forests is a long-term complex process, essential for developing effective forest management strategies to meet the local industry’s needs. The most important intermediate treatment is thinning, where trees can be left or removed according to different strategies, such as selective or systematic cutting that regulates the available natural resources in the managed stand (light, water, nutrients) and influences tree growth and intra-tree competition (Cruz et al. 2018). It is a challenge to establish a full-cycle prescription of interventions to optimize forest production (Allen 2008), as long-term research is necessary to define the best practical solutions for managers.

In Tierra del Fuego, Nothofagus pumilio (Poepp. et Endl.) Krasser (lenga) forests are the main resource for timber in the sawmill industry (Martínez Pastur et al. 2009), and are harvested and managed through various practices. From the 1960s to the present day, several large-scale trials have been established (Alonso et al. 1968), promoting the implementation of new silvicultural strategies (Gea-Izquierdo et al. 2004) including selective and clear-cuttings (Mutarelli and Orfila 1973), shelterwood cuts (Martínez Pastur et al. 2000) and variable retention cuts (Martínez Pastur et al. 2009). These systems open the canopy and stimulate the establishment of natural regeneration (Martínez Pastur et al. 2009) by modifying microclimate and natural cycles (Martínez Pastur et al. 2009, 2013, Soto et al. 2019). Therefore, established regeneration leads to an even-aged secondary forest with high density that is generally unmanaged (Gea-Izquierdo et al. 2004) and follows a natural dynamic with a rotation length that varies from 110 to 170 years depending on site quality (Martínez Pastur et al. 2004). Aggressive self-thinning occurs, making implementation of intermediate treatments (intermediate cuts - thinning) necessary to open the canopy and advantageously transform the secondary forest into high-productivity managed forests, reducing the time needed to produce larger and higher priced logs. In effect, growth is concentrated in selected trees, and an intermediate yield is harvested, considering that it would otherwise be lost due to the high intra-tree competition and the consequent natural mortality (Cruz et al. 2018).

Most of the harvested lenga forests were not managed after logging (Gea-Izquierdo et al. 2004). However, several theoretical proposals (e.g.Schmidt and Urzúa 1982, Uriarte and Groose 1991) and a few long-term studies were carried out to improve and manage lenga secondary forests (Martínez Pastur et al. 2001, 2013, Contreras 2004, Mundo et al. 2020). Several advantages and trade-offs were identified for thinning practices. For example, high-intensity selective thinning improves tree quality and growth, reduces forest longevity, but encourages the growth of lateral branches. Meanwhile, low intensity thinning reduces the negative visual impact, increases stability, and leaves a branch-free stem, but increases operation costs (Peri et al. 2013). Many authors described outputs in the short-term (Contreras 2004), but more long-term studies that quantify the effect of different intensities and thinning strategies are needed, such as evaluation of the periodicity between interventions and the effects of combining different thinning strategies (selective and systematic cuttings).

Radial growth response was the preferred proxy for managers to define levels of thinning and annual growth of trees through measurement of rings. Using this methodology, effects of thinning on annual radial growth have been reported for different species around the world (Misson et al. 2003, Pérez-de-Lis et al. 2011). In Tierra del Fuego, Franco et al. (2019) analyzed the response to different combinations of thinning in secondary forests of Nothofagus betuloides (Mirb.) Oerst. (guindo, coihue de Magallanes), and Mundo et al. (2020) researched the effects of two thinning levels in secondary lenga forests.

This study presents cases of long-term thinning strategies for secondary forests, however little is known about their effects nor the best recommendations for its application. Therefore, the objective was to determine and quantify the radial growth response to long-term thinning strategies in lenga forests in Tierra del Fuego and provide tools for the development of intermediate treatments in the long term. We sought to answer the following questions: (i) What is the response of trees’ radial growth to each combination and intervention strategy (selective or systematic cuttings) as a function of years-after-thinning (YAT)?; (ii) What is the growth response time of the remaining trees as a function of applied thinning intensity?; (iii) Which is the recommended strategy for secondary lenga forests based on the studied treatments? It is hypothesized that the effect of thinning on growth will be positive and highly correlated, and that response time and annual growth rates will be proportionally higher according to thinning intensity. Furthermore, the impact of different thinning strategies on growth as well as potential costs must be evaluated, as selective cuts can be more effective at stimulating individual growth but more expensive than systematic cutting.

METHODS

History of management of the study area

The study was carried out in a long-term permanent plot of the PEBANPA network (Parcelas de Ecología y Biodiversidad de Ambientes Naturales en Patagonia Austral, INTA-UNPA-CONICET, Argentina) (Peri et al. 2016). The study area is located in Tierra del Fuego province, 14 km from the town of Tolhuin, in the Aguas Blancas forests (54° 36' 35" S, 67° 15' 43" W) which is characterized by a subantarctic oceanic climate (Cfc by Köppen), an elevation of 200 meters above sea level, and the presence of inceptisols. The studied stands are pure lenga forests of medium site quality (Martínez Pastur et al. 1997) with dominant heights fluctuating between 22 and 26 m at maturity, which is reached around 200 - 250 years old (Schmidt and Urzúa 1982).

Site history

In the 1960s, the University of Buenos Aires (UBA) began a series of research studies in Andean-Patagonian forests with the purpose of evaluating the necessary elements to ensure the rational management of these native forests. Through these studies, they aimed to establish the most convenient silvicultural treatments and the best timing and intensity of the interventions in order to grow the productive stands to harvesting age with the maximum possible value according to the technological criteria to balance growth and quality (Mutarelli and Orfila 1973). These investigations were supported by the National Forestry Administration (ANB) and later by the National Forestry Institute (IFONA), and were gradually abandoned, as the last monitoring and trials were carried out in the mid-1980s.

Between 1965 and 1967, the ANB and UBA established one permanent plot at the Aguas Blancas forests in 16 ha of unmanaged lenga forests with high homogeneity in terms of composition and structure of vegetation, soil type, and topography. This forest was harvested for timber purposes using clear cuts in alternating NE-SW oriented strips, 35 meters wide and 400 meters long (Cozzo et al. 1967). Within the harvested strips, all the woody material was removed (> 5 cm), allowing the establishment of natural abundant regeneration of the harvested areas. During the following years, the remnant strips were partially harvested (selective cuts), and then the entire area was abandoned and no management was conducted for 15 years. In 1984, IFONA carried out different commercial thinning systems on harvested strips. The regeneration presented forests of approximately 2 meters of height, 2.5 cm DBH, and a tree density of 80,000 - 100,000 ind ha-1 (Mundo et al. 2020). These start-up trials were applied at one per strip without replications (selective or systematic cuttings): (i) heavy thinning from below (HT), (ii) light thinning from below (LT), (iii) extraction of competitors through selection thinning (S), (iv) 2 meter strips parallel to the mother strips (longitudinal strips, LS), (v) 2 meters strips perpendicular to the mother strips (transverse strips, TS), and (vi) checkerboards (with parallel and perpendicular strips every 2 meters, D) (figure 1). Control areas (C) without cuttings were left in every strip.

Figure 1 . Start-up thinning treatments (1984) in Aguas Blancas permanent plots. A = control, B = heavy thinning, C = light thinning, D = selection, E = selection strip by competitors, F = transversal strips, G = longitudinal strips, H = checkerboard. Scale in meters. Tratamientos de la primera intervención (1984) en las parcelas permanentes de Aguas Blancas. A = control, B= raleo fuerte, C = raleo suave, D = selección, E = extracción de competidores, F = franjas transversales, G = franjas longitudinales, H = damero. Escala en metros. 

New monitoring and management of the permanent plot was conducted starting in 1999, led by Centro Austral de Investigaciones Científicas (CADIC CONICET), Facultad de Ciencias Agrarias y Forestales (UNLP), the Forest Office (Government of Tierra del Fuego) and private companies. New trials and cuttings for thinning were applied in 1999 - 2000, 2010, and 2020 in the harvested strips. During 1999 - 2000, new subtreatments were installed without replications of the analyzed treatments, and one new treatment was established in an undisturbed sector. In the start-up trials, the areas were subdivided into three sectors, where subtreatments were defined during the second thinning: (i) heavy thinning (> 50 % basal area was removed), (ii) light thinning (< 50 % basal area was removed), and (iii) a thinning control, where no new cuttings were applied on the areas thinned in 1985. Both thinning levels removed individuals of low crown class with poor health and deformations in the stem, prioritizing an equivalent distance among remnant individuals. The objective of these treatments was to promote the growth of remnant trees by increasing rainfall and temperature at the understory level, as well as reducing competition for nutrients.

The new trial established in 2000 defined a new thinning strategy based on selective cuttings in two stages (T2000). For this purpose, target trees were selected across the plot (3 x 3 m; 1,111 trees ha-1). During the first stage all the crown competitors of the target trees were removed (year 2000), and in the second stage the remaining non-target trees were removed (year 2010). The target trees were dominant and presented good health and forest stem quality. The objective of this treatment was to promote the crown development of the target trees before the final cut, in which the canopy was abruptly opened. In 2010, the third intervention was conducted, where all subtreatments were thinned in order to homogenize remnant forest structure. Finally, in 2021, the fourth intervention was implemented, applying heavy thinning (50 - 60 % of the BA was removed), removing individuals with poor health and deformations in the stem, and prioritizing an equivalent distance among remnant individuals (not analyzed in this study). This thinning homogenizes the area without consideration of the previously established subtreatments of each harvested strip. A total of 19 thinning strategies were considered. Evolution over years of treatments and its subdivisions are presented in table 1, representing thinning intensities as compared to BA of the control. Combinations with their abbreviations are presented in annex 1.

Table 1 Classification of treatments and subtreatments. Percentages represent the proportion of BA remnants compared to C. Clasificación de los tratamientos y subtratamientos. Los porcentajes representan la proporción de AB remanente en comparación con C.  

Where: T17= thinning 17 YAH (first, 1984), T33= thinning 33 YAH (second, 1999), T44= thinning 44 YAH (third, 2010), T54= thinning 54 YAH (fourth, 2021). H= heavy thinning, less than 50% of BA of Control. L= light thinning, more than 50% of BA of Control. No= no thinning applied. Row values of C represent BA (m2.ha-1). PRE= values previous to thinning. POST= values post thinning. C= control, HT= heavy thinning, LT= light thinning, S= selection, TS= transversal strips, LS= longitudinal strips, D= checkerboard, T2000= thinning year 2000.

Sampling and dendrochronological analysis

Forest structure was measured for all treatments prior to cuttings and samplings. We used the Bitterlich point sampling method with a forest caliper and a Criterion RD-1000 (Laser Technology, US), applying a K coefficient between 1 and 6. We obtained: (i) height (H, m), (ii) basal area (AB, m2 ha-1), (iii) tree density (individuals ha-1), and (iv) average diameter at breast height (DBH, cm). After cutting, slides were obtained from harvested trees during thinning in 2021. A total of 100 slices at DBH were collected considering trees of good quality (straight shafts without bifurcations), DBH of 17 - 23 cm and dominant heights of up to 16 meters (6 treatments x 3 subtreatments x 5 replicates + 5 replicates of C + 5 replicates of T2000), resulting in the 19 different thinning strategies and combinations of treatments across the years reported in table 1. The slides were air-dried and polished in order to obtain smooth surfaces to facilitate the identification of growth rings under magnification (10x - 50x).

Cross sections were processed using standard dendrochronological techniques (Stokes and Smiley 1968) and two perpendicular radii were measured per sample. Growth rings were dated according to the Schulman convention (Schulman 1956) for the southern hemisphere, and the year of each growth ring was assigned according to the year it began to form. Ring widths were measured with a digital caliper to the nearest 0.01 mm. COFECHA software was used to evaluate the quality of the dating and to calculate Mean Intercorrelation (MI) and Mean Sensitivity (MS). For analyses, average ring width (RW), basal area increment (BAI) assuming circular sections, and periodic annual increment (PAI) were calculated, using the average of 5-year periods before and after the interventions, grouping dendrochronological series in treatments and subtreatments. The difference between the means of PAIpre and PAIpost was calculated by statistical tests to compare growth before and after the interventions. Thus, the periods were established as such: for the 1984 intervention (first one, 17 YAH), pre (1980 - 1984), post (1985 - 1989); for 1999 (second one, 33 YAH), pre (1995 - 1999), post (2000 - 2004); for 2010 (third one, 44 YAH), pre (2006 - 2010), post (2011 - 2015). The year 2020 was discarded because cuts were conducted during ring growth. The Kruskal-Wallis non-parametric test was applied to detect differences between means with a significance level of P < 0.05, and the Mann-Whitney U test was applied for comparison between means. Statistical analyses were performed using STATGRAPHICS Centurion XVI.I software (USA, StatPoint Technologies, Inc). Further, the response time of the interventions was estimated as the number of years with RW values greater than the average of the 5 years prior to thinning.

RESULTS

Regarding forest structure in 2021 (table 2), the lowest average DBH value corresponded to C (11.2 cm), but with the highest basal area (52.2 m2 ha-1) and density (7,420 ind ha-1). The treatments presented average DBH between 18.3 and 28.1 cm, with the highest value corresponding to LS followed by D (21.9 cm) and the lowest to S. Basal area fluctuated between 30.4 and 34.9 m2 ha-1 and density between 919 and 1,239 ind ha-1 in a treatments thinned. Regarding T2000, and considering that its structure was the result of two thinning treatments, it reached lower values of basal area (28.0 m2 ha-1) and average DBH (18.4 cm) than the other treatments with three interventions. Heights fluctuated between 13.83 and 16.46 m, where higher values corresponded to LS, LT and D. The highest gross BAI corresponded to HTHT and DHT, being the subtreatments with the biggest trees (so highest response to thinning). No mortality was reported in thinned sites, but it was reported in control sites due to natural stand dynamics (self-thinning, 6 % annually on average). Exceptional cases of wind throw were reported. The slides sampled presented a diameter of 22.5 ± 6.3 cm, considering the different sizes of samplings.

Table 2 Forest structure in 2021 for subtreatments. 

Where: Treatment (Tr), basal area (BA, m2 ha1), density (DEN, ind ha-1), diameter at breast height -130cm- (DBH, cm), quadratic mean diameter (QMD, cm), Gross BAI (sum of BAI from 45 to 54 YAH -after 3rd thinning-, cm2).

The dendrochronological series showed a year of onset that varied between 1960 and 1977, considering the initial age at DBH, and indicated that several individuals were probably part of the seedling bank prior to clear-cutting. Likewise, there was evidence of an established cohort after harvesting in all treatments. Regarding the common variability between individual series depending on treatment, the average MI prior to the first intervention was relatively low (0.23 across all series, varying between 0.23 and 0.53 if divided by treatments), increasing to 0.52 after the first intervention, 0.62 after the second, and 0.71 after the third. The number of series varied among treatments according to the age of the trees studied (table 3).

Table 3 Ring width statistics for periods between interventions. Estadísticas de ancho de anillos para los períodos entre intervenciones.  

Where: Treatment (Tr), First year (Fy), mean intercorrelation (MI), mean sensitivity (MS), number of series (NS), subtreatment (Sub-T). C= control, HT= heavy thinning, LT= light thinning, S= selection, TS= transversal strips, LS= longitudinal strips, D= checkerboard, T2000= selection strip by competitors.

The analysis of the effects of thinning using δPAI comparisons showed positive effects on radial tree growth for the 3 interventions. The greatest effects were observed for the first and third thinning treatments, with values close to 0 in δPAIRW for the second thinning, and negative values for the control. Similarly, significant differences were observed when treatments were compared between the 3 interventions, finding for all cases that the control presented the lowest values, negative for the second and third interventions (figure 2). For the first intervention, the treatment with the greatest effect corresponds to HT, showing significant differences with the rest of the treatments when δPAIRW (P < 0.001) and δPAIBAI (P < 0.001) were compared. Analyzing the results of δPAIBAI for the second intervention, the different thinning treatments did not show significant differences. However, when analyzing the third intervention, significant differences between treatments both for δPAIRW (P < 0.001) and δPAIBAI (P < 0.001) were found, where the largest effects correspond to T2000, HT and systematic treatments.

Figure 2 . Differences in periodic annual increment (δPAI) for 5 years before and after 3 interventions for ring width (RW, left panels) and basal area increment (BAI, right panels). A (A1 and A2): for first thinning (1980 - 1984) - (1985 - 1989). B (B1 and B2): for second thinning (1995 - 1999) - (2000 - 2004). C (C1 and C2): for third thinning (2006 - 2010) - (2011 - 2015). Bars represent mean values, with their corresponding standard errors of the mean. Different letters show significant differences by the Mann-Whitney U test (P < 0.05 based on Kruskal-Wallis test). C = control, HT = heavy thinning, LT = light thinning, S = selection, TS = transversal strips, LS = longitudinal strips, D = checkerboard, T2000 = selection strip by competitors. Diferencias en el incremento periódico anual (δPAI) durante 5 años antes y después de 3 intervenciones para anchos de anillo (RW, paneles de la izquierda) y el incremento del área basal (BAI, paneles de la derecha). A (A1 y A2): para 1ra intervención (1980 - 1984) - (1985 - 1989). B (B1 y B2): para 2da intervención (1995 - 1999) - (2000 - 2004). C (C1 y C2): para 3ra intervención (2006 - 2010) - (2011 - 2015). Las barras representan la diferencia de medias, con su correspondiente error estándar de la media. Letras diferentes indican diferencias significativas de acuerdo a Mann-Whitney U test (P < 0,05 según la prueba de Kruskal-Wallis). C = control, RF = raleo fuerte, RS = raleo suave, S = selección, TS = franjas transversales, LS = franjas longitudinales, D = damero, T2000 = extracción de competidores. 

The evolution of the effects after the first thinning (figure 3, 17 YAH) is consistent with those previously mentioned, showing similar response patterns for all treatments. HT was the treatment with the highest response and growth for RW and IAB, considering also that at the time of thinning it started from lower values than the rest of the treatments. This was the only one that differed from systematic treatments. When compared with C, all treatments followed similar patterns, and after 25 YAH, RW values kept constant, but with higher rates (excepting for S).

Figure 3 . Mean RW (A) and BAI (B) mean series for treatments for period 3 to 33 YAH. The intervention year is marked with dashed lines (17 YAH, 1st intervention). HT = heavy thinning, LT = light thinning, S = selection, TS = transversal strips, LS = longitudinal strips, D = checkerboard, C = control. Series medias de RW (A) y BAI (B) para los tratamientos del período 3 a 33 YAH. El año de intervención está marcado con líneas discontinuas (1984, 1ra intervención). HT = raleo fuerte, LT = raleo suave, S = selección, TF = franjas transversales, LF = franjas longitudinales, D = damero, C = control. 

When subtreatments (2nd thinning) were applied, the responses were differential according to the intensity of thinning applied (figure 4). First, all thinned subtreatments showed positive effects and differed in C growth values. For heavy thinning, the highest response was for LT (subT = LTHT), followed by HT (subT = HTHT). Systematic treatments showed similar patterns, with D (subT = DHT) showing the greatest growth rates, not only after the second thinning, but also after the third one. For light thinning, systematic treatments differed from non-systematic ones in higher growth, except HT (subT = HTLT). For controls, values remained constant until the application of the third thinning (44 YAH), where all subtreatments showed similar effects. However, when compared with subtreatments that were applied in 1999, RW and IAB values were lower for controls, highlighting the successive effect of thinning application. Finally, T2000 highlights a behavior similar to C, modifying this pattern considerably at the time of the second intervention (1999), reaching similar rates to the other subtreatments.

Figure 4 . RW (left) and BAI (right) mean series for subtreatments for period 33 to 53 YAH, classified as heavy thinning (A1 and A2), light thinning (B1 and B2) and control (C1 and C2). Intervention years are marked with dashed lines (33 YAH, 2nd intervention; 44 YAH, 3rd intervention). HT = heavy thinning, LT = light thinning, S = selection, TF = transversal strips, LS = longitudinal strips, D = checkerboard, T2000 = selection strip by competitors, C = control. Series mediasde RW (izquierda) y BAI (derecha) para los subtratamientos para el periodo 33 a 53 YAH, clasificado en raleo fuerte (A1 y A2), raleo suave (B1 y B2) y control (C1 y C2). Los años de intervención están marcados con líneas discontinuas (1999, 2da intervención; 2010, 3ra intervención). HT = raleo fuerte, LT = raleo suave, S = selección, TF = franjas transversales, LF = franjas longitudinales, D = damero, C = control, T2000 = extracción de competidores. 

The response time (lag) of the positive effects of the first intervention varied between 1 and 13 YAT, with variability according to the treatment applied and the % BA removed. In this case, the greatest effect was observed in HT, followed by D and LT. Regarding the effects of the second intervention, the application of subtreatments highlighted the difference between controls and interventions. In the controls, growth rates remained below the pre-intervention average. In contrast, heavy and light thinning for each treatment had positive effects on radial tree growth in all cases, with the year of greatest growth ranging from the third to the seventh year after thinning, with lags of 6 to 10 YAT (table 4). For all treatments, the average ring widths post-intervention were larger than before the intervention, except for HT pruned. Finally, for the third intervention, where all the subtreatments were applied, the effects were positive in all cases, and the largest growth occurred between the second and third year after thinning, with a lag between 7 and 10 YAT, considering that at the time of the last intervention (4th thinning, 2021) the trees still presented higher average rates than those prior to the intervention. In the case of C, a drop in average growth rates was observed throughout the entire period studied.

Table 4 Treatment (Tr), intensity of thinning represented as percentage of basal area thinned (% AB), average of ring width (mm) 5 years before thinning (XRW pre), average of ring width (mm) 5 years after thinning (XRW post), year of highest growth and value (Year (Max)), response time (Lag), subtreatment (Sub-T). Tratamiento (Tr), intensidad del raleo representado como porcentaje del área basal raleado (% AB), media de anchos de los anillos (mm) 5 años antes del raleo (XRW pre), media de anchos de los anillos (mm) 5 años después del raleo (XRW post), año de mayor crecimiento y valor (Año (Max)), tiempo de respuesta (Lag), subtratamiento (Sub-T).  

DISCUSSION

As expected, the results of the present work show that the effects of thinning are positive, considering that they generate more growth in the remaining forest stand when compared with control sites. Additionally, the evaluation of the effects by statistical analysis and interpretation of graphs seems to be more sensitive when BAI is applied. The present findings are consistent with those of Mundo et al. (2020), who conducted a study at the same site but applied only HT and its subtreatments. In the present work, the study was expanded with the incorporation of a new intervention and the remaining treatments, comparing systematic and non-systematic interventions. The mean intercorrelation values obtained are similar to those presented by Mundo et al. (2020), highlighting the increase of these values as the interventions elapse. These values indicate that thinning produces a homogeneous growth effect among the remaining trees in the stand, a condition that can be modified when competition for resources increases.

Analysis of the effects of thinning shows that more intense interventions generate more annual radial growth and a longer response period. Similar results were found by Pérez-de-Lis et al. (2011) where the effects of heavy thinning were greater than those of light thinning for Pinus canariensis Sweet ex Spreng. This is due to fewer individual trees that are able to take advantage of more resources for a longer period (McDowell et al. 2003). The time limit seems to be the moment of crown closure, when greater inter-tree competition for light is triggered, and annual growth rates return to pre-intervention values. Annual growth surveys could allow for determination of opportune moments for interventions, and the results of the present work could be used as a reference, considering the varying situations studied. However, the intensity of thinning can also have consequences on the stability of the remaining canopy (Peri et al. 2002) and the excessive formation of lateral branches (Martínez Pastur et al. 2001).

When the treatments are analyzed separately, a gradual increase in growth rates is observed, reaching a peak that varies between year 2 and 6 after the interventions. Similar patterns of the effects of thinning were found by Franco et al. (2019) for Nothofagus betuloides. It is a challenge to determine whether it is best to apply thinning once these maximum values are reached, or at the time when previous values are recovered. Economic aspects and social context that are not considered in this work should be discussed.

It may be considered that the second thinning should have been implemented with greater intensity to observe more effects on the trees. Thus, it is estimated that four thinning treatments are sufficient to complete the silvicultural system in a rotation of approximately 80 years, maintaining radial growth rates of ≈ 2 mm per year (finally, sawn logs with a diameter of 40 cm, forest density of ≈ 600 ind ha-1). These intermediate treatments are suitable for application after logging (in this case, 17 YAH, but it could be years later, depending on site quality and height of saplings), both for shelterwood cuts and retention strategies. Furthermore, thinning generates products (smaller than timber products) that can potentially be used by farmers without having to wait for the final felling year. Considering differences between systematic and non-systematic interventions, a primary systematic intervention (scripts or checkerboarded) seems to be successful in terms of logistics and effects on growth rates, which are similar to non-systematic treatments (except for HT), as evidenced in figure 2. Pruning does not make a big difference in radial growth after the second intervention, but it does after the third one, comparable with results found by Mundo (2020), due to more carbon accumulation in the stem and not in the branches. Therefore, other interventions should be conducted, considering that the effects of thinning decay and growth rates recover to previous values 6 - 10 years after thinning. It is then appropriate to apply a second thinning treatment in which about 50 % of the remaining basal area is removed from those that were thinned systematically in the first intervention. In this way, canopy release is achieved, which was not completed during the first intervention. If second thinning is not applied, effects on systematic treatments will be lower than non-systematic ones after third interventions, as evidenced in figure 4. This combination of interventions generates similar effects at the time of a third thinning treatment, when compared to originally non-systematic thinning. In addition, in the year 2020 (54 YAH), these combinations reached higher average DBH values. The threshold of interventions results in a combination of periods between thinning and intensities applied, reaching a final density before shelterwood cuts. In summary, a combination similar to DHT and DLT appears to be the most convenient to apply, considering both growth rates and operability. Considering that the results of this study are from experimental plots, it is expected to be useful as a tool for real application on a larger scale, taking advantage of the abundance of secondary forest resources available in the province of Tierra del Fuego.

Finally, repeated situations in all the subtreatments were observed in 2003 (31 YAH), 2011 (45 YAH) and 2018 (52 YAH), moments in which the growth rates decreased noticeably and punctually, recovering previous growth values in the following years. This can be explained by years with adverse climatic conditions (e.g. low rainfall) (Franco et al. 2019) or mast years where photo assimilates are mostly used to synthesize reproductive organs instead of carbon fixation in growth rings (Hadad et al. 2021). Although it has been shown that the application of intermediate treatments generates effects on ring widths, annual variation due to climatic variations cannot be eliminated by these silvicultural practices (Spiecker 2002). Indeed, these variations must be considered to characterize the effectiveness of thinning after interventions.

CONCLUSIONS

Results show that thinning generates positive effects on the radial growth of remnant trees in secondary forests of Nothofagus pumilio. These interventions produce releases of remnant trees in all treatments, and the effects are quantified as compared to control sites. Indeed, the application of dendrochronological techniques makes it possible to demonstrate the periodic impact of thinning. According to analysis of these effects, non-systematic interventions did not differ from systematic ones, so the application of the latter is considered more convenient, though it is also necessary to apply others to maintain higher growth rates than control sites. Annual growth surveys could allow for determination of opportune moments for interventions, and the values of the present work could be considered as a reference, taking into account the combined effect of the situations studied. Finally, these results can contribute to better silvicultural decision-making for sustainable forest management in secondary forests, improving the management of resources and their abundance. It is important to continue research in such a valuable sample site, as it could be a model for application in other sites and Nothofagus forests.

ACKNOWLEDGMENTS

To the researchers, technicians, students, and landowners (ranch and sawmill companies) that support this research. It was impossible to obtain this invaluable data without their assistance

REFERENCES

Allen HL. 2008. Silvicultural Treatments to Enhance Productivity. In Evans J ed. The Forests Handbook vol 2: Applying Forest Science for Sustainable Management. New Jersey, USA. Blackwell Science Ltd. p 129-139. DOI: https://doi.org/10.1002/9780470757079.ch6Links ]

Contreras C. 2004. Evaluación del efecto de un raleo selectivo sobre el crecimiento en un renoval de Nothofagus pumilio (Poepp. et Endl.) Krasser, en Tierra del Fuego, XII Región. Forestry Engineer Thesis. Santiago, Chile. Escuela de Ciencias Forestales. Universidad de Chile. 49 p. [ Links ]

Cozzo D, E Mutarelli, E Orfila. 1967. Plan de Investigaciones silviculturales y dasonómicas necesarias para la organización económica de los bosques subantárticos argentinos. Segundo Informe. Plan nº 118. Buenos Aires, Argentina. Convenio Cátedra de Dasonomía - UBA y CAFPTA. 260 p. [ Links ]

Cruz G, FA Rodriguez, PA Tapia, A Patricio, HE Bown. 2018. Growth responses after a future crop tree thinning and a thinning from below in a second-growth Nothofagus pumilio forest in Tierra del Fuego, Chile. Ciencia e Investigación Agraria 45: 263-276. DOI: https://dx.doi.org/10.7764/rcia.v45i3.1856Links ]

Franco MG, IA Mundo, G Martínez Pastur, MD Barrera. 2019. Radial growth responses to thinning and climate in native Nothofagus betuloides forests in Tierra del Fuego, Argentina. Dendrochronologia 57: e125625. DOI: https://doi.org/10.1016/j.dendro.2019.125625 Links ]

Gea-Izquierdo G, G Martínez Pastur, JM Cellini, MV Lencinas. 2004. Forty years of silvicultural management in southern Nothofagus pumilio (Poepp. et Endl.) Krasser primary forests. Forest Ecology and Management 201: 335-347. DOI: https://doi.org/10.1016/j.foreco.2004.07.015Links ]

Hadad M, FA Roig, JG Arco Molina, A Hacket-Pain. 2021. Growth of male and female Araucaria araucana trees respond differently to regional mast events, creating sex-specific patterns in their tree-ring chronologies. Ecological Indicators 122(0). 107245. DOI: https://doi.org/10.1016/j.ecolind.2020.107245Links ]

Martínez Pastur G, PL Peri, R Vukasovic, S Vaccaro, V Piriz Carrillo. 1997. Site index equation for Nothofagus pumilio Patagonian forest. Phyton 6(1/2): 55-60. URL: https://agris.fao.org/agris-search/search.do?recordID=US201302880251Links ]

Martínez Pastur GJ, JM Cellini, PL Peri, R Vukasovic, C Fernández. 2000. Timber production of Nothofagus pumilio forests by a shelterwood system in Tierra del Fuego (Argentina). Forest Ecology and Management 134(1-3): 153-162. DOI: https://doi.org/10.1016/S0378-1127(99)00253-4Links ]

Martínez Pastur G, JM Cellini, MV Lencinas, R Vukasovic, R Vicente, F Bertolami, J Giunchi. 2001. Modificación del crecimiento y de la calidad de fustes en un raleo fuerte de un rodal en fase de crecimiento óptimo inicial de Nothofagus pumilio (Poepp. et Endl.). Ecología Austral 11: 95-104.URL: http://ojs.ecologiaaustral.com.ar/index.php/Ecologia_Austral/article/view/1574Links ]

Martínez Pastur GJ, MV Lencinas, R Vukasovic, PL Peri, B Díaz, JM Cellini. 2004. Turno de corta y posibilidad de los bosques de lenga (Nothofagus pumilio) en Tierra del Fuego (Argentina). Bosque 25(1): 29-42. DOI: http://dx.doi.org/10.4067/S0717-92002004000100003 Links ]

Martínez Pastur GJ, MV Lencinas, JM Cellini. 2009. Timber management with variable retention in Nothofagus pumilio forests of Southern Patagonia. Forest Ecology and Management 258(4): 436-443. DOI: https://doi.org/10.1016/j.foreco.2009.01.048Links ]

McDowell N, JR Brooks, SA Fitzgerald, BJ Bond. 2003. Carbon isotope discrimination and growth response of old Pinus ponderosa trees to stand density reductions. Plant, Cell & Environment 26(4): 631-644. DOI: https://doi.org/10.1046/j.1365-3040.2003.00999.xLinks ]

Misson L, C Vincke, F Devillez. 2003. Frequency responses of radial growth series after different thinning intensities in Norway spruce (Picea abies (L.) Karst.) stands. Forest Ecology and Management 177: 51-63. DOI: https://doi.org/10.1016/S0378-1127(02)00324-9Links ]

Mundo IA, DA Palazzini, JA Barotto, GJ Martínez Pastur, MD Barrera. 2020. Dendroecology applied to silvicultural management in the Southern Patagonian forests: A case of study from an experimental forest in Tierra del Fuego, Argentina. In Pompa-García M, JJ Camarero eds. Latin American Dendroecology: Combining tree-ring sciences and ecology in a megadiverse territory. Cham, Switzerland. Springer. p. 317-330. DOI: https://doi.org/10.1007/978-3-030-36930-9_14Links ]

Mutarelli EJ, E Orfila. 1973. Algunos resultados de las investigaciones de manejo silvicultural que se realizan en los bosques andino-patagónicos de Argentina. Revista Forestal Argentina 17(3): 69-75. [ Links ]

Pérez-de-Lis G, I García-González, V Rozas. 2011. Effects of thinning intensity on radial growth patterns and temperature sensitivity in Pinus canariensis afforestations on Tenerife Island, Spain. Annals of Forest Science 68: 1093-1104. DOI: https://doi.org/10.1007/s13595-011-0125-8 Links ]

Peri PL, G Martínez Pastur, R Vukasovic, B Díaz, MV Lencinas, JM Cellini. 2002. Thinning schedules to reduce risk of windthrow in Nothofagus pumilio forests of Patagonia, Argentina. Bosque 23(2): 19-28. DOI: https://dx.doi.org/10.4067/S0717-92002002000200003Links ]

Peri PL, GM Martínez Pastur, L Monelos. 2013. Natural dynamics and thinning response of young lenga (Nothofagus pumilio) trees in secondary forests of Southern Patagonia. Bosque 34(3): 273-279. DOI: http://dx.doi.org/10.4067/S0717-92002013000300003 Links ]

Peri PL, MV Lencinas, J Bousson, R Lasagno, R Soler, HA Bahamonde, G Martínez Pastur. 2016. Biodiversity and ecological long-term plots in southern Patagonia to support sustainable land management: the case of PEBANPA network. Journal of Nature Conservation 34: 51-64. DOI: https://doi.org/10.1016/j.jnc.2016.09.003Links ]

Schulman E. 1956. Dendroclimatic changes in semiarid America. Tucson, USA. University of Arizona Press. 142 p. [ Links ]

Spiecker H. 2002. Tree rings and forest management in Europe. Dendrochronologia 20(1-2): 191-202. DOI: https://doi.org/10.1078/1125-7865-00016 Links ]

Soto DP, JK Puettmann, C Fuentes, DF Jacobs. 2019. Regeneration niches in Nothofagus-dominated old-growth forests after partial disturbance: Insights to overcome arrested succession. Forest Ecology and Management 445: 26-36. https://doi.org/10.1016/j.foreco.2019.05.004Links ]

Stokes M, T Smiley. 1968. An introduction to tree-ring dating. Chicago, USA. University of Chicago Press. 73 p. [ Links ]

Schmidt H, A Urzúa. 1982. Transformación y manejo de los bosques de lenga en Magallanes. Ciencias Agrícolas Nº 11. Santiago, Chile. Facultad de Ciencias Agrarias, Veterinarias y Forestales, Universidad de Chile. 62 p. [ Links ]

Uriarte A, H Grosse. 1991. Informe ténico 126: Los bosques de lenga. Una orientación para su uso y manejo; recopilación bibliográfica. Concepción, Chile. CORFO-INFOR. 103 p. DOI: https://doi.org/10.52904/20.500.12220/6608Links ]

1AUTHOR CONTRIBUTIONS Guillermo Martinez Pastur, Juan Manuel Cellini and Pablo Luis Peri conceived and designed the experiments; Guillermo Martinez Pastur, María V Lencias, Dardo Paredes, Jimena E Chavez, Fidel A Roig and Julián Rodríguez-Souilla collaborated in fieldwork and data analysis, and all the authors collaborated in writing the manuscript. All authors have read and agreed to the published version of the manuscript.

FUNDING This research was funded by PDTS-0398 (2020-2023) “Manejo sostenible de los bosques de Nothofagus y ambientes naturales de Tierra del Fuego: Compatibilizando la producción y la conservación de la biodiversidad” of MINCyT (Argentina)

Annex 1. Combinations: C = control, T2000 = thinning year 2000, HTC= heavy thinning - control, HTHT= heavy thinning- heavy thinning, HTLTP= heavy thinning - light thinning - pruned, LTC= light thinning - control, LTLT= light thinning - light thinning, LTHT= light thinning - heavy thinning, SC= selection - control, SLT= selection - light thinning, SHT= selection - heavy thinning, LSC= longitudinal strips - control, LSLT = longitudinal strips - light thinning, LSHT= longitudinal strips - heavy thinning, TSC= transversal strips - control, TSLT= transversal strips - light thinning, TSHT= transversal strips - heavy thinning, DC= checkerboard - control, DLT= checkerboard - light thinning, DHT= checkerboard - heavy thinning. Combinaciones: C= control, T2000= extracción de competidores, HTC= raleo fuerte - control, HTHT= raleo fuerte - raleo fuerte, HTLTP= raleo fuerte - raleo suave - podado, LTC= raleo suave - control, LTLT= raleo suave - raleo suave, LTHT= raleo suave - raleo fuerte, SC= selección - control, SLT= selección - raleo suave, SHT= selección - raleo fuerte, LSC= fajas longitudinales - control, LSLT= fajas longitudinales - raleo suave, LSHT= fajas longitudinales - raleo fuerte, TSC= fajas transversales - control, TSLT= fajas transversales - raleo suave, TSHT= fajas transversales - raleo fuerte, DC= damero - control, DLT= damero - raleo suave, DHT= damero - raleo fuerte

Received: August 03, 2022; Accepted: January 30, 2023

* Corresponding author: j.rodriguez@conicet.gov.ar

CONFLICTS OF INTEREST The authors declare no conflict of interest

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License