Spatiotemporal pattern and driving factors of leaf area index in a karst old-growth forest in southwest China

Abstract

Leaf area index (LAI) is a critical parameter for characterizing canopy structures. The spatiotemporal patterns of LAI reflect the multi-scale adaptive mechanisms of forests, spanning from individual leaf expansion and population structure optimization to coordinated ecosystem functionality. However, research on the spatiotemporal dynamics of LAI remains limited, particularly in karst primary forests that feature complex community structures, high rock exposure rates, and heterogeneous soil distribution patterns. Our research focused on a 1.28 ha karst forest plot in southwest China, where we quantified LAI seasonality using the LAI-2200 analyzer during four critical phenophases: January (dormant period), April (growth period), July (maturation period) and October (senescence period). Structural equation modeling (SEM), variance partitioning and hierarchical partitioning analysis were employed to quantify the relative contributions of biological, abiotic and spatial factors (represent community spatial connectivity shaped by ecological processes including seed dispersal, resource distribution and competition) to LAI variation. The results revealed distinct seasonal fluctuations in LAI, with significant spatial autocorrelation occurring at characteristic scales of 11.4 m (dormant period), 20.1 m (growth period), 66.3 m (maturation period), and 19.2 m (senescence period), respectively. Spatial factors were the dominant factors in explaining LAI variation across growth, maturation and maturation periods, accounting for 32% to 80% of the explained variance. During the dormancy period, however, the proportion of deciduous trees was the most important factor (32%). The spatial factors exhibited a positive correlation with LAI during the growth and dormancy periods but showed an inverse relationship in other periods. Biological drivers differentially regulated LAI across growth, senescence and dormancy period. Tree height and deciduous species proportion reduced LAI, whereas crown width, density of small trees and stand density increased it. Potassium (inhibiting) and phosphorus (promoting) elements were the primary soil elements influencing LAI. Topographic factors (altitude, slope position and slope aspect) affected LAI through direct, indirect, or combined effects, with slope aspect being the dominant topographic factor. In summary, spatial distribution of LAI exhibited seasonal dependence, which was influenced by spatial, stand, soil and topographic factors, etc. This study identified the primary drivers of LAI in karst evergreen and deciduous broadleaf forests during different seasons, and simultaneously demonstrated the importance of long-term forest monitoring and targeted management strategies, while highlighting the necessity of incorporating spatial factors into LAI estimation and modeling frameworks.

Keywords

Karst forest / Leaf area index / Spatial autocorrelation / Structural equation modeling

References

[1] Ali A, Lin SL, He JK, Kong FM, Yu JH, Jiang HS. Big-sized trees overrule remaining trees’ attributes and species richness as determinants of aboveground biomass in tropical forests. Glob Change Biol. 2019, 25(8): 2810-2824.

[2] Bequet R, Campioli M, Kint V, Muys B, Bogaert J, Ceulemans R. Spatial variability of leaf area index in homogeneous forests relates to local variation in tree characteristics. For Sci. 2012, 58(6): 633-640.

[3] Bequet R, Kint V, Campioli M, Vansteenkiste D, Muys B, Ceulemans R. Influence of stand, site and meteorological variables on the maximum leaf area index of beech, oak and Scots pine. Eur J Forest Res. 2012, 131(2): 283-295.

[4] Bequet R (2011) Environmental determinants of the temporal and spatial variability in leaf area index of Fagus sylvatica L., Quercus robur L. and Pinus sylvestris L.

[5] Blanken PD, Black TA. The canopy conductance of a boreal aspen forest, Prince Albert National Park, Canada. Hydrol Process. 2004, 18(9): 1561-1578.

[6] Brinkhoff R, Mendham D, Hunt M, Unwin G, Hovenden M. Nitrogen fertiliser only increases leaf area in the lower crown of mid-rotation Eucalyptus nitens plantations. For Ecol Manag. 2022, 508. 120048

[7] Chen JM. Optically-based methods for measuring seasonal variation of leaf area index in boreal conifer stands. Agric for Meteorol. 1996, 80(2–4): 135-163.

[8] Chen JM, Black TA. Defining leaf area index for non-flat leaves. Plant Cell Environ. 1992, 15(4): 421-429.

[9] Chen JM, Cihlar J. Plant canopy gap-size analysis theory for improving optical measurements of leaf-area index. Appl Opt. 1995, 34(27): 6211-6222.

[10] Chen YF, Yu FH, Dong M. Scale-dependent spatial heterogeneity of vegetation in Mu Us sandy land, a semi-arid area of China. Plant Ecol. 2002, 162(1): 135-142.

