Skip to main content
Download PDF

Carbon fluxes resulting from land-use changes in the Tamaulipan thornscrub of northeastern Mexico

Carbon Balance and Management volume 3, Article number: 6 (2008) Cite this article

Abstract

Information on carbon stock and flux resulting from land-use changes in subtropical, semi-arid ecosystems are important to understand global carbon flux, yet little data is available. In the Tamaulipan thornscrub forests of northeastern Mexico, biomass components of standing vegetation were estimated from 56 quadrats (200 m2 each). Regional land-use changes and present forest cover, as well as estimates of soil organic carbon from chronosequences, were used to predict carbon stocks and fluxes in this ecosystem.

For the period of 1980–1996, the Tamaulipan thornscrub is presenting an annual deforestation rate of 2.27% indicating that approximately 600 km2 of this plant community are lost every year and that 60% of the original Mexican Tamaulipan thornscrub vegetation has been lost since the 1950's. On the other hand, intensive agriculture, including introduced grasslands increased (4,000 km2) from 32 to 42% of the total studied area, largely at the expense of the Tamaulipan thornscrub forests. Land-use changes from Tamaulipan thornscrub forest to agriculture contribute 2.2 Tg to current annual carbon emissions and standing biomass averages 0.24 ± 0.06 Tg, root biomass averages 0.17 ± 0.03 Tg, and soil organic carbon averages 1.80 ± 0.27 Tg. Land-use changes from 1950 to 2000 accounted for Carbon emissions of the order of 180.1 Tg. Projected land-use changes will likely contribute to an additional carbon flux of 98.0 Tg by the year 2100. Practices to conserve sequester, and transfer carbon stocks in semi-arid ecosystems are discussed as a means to reduce carbon flux from deforestation practices.

Background

Carbon dioxide concentrations in the atmosphere have increased by nearly 30% since the beginning of the Industrial Revolution [1–3]. The climatic and other environmental consequences of increasing carbon dioxide concentrations in the atmosphere have been widely discussed [4]. Carbon is released to the atmosphere by a variety of activities including the combustion of fossil fuels, deforestation, land-use changes, and cement production [5].

Deforestation and other land-use changes have contributed on the average 28% of the total global carbon dioxide emissions for the period between 1980 and 1998 [6]. Forest ecosystems cover more than 4.1 × 109 hectares of the Earth's land area [7] and contain about 1146 Pg C, with approximately 378 Pg in vegetation and the remaining 768 Pg in soils. Assessments of CO2 emissions from deforestation practices have been conducted at global and regional scales [e.g., [5, 7–10]], at national scales [11, 12], and at local levels [13–15]. Most studies at local or regional scales have focused on tropical regions, but little is known about carbon emissions of subtropical, semi-arid ecosystems. Masera et al., [11] suggested that conversion of forest lands to other uses can be contributing as much as 41% of the CO2 for Mexico and for the globe it is approximately 25% [5]. More realistic assessments must be conducted at local scales to improve understanding and to evaluate the feasibility of conducting practices leading to increasing carbon stocks to reduce this flux.

The aim of this study is to estimate the carbon flux associated with land-use change in the Tamaulipan thornscrub of northeastern Mexico as an example of the contribution to carbon flux of subtropical, semi-arid ecosystems. Management options for conserving stocks and sequestering carbon are discussed for the Tamaulipan thornscrub of northeastern Mexico with the possibility of extending these options to other areas.

Results

Estimated aboveground and root biomass components

Observed aboveground biomass components measured during the period of 1999–2002 are presented in Table 1. The spatial variation of biomass components was quite large, and most variation was observed between Mexican states. The semi-arid northwestern portion had average and confidence intervals (P = 95%) of biomass components of 12.93 ± 3.24 Mg ha-1, 10.17 ± 1.82 Mg ha-1, and 23.10 ± 5.06 Mg ha-1 for aboveground, root, and total standing biomass, respectively. The other two locations had 48.40 ± 7.40 Mg ha-1, 30.04 ± 4.15 Mg ha-1, and 78.44 ± 11.54 Mg ha-1 for aboveground, root, and total standing biomass, respectively. Total weighted root biomass estimates approximated 15.62 Mg ha-1 and accounted for 64% and 39% of the total aboveground and total standing biomass, respectively. That is, weighted carbon stocks in biomass components of Tamaulipan thornscrub are 60% in aboveground biomass (11.35 Mg ha-1) and 40% in root biomass (7.81 Mg ha-1).

Table 1 Carbon estimates for the period 1999–2002 in plant compartments and soils in the Tamaulipan thornscrub ecosystem of three Mexican States located within northeastern Mexico.
Full size table

Estimated aboveground and root carbon components in fallow lands

In fallow lands, carbon stocks in aboveground biomass reach a weighted average of 0.44, 2.11, and 5.25 Mg ha-1 for sites with 5, 10, and 15 years of abandonment, respectively. The mean periodic and mean annual biomass increment is: 0.61 and 0.24; 1.36 and 0.57; and 2.21 and 0.95 Mg ha-1 year-1, for sites with 5, 10, and 15 years of abandonment, respectively. That is, in 15 years of fallow, carbon stocks in aboveground biomass attain less than 50% of the carbon stocks in the original Tamaulipan thornscrub. Carbon stocks in root biomass attain a weighted average of 5.84 Mg ha-1 in sites with 15 years of abandonment.

Soil carbon stocks in the Tamaulipan thornscrub forest

The soil contains the largest C pool in the Tamaulipan thornscrub ecosystem. It has a weighted average of approximately 136.5 Mg ha-1, or 88% of the total organic carbon in the Tamaulipan thornscrub of northeastern Mexico. The rates of soil carbon accumulation during the establishment of Tamaulipan thornscrub in aboveground and root biomass and soils is approximately 5.25, 3.28, and 7.72 Mg ha-1 in 15 years of site abandonment from agricultural practices. These values translated to annual rates per m-2 approximates to 35, 22, and 51 g m-2 y-1 and they are consistent with values reported by Post and Kwon [16] for several semi-arid ecosystems.

Land-use changes in the Tamaulipan thornscrub of northeastern Mexico

The land-use changes from 1980 to 1996 observed in a 30,000 km2 block placed in the center of the distribution of the Tamaulipan thornscrub of northeastern Mexico are large (Table 2). The disappearance of Tamaulipan thornscrub forest correlates to the appearance of agriculture and grasslands in the region (Figure 1). The proportion of the total studied area (30,000 km2) covered by native Tamaulipan thornscrub forests estimated from the 1980 data was 42% (12 282 km2) and only 27% remained intact for 1996 (4 905 km2). Intensive agriculture, including induced grasslands increased from 32 to 42% of the total studied area. This relationship implies that the significant increase in area converted to mainly agricultural fields (3,000 km2) and grasslands (1,000 km2) occurred largely at the expense of the Tamaulipan thornscrub forests. Plant communities dominated by Prosopis and Acacia forests and degraded forms of scrublands accounted for by the rest of the Tamaulipan thornscrub area lost during the studied period. Areas with gentle slopes, characterized by deep dark soils, are the most affected by shifting cultivation practices and intensive agriculture. Fallow lands, on the other side, increased at the expense of Mezquitales (1 000 km2) and abandoned agriculture and grasslands (1 500 km2). The annual rate of abandoned, fallow lands approximates to 2.33%, and surprisingly, this figure is quite similar to the annual deforestation rate. That is, for the study period, the area cleared is approximately the same as the area abandoned from other land uses which may eventually be converted to Tamaulipan thornscrub forests.

Table 2 Lost area and land-cover change of several plant communities in the Tamaulipan thornscrub ecosystem of northeastern, Mexico for the period of 1980–1996.
Full size table
Figure 1
figure 1

The distribution range of plant cover in northeastern Mexico, including Tamaulipan thornscrub, agriculture and grasslands.

Full size image

The annual rate of clearing Tamaulipan thornscrub is approximately 2.37%. That is, according to the last forest inventory, this ecosystem covers an approximate area of 32,188 km2, and it is disappearing in the natural state at an annual average rate of 600 km2. This figure is consistent with land use changes reported in other studies for the region (Table 3). SARH [17] reported deforestation of 3350 km2 resulting from a combination of the opening of the irrigation districts 'Bajo Rio Bravo' (2400 km2), 'Bajo Rio San Juan' (800 km2), and 'Las Lajas' (150 km2), which occurred because of the construction of the International reservoirs 'La Amistad' and 'Falcon' and the Mexican reservoir 'El AzĂºcar' between the 1950's and the 1970's. The state government of Nuevo Leon reported deforestation figures of 1800 km2 between 1980 and 1986 [18]. Maldonado [19] found that 11% of the plains of the Gulf of Mexico of the State of Nuevo Leon were deforested between 1976 and 1986. Treviño et al., [20] reported for the municipality of Linares annual deforestation rates of 2.2% for the period of 1973 to 1994.

Table 3 Land-use changes from the 1950's to the 2000's in the Tamaulipan thornscrub of northeastern Mexico and southeastern USA.
Full size table

Considering the average annual deforestation rate (1.85%) estimated from a) this study as 2.37%, b) by Treviño et al., [20] for the period of 1973–1994 as 2.2%, and c) by Maldonado [19] for the period of 1976–1986 as 1.1%, and current area reported by Palacios-Prieto [21], the Tamaulipan thornscrub of northeastern Mexico would have covered an area of 80,490 km2 in 1950. This estimate compares to the deforestation area reported by the state of Nuevo Leon (1800 km2) for the period of 1980–1986 [18] and the sum all 6-year periods since 1950 and for the other two states (Tamaulipas and Coahuila). Considering the average rate of deforestation of 1.85%, by the year 2000 there was 60% less area covered by Tamaulipan thornscrub in northeastern Mexico. The World Wildlife Fund [22] estimated that over 25% (8 047 km2) of this ecoregion will disappear due to changes in land-use by the year 2020 since this rate will reduce the area covered by Tamaulipan thornscrub by 31% of the area reported by Palacios-Prieto [21]. Therefore, under this deforestation scenario, the area covered by this ecosystem would be 28% (22 156 km2) and 6% (4 974 km2) of the area covered in 1950 for the years 2020 and 2100, respectively.

Carbon fluxes by land-use changes

The carbon released in this ecosystem due to land-use changes was 180.1 Tg for the period between 1950–2000. Under the rate of the deforestation scenario, carbon fluxes for the period of 2000–2100 would be 98.0 Tg (Figure 2). For the year 2000, in general, clearing Tamaulipan thornscrub for agriculture is releasing carbon at an annual rate of 2.19 Tg C y-1 of which aboveground biomass and root biomass accounted for by 10 and 8% and soils explained the rest 88%. Shifting land use from agriculture to Tamaulipan thornscrub increases carbon stocks at an annual rate of 0.74 Tg C y-1 of which aboveground and root biomass explained 16 and 9% and soil the remaining 76%. Therefore, there is a net release of carbon stock of 2.19 Tg C y-1 of which soil looses approximately 1.80, aboveground biomass 0.22 and roots 0.17 Tg C y-1. The conversion of native forests and grasslands to farms has released globally about 100 Pg of carbon over the past 150 years [23]. During the 150 years scenario studied in this report (1950–2100), this ecosystem would likely release 0.278 Pg C.

Figure 2
figure 2

Land-use change and carbon flux resulting from land-use change in the Tamaulipan thornscrub of northeastern Mexico.

Full size image

Conclusion

Estimated aboveground and root biomass components

Biomass estimates are consistent with other biomass measurements and estimates in the Tamaulipan thornscrub of south central Nuevo Leon. NĂ¡var et al. [24] measured for the subtropical Tamaulipan thornscrub from 25 m2 ha plots an average aboveground biomass of 44.4 Mg ha-1. Heiseke [25] reported biomass stock figures between 34.21 and 62.70 Mg ha-1 in the plains and between 26.06 and 37.70 Mg ha-1 in the hills and concluded that this biome has a maximum standing biomass of 50 Mg ha-1. Carstens [26] measured in this plant community between 35 and 47 Mg ha-1. Cairns et al. [14] recorded slightly smaller biomass figures for chaparrales and xerophytic vegetation of southern Mexico (24 and 28 Mg ha-1, respectively). Castellanos et al. [27] reported larger biomass figures for dry tropical forests of western Mexico (73 Mg ha-1). Root biomass estimates are consistent with estimates in cloud forest ecosystems, but they are smaller than reported figures for tropical ecosystems [15, 27].

Total soil carbon stock per hectare in this ecosystem (138 Mg ha-1) is within the range of measurements and estimates reported by De Jong et al., [15] for secondary shrubs (average of 140 Mg ha-1), for pastures (120 Mg ha-1), and for cultivated lands (140 Mg ha-1) of the selva Lacandona in Chiapas Mexico. The Tamaulipan thornscrub has most of its carbon pool allocated in organic soil matter and aboveground biomass (weighted average of 136.46 Mg ha-1 and 11.35 Mg ha-1, respectively). The ratio of standing biomass carbon/soil organic carbon approximates to 7% and it is consistent with the ratio reported by Watson et al. [6] for tropical savannas. De Jong et al., [13] summarized several studies on tropical ecosystems and observed that on the average the soil and standing biomass compartments are balanced with 120 Mg ha-1. Watson et al., [6] also reported similar carbon stocks in soils and standing biomass for tropical forests.

Rates of deforestation in the Tamaulipan thornscrub ecosystem

Deforestation figures reported for the Tamaulipan thornscrub are also within the average recorded for several studies. For Mexico, for the period of 1993–2000, Semarnat [28] recorded deforestation rates of 3.93%. For the tropical rainforest of 'Los Tuxtlas' for the period of 1967–1986, Dirzo and Garcia [29] reported 4.3%. Cairns et al. [30], Cairns et al., [14], De Jong et al., [15], and Ochoa-Gaona [31] recorded similar deforestation figures for tropical forests of southern Mexico. Dirzo and García [29] also reported similar deforestation rates for tropical dry forests of western Mexico.

Land-management options to conserve carbon sequestration in the Tamaulipan thornscrub ecosystem

This study suggests that eliminating land-use changes in the Tamaulipan thornscrub would conserve carbon stocks by 2.19 Tg annually. In addition, this action would provide other environmental benefits such as habitat conservation for diversity [22], soil and water conservation [32], buffering local microclimate and global climatic changes, and other additional benefits.

Deforestation promotes soil carbon emissions of the order of 1.8 Tg yr-1 (82% of total emissions). Removals of standing aboveground and root biomass account for the remaining 18%. Conservation of SOC must be the primary concern and can be approached in several ways: (1) conserving the natural ecosystem, (2) protecting the soil carbon, and (3) sequestering soil carbon that has been lost.

The conservation of the natural ecosystem is of critical local, national, and global importance. The promotion of new policies on land change and management can reduce carbon emissions in the agricultural sector. For example, modifying the forestry law by reducing the scale of land-use change or increasing the time allotted for abandoned lands to recover the original plant community would cause less private and communal lands to be cleared. Conservation of carbon stocks under the new policy scenario is a function of the area allowed to be cleared and the time to return to forests. For example, at the present time, if fallow last 5 or 10 years, SOC sequestration would be reduced to 0.07 and 0.01 Tg C y-1, respectively. On the other hand, if fallow lasts for 25 years, SOC sequestration would be 0.52 Tg C y-1.

In recently cleared lands, conservation of the SOC must be the primary concern because organic carbon in total aboveground biomass would only be released, accounting for 10% of the total organic carbon stocks in the system. Currently, approximately 600 km2 year-1 are cleared for farming in this ecosystem and lose much soil organic carbon (65%) in 5 years [33]. Since most of the newly cleared lands are used for rain-fed agriculture, a feasible way to preserve carbon is to implement non-till, minimum tillage or tillage conservation practices at the beginning of farming. Of the area cleared each year, approximately 20 Tg C could be conserved in 10 years of no-till conservation practices. However, the expense of high fertilizer and pesticide expenditures used to control weeds may limit the ability of farmers to carry out these practices unless a carbon credit or payment to farmers is associated with the carbon sink in cropland agriculture.

Carbon can be sequestered in soils that have been under agricultural use for more than 15 years by promoting no-till or conservation practices. No-till and conservation tillage practices associated to efficient irrigation, fertilizer, and pesticide applications can increase SOC [34–38]. The increase can occur via several means: improving yields and subsequent organic matter additions to the soil or reductions in the rate of soil organic carbon loss [39], decreasing carbon dioxide emissions caused by the reduction of fossil-fuel use by tractors and other tillage equipment, and reducing the area cleared by increasing crop yield. Since the rate of carbon sequestration can be doubled in soil with depleted organic carbon [40], conservation tillage practices can be effectively carried out in the irrigation districts 'Bajo Bravo', 'Bajo San Juan', and 'Las Lajas', where agriculture is highly mechanized, irrigated, and with high inputs of fertilizers, pesticides, and herbicides. There are approximately 3500 km2 available for this purpose within the distribution range of the Tamaulipan thornscrub in northeastern Mexico. The rate of carbon sequestration and an economic analysis of carbon flux in the irrigation districts require further study. The former is a function of the observed turnover time and the amount of active carbon in the system [40]. Irrigation schedules must be planned properly since leaching of humic substances can be a factor of reducing carbon sequestration in soils.

The potential of carbon sequestration by shifting from conventional tillage to no-tillage in agricultural soils is dependent on several factors [41]. In the US, carbon can be sequestered on average in the range of 337 ± 108 Kg C ha-1 yr-1 to a depth of 30 cm., and these rates vary from 300–600 kg C ha-1 yr-1 in the Great Plains of US, to 100–500 kg C ha-1 yr-1 in the Canadian prairie region [37]. This rate of carbon sequestration can continue for 20 years, and then it will decline [39]. Therefore, conversion from conventional tillage to no-till farming of the 3,500 km2 of the irrigation districts could potentially sequester 2.4 Tg C in 20 years. Information on C sequestration by tillage operations is scarce in northeastern Mexico, however, our data on abandoned agricultural lands suggest that C sequestration in fallow lands could be on the order of 1000 kg ha-1 yr-1 when organic carbon in roots and above ground biomass is included in the estimate.

Economic incentives to promote carbon-conservation and sequestration practices in native forests and soils are not presently recognized by the United Nations Framework Convention on Climate Change. However, carbon sequestration in agricultural ecosystems soils may later be added to the Kyoto Protocol or other protocols that may develop as a response to climate change, and then carbon credits or payments coming from international funds would be available to farmers. In the meantime, Mexican forestry law is promoting local markets for the environmental benefits provided by forest ecosystems.

The deforestation is a cause of concern since world wide it is close to 13 M ha. Gurney and Raymond [42] proposed the approach called the 'Preservation Pathway' that combines the desire for forest preservation with the need to reduce emissions associated with forest loss by focusing on the relative rate of change of forest cover as a criteria by which countries gain access to trading preserved forest carbon stocks. To reduce deforestation in poor countries, Kindermann et al., [43] advanced a combination of incentives and taxes as financial mechanisms. The Preservation Pathway with a combination of incentives and taxes can be downscaled at the ecosystem level, since tropical countries and tropical ecosystems show the largest potential for carbon sequestration despite projections that most of the agricultural expansion will be in these regions [44].

Conclusion

Semi-arid and subtropical semi-arid ecosystems comprise 28% of the globe [45]. Land-use changes in the Tamaulipan thornscrub ecosystem of northeastern Mexico have instigating changes in carbon stocks. Carbon stocks were observed primarily in soils (82%), standing aboveground biomass (10%), and roots (8%). Currently land-use changes have contributed to approximately 2.19 Tg C yr-1 emissions by clearing approximately 300 km2 yr-1 and abandoning approximately 300 km yr-1 to fallow lands. Several practices are discussed and recommended to conserve and sequester carbon in soils. This information is critical for sustainable management and understanding of several of the environmental services provided by the Tamaulipan thornscrub of northeastern Mexico. This type of analysis would be useful to conduct in other regions of subtropical, semi-arid ecosystems when information on carbon flux and land-use change are poor and where there is opportunity for improved carbon storage by means of changes in land-management practices.

Methods

The Tamaulipan thornscrub forest of northeastern Mexico

The Tamaulipan thornscrub forest is a subtropical, semi-arid ecosystem that distributes in the Mexican states of Coahuila, Nuevo Leon, and Tamaulipas, covering a total area of 32,188 km2 according to the 2000 forest inventory by Palacios-Prieto [21]. See Figure 1. This ecosystem is limited to the northwest by the Chihuahuan Desert, to the west by the Sierra Madre Oriental mountain range, and to the south by the tropical rainforest of the Sierra Azul mountain range [46]. Native thornscrub forests historically dominate land cover in this region. However, shifting cultivation is a typical practice in the community-based land tenure system, ejido, which is rapidly reducing the area of the thornscrub ecosystem. Like other subtropical thornscrub forests in the world, the ecosystem has been extensively used as pastureland for the last 350 years [47] and, currently, is used for fuel, timber for fencing and construction, food, and medicinal plants [48, 49]. Even though selective harvesting is common, it is lightly carried out and does not have a noticeable impact on the Tamaulipan thornscrub structure. Grazing is the major disturbance that promotes important structural changes in pasture lands [50, 51].

In northeastern Mexico, the Tamaulipan thornscrub is a dense and diverse ecosystem. It contains at least four endemic genera of woody plants, which indicate a significantly high diversity of woody plant species [52]. It is also rich in terms of native cacti, as well as endangered plant species [53]. Spiny shrubs and low trees dominate the thornscrub, but grasses, forbs, and succulents are also prominent [54]. Romero [55] and Manzano and NĂ¡var [50] recorded 22 shrub species in 0.1 ha plots and more than 5000 shrubs per ha in 0.025 ha plots. Medium and small shrubs (less than 10 m in height) are common life forms, and they are disappearing because of selective harvesting for fuel wood and timber extraction and a reduction in the cycle of shifting cultivation practices in the area.

The southern part of the ecosystem is characterized by a moist, subtropical climate typical of southeastern Nuevo Leon and western Tamaulipas, while the northern bordering region is characterized by a semi-arid climate. Average annual precipitation ranges from 400 to 500 mm in the northern part of the three-bordering Mexican states; to 700 to 800 mm at the south central region, to 1000 mm at the piedmont of the Sierra Madre Oriental mountain range. Pan evaporation is less variable than annual precipitation and approximates 2000 mm in the plains of the northern Gulf of Mexico.

Litosols and rendzins dominate the hilly slopes at the piedmont of the Sierra Madre Oriental mountain range and the smaller isolated mesetas of the plains of the northern Gulf of Mexico. Yermosols and xerosols are widely distributed in the semi-arid western and northwestern region, and deep Vertisols dominate the lowlands of the plains of the northern Gulf of Mexico.

Carbon stocks and fluxes

Carbon stocks in the Tamaulipan thornscrub ecosystem were estimated from biomass components in (a) aboveground vegetation, (b) roots, and (c) soils. The soil carbon of the litter layer was not included in the measurements since late field observations indicate that it may account for less 1 Mg ha-1 and have large spatial variation.

Within the Tamaulipan thornscrub of northeastern Mexico, 56 plots (each 5 m Ă— 40 m) were used to estimate aboveground biomass to derive carbon stocks in three locations. Quadrats were systematically placed following the distribution of the thornscrub from the northwestern to the southeastern portion: (a) the northwestern (NW) portion, covering the northern part of the states of Coahuila, Nuevo Leon, and Tamaulipas (21 plots), (b) the south-central (SC) region of Nuevo Leon (20 plots), and (c) the piedmont of the eastern Sierra Madre (SM) mountain range in the state of Nuevo Leon and Tamaulipas (15 plots). At each location, quadrats were stratified randomly placed to account for typical conditions and major sources of variation in the physical characteristics of the environment. All woody shrubs were measured for basal diameter (d), top height (h), and vertical projection of canopy cover (ct) by species (s). These data provided information to estimate aboveground biomass by applying 17 allometric equations derived by NĂ¡var et al. [56] to estimate total tree biomass.

There is a smooth transition of Tamaulipan thornscrub from subtropical to semiarid climates. Therefore, the area covered by this ecosystem in the semi-arid portion of Nuevo Leon and Tamaulipas is not well known. Therefore, average biomass densities for quadrats measured in the semi-arid (NW) and subtropical (SC and SM) plots were utilized to estimate biomass for these states. The NW portion is characterized by P. glandulosa, the SC by D. texana, P. pallens, and A. berlandieri, and the SM by C. boissierii, B. myricaefolia, and P. pallens as the most important plant species. This approach assumes that 50% of the area is dominated by P. glandulosa (semi-arid Tamaulipan thornscrub) and the remaining 50% by D. texana, C. boissierii and P. pallens (subtropical Tamaulipan thornscrub).

An independent equation derived by NĂ¡var et al., [57] was used to estimate biomass of roots with diameters larger than 1 mm. The equation uses total aboveground biomass as an independent variable and estimates root biomass in the A and B soil horizons. The same weighting procedure to calculate biomass density for each of the Mexican states was used to estimate root biomass.

Data on soil organic carbon (SOC) for the SC and SM locations were provided by Bravo [33] for undisturbed thornscrub forests, abandoned fallow lands, and for agricultural soils cultivated for 15 years. Undisturbed thornscrub forests and agriculture soils have approximately 43 and 15 g kg-1, respectively, in the top soil horizon, and SOC decays hyperbolically as a function of time and soil depth [58] as has been observed by Post et al., [59] in other parts of the world. The depth of tillage operations in the lowland plains of northern Mexico rarely exceeds 50 cm. Therefore, computations of SOC stocks and in our releases during cultivation and gains in fallow lands considered this soil depth as the most dynamic SOC stock in the soil system. Estimates of soil mass in the soil profile (50 cm of soil depth) were calculated according to:

M = Pb*V,

Where: M = soil mass (Mg ha-1), Pb = soil bulk density (Mg m-3); V = soil volume (m3 ha-1). These equations were empirically derived from data reported by Bravo [58] for SOC, which are consistent with the mathematical functions reported by Post et al., [59] and for soil bulk density from NĂ¡var et al., [60]. The equations are:

S O C ( t ) = 14.193 + 26.60 * 1.81 1.81 * t P b ( z ) = 0.9223 Z 0.063 S O C ( z ) = 40.80 * E x p − 0.0089 Z MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqbaeaabmqaaaqaaiabdofatjabd+eapjabdoeadjabcIcaOiabdsha0jabcMcaPiabg2da9iabigdaXiabisda0iabc6caUiabigdaXiabiMda5iabiodaZiabgUcaRKqbaoaalaaabaGaeGOmaiJaeGOnayJaeiOla4IaeGOnayJaeGimaaJaeiOkaOIaeGymaeJaeiOla4IaeGioaGJaeGymaedabaGaeGymaeJaeiOla4IaeGioaGJaeGymaeJaeiOkaOIaemiDaqhaaaGcbaGaemiuaaLaemOyaiMaeiikaGIaemOEaONaeiykaKIaeyypa0JaeGimaaJaeiOla4IaeGyoaKJaeGOmaiJaeGOmaiJaeG4mamJaemOwaO1aaWbaaSqabeaacqaIWaamcqGGUaGlcqaIWaamcqaI2aGncqaIZaWmaaaakeaacqWGtbWucqWGpbWtcqWGdbWqcqGGOaakcqWG6bGEcqGGPaqkcqGH9aqpcqaI0aancqaIWaamcqGGUaGlcqaI4aaocqaIWaamcqGGQaGkcqWGfbqrcqWG4baEcqWGWbaCdaahaaWcbeqaaiabgkHiTiabicdaWiabc6caUiabicdaWiabicdaWiabiIda4iabiMda5iabdQfaAbaaaaaaaa@75F5@

Where: t = time (years), Z = soil depth (cm). Note that 40.80 g kg-1 is the initial SOC content at the top soil surface for the Tamaulipan thornscrub. When Tamaulipan thornscrub forests are cleared, 40.80 g kg-1 decays hyperbolically as a function of time and exponentially as a function of soil depth.

Changes in SOC recorded during fallow are also available for these computations. The rate of increase is linear, as reported by Bravo [58]:

SOC(t) = 14.5 + 0.3277* fallow·length

The rate of carbon increase, in aboveground biomass, AC(t), and root biomass, RC(t), were derived from choronosequence biomass data collected in abandoned fallow lands. The equation used to estimate root biomass from aboveground biomass was reported previously by NĂ¡var et al. [61]. Biomass growth reported from chronosequences is consistent with periodic and mean annual biomass increments recorded by NĂ¡var et al., [24]:

AC(t) = 0.0164t2.2422

RC(t) = 1.465 + 0.56 * [0.0164t2.2422]

The power equation to estimate aboveground biomass growth works well for fallow lands with less than 20 years of abandonment. Data from fallow lands with longer periods of time since abandonment is required to fit the classical s-shape function to biomass growth as a function of time.

Soil carbon data for the semi-arid NW region was not available (SOCnw); hence, it was calculated by multiplying the average SOC for the SC and SM (SOCsc and SOCsm) by the ratio of average total aboveground biomass for the NW (TAGBnw) and the average aboveground biomass for SC and SM (TAGBsc and TAGBsm locations). That is, (SOCnw/SOCsc, sm) = (TAGBnw/TAGBsc, sm). The justification for this equation is that plant cover controls the soil organic carbon content.

Land-use changes

The spatial changes in plant cover through time were measured using Landsat MSS images from February 1980 and March 1996. These data were geometrically corrected, and an unsupervised classification was used to define spectral classes based on a maximum -likelihood clustering algorithm. Ten different spectral classes were obtained, identified, and confirmed during field work. Vegetation maps for 1980 and 1996 were overlaid in order to estimate the fate of land cover with special emphasis on Tamaulipan thornscrub ecosystem, agricultural uses, and grasslands.

Additional data on land-use changes were obtained from several sources [17–20] for different periods and locations encompassing the Tamaulipan thornscrub of northeastern Mexico. An average deforestation rate was estimated for the area covered by Tamaulipan thornscrub using land-use changes observed in this study and reported land-use changes, consistent with the total area covered by Tamaulipan thornscrub in Mexico and US for the 1950's [62, 63] and 2000 [31]. An average deforestation rate has the advantage that it is consistent with other regional land-use change estimates and smoothes the short spatial and temporal changes, but it has the disadvantage in that it is fixed over time, since land-cover change is a disjointed process, with irregular periods of rapid change [64].

Procedure for estimating carbon stocks and fluxes

The procedure for estimating carbon stocks and fluxes employed five steps: 1) an estimation was made of the average annual deforestation rate; 2) projections of the area covered by Tamaulipan thornscrub was done using the average annual deforestation rate and as a baseline the area recorded from the last forest inventory of 2000 [21]; 3) the rate of abandoned, fallow lands was estimated from other plant communities (grasslands, agriculture, scrublands) shifting to Tamaulipan thornscrub as a function of the area covered by Tamaulipan thornscrub and subtracted from the projected area covered by Tamaulipan thornscrub; 4) carbon stocks in aboveground and root biomass were derived from biomass estimates using allometric equations and employing a biomass to carbon content factor of 0.50; and 5) SOC mass estimates were based on projections of soil bulk density with soil depth, soil mass per hectare basis, the annual rate of SOC depletion by shifting cultivation using the hyperbolic equation, the annual SOC increment in abandoned, fallow, agricultural lands.

The major assumptions of this methodology are: a) the annual rate of deforestation and abandoned, fallow lands is constant over time, b) there is a continuous shifting cultivation cycle of clearing, agricultural use, and fallow lasting 15 years; c) there is a mosaic of farms with different times of agricultural use and fallow; d) in agricultural lands SOC depletes in time to a minimum over 15 years as specified by equation [2] and in fallow lands SOC recovers in 15 years according to equation [5]; aboveground and root biomass recover during fallow according to equations [6] and [7], respectively; e) the carbon content in biomass is 0.5; d) the soil bulk density function with soil depth remains constant in Tamaulipan thornscrub, agriculture, and fallow lands. This methodology is quite similar to other procedures widely used to estimate carbon stock and fluxes [13–15, 29].

References

  1. Vitousek PM, Mooney HA, Lubachenko J, Melillo JM: Human domination of Earth's ecosystems. Science 1997, 277: 494–499. 10.1126/science.277.5325.494

    Article  CAS  Google Scholar 

  2. Melillo JM, Houghton RA, Kicklighter DW, McGuire AD: Tropical deforestation and the global carbon budget. Annual Review of Energy and the Environment 1996, 21: 293–310. 10.1146/annurev.energy.21.1.293

    Article  Google Scholar 

  3. Brown S: Los bosques y el cambio climatico: el papel delos terrenos forestales como sumideros de carbono. Actas del XI Congreso Mundial Forestal: Recursos Forestales y Arboles. Antalya Turkia 1997, I: 13–22. 13–22 October of 1997 13-22 October of 1997

    Google Scholar 

  4. IPCC: Climate Change: Impacts, adaptation, and vulnerability. In Summary for policy makers. Cambridge University Press. Cambridge, UK; 2007.

    Google Scholar 

  5. Houghton RA: Land-use cange and the carbon cycle. Global Change Biology 1995, 275–287. 10.1111/j.1365-2486.1995.tb00026.x

    Google Scholar 

  6. Watson RT, Noble IR, Bolin B, Ravindranath NH, Verardo DJ, Dokken DJ: Land use, land-use change, and forestry. In A Special Report. Intergovernmental Panel on Climate Change. Cambridge University Press. Cambridge, UK; 2000.

    Google Scholar 

  7. Dixon RK, Brown S, Houghton RA, Solomon AM, Trexler MC, Wisniewski J: Carbon pools and flux of global forest ecosystems. Science 1994, 263: 185–190. 10.1126/science.263.5144.185

    Article  CAS  Google Scholar 

  8. Houghton RA: Tropical deforestation and carbon dioxide. Climate Change 1991, 19: 99–118. 10.1007/BF00142217

    Article  CAS  Google Scholar 

  9. Schimel DS: Terrestrial ecosystems and the carbon cycle. Global Change Biology 1995, 1: 77–91. 10.1111/j.1365-2486.1995.tb00008.x

    Article  Google Scholar 

  10. Grace PR, Post WM, Hennessy K: The potential impact of climate change on Australia's soil carbon resources. Carbon Balance and Management 2006, 1: 1–14. 10.1186/1750-0680-1-14

    Article  Google Scholar 

  11. Masera OR, Bellon MR, Segura G: Forest management options for sequestering carbon in Mexico. Biomass & Bioenergy 1995, 8: 357–367. 10.1016/0961-9534(95)00028-3

    Article  Google Scholar 

  12. Fang JY, Wang GG, Liu GH, Xu SL: Forest biomass of China: an estimate based on the biomass-volume relationship. Ecological Applications 1998, 8: 1084–1091.

    Google Scholar 

  13. De Jong BHJ, Cairns MA, Haggerty PK, Ramirez-Marcial N, Ochoa-Gaona S, Mendoza-Vega J, Gonzalez-Espinosa M, March-Mifsut I: Land-use change carbon flux between 1970's and 1990's in central highlands of Chiapas, Mexico. Environmental Management 1999, 23: 373–385. 10.1007/s002679900193

    Article  Google Scholar 

  14. Cairns MA, Haggerty PK, Alvarez R, De Jong BHJ, Olmsted I: Tropical Mexico's recent land-use change: a region's contribution to global carbon cycle. Ecological Applications 2000, 10: 1426–1441. 10.1890/1051-0761(2000)010[1426:TMSRLU]2.0.CO;2

    Article  Google Scholar 

  15. De Jong BHJ, Ochoa-Gaona S, Castillo-Santiago MA, Ramirez-Marcial N, Cairns MA: Carbon flux and patterns of land-use/land-cover change in the Selva Lacandona, Mexico. Ambio 2000, 29: 504–511. 10.1639/0044-7447(2000)029[0504:CFAPOL]2.0.CO;2

    Article  Google Scholar 

  16. Post WM, Kwon CK: Soil carbon sequestration and land-use change: processes and potential. Global Change Biology 2000, 6: 317–327. 10.1046/j.1365-2486.2000.00308.x

    Article  Google Scholar 

  17. SARH: Secretaria de Agricultura y Recursos Hidraulicos: Rehabilitacion de los distritos de riego Bajo Rio Bravo (25) y del Rio San Juan (26), Tamaulipas. Resumen. In Infraestructura, Recursos y Servicios. SARH. Mexico. D.F; 1980.

    Google Scholar 

  18. Proderleon: Promotora del Desarrollo Rural de Nuevo Leon: Estadisticas de cambio de uso de suelo 1980–1986. In Gobierno del Estado de Nuevo Leon. Monterrey, N.L., Mexico; 1990:1980–1986.

    Google Scholar 

  19. Maldonado A: DesertificaciĂ³n en la subcuenca del Rio Limon, Linares, N.L. Tesis Profesional de Licenciatura. Facultad de Ciencias Forestales, UANL. Linares, N.L., Mexico 1992.

    Google Scholar 

  20. Treviño E, Akca A, NĂ¡var J, Jimenez J, Aguirre O: Detection of land use change by satellite imagery in the municipality of Linares, Nuevo Leon, Mexico. In Proceedings of the Fifth International Conference on Desert Development. Volume I. Texas Tech University; August 12 – 17, 1996 August 12 – 17, 1996

    Google Scholar 

  21. Palacios-Prieto E: La condiciĂ³n actual de los recursos forestales en Mexico. Resultados del Inventario Nacional Forestal 2000. Investigaciones Geograficas 2000, 43: 183–203.

    Google Scholar 

  22. WWF: Ecoregions of the World. Neartic Desert and Shrublands. Tamaulipan Mezquital (NA1312) 2001.

    Google Scholar 

  23. Siegenthaler U, Sarmiento JL: Atmospheric carbon dioxide and the ocean. Nature 1993, 365: 119–125. 10.1038/365119a0

    Article  CAS  Google Scholar 

  24. NĂ¡var J, NĂ¡jera J, Jurado E: Preliminary estimates of biomass growth in the Tamaulipan thornscrub of northeastern Mexico. Journal of Arid Environments 2001, 47: 281–290. 10.1006/jare.2000.0708

    Article  Google Scholar 

  25. Heiseke D: El matorral como recurso forestal. Reporte cientĂ­fico NĂºmero 1. Facultad de Ciencias Forestales UANL, Linares, N. L., MĂ©xico 1985, 31.

    Google Scholar 

  26. Carstens A: Gottinger Beitrage zur Land – und Forstwirtschaft in der Tropen Subtropen. In Heft 40. University of Gottingen. Gottingen, Germany; 1987.

    Google Scholar 

  27. Castellanos J, Maass M, Kummerow J: Root biomass of a dry deciduous tropical forest in Mexico. Plant and Soil 1991, 131: 225–228. 10.1007/BF00009452

    Article  Google Scholar 

  28. SEMARNAT: Secretaria del Medio Ambiente y Recursos Naturales: Estadisticas de Deforestacion en Mexico: 1993–2000. Mexico, D.F. 2001.

    Google Scholar 

  29. Dirzo R, Garcia C: Rates of deforestation in Los Tuxtlas, a netropical area in southeast Mexico. Conservation Biology 1991, 6: 84–90. 10.1046/j.1523-1739.1992.610084.x

    Article  Google Scholar 

  30. Cairns MA, Dirzo R, Zadroga F: Forests of Mexico: a diminishing resource? Journal of Forestry 1995, 93: 21–24.

    Google Scholar 

  31. Ochoa-Gaona S: Traditional land-use systems and patterns of forest fragmentation in the highlands of Chiapas, Mexico. Environmental Management 2001, 27: 571–586. 10.1007/s002670010171

    Article  CAS  Google Scholar 

  32. NĂ¡var J, Synnott TJ: Soil erosion and land use in northeastern Mexico. Terra 2000, 8: 247–253.

    Google Scholar 

  33. Bravo MRG: Distribucion de la materia organica del suelo en ecosistemas naturales e inducidos en el estado de Nuevo Leon, Mexico. Tesis Profesional de Licenciatura. Facultad de Ciencias Forestales. Universidad Autonoma de Nuevo Leon. Linares, N.L., Mexico; 1999.

    Google Scholar 

  34. Boyd R, Uri ND: A note on the use of conservation practices in US agriculture. Environmental Monitoring and Assessment 2001, 72: 141–178. 10.1023/A:1012088919587

    Article  CAS  Google Scholar 

  35. Campell CA, Zentner RP, Selles F, Liang BC, Blomert B: Evaluation of a simple model to describe carbon accumulation in Brown Chernozem under varying fallow frequency. Canadian Journal of Soil Science 2001, 81: 383–394.

    Article  Google Scholar 

  36. Hao X, Chang C, Lindwall CW: Tillage and crop sequences on organic carbon and total nitrogen content in an irrigated Alberta soil. Soil & Tillage Research 2001, 62: 167–169. 10.1016/S0167-1987(01)00222-7

    Article  Google Scholar 

  37. Follett RF: Soil management concepts and carbon sequestration zin cropland soils. Soil & Tillage Research 2001, 61: 77–92. 10.1016/S0167-1987(01)00180-5

    Article  Google Scholar 

  38. Salinas-Garcia JR, Baez-Gonzalez AD, Tiscareño-Lopez M, Rosales-Robledo E: Residue removal and tillage interaction affects on soil properties under rain-fed corn production in Central Mexico. Soil and Tillage Research 2001, 59: 67–79. 10.1016/S0167-1987(00)00187-2

    Article  Google Scholar 

  39. Lal R, Follett RF, Kimble J, Cole CV: Managing US cropland to sequester carbon in soils. Journal of Soil and Water Conservation 1999, 54: 374–381.

    Google Scholar 

  40. Harrison KG, Post WM: Soil carbon turnover in a recovering temperate forest. Global Biogeochemical Cycles 1995, 9: 449–454. 10.1029/95GB02380

    Article  CAS  Google Scholar 

  41. Yang XM, Kay BD: Impacts of tillage practices on total, loose- and occluded-particulate, and humified organic carbon fractions in soils within a field in southern Ontario. Canadian Journal of Soil Science 2001, 81: 149–156.

    Article  CAS  Google Scholar 

  42. Gurney KR, Raymond L: Targeting deforestation rates in climate change policy: a preservation patway approach. Carbon Balance and Management 2008, 3: 2. 10.1186/1750-0680-3-2

    Article  Google Scholar 

  43. Kindermann JE, Obersteiner M, Rametsteiner E, McCallum I: Predicting the deforestation-trend under different carbon-prices. Carbon Balance and Management 2006, 1: 1–15. 10.1186/1750-0680-1-15

    Article  Google Scholar 

  44. Minnen GJV, Strengers BJ, Eickhout B, Swart RJ, Leemans R: Quantifyimg the effectiveness of climate change mitigation through forest plantations and carbon sequestration with an integrated land-use model. Carbon Balance and Management 2008, 3: 3. 10.1186/1750-0680-3-3

    Article  Google Scholar 

  45. UNSO-UNDP: Office to combat desertification and drought. Aridity zones and dryland population. An assessment of population levels in the worlds drylands. New York. USA. 1997, 32.

    Google Scholar 

  46. GonzĂ¡lez M: El lĂ­mite sur de la Provincia BiĂ³tica Tamaulipeca. II Simposio Internacional sobre la Provincia BiĂ³tica Tamaulipeca UAT-UNAM 1985, 7.

    Google Scholar 

  47. Del Hoyo E: Historia del Reino de Nuevo Leon (1577–1723). Editorial Libros de Mexico. Mexico, D.F., Mexico 2nd edition. 1979.

    Google Scholar 

  48. Reid N, Smith DMS, Marroquin J, Beyer-Munzel P: Floristic and structural variation in the Tamaulipan thornscrub, northeastern Mexico. Journal of Vegetation Science 1990, 1: 529–538. 10.2307/3235787

    Article  Google Scholar 

  49. Reid N, Marroquin J, Beyer-Munzel P: Utilization of shrubs and trees for browse, fuelwood and timber in the Tamaulipan thornscrub, northeastern México. Forest Ecology and Management 1990, 36: 61–79. 10.1016/0378-1127(90)90064-I

    Article  Google Scholar 

  50. Manzano MG, NĂ¡var J: Processes of desertification by goats overgrazing in the tamaulipan thornscrub (matorral) in north-eastern Mexico. Journal of Arid Environments 2000, 44: 1–17. 10.1006/jare.1999.0577

    Article  Google Scholar 

  51. Manzano MG, NĂ¡var J, Pando M, Martinez A: Overgrazing and desertification in northern Mexico: highlights on northern region. Annals of Arid Zone 2000,39(3):285–304.

    Google Scholar 

  52. Rzedowski J: Vegetacion de Mexico. Limusa, Mexico 1978, 432.

    Google Scholar 

  53. Suzan UK, Malda G: Monitoreo demografico en 5 cactaceas en peligro de extincion en Tamaulipas. XI Congreso Mexicano de Botanica. Sociedad Botancia de Mexico. Oaxtepec Morelos 1990, 225–228.

    Google Scholar 

  54. Crosswhite FS: Dry country plants of the south Texas plains. Desert Plants 1980, 2: 141–179.

    Google Scholar 

  55. Romero FG: CaracterizaciĂ³n ecolĂ³gica y definiciĂ³n deesquemas de muestreo en el matorral espinoso tamaulipeco del nordeste de MĂ©xico. Tesis Profesional de MaestrĂ­a en Ciencias. Facultad de Ciencias Forestale s. UANL. Linares, N.L., MĂ©xico 1999.

    Google Scholar 

  56. NĂ¡var J, MĂ©ndez E, Graciano J, Dale V, Parresol B: Biomass equations for shrub species of Tamaulipan thornscrub of northeastern Mexico. Journal of Arid Environments 2004, 59: 657–674. 10.1016/j.jaridenv.2004.02.010

    Article  Google Scholar 

  57. NĂ¡var J, NĂ¡jera , Jurado E: Biomass estimation equations in the Tamaulipan thornscrub of northeastern MEXICO. Journal of Arid Environments 2002, 52: 167–179. 10.1006/jare.2001.0819

    Article  Google Scholar 

  58. Bravo MRG: Degradation and natural rehabilitation of Vertisols of northeastern Mexico. In M.Sc. For. Thesis. Faculty of Forestry. University of Toronto. Notario, Canada; 2002.

    Google Scholar 

  59. Post WM, Izaurralde RC, Mann LK, Bliss N: Monitoring and verifying changes of organic carbon in soil. Climate Change 2002, 51: 73–99. 10.1023/A:1017514802028

    Article  Google Scholar 

  60. NĂ¡var J, Synnott T, Pando M, DomĂ­nguez P, Jurado E, Trevino E, Manzano M: Land degradation by land use changes in northeastern Mexico. Proceedings of the Fifth International Conference on Desert Development. Volume I. Texas Tech University; 1996. August 12 – 17, 1996

    Google Scholar 

  61. NĂ¡var J, MĂ©ndez E, Dale VH: Estimating stand biomass in the Tamaulipan thornscrub of northeastern Mexico. Annals of Forest Science 2002, 59: 813–821. 10.1051/forest:2002079

    Article  Google Scholar 

  62. Johnston MC: Past and present grasslands of southern Texas and northeastern Mexico. Ecology 1963, 44: 456–466. 10.2307/1932524

    Article  Google Scholar 

  63. Udvardy MDF: A classification of the biogeographical provinces of the world. IUCN, Morges. Switzerland. Ocas Paper No 18 1975, 48.

    Google Scholar 

  64. Geist HJ, Lambin EF: Dynamic causal patterns of desertification. Bio Science 2005, 12: 817–829.

    Google Scholar 

Download references

Acknowledgements

This report was written while Dr. NĂ¡var was on sabbatical leave at the Environmental Sciences Division of the Oak Ridge National Laboratory. Oak Ridge National Laboratory is managed by UT-Battelle, LLC, for the U.S. Department of Energy under contract DE-AC05-00OR22725.

Author information

Authors and Affiliations

  1. Professor, CIIDIR-IPN Unidad Durango, Sigma s/n. Fracc. 20 de Noviembre II. Durango, Dgo., Mexico, CP 34220

    Jose de Jesus NĂ¡var-Chaidez

Authors
  1. Jose de Jesus NĂ¡var-Chaidez
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding author

Correspondence to Jose de Jesus NĂ¡var-Chaidez.

Additional information

Competing interests

The author declares that he has no competing interests.

Authors' contributions

JN contributed entirely to this manuscript.

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 2.0 International License (https://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

NĂ¡var-Chaidez, J.d.J. Carbon fluxes resulting from land-use changes in the Tamaulipan thornscrub of northeastern Mexico. Carbon Balance Manage 3, 6 (2008). https://doi.org/10.1186/1750-0680-3-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1750-0680-3-6

Keywords

Carbon Balance and Management

ISSN: 1750-0680

Contact us