Driving factors of tree biomass and soil carbon pool in xerophytic forests of northeastern Argentina

Background The conversion of forests into agricultural lands can be a threat because the forests carbon stored could be a source of emissions. The capacity to improve the predictions on the consequences of land use change depends on the identification of factors that influence carbon pools. We investigated the key driving factors of tree biomass and soil carbon pools in xerophytic forests in northeastern Argentina. Based on analyses of forest structure variables and abiotic factors (topography and soil properties) from 18 mature forests, we evaluated carbon pools using uni-and multivariate (redundancy analysis) methods. Results The total carbon pool was estimated at 102.4 ± 24.0 Mg ha −1 . Soil organic carbon storage is the single largest carbon pool relative to tree biomass, representing 73.1% of total carbon. Tree canopy cover and basal area were positively correlated with biomass carbon pool ( r = 0.77 and r = 0.73, p < 0.001, respectively), proving to be significant drivers of carbon storage in this compartment. Slope, soil clay content and cation-exchange capacity had a better explanation for the variability in soil carbon pools, and all showed significant positive correlations with soil carbon pools ( r = 0.64, 0.60 and 0.50; p < 0.05 , respectively). The vertisols showed a 27.8% higher soil carbon stock than alfisols. Conclusions The relevance of our study stems from a dearth of information on carbon pools and their drivers in xerophytic forests, and in particular, the importance of this ecosystems’ type for Argentina, because they cover 81.9% of native forest area. Basal area and tree canopy cover exert a strong effect on the carbon pool in tree biomass but not in the soil. The results suggests that there is a potentially major SOC accumulation in forests located in slightly sloping areas and soils with higher topsoil clay content, such as vertisols. This could provide an important reference for implementing forestry carbon sink projects.


Introduction
In the context of weather patterns changes due to the increase in greenhouse gases (GHGs) concentration in the atmosphere, the sustainable management and conservation of forests are key to ensure the provision of the service they provide as carbon (C) sinks (Salas-Macías et al. 2020).Carbon sequestered in forest equate to ~ 45% of terrestrial carbon stocks (Rodriguez-Veiga et al. 2017).However, warnings are currently being raised about the reduction of capacity forests to supply services of carbon capture and storage due to several anthropogenic disturbances such as deforestation associated to land use change, fires, selective logging and unplanned grazing management that lead to forest degradation (Baccini et al. 2017;Ordway and Asner 2020;Hanan et al. 2021).
The soil and tree biomass are the largest carbon reservoirs of forests (González Roglich et al. 2014;Motharfard et al. 2019).Soil organic carbon (SOC) stocks exceed those in plant biomass in most climatic regions and represent between 30% (tropical forests) and 70% (boreal forests) of total stock (Pan et al. 2011).Despite, there currently remains uncertainty on the SOC stocks, and these have received limited attention from policies of climate change mitigation (Scharlemann et al. 2014).
At a global level, regional climate conditions (e.g., temperature and precipitation) are key drivers of soil carbon storage; however, climate variability at local scales is frequently small, so climatic control of carbon soil storage at a local level is less relevant in contrast to other factors.C input and C stabilization within the soil are most important aspects of SOC storage (Wiesmeier et al. 2019).Several studies in various climatic zones indicated that vegetation type affects SOC stock by controlling both the input and decomposition of C.Even if, significantly different SOC stocks were found between plant functional types due to different C allocation patterns (Jobbágy and Jackson 2000), in forests there is no consensus on the effect of tree species on SOC storage.For example, significant differences in C stocks between conifer and broadleaf forests were reported by Schulp et al. (2008a, b) while some studies indicate that species-induced differences in C contents are easily and relatively soon seen in forest floors but much less so in the mineral soil (Dijkstra and Fitzhugh 2003;Vesterdal et al. 2008).
Soil texture is probably one of the most important factors controlling SOC stabilization; therefore, the silt and clay content represents a key indicator (Wiesmeier et al. 2019).Although several studies on the key drivers of carbon stocks can be found, the wide heterogeneity of forest ecosystems makes it difficult to understand interactions between these factors (Ordway and Asner 2020;Ahirwal et al. 2021;Saimun et al. 2021).In addition, drivers of carbon pools in xerophytic forests have been comparatively less studied than in tropical forests.The xerophytic forests cover the largest native forest area in Argentina (81.9%).Although the influence of biotic and abiotic factors on tree biomass and soil C pools is widely recognized, only a few studies have been limited to specific xerophytic ecosystems, as Caldenal, and have analyzed in detail these influences (González Roglich et al. 2014); thus our understanding of drivers of C stocks in xerophytic ecosystems is still limited.Therefore, it is relevant to obtain information on carbon stocks in this type of ecosystems and the factors that influence their dynamics as input to improve the models that simulate carbon storage and to understand the consequences of land use change and forest degradation.
The research presented in this paper aimed to find out the main driving factors of tree biomass and soil carbon pool in xerophytic forests in northeastern Argentina.We defined the following specific objectives: (i) to determine the carbon pool compartmentalization of tree biomass and soil; (ii) to identify the effects of forest structure, topographic factors, and soil properties on biomass and soil carbon pools; and (iii) to assess differences in carbon pools between forests developing on different soil types.In this study, we hypothesize that: (i) forest structure has a stronger effect on carbon pool in tree biomass than in soil, and (ii) forests developing on diverse soil types have different carbon stocks both in biomass and soil.

Study area
The study was conducted in the Feliciano stream basin (Entre Ríos Province, northeastern Argentina).The basin is located between 30° 15ʹ to 31° 06ʹ SL and 58° 20ʹ to 59° 36ʹ WL.The total area of the basin (Fig. 1) is covered with about 70% (589,800 ha) of forests.The study area corresponds to the Espinal phyto-geographical province (Cabrera 1976).The Espinal woodlands have undergone rapid deforestation for the advance of the agricultural frontier, due to a combination of favorable conditions for agribusiness oriented toward global markets (Bouza et al. 2015;Cherlet et al. 2018).As well, the remaining forests are currently showing signs of degradation due to anthropogenic management for economic purposes (selective logging and unplanned grazing management).
The climate is subtropical and humid (1000 mm year −1 ); the mean annual temperature is 18.5 °C and precipitation is concentrated in summer (37.2%) and spring (32.5%), with a dry season during the winter months (11.2%).The soils at the study area are mainly vertisols and alfisols (Plan Mapa de Suelos 1990;Soil Survey Staff 2014).These soils are a specific soil type characterized by the dominant presence of smectite clays, displaying a strong shrinking-swelling capacity under repeated dry-wet seasonal cycles (Schaetzl and Anderson 2005;Wilson et al. 2013).Vertisols and alfisols developed from same parent material (Hernandarias Formation, e.g., calcareous silt material, high clay contents and carbonates presence) but they occupy different positions of the landscape.Alfisols developed mainly in higher and flat landscape positions and slightly undulating areas, with risk of waterlogging; the surface horizon has a depth of approximately 6 to 10 cm, with a low level of organic matter and clay due to leaching.The vertisols are mostly located in slightly sloping terrain (0.5-2.5%) and they are often darker, which are generally attributed to high biological activity and high organic matter content; these soils have clay contents from 35 to 50% (Plan Mapa de Suelos 1990).

Sampling sites
A re-cognisance study was conducted in the xerophytic forest ecosystem of Espinal to gather baseline information and identify possible sampling sites.Based on the forest structure, we selected 18 mature forests corresponding to "ñandubay forests" according to the First National Inventory of Native Forests (SAyDS 2007)."Ñandubay forest" is xerophytic forest representative of Espinal where P. affinis or V. caven presents a basal area higher than 2 m 2 ha −1 and conforms a tree stratum of variable density.In each forest, a 1000m 2 circular shape plot was located randomly.The minimum distance between sites was 2 km.At each site, we collected data on 11 variables: three were gathered during fieldwork (elevation, slope, tree canopy cover), three were determined through post-fieldwork processing (tree density, basal area, relative dominance of P. nigra), and finally, the other variables (soil type, clay content, soil pH, soil C:N ratio and cation-exchange capacity-CEC) were obtained from Soil Charts of Entre Rios Province, Argentina (Plan Mapa de Suelos 1990).

Field data collection and measurements
Geographic coordinates and elevation data were recorded at the center of each plot with Global Positioning System (GPS).Slopes (in percentage) were measured using a clinometer, the Suunto PM-5.At each plot, all trees with a basal diameter (0.3 m aboveground) ≥ 10 cm were identified at species level and measured for diameter and height.The trunk diameter was determined with a Mantax caliper and the height with a telescopic rod.In C. ehrenbergiana individuals, the crown area was also determined to estimate the biomass from this variable with allometric function from Conti et al. (2013).For the estimation of tree canopy cover, five rectangular shape plots of 20 m 2 were located in the middle of the 1000-m 2 circular shape plot and at their N, S, E, and W extremes to reduce the heterogeneity.In each rectangular shape plot, the canopy cover was estimated through ocular estimation (Bunnell and Vales 1990), as the proportion of the forest floor covered by the vertical projection of the tree crowns (Jennings et al. 1999).Relative dominance (RD) was calculated as the proportion of the total basal area of each species.
The average tree density, total basal area, tree canopy cover, and the relative dominance of species within the sampling plots are presented in Table 1.The total number of stems (≥ 10 cm basal diameter) inventoried per plot fluctuated from 130 per hectare to 590 trees per hectare.In the forest sites evaluated, only three species (P.nigra, P. affinis and V. caven) accounted for 94 ± 8.3% of the total basal area, while C. ehrenbergiana, A. quebracho-blanco, P. alba, and M. cisplatensis had relative dominance values of less than 5%.
To determine soil organic carbon (SOC) content, ten soil samples were collected from each extreme and center of the circular plots.Samples were taken at two depths in topsoil (0-10 cm and 10-30 cm) according to IPCC Guidelines (2006).A soil auger was used.The samples for each depth were mixed to obtain a composite soil sample for each plot.The samples were sent to the soil laboratory of the National University of Entre Ríos for further analyses.In addition to changes in the carbon concentration of soil, changes in soil bulk density (BD) should also be considered.The litter is controlled mostly by physical processes like swelling and shrinking (Wendt and Hauser 2013).In vertisol soils of the study area, Wilson et al. (2013) determined volumetric changes greater than 19% as a function of soil moisture content variation.Therefore, it is necessary to define an equivalent mass of soil from a standardized BD value (SBD) to calculate the carbon pool on an area basis (Mg ha −1 ).For this purpose, soil samples were collected to quantify the soil moisture content (%) and BD (Mg m −3 ).Three samples at each depth (0-10 cm and 10-30 cm) from different points of the circular plot using a cylinder were taken.

Soil carbon pool determination
Soil moisture content was measured through the ovendry method at 105 °C (O'Kelly 2004), and BD was determined by the cylinder method (Forsythe 1975).
The available water range is the difference between soil water content at field capacity (FC) and soil water near the permanent wilting point (WP).In the study area, soils with a dominant presence of smectite clays (vertisols and alfisols) showed a linear relationship between soil moisture content and BD within the available water range (Wilson and Cerana 2004).Therefore, a linear equation for each soil layer was generated (Eqs. 1 and 2) from the BD and soil moisture content data obtained in this work.The available water range in the 0-10 cm soil layer was between 0.20 g g −1 (WP) and 0.35 g g −1 (FC) (Wilson et al. 2013), while in the 10-30 cm layer, it was between 0.24 g g −1 (WP) and 0.40 g g −1 (FC) (Wilson, pers. comm).
where SBD: soil standardized bulk density; x: soil water content (%).We proposed to consider as SBD the BD value with soil water content corresponding to 50% of available water.Using Eqs. 1 and 2, SBD values determined were 1.17 Mg m −3 and 1.24 Mg m −3 for 0-10 cm and 10-30 cm soil layers, respectively.From these SBD values, the equivalent soil mass in the 0-10 cm layer was estimated at 1166.80 Mg ha −1 whereas, in the 10-30 cm layer, the equivalent soil mass was 2484.20 Mg ha −1 .The procedure described requires no bulk density sampling and thus avoids imprecision in its determination that occurs through volumetric changes under repeated drywet seasonal cycles.
For SOC determination, samples were oven-dried at temperatures below 40 °C and sieved (2 mm).SOC ( 1) (2) SBD (10-30) = −0.005601x+ 1.416 R 2 = 0.30; p < 0.007 content was determined by Walkley and Black method (1934) and was expressed as mass of carbon per unit mass of soil (g kg −1 ) or percentage.For the calculation of the soil carbon pool on an area basis (Mg ha −1 ), SBD and soil depth data were used, as follows (Eq.3): where SBD: soil standardized bulk density; SOC: soil organic carbon.
Soil carbon pool at 30 cm was obtained from the sum of the carbon pool at two depths in topsoil (0-10 cm and 10-30 cm).

Tree aboveground biomass (AGB) and carbon pool estimation
To estimate tree AGB, we used functions developed for xerophytic forests of Espinal and Arid Chaco regions (Table 2).For species without models (P.alba and M. cisplatensis), we applied the general Eq.4: where AGB: aboveground biomass; bd: basal diameter trunk; h: tree high; wd: wood density.We used the wood density data from INTI-CITEMA (2007).
We used an average value of C fraction (0.475) for the Espinal dominant species (Sione et al. 2021).Finally, at each site, the total C pool was calculated from the sum of tree biomass and soil C pools.

Statistical analysis
Redundancy analysis (RDA) was developed to assess how much of the variation in one set of variables (response (3) variables: tree biomass and soil C pools) can be explained by the variation in another set of variables (explanatory variables: tree density, basal area, tree canopy cover, relative dominance of P. nigra, elevation, slope, soil clay content, soil CEC, soil pH, and soil C:N ratio).RDA, a direct extension of multiple regression in that it models the effect of an explanatory matrix X on a response matrix Y, was carried out using R software (vegan package).Pearson's correlation was calculated to assess the relationships between dependent variables (tree biomass carbon, soil carbon, and total carbon stock) and independent variables (forest structure data and abiotic factors).To show the variability of tree biomass and soil carbon pools for each soil type, the coefficient of variation (CV) was calculated (Guerra Dávila 2014).We performed oneway analysis of variance (ANOVA) of tree biomass and soil pools for different types of soils.For these ANOVAs, we employed Tukey's test at p < 0.05 to separate means.All statistical analyses were performed using the RStudio software 2023.06.0.

Tree biomass and soil carbon pools
The total carbon pool in the xerophytic forests in study sites was estimated at 102.4 ± 24.0 Mg ha −1 .The carbon stored in the first 30 cm depth of soil accounts for 73.1 ± 6.6% of total carbon per hectare.Tree biomass carbon pool per unit area varied widely, ranging from a low of 11.0 Mg ha −1 to a high of 44.7 Mg ha −1 .The soil carbon pool showed a significant positive correlation (R 2 = 0.25; p < 0.05) with tree biomass carbon (Fig. 2).

Comparison of the carbon pool among forests developing on different soil types
There were large differences in carbon pools among soil orders (Fig. 3 and Additional file 1: Table S1).Vertisols had the greatest soil carbon pool (86.1 Mg ha −1 ).
A differential behavior was observed in the variability of the soil carbon pool, resulting in a higher CV in vertisols (24.5%) than in alfisols (16.8%).The tree biomass carbon pool in vertisols also resulted in higher levels than in alfisols (Fig. 3); however, these differences were not significant (p > 0.05).The variation coefficients in the mentioned carbon pool were estimated at 39.1% and 30.8% for vertisols and alfisols, respectively.

Driving factors of carbon pools (tree biomass and soil pools)
Both forest structure variables and abiotic factors showed significant associations with carbon pools in the different compartments (tree biomass and soil pool) at a plot level (R 2 = 0.86; adj.-R 2 = 0.66; F = 4.30; p = 0.015).Redundancy analysis (RDA) showed that basal area, canopy cover, relative dominance of P. nigra, tree density, slope, soil clay content, soil CEC, elevation, soil pH, and soil C:N ratio were important environmental factors influencing tree biomass and soil carbon pool.The results of the RDA test indicated that the RDA biaxial explained 86.0% of the variation in carbon pools.Axis 1 of the RDA plot explained roughly 73.9% of the variation, while axis 2 explained a further 12.0% (Table 3 and Fig. 4).Among the environmental factors affecting soil carbon pools, slope, clay content and soil CEC had a better explanation for the variability in soil carbon pools.Among them, slope explained the highest degree, and all showed significant (p < 0.05) positive correlations with soil carbon pools (Figs. 4 and 5).Tree biomass carbon pool was found highly significantly positively (p < 0.001) correlated with cover canopy and basal area (Figs. 4 and 5).

Discussion
The soil and vegetation have complex interrelations because they develop together over a long period of time (Ahirwal et al. 2021).Soil organic carbon storage is controlled by the balance of carbon inputs from vegetation and outputs through decomposition; so amount organic carbon in the soil affects and is affected by vegetation (Schlesinger 1977).We found that xerophytic forests of northeastern Argentina contain a total carbon pool (tree biomass + soil for the first 30 cm depth) of 102.4 Mg ha −1 .These estimates are within the range of earlier observations of the comparable ecosystems considering aboveground biomass and upper soil profiles (Fig. 6).
The estimated carbon pool in our work is within the range of estimates for seasonally dry Chaco forests, including all aboveground compartments and organic soil profiles up to 1 m depth (Gasparri et al. 2008), and it is similar to the closed semiarid Caldenal forests, considering woody and herbaceous, above-and belowground litter, and deeper organic soil profiles up to 2 m depth (González Roglich et al. 2014).Other work, including carbon in aboveground plant standing biomass and organic soil profiles up to 2 m depth, resulted in 121.3 Mg C ha −1 of total carbon storage (Conti et al. 2014), similar to the ones determined in our study.
According to previous studies, we found that carbon stored in soils to be the single largest carbon pool (74.7 Mg ha −1 ), representing 73.1% of total carbon per hectare, with more than half the sites falling between 72.1 and 85.7% of total pool.In comparable forests (Fig. 6), such as the Caldenal in semiarid Espinal, and semiarid Chaco forests, between 53 and 63% of total ecosystem carbon stock was in the soil, respectively (Bonino 2006;  Studies carried out in different forest regions of Argentina reported a wide range of tree biomass carbon pools among forest types.For instance, the xerophytic forest of Espinal region stored between 38.9 and 43.9 Mg ha −1 of carbon (Sione et al. 2021 andLedesma et al. 2021, respectively); the semiarid Chaco forests sequestered up to 43.0 Mg ha −1 (Gasparri and Manghi 2004;Conti et al. 2014;Loto 2021), while a carbon stock of 270.1 Mg ha −1 was determined in the aboveground biomass of Andean-Patagonian forests (Gasparri and Manghi 2004).The average carbon value estimated in our work was 27.8 Mg ha −1 , with half the sites being between 20 and 35 Mg ha −1 .These values showed moderate variability (34%).The remaining forests of the study area are heterogeneous in coverage and degradation level due to several anthropogenic disturbances such as selective logging, fires, and unplanned grazing management.Previous studies have found high variability both in biomass carbon pools and in soil organic carbon pools of degraded forests.In the case of Chaco forests, for example, the forest was degraded, its biomass and carbon decreased significantly because of the reduction in density and size of individuals; as a result, total biomass exhibited significant variations among the communities.However, a moderate decrease in soil carbon content with respect to vegetation degradation was observed (Bonino 2006).
Biomass carbon storage depends on many factors, such as the type, density, and diversity of the vegetation.On the one hand, a strong positive correlation between carbon stock in trees was observed, and on the other tree density, size, and height (Aryal et al. 2018;Dimobe et al. 2019).In relation to previous studies (González Roglich et al. 2014;Sheikh et al. 2020;Osei et al. 2022), we found that tree canopy cover and basal area proved to be significant drivers of biomass carbon storage, as well as the total carbon pool stocks of forests.Instead, we did not find a relationship between tree density and these pools.The size of tree species and canopy cover explained more of the variation observed in biomass carbon stocks than in tree density (Gebeyehu et al. 2019).Burrows et al. (2000) indicate that basal area is usually used to estimate biomass carbon pools since it shows a significant relationship with this reservoir, integrates both the number and size of trees, and is easily measurable.Osei et al. (2022) found that basal area explained > 88.0% variability in aboveground woody carbon stocks.The significant importance of tree canopy cover on biomass carbon pool could be attributed to the physiognomy of dominant species of Espinal forests.For instance, P. nigra and P. affinis individuals invest most of their resources in the development of crowns rather than stems (Sione et al. 2019(Sione et al. , 2020a)), therefore, high tree canopy cover values are associated with higher carbon pools.
An effect of P. nigra relative dominance on biomass carbon storage (r = 0.43) was observed, although it was not significant.Earlier findings revealed that P. nigra was associated with the highest biomass and carbon stocks per unit or basal area (Sione et al. 2020a, b).In addition, at the individual tree level, this species stores more carbon per unit of biomass (48%) than other dominant species of forests evaluated (Sione et al. 2021).For instance, applying existing biomass models allows us to know that the aboveground biomass carbon storage in P. nigra individual (with basal diameter of 30 cm) was 45% higher than for P. affinis individual (same diameter), even when the latter has higher specific wood density.This implies that forests could have comparable basal area per hectare but higher biomass carbon storage in those with higher relative dominance of these species.
Neither elevation nor slope was significantly associated with carbon pools of tree biomass in our study.Existing literature is not consistent in showing the relationship between topographic factors and carbon stocks: some sources provide evidence that carbon stocks increase gradually with increasing slope and elevation (Gebeyehu et al. 2019;Pragasan 2020;Saimun et al. 2021), others indicate that stocks are reduced (Bohara et al. 2021) or no relation is shown (Salas Macías et al. 2020).Our results showed that slope, soil clay content, and CEC are the principal factors controlling soil carbon pool, as well as total carbon pool stocks of forests, although the correlation between CEC and total pool (r = 0.46) was not significant.The strong influence of clay content on soil carbon storage is consistent with earlier observations in this same region (González Roglich et al. 2014;SAGyP 2022), reinforcing the previous findings in different ecosystems worldwide.Several authors reported that soil texture plays an important role, as clay content increases, carbon outputs decrease because clay has a stabilizing and protecting impact on soil organic carbon (Paul 1984;Jobbágy and Jackson 2000).A high clay content may cause to clay surfaces adsorb more organic C molecules due to the larger surface area and the presence of polyvalent cation forming organic-mineral complexes to control the protection of SOC from microbial and enzymatic decomposition, increasing SOC storage (Zaffar and Lu 2015;Zhong et al. 2018).It is commonly accepted that SOC binding potential is higher in soils dominated by 2:1 clay minerals (such as montmorillonite, clay type dominant in the soils of study area) than in soils dominated by 1:1 clay minerals.The higher cation-exchange capacity and larger specific surface of 2:1 clay minerals, which would promote organic-mineral complexion and their capacity of these minerals to promote aggregation could explain the higher SOC stabilizing capacity of soils dominated by this type clays (Greenland 1965;Kleber et al. 2007;Fernández-Ugalde et al. 2013).Besides, we did not find a significant relationship between elevation and soil carbon stock.A review of drivers of soil organic carbon at various scales (Wiesmeier et al. 2019) indicates that there is no general agreement on the importance of topographical features for soil carbon storage, due to various scales of observation.Local terrain attributes such as slope and curvature appear to influence soil carbon storage only at small spatial scales, while landscape position may become an important key at larger scales.These authors indicated that the importance of topographical features for soil carbon storage is due to their contribution to erosional processes because they control water flow paths, water accumulation, and discharge.Sites with steep slopes and convex curvatures lead to higher water discharge, whereas low inclination and concave curvatures favor water accumulation.Fritsch and Fizpatrick (1994) propose a link between waterlogging and organic matter and clay migration down to the lower horizons.Therefore, SOC and clay loss occurs in surface horizons, which leads to poor soil aggregate stability (Indelángelo et al. 2007;Gabioud et al. 2011).We found that forests developing on diverse soil types have different carbon stocks in soil but not in biomass, partially reinforcing our second hypothesis.The vertisols showed 27.8% higher soil carbon stock than that in alfisols.The observed differences are probably explained by differences in slope, clay content, and CEC.These two soil types, developed from the same parent material, occupy different positions in the landscape.Alfisols develop mainly in higher and flat landscape positions and slightly undulating areas, with a risk of waterlogging while vertisols are mostly located in slightly sloping terrain.However, as was reported by Wiesmeier et al. (2019), different soil carbon pools among soil types could be attributed to the different slopes rather than to the position in the landscape.In our work, 61% of the plots were located on alfisols and 39% on vertisols.All plots located on alfisols were located on slopes less than 1%, whereas around 64% were on areas with higher elevations than 65 m.In vertisols, the elevation of the plots was lower than 65 m and the slope varied between 1 and 2%.In addition, vertisols showed average clay contents of 327.6 g kg −1 , which are 16% higher than those determined in alfisols (Additional file 1: Table S1).These outputs indicate that soil type is not independent controlling driving factor but integrates a number of parent materials and topographical situation-related variables, such as soil texture and moisture contents, which directly affect the potential of soils to store C (Wiesmeier et al. 2019).
Previous literature showed controversial results on the importance of slope as a key driver of soil carbon storage.In the semiarid savannahs of Argentina, no statistically significant relationship between slope and soil carbon pool was found (Gonzalez Roglich et al. 2014).However, some studies determined that a significant proportion of soil carbon storage variability was explained by slope gradient (Chaplot et al. 2010) while other studies found that sites on steeper inclination have been found to contain less SOC (Ritchie et al. 2007).Even when our work showed that forest structure variables have no significant effect on soil carbon storage, tree canopy cover was the variable of higher correlation (r = 0.44), although it was not significant.Consistent with our results, earlier studies found that stand density (expressed as basal area per hectare) had no significant effect on SOC storage in the whole soil profile (Osei et al. 2022).
There is no consensus on the effect of tree species on soil carbon storage in forests (Wiesmeier et al. 2019).For example, tree nitrogen-fixing species have larger effects on forest soils than other species, which include consistent increases in soil organic carbon; this could be attribute to their low C:N ratios, its therefore their relatively fast decomposition rates (Mendoza et al. 2014).On the other hand, soil storage in coniferous forests may be relatively high due to an accumulation of acidic litter in the organic layer (Schulp et al. 2008a, b).However, other studies did not showed significant differences between soil carbon pools among forest types (Vesterdal et al. 2008).Both soil pH and soil C:N ratio were negatively associated with carbon pools in tree biomass and soil, although these correlations were not statistically significant.In dry forests, Gebeyehu et al. (2019) found that moderate levels of soil pH had a positive effect on soil carbon storage, but extreme levels (either low or high) of soil pH caused a decrease in this stock.We found that soil pH values ranged from 5.2 to 7.2, with 89% of sites falling below 6.3 (Additional file 1: Table S1).Soils with these pH values are considered moderately to slightly acid, characteristic of humid regions.
Early findings revealed that soils with C:N ratio in the range of 10 to 12 contribute to slow organic matter decomposition as they require an excess of microbial activities (Landon 1991;Gebeyehu et al. 2019).The 72.2% of sites evaluated in our work corresponded to this range, suggesting slow organic matter decomposition.The main limitation of our study is that we have only quantified carbon storage in aboveground tree biomass and soil in first 30 cm depth.Its therefore, data obtained are likely to be underestimates of carbon stocks in forest of region, carbon emissions from land use change and their CO 2 mitigation potential, as we have not included carbon storage in shrubs, non-woody vegetation and coarse woody debris, belowground biomass, and the carbon stock in the deepest soil layers.

Conclusions
This study contributes to understanding of the key driving factors of carbon pools in xerophytic forests in northeastern Argentina.Its relevance stems from a dearth of information on carbon stocks in this forest type relative to tropical forests, and the importance of xerophytic forests for Argentina, because they cover the largest native forest area (81.9%).The results suggests that there is a potentially major accumulation of SOC in forests located in slightly sloping areas and soils with higher topsoil clay content.These characteristics determinate to the fact that SOC storage in forests developing on vertisols were higher compared to alfisols.
We found that basal area and tree canopy cover exert a strong effect on the carbon pool in tree biomass but not in the soil, being consistent with our first hypothesis.We also showed that canopy cover is a driving factor of total carbon pool.This suggests that actions to manage tree species can affect the variability of carbon stocks.Management practices that maintain forest cover are essential to conserving carbon in tree biomass.Furthermore, the high association between canopy cover and carbon stock in tree biomass raises challenges for the application of remote sensing to the quantification of carbon stored in tree biomass due to canopy cover being one of the variables that can most easily be measured using this tool.Soil organic carbon stock is the single largest carbon pool in xerophytic forests in northeastern Argentina, accounting for more than 70% of the total pool.These findings reinforce the need to provide certainty on estimates both on soil carbon stocks and emissions from soils, as well as to focus on this carbon pool in climate change mitigation policies.The slope and clay content could provide an important reference for implementing forestry that analyze carbon sink projects.In soils that have higher carbon storage potential, such as vertisols, forests can offer major opportunities for increased carbon stock through sustainable management practices.The information presented here may be especially useful to improve carbon storage models and predictions on consequences of climate change, land use change, and forest degradation, as well as for regional policy development on emissions reduction by deforestation and forest degradation (e.g., REDD+).

Note
A disintegration of the genus Prosopis L. has been recently proposed by Hughes et al. (2022) based on molecular phylogeny, and new nomenclatural combination for American species has been published under genus Neltuma.Nevertheless, the classification of Prosopis by Burkart's (1976) is maintained here due to nomenclatural stability, prevailing current use and disadvantageous changes (Articles 14, recommendation 14A.1;Article 34, recommendation 34A.1;Article 56 recom-mendation56A, ICBN (Turland et al. 2018), until acceptance of nomenclatural changes by the International Committee of Nomenclature.

Fig. 1
Fig. 1 Study area (Espinal phyto-geographical province, northeastern Argentina): A inbox, location of the study area in South America (modified from MAyDS 2020); B study area map (Stream Feliciano Basin, Entre Ríos Province) showing the distribution of the three main forest types (ñandubay-espinillo forests, dark green-shaded area; ñandubay-espinillo forests which included other species, green-shaded area; and ñandubay open forests, light green-shaded area), and croplands, red-shaded area; C and D typical Espinal forest landscape, where this study was conducted

Fig. 4
Fig. 4 Redundancy analysis (RDA).Explanatory variables (forest structure and abiotic factors) were indicated by arrows.Arrows pointing toward a response variable (tree biomass and soil carbon pool) indicate a high positive correlation.Arrow lengths indicate the strength of the relationship between carbon pools and explanatory variables.Forest structure data (density, basal area, canopy cover, relative dominance of P. nigra) and abiotic factors (elevation, slope, soil clay content, soil CEC, soil pH, and soil C:N ratio) explained 84.0% of the total variance in carbon pools (tree biomass and soil), with 73.9% of the variance on RDA axis 1 and 12.0% on RDA axis 2. pH soil pH, C:N soil C:N ratio, CEC cation-exchange capacity, CLC clay content, SLO slope, ELE elevation, PND relative dominance of P. nigra, TCC tree canopy cover, TD tree density and BA basal area

Fig. 5 Fig. 6
Fig. 5 Correlation coefficient (Pearson's r) of the relationship among carbon pools (tree biomass carbon, soil carbon and total carbon) and different factors.pH soil pH, C:N soil C:N ratio, CEC cation-exchange capacity, CLC clay content, SLO slope, ELE elevation, PND relative dominance of P. nigra, TCC tree canopy cover, TD tree density and BA basal area.*significant (p < 0.05), and **highly significant (p < 0.01)

Table 3
Redundancy analysis (RDA) and biplot scores for explanatory variables pH soil pH, C:N soil C:N ratio, CEC cation-exchange capacity, CLC clay content, SLO slope, ELE elevation, PND relative dominance of P. nigra, TCC tree canopy cover, TD tree density and BA = basal area