[11] Dai FQ, Zhou QG, Lv ZQ, Wang XM, Liu GC. Spatial prediction of soil organic matter content integrating artificial neural network and ordinary Kriging in Tibetan Plateau. Ecol Indic. 2014, 45: 184-194.

[12] de Mattos EM, Binkley D, Campoe OC, Alvares CA, Stape JL. Variation in canopy structure, leaf area, light interception and light use efficiency among Eucalyptus clones. For Ecol Manage. 2020, 463. 118038

[13] Dray S, Legendre P, Peres-Neto PR. Spatial modelling: A comprehensive framework for principal coordinate analysis of neighbour matrices (PCNM). Ecol Model. 2006, 196(3–4): 483-493.

[14] Dutilleul P, Legendre P. Spatial heterogeneity against heteroscedasticity: An ecological paradigm versus a statistical concept. Oikos. 1993, 66(1): 152.

[15] Ewers BE, Pendall E. Spatial patterns in leaf area and plant functional type cover across chronosequences of sagebrush ecosystems. Plant Ecol. 2007, 194167-83.

[16] Fitzjarrald DR, Acevedo OC, Moore KE. Climatic consequences of leaf presence in the Eastern United States. J Climate. 2001, 14(4): 598-614.

[17] Gazol A, Ibáñez R. Plant species composition in a temperate forest: Multi-scale patterns and determinants. Acta Oecol. 2010, 36(6): 634-644.

[18] Gumede N, Mutanga O, Sibanda M. Mapping leaf area index of the Yellowwood tree species in an Afromontane mistbelt forest of Southern Africa using topographic variables. Remote Sens Appl Soc Environ. 2022, 27. 100778

[19] Hardwick SR, Toumi R, Pfeifer M, Turner EC, Nilus R, Ewers RM. The relationship between leaf area index and microclimate in tropical forest and oil palm plantation: Forest disturbance drives changes in microclimate. Agric for Meteor. 2015, 201: 187-195.

[20] Jin P, Xu M, Yang QP, Zhang J. The influence of stand composition and season on canopy structure and understory light environment in different subtropical montane Pinus massoniana forests. PeerJ. 2024, 12. e17067

[21] Jones MM, Tuomisto H, Borcard D, Legendre P, Clark DB, Olivas PC. Explaining variation in tropical plant community composition: Influence of environmental and spatial data quality. Oecologia. 2008, 155(3): 593-604.

[22] Joshi RK, Garkoti SC. Litter dynamics, leaf area index and forest floor respiration as indicators for understanding the role of Nepalese alder in white oak forests in central Himalaya, India. Ecol Indic. 2020, 111. 106065

[23] Konôpka B, Pajtík J, Šebeň V, Surový P, Merganičová K. Biomass allocation into woody parts and foliage in young common aspen (Populus tremula L.)—trees and a stand-level study in the western carpathians. Forests. 2020, 11(4. 464

[24] Lai JS, Zou Y, Zhang S, Zhang XG, Mao LF. Glmm.hp: an R package for computing individual effect of predictors in generalized linear mixed models. J Plant Ecol. 2022, 156): 1302-1307.

[25] Le Bagousse-Pinguet Y, Gross N, Maestre FT, Maire V, de Bello F, Fonseca CR, Kattge J, Valencia E, Leps J, Liancourt P. Testing the environmental filtering concept in global drylands. J Ecol. 2017, 105(4): 1058-1069.

[26] Leblanc SG, Chen JM. A practical scheme for correcting multiple scattering effects on optical LAI measurements. Agric for Meteor. 2001, 110(2): 125-139.

[27] Leblanc SG, Chen JM, Fernandes R, Deering DW, Conley A. Methodology comparison for canopy structure parameters extraction from digital hemispherical photography in boreal forests. Agric for Meteor. 2005, 129(3–4): 187-207.

[28] Legendre P, Mi XC, Ren HB, Ma KP, Yu MJ, Sun IF, He FL. Partitioning beta diversity in a subtropical broad-leaved forest of China. Ecology. 2009, 90(3): 663-674.

[29] Legendre P, Legendre L (2012) Numerical ecology: Elsevier.

[30] Lewis RJ, Pakeman RJ, Marrs RH. Identifying the multi-scale spatial structure of plant community determinants of an important national resource. J Veg Sci. 2014, 25(1): 184-197.

[31] Li H, Reynolds JF. On definition and quantification of heterogeneity. Oikos. 1995, 73(2): 280.

[32] Li HB, Wang ZQ, Wang QC. Theory and methodology of spatial heterogeneity quantification. Chin J Appl Ecol. 1998, 96651-657(in Chinese)

[33] Li X, Ramos Aguila LC, Luo JC, Liu Y, Wu T, Lie ZY, Liu XJ, Cheng Y, Jiang F, Liu JX. Carbon storage capacity of Castanopsis hystrix plantations at different stand–ages in South China. Sci Total Environ. 2023, 894. 164974

[34] Liu ZL, Jiang F, Zhu Y, Li FR, Jin GZ. Spatial heterogeneity of leaf area index in a temperate old-growth forest: Spatial autocorrelation dominates over biotic and abiotic factors. Sci Total Environ. 2018, 634: 287-295.

[35] Liu ZL, Li B, Jin GZ. Scale-dependent changes in the contributions of biotic and abiotic factors to leaf area index in a natural forest in Northeast China. For Ecol Manag. 2021, 479. 118540

[36] Lysenko SA. Forward and backward linkages between land surface temperature and leaf area index for the summer in Belarus. Izv Atmos Ocean Phys. 2023, 59(9): 1137-1149.

[37] Ma JZ, Niklas KJ, Liu LY, Fang ZD, Li YR, Shi PJ. Tree size influences leaf shape but does not affect the proportional relationship between leaf area and the product of length and width. Front Plant Sci. 2022, 13. 850203

[38] Madhumali C, Wahala S, Sanjeewani N, Samarasinghe D, De Costa J. Influence of climate, soil and vegetation diversity on the leaf area index as measured by canopy hemispherical photography in tropical rainforests across a wide elevational gradient. Geomatica. 2024, 76(1. 100001

[39] Mafi-Gholami D, Zenner EK, Jaafari A, Ward RD. Modeling multi-decadal mangrove leaf area index in response to drought along the semi-arid southern coasts of Iran. Sci Total Environ. 2019, 656: 1326-1336.

[40] Majasalmi T, Stenberg P, Rautiainen M. Comparison of ground and satellite-based methods for estimating stand-level fPAR in a boreal forest. Agric for Meteorol. 2017, 232: 422-432.

[41] Malhado ACM, Costa MH, de Lima FZ, Portilho KC, Figueiredo DN. Seasonal leaf dynamics in an Amazonian tropical forest. For Ecol Manag. 2009, 2587): 1161-1165.

[42] Marenco RA, Antezana-Vera SA, Dias DP, Cândido LA. Litter production and leaf nutrient concentration and remobilization in response to climate seasonality in the central Amazon. J for Res. 2024, 35. 54

[43] Medhurst JL, Beadle CL. Crown structure and leaf area index development in thinned and unthinned Eucalyptus nitens plantations. Tree Physiol. 2001, 21(12–13): 989-999.

[44] Méndez-Toribio M, Meave JA, Zermeño-Hernández I, Ibarra-Manríquez G. Effects of slope aspect and topographic position on environmental variables, disturbance regime and tree community attributes in a seasonal tropical dry forest. J Veg Sci. 2016, 27(6): 1094-1103.

[45] Moser-Reischl A, Franceschi E, Rahman MA, Rodrigues-Leite J, Pretzsch H, Pauleit S, Rötzer T. Spatial and temporal dynamics of the leaf area index (LAI) of selected tree species in urban environments. Urban for Urban Green. 2025, 107. 128795

[46] Naithani KJ, Baldwin DC, Gaines KP, Lin H, Eissenstat DM. Spatial distribution of tree species governs the spatio-temporal interaction of leaf area index and soil moisture across a forested landscape. PLoS ONE. 2013, 8(3. e58704

[47] Nunes FP, Castro SAB, Lemos-Filho JP, Garcia QS. Correction to: canopy dynamics in a tropical Atlantic forest mosaic in southeastern Brazil. Braz J Bot. 2020, 433): 663.

[48] Oksanen J (2013) Multivariate analysis of ecological communities in R: vegan tutorial. r package version. 1.

[49] Peng WX, Zhu YF, Song M, Du H, Song TQ, Zeng FP, Zhang F, Wang KL, Luo YQ, Lan X, Zhang JY. The spatial distribution and drivers of soil microbial richness and diversity in a karst broadleaf forest. For Ecol Manag. 2019, 449. 117241

[50] Qi YJ, Zhang GQ, Luo GL, Yang TL, Wu QC. Community-level consequences of harsh environmental constraints based on spatial patterns analysis in karst primary forest of Southwest China. For Ecol Manag. 2021, 488. 119021

[51] Ramette A, Tiedje JM, Boetius A. Impact of space, time and complex environments on microbial communities. Clin Microbiol Infect. 2009, 15: 60-62.

[52] Rocha AD, Groen TA, Skidmore AK, Darvishzadeh R, Willemen L. Machine learning using hyperspectral data inaccurately predicts plant traits under spatial dependency. Remote Sens. 2018, 1081263.

[53] Rogers C, Chen JM, Croft H, Gonsamo A, Luo XZ, Bartlett P, Staebler RM. Daily leaf area index from photosynthetically active radiation for long term records of canopy structure and leaf phenology. Agric for Meteor. 2021, 304-305. 108407

[54] Rosseel Y. lavaan: AnRPackage for structural equation modeling. J Stat Soft. 2012, 4821-36.

[55] Rossi RE, Mulla DJ, Journel AG, Franz EH. Geostatistical tools for modeling and interpreting ecological spatial dependence. Ecol Monogr. 1992, 62(2): 277-314.

[56] Sarker SK, Rashid S, Sharmin M, Haque MM, Sonet SS, Nur-Un-Nabi M. Environmental correlates of vegetation distribution in tropical Juri forest. Bangladesh Trop Ecol. 2014, 55(2): 177-193

[57] Sarma D, Burman PKD, Chakraborty S, Gogoi N, Bora A, Metya A, Datye A, Murkute C, Karipot A. Quantifying the net ecosystem exchange at a semi-deciduous forest in Northeast India from intra-seasonal to the seasonal time scale. Agric for Meteor. 2022, 314. 108786

[58] Schleppi P, Thimonier A, Walthert L. Estimating leaf area index of mature temperate forests using regressions on site and vegetation data. For Ecol Manag. 2011, 261(3): 601-610.

[59] Shang X, He ZQ, Chen WB, He L, Yang H. Changes and response mechanisms of leaf area index and evapotranspiration in the typical natural landscapes of the Loess Plateau in northern Shaanxi of China under the human intervention. Ecol Indic. 2023, 154. 110517

[60] Spadavecchia L, Williams M, Bell R, Stoy PC, Huntley B, Van Wijk MT. Topographic controls on the leaf area index and plant functional type of a tundra ecosystem. J Ecol. 2008, 96(6): 1238-1251.

[61] Tanioka Y, Ida H, Hirota M. Relationship between canopy structure and community structure of the understory trees in a beech forest in Japan. Forests. 2022, 13(4. 494

[62] Vanhala P, Karhu K, Tuomi M, Björklöf K, Fritze H, Hyvärinen H, Liski J. Transplantation of organic surface horizons of boreal soils into warmer regions alters microbiology but not the temperature sensitivity of decomposition: warming, soil biology and decomposition. Glob Change Biol. 2011, 171): 538-550.

[63] Wang HQ, Sun BP, Yu XX, Xin ZB, Jia GD. The driver-pattern-effect connection of vegetation dynamics in the transition area between semi-arid and semi-humid northern China. Catena<. 2020, 194. 104713

[64] Warren CR, Adams MA, Chen ZL. Is photosynthesis related to concentrations of nitrogen and Rubisco in leaves of Australian native plants?. Aust J Plant Physiol. 2000, 27(5): 407-416.

[65] Xu YZ, Franklin SB, Wang QG, Shi Z, Luo YQ, Lu ZJ, Zhang JX, Qiao XJ, Jiang MX. Topographic and biotic factors determine forest biomass spatial distribution in a subtropical mountain moist forest. For Ecol Manag. 2015, 357: 95-103.

[66] Yuan ZQ, Gazol A, Wang XG, Lin F, Ye J, Bai XJ, Li BH, Hao ZQ. Scale specific determinants of tree diversity in an old growth temperate forest in China. Basic Appl Ecol. 2011, 12(6): 488-495.

[67] Zhou SL, Zhang WC, Wang SH, Zhang B, Xu Q. Spatial–temporal vegetation dynamics and their relationships with climatic, anthropogenic, and hydrological factors in the Amur River Basin. Remote Sens. 2021, 13(4. 684

[68] Zhou GH, Long FY, Zu L, Jarvie S, Peng Y, Zang LP, Chen DM, Zhang GQ, Sui MZ, He YJ, Liu QF. Stand spatial structure and microbial diversity are key drivers of soil multifunctionality during secondary succession in degraded karst forests. Sci Total Environ. 2024, 937. 173504

[69] Zhu WJ, Xiang WH, Pan Q, Zeng YL, Ouyang S, Lei PF, Deng XW, Fang X, Peng CH. Spatial and seasonal variations of leaf area index (LAI) in subtropicalsecondary forests related to floristic composition and stand characters. Biogeosciences. 2016, 1312): 3819-3831.

RIGHTS & PERMISSIONS

Northeast Forestry University