Revista Internacional de Contaminación Ambiental, vol.. 32, no. 2, 2016
Universidad Nacional Autónoma de México

Available in:

Received: 15 February 2015

Accepted: 15 November 2015

Abstract: Sugarcane cultivation is an agricultural activity of worldwide importance. This crop has been cultivated in Mexico for centuries, impacting important productive areas like the study region called Huasteca Potosina which is located at the Central-East part of Mexico. The relationship between soil conditions, weather and production per hectare allows identifying the edaphological and climatic conditions (aptitude levels) for cultivating sugarcane in the study area. The objective of this research work is to analyze the relation between the cultivated hectares and production of sugarcane. Likewise, to explain its behavior and interpret the contribution of N and C released to the atmosphere for burning sugar cane in association with events and meteorological parameters involved in the flux of water between soil and atmosphere such as drought, evapotranspiration and rainfall. The methods and indicators for burning waste provided by the Intergovernmental Panel on Climate Change (IPCC), the value of potential evapotranspiration as Thornthwaite and drought severity index were applied in order to understand the correlations between the drought severity index and potential evapotranspiration, and between C and N release and potential rainfall availability. The results show the increment of harvested hectares and the loss of water from the soil by the increase of periods of drought events and evapotranspiration. The volatilized nutrients in soil varied from 1.32 x 105 to 2.17 x 105 t for C and from 1.32 x 103 t to 2.17 x 103 t for N during the burning of sugarcane, affecting production levels with values fluctuating between 38 t/ha to 77 t/ha for the term of 1990-2010.

Keywords: soil quality, meteorological events, waste burning, water flux.

Resumen: El cultivo de caña de azúcar es una actividad agrícola de importancia mundial llevada a cabo en México durante siglos, impactando importantes áreas productivas como la región de estudio llamada Huasteca Potosina, localizada al este central de México. La relación entre las condiciones del suelo, clima y producción por hectárea permite identificar las condiciones edafológicas y climáticas (niveles de aptitud) para cultivar la caña de azúcar en el área de estudio. Bajo este contexto, el objetivo fue analizar la relación entre las hectáreas cultivadas y la producción de caña de azúcar y explicar su comportamiento al interpretar las contribuciones del N y C liberados a la atmósfera por la quema de caña asociados a eventos y parámetros meteorológicos involucrados en el flujo de agua entre el suelo y la atmósfera como sequía, evapotranspiración y lluvia. Los métodos e indicadores para la quema de residuos provistos por el Panel Intergubernamental del Cambio Climático (IPCC. por sus siglas en inglés), el valor de la evapotranspiración potencial por Thornthwaite y el índice de severidad de la sequía fueron aplicados para comprender las correlaciones entre el índice de severidad de la sequía y la evapotranspiración potencial, así como entre el C y el N liberados y la disponibilidad potencial de lluvia. Los resultados muestran el incremento en las hectáreas cosechadas y la pérdida de agua del suelo por el aumento de los periodos de sequía y la evapotranspiración. Los nutrientes del suelo volatilizados, varían de 1.32 x 105 ton a 2.17 x 105 ton para C y 1.32 x 103 ton a 2.17 x 103 ton para N durante la quema de caña de azúcar, lo que afecta sus niveles de producción con valores que fluctúan entre 38 t/ha y 77 t/ha durante el periodo 1990-2010.

Palabras clave: calidad del suelo, eventos meteorológicos, quema de residuos, flujo de agua.


The sugarcane crop is easily adaptable to tropical and subtropical climates and is of economic importance in more than 100 countries around the world (Rein et al. 2011). Sugarcane cultivation has caused enormous changes in land use with the removal of tropical forest and other vegetation types (Scherr and Sthapit 2009). Some practices related to sugarcane such as burning bring with them the further decline of air and soil quality and affect the regional climate due to variations of soil temperature (Sandhu et al.2013), which can be reflected for example in the heat flux and moisture transport at this scale. The common practice of sugarcane burning methods, before and after the harvest, is an anthropogenic source of emissions of greenhouse gases (GHG) and aerosol to the atmosphere that may also cause respiratory toxicity (França et al. 2012, Mazzoli-Rocha et al.2014). Moreover, the procedure includes the strippage of weed from soil, nutrients uptake, and control of pests among others (Hernández-Valencia and López-Hernández 2002). The sugarcane burning can increase its efficiency if the moisture levels in the fuel are poor (Chen et al. 2010, França et al. 2012).

Soil quality depends on its capability to accept, store and recycle nutrients and water resources in maintaining production yields without affecting the environment (Haynes 1997), but the sugarcane burning can affect several of these factors. Particularly, the process of combustion during vegetation burning contributes to reduce the organic matter in the soil due to the release of nutrients such as carbon and nitrogen. This can be detrimental to crops if the nutrients do not recuperate through natural mechanisms such as rain, atmospheric fixing, or by using nitrogen fertilizers (Hernández-Valencia and López-Hernández 2002).

Potential evapotranspiration (PET) is an indicator of the maximum evaporation from the soil and plant transpiration and depends on weather conditions of the study region where factors such as temperature, solar radiation, and water availability among others have an important role. There are several methods to calculate PET (Tukimat et al. 2012, De Carvalho et al. 2013), where a variety of parameters are used to analyze the atmospheric forcing (Bautista et al. 2009). Comparing the change in the water vapor concentration that occurs in the air with precipitation levels, the water flux between the atmosphere and the soil surface, it is used (Swenson and Wahr 2005) to know the potential rainfall availability at the surface (SENASICA 2010). In humid regions, such as in the Huasteca region in Mexico, land use change from natural vegetation to cropping and grazing occurs (Reyes et al. 2005). Thes echanges may cause variations in humidity and precipitation levels, affecting negatively the availability of water resources, leading to an increase in drought periods (Algara-Siller and Contreras-Servin 2009).

There are several ways to quantify drought periods, such as the Palmer drought severity index (PDSI), the drought severity index (DSI), and the standardized precipitation index (SPI) among others. To calculate the PDSI it is necessary to know the precipitation, the temperature and the soil moisture (Dai et al. 2004). However, in some regions such as San Luis Potosí State, Mexico, this last parameter is not measured by weather stations. On the other hand, the DSI method is defined as a statistical analysis based on the deficit of precipitation which is commonly calculated with data usually registered by regular meteorological stations. Since it does not depend on the type of climate, it can be applied to any region (Hernández and Valdez 2004, Algara-Siller and Contreras-Servin 2009).

Methane (CH4) and nitrous oxide (N2O) are released to the atmosphere by agricultural activities significantly contributing to total emissions of GHG and to radiative forcing due to their effectiveness to absorb heat compared to carbon dioxide (CO2) (Scheele and Kruger 2006). The IPCC developed a series of methodologies to quantify emissions. Inventories of GHG in Mexico have included agricultural waste burning, dominantly of sugarcane crops (Ordóñez et al. 2012). However, these inventories are still deficient in variables that have an impact on the loss of soil nutrients. It is also necessary to investigate the fraction of carbon in organic waste to obtain national emissions factors (Ordóñez and Hernández 2006) and the annual amount of organic matter andcarbon in different types of soil (Sánchez-Hernández et al. 2011), among other issues. The calculation of C and N released during sugarcane burning depends on parameters such as fraction of burned area, aboveground biomass density and aboveground burned biomass, among others. These values can be obtained specifically for the study region (Tier 2) or by using the IPCC proposed standard values (Tier 1) (IPCC 1996).

Mexico is the seventh-largest world producer of raw material for the sugar industry. Sugarcane fields cover up 790 650 ha supplying raw material to 57 sugar mills, located in 15 states and 6 regions of Mexico (Fig. 1 and Fig. 2). Sugarcane represents 3 % of total Mexican crop area (Aguilar-Rivera et al. 2012, Aguilar-Rivera et al. 2013, CONADESUCA 2014). The conventional harvest system in Mexico starts when the sugarcane is burned before manual cutting is carried out in order to eliminate dry foliage and to facilitate the harvest process. After removing the cane, there is a waiting period until the remaining foliage is completely dry. Finally, a second burning is applied. The pollutants that are emitted to the atmosphere are nitrogen monoxide, carbon monoxide, sulfur oxide, hydrocarbons, carbon dioxide and other substances (Vilaboa and Barroso 2013, Sentíes-Herrera et al. 2014).

Fig. 1
National distribution of sugarcane production

Sugarcane fields shown in the centralpicture and sugar mills shown in the East zone and in the West zone pictures (author’sown production)

The Mexican government has implemented efforts to decrease the burning of sugarcane through the National Program of the Sugarcane Agroindustry (PRONAC) (SEGOB 2014), in which, for example, for the period 2006-2012 the target was to reduce the cane burning from 90 % to 70 % (SAGARPA 2007, FIRCO 2012). Particularly, in the cane region of San Luis Potosí State, green harvest is practiced in 5 % of the total cane production. Additionally, there are experimental cane fields in which the government implement strategies that contribute to the gradual recovery of ecological balance (CNIDTACHP 2014).

On the other hand the use of edaphological mapping and satellite imagery allow identifying spatial attributes of soil fertility for sugarcane crop. The method is based on physical, chemical and biological soil properties affecting the production levels in a specific sugarcane region (Ortega and Corvalán 2007, Quintero-Durán 2008, Geissen et al. 2009, Aguilar-Rivera et al. 2010, Aguilar-Rivera et al. 2013). The soil characteristics are better classified when taking into account the weather conditions that interact directly with sugarcane crops, showing a robust basis for a classification on edaphological and climatic coditions under wich the cane cultivation takes place (Subiros 2002, Sánchez-Hernández et al. 2011), thus the term “aptitude levels” is adopted from now on.

According to previous studies, during the period 1990-2010 while the field performance (productionyield) decreased, the harvested surface increased substantially (Aguilar-Rivera et al., 2012, Aguilar-Riveraet al. 2013). The hypothesis is that changes in theland use in Huasteca Potosina and in the water flux between soil and atmosphere have contributed to the increase in the release of C and N to the atmosphere.The objective of this research work is to analyze the relation between the cultivated hectares and production of sugarcane. Likewise, explain its behavior and interpret the contributions of N and C released to the atmosphere for burning sugar cane in association with events and meteorological parameters involved in the flux of water between soil and atmosphere such as drought, evapotranspiration and rainfall.


Sugarcane production conditions at the HuastecaPotosina region

The Huasteca region is located in the Central-East part of Mexico covering around of 65 675.85 km2 (CONABIO 2012). An important area of this region is called Huasteca Potosina, located in the eastern part of San Luis Potosí State (Fig. 2). This region shows year-long rainfall conditions most of the year where its maximum precipitation values are in the months of September and July and range from 200 to 500 mm. Likewise, highest monthly mean temperatures varies between 22 ºC and more than 40 ºC and the coldest months have temperatures fluctuating from –3º C to 18 ºC (Pineda et al. 2006). The predominant type of soil is argillaceous Pelic Vertisol (Aguilar-Rivera et al. 2012). This kind of soil reaches deep in the ground, and develops within dry to sub-humid climates. Although it is not considered as a fertil esoil, it has been adapted to sugarcane fields with an agricultural irrigation system (SEMARNAT 2013).

Fig. 2
Study area. Meteorological stations as black dots (author’sown production)

The plant water availability depends largely on soil-atmosphere interaction. In Huasteca Potosina, the study of these elements is carried out by the use of the drought severity index (Hernández and Valdez 2004, Algara-Siller and Contreras-Servin 2009), precipitation (P), potential evapotranspiration (PET) and the relation between the latter two parameters P-PET to investigate the net flux of water between the soil and the atmosphere (Swenson and Wahr 2005, SENASICA 2010).

The soil parameters based on physical, chemical and biological properties and their relationship with those weather conditions that contribute with sugarcane production levels are shown in table I and they are also used as a starting point to analyze the stages under which carbon (C) and nitrogen (N) are released in the sugarcane burning. This approach includes the classification of production per hectare based on aptitude levels (Table I), the relationship analysis between DSI and PET and the release of N and C during the burning process related to the availability of water in the ground (P-PET).


Sugarcane production data and harvested area used to classify aptitude levels for sugarcane crop for the period 1990-2010 were obtained from the Agro-food system reference information software (SIACON 2010) administrated by the Servicio de Información Agroalimentaria y Pesquera (SIAP) (Agro-food and fishery information system) from the Mexican government. Weather data, specifically precipitation and maximum and minimum temperatures were collected from the Base de Datos Climatológica Nacional (CLICOM) (National Climatological Database) (CLICOM 2014) whose data are administrated by the Servicio Meteorológico Nacional (NMS) (The National Meteogological Service).

Thornthwaite potential evapotranspiration anddrought severity index

Potential evapotranspiration represents the maximum possible evapotranspiration under existing circumstances. There are few models for its quantification (Tukimat et al. 2012, De Carvalho et al. 2013). In recent years Thornthwaite method has been utilized to calculate potential evapotranspiration, particularly for humidity balance studies (SEMARNAT 2009, Rivas et al. 2011). Moreover, it has shown great reliability during rainy periods, disregarding any adjustment (Bautista et al. 2009, Gutiérrez et al. 2013). Based on this evidence, due to the humid climate of the Huasteca (Pineda et al. 2006) and data availability, the calculation of potential evapotranspiration parameter by this method is the selected choice for the present study.


where Hs is maximum number of sunlight hours, dm is number of days per month and ET is evapotranspiration without corrections, which in turn is rewritten as


where the constant C=16, tmed is monthly mean temperature, and I is the annual heat index expressed as I = Σ i where


Finally, a is an exponent on a logarithmic basis for the annual index, where both are related by a third polynomial order


In the same way, the precipitation (P) minus the potential evapotranspiration (PET) allows the calculation of the net flux of water between soil and atmosphere or the potential rainfall availability, namely P-PET (Swenson and Wahr 2005, SENASICA 2010).

The drought severity index (DSI) shows the percentage of rain scarcity in a year with respect to average rain values and for every meteorological station. It is calculated from the following statistical analysis (Hernández and Valdez 2004, Algara-Siller and Contreras-Servin 2009)


where ΣX – is the sum of the accumulated monthly precipitation (ΣX) for the year calculated minus the annual average precipitation ( ) and where n represents the number of weather stations. The yearly mean DSI considering all stations is:


Drought severity index comprises seven levels as shown in table II. These levels show that the DSI is strong when the negative value increases, while the strength of DSI decreases when it is close to zero.

Calculation of carbon and nitrogen release

The parameters based on IPCC guides to calculate the release of C and N values during the sugarcane burning, correspond to the information available for San Luis Potosí state. The amount of waste and the amount of dry residue described in equations (7) and (8), the total burned biomass (equation (9)), and C and N release (equations (10) and (11)) are obtained by the IPCC (2000) methods. Some variable values used in the equations correspond to the geographical region, which contains the study area.

After reviewing previous studies, it was found that the value of 0.2230 corresponding to fraction of C in waste (FCW) (Table III) obtained from Mexican reference (Ordóñez and Hernández 2006) represents 50 % of the IPCC value. This value has not been verified for sugarcane crops in the region of Huasteca Potosina. Likewise, there are other parameters such as fraction burned in fields (FBF) and fraction oxidized (FO) that have not been calculated for Mexico. These are the main reasons why in this study the reference data to calculate the C and N release during the burning of sugarcane were obtained using the values of the IPCC methods.

The Mexican reference values of nitrogen-carbon ratio (NCR) are obtained from the inverse of C/N, shown in table I as a first approach, while the value used in this research was obtained from IPCC methods (Table III; IPCC 1996). Also, the AP variable in equation (7) corresponds to sugarcane annual production for Huasteca Potosina. Finally, the values of residues-crop ratio (RCR) and dry matter fraction are proposed by Ordóñez and Hernández (2006).







Annual aptitude level classification for sugarcane crop

The relationship between extension of harvested sugarcane, its production per hectare considering both irrigated and rain-fed, and its relationship to aptitude levels (Table I) is shown in figure 3. The mean level was between 55 and 80 t/ha.

Fig. 3
Aptitude levels for sugarcane for the period 1990-2010 based on production levels given in t/ha

The correlation coefficient between the harvested area and the production per hectare is -0.1

Author’s own production

The slope in the production per hectare is negative and the trend of the crop surface is towards growth, hence the correlation coefficient between both parameters is very small (-0.1) because the aptitude levels of sugarcane keep a relation with specific weather and soil conditions. For example, relating the production per hectare (Fig. 3) to aptitude levels during the study period, the DSI presented a change from slight to severe (Table I). Also, the C/N relation ranged from 12 to 30 affecting the annual production per hectare of this crop (Table I). Thus, an increase in cultivated hectares does not imply growth in sugarcane production (Fig. 3) because of the influence of soil and climate parameters such as PET, DSI, PRA (potential rainfall availability) and the release to the atmosphere of C and N as will be discussed further.

Potential evapotranspiration (PET) and drought severity index (DSI) analysis

While the PET tendency increased in the 1990- 2010 period, the DSI remained practically constant with a very strong severity index ranging between –0.4 and –0.5 (Table II) and consequently their correlation coefficient was 0.6 (Fig. 4). However, for the first eight years it reached a value of 0.8 due to an increase of DSI and a decrease of PET because of deficit precipitation and decrease of water flux to the atmosphere respectively. For the period 1998-2010 this correlation decreased to 0.45 due to the increase of both parameters. This phenomenon is mainly present in 1999 where both parameters had contrast high values.


Fig. 4
Comparison between DSI (drought severity index) and PET (potential evapotranspiration) for the 1990-2010 period

The correlation coefficient between both parameters is 0.6

Author’s own production

In contrast, the annual variation levels for the DSI ranged from slight (-0.2 to -0.35) to severe (–0.5 to –0.6) drought conditions as observed for the years 1992 and 1999 respectively (Fig. 4). These calculated values from the DSI method are consistent with the information obtained from yearly production per hectare and aptitude levels of sugarcane (Table I and Fig. 3).

As mentioned above, the rise of evaporation and transpiration from soil and plants respectively increase the potential evapotranspiration. Also, if both phenomena are combined with an increase of precipitation levels as occurred from 1990 to 1992, then the DSI values decrease. However, in 1993 and 1994 the precipitation levels decreased by 60 % compared to 1991 (i.e. from 1529 mm to 912 mm) helping to increase the DSI and reducing the PET (Fig. 4) due to plants and soil not having the same water availability as in previous years.

The DSI indicated a severe drought event (–0.5) in 1995 and 1997 where lowest PET values were registered with 261.9 and 261.8 mm respectively (Fig. 4). Meanwhile, the same years presented low accumulated precipitations of 879 and 905 mm. These values are in contrast with the highest precipitation of 1529 mm in 1991, 1516 mm in 1993, 1411 mm in 2007 and 1464 mm in 2008, where there was a greater availability of water.

Year 1996 reveals that a severe drought condition of –0.5 is accompanied by a high PET of 271 mm. This value was derived from the second maximum temperature (30 ºC) in the year of the study period. Likewise, the highest mean temperature of 31.8 ºC corresponds to 1998 which is reflected in an increase of the PET value (266 mm) with respect to 1997.

The strongest increment of DSI, up to a severe drought event of -0.68 is found in 1999 and is associated with an increase of PET. The period 2000 to 2008 maintained values for PET above the mean of 271.5 mm while the DSI value oscillated from strong to very strong levels. During the last two years of the study, the DSI varied from very strong to severe drought condition in chronological order.

Carbon and nitrogen release and their relation with potential rainfall availability (PRA)

The water flux between soil and atmosphere can be reflected in the analysis of the precipitation (P) minus the PET (P-PET) and is termed potential rainfall availability (PRA) (SENASICA 2010).

Both C and N release showed levels above its annual average, 17.8 ×104 t and 17.8 × 102 t respectively when the PRA (P-PET) was below annual average value (10.8 × 102 mm) from 1994 (Fig. 5). In general, a PET increase and its derived PRA (P-PET) decrease provoked an increment on the release of both gases for the period 1990-2010 (Fig. 5a and 4b). The linear tendency shown in figure 5a and 5b revealed a positive slope for C and N release and a negative slope for P-PET. This behavior is due first to the decrease of precipitation levels from 1994 to 2006, mainly below 1000 mm, and secondly to the increase of PET (> 270 mm) during the period.

Fig. 5
Annual averages for C (a) and N (b) release during sugarcane crop burning and precipitation (P) minus potential evapotranspiration (PET) (P-PET), which is the potential rainfall availability (PRA)
Author’s own production


The variability in PET and DSI shown in figure 4 are related with the increase or decrease occurred in the release of C and N (Fig. 5a and 4b) during the sugarcane burning due to water availability while the combustion process occurs. However, the increase in cultivated areas also favored the rise of soil evaporation and plant transpiration, reflected in PET parameter (Fig. 4) and the growth of burned areas as well (in the period 2000-2008 and in year 2010).

The water flux between soil and atmosphere is studied with the PRA through the P-PET. But, the correlation analysis between these parameters with the C and N released let us know that for both cases (Fig. 5a and 4b) it was close to zero (–0.28), mainly during the 13 years of the study period (1994-2006), where P-PET was below their average because there was PET increase with respect to precipitation levels. Likewise, for the period 1990 to 1993 the water availability levels reduced the C and N emitted during sugarcane burning, but during the period 2007-2010 the increase of burned areas did not allow these nutrients to reduce their emission rate, although the precipitation levels grew compared to the first years.

The differences in the order of magnitude between C and N release in IPCC and Mexico’s results are indicated in table III. The correlation coefficient in the C release for both analyses (with Mexican reference values (*) and IPCC reference values) is practically equal to 1 because the fraction of carbon in waste (FCW) was constant in every yearly calculus (0.2230 for Mexican values and 0.4325 for IPCC; Fig. 6a). Nevertheless, for the N the value is 0.8 due to variations of nitrogen-carbon ratio (NCR) depending on annual production per hectare levels (tables I and III; Fig. 6b).


Fig. 6
Carbon (a) and nitrogen (b) releases

Calculated according to Intergovernmental Panel on Climate Change and Mexico’s own (*) reference data respectively

Author’s own production

According to sugarcane annual production per hectare (Fig. 2), San Luis Potosí keeps a medium aptitude level (Table I) for most of the year. Although, there is insufficient information to analyze the annual levels of C and N during the study period. The only year in which data exist for N levels in the sugarcane zone (Fig. 3) is 2009, where N availability was less than 20 kg/ha (SAGARPA 2009a, 2009b, 2009c) placing the aptitude level of sugarcane as “not suitable” (Table I). However, this condition does not match with the aptitude level obtained based on the annual production because it was superior (medium level). This means that the floor area of sugarcane is being forced to support the cane crop, affecting the soil health.

The increase of potential evapotranspiration and the decrease of rainfall in the period 2000-2008 and in year 2010 were not favorable for production because the net water flux is principally evaporated to the atmosphere (P-PET) and the sugarcane must extract water at a rate fast enough not to lose it during this process. Also, the decrease in soil water favors greater loss of C and N because of the increase in efficiency in the combustion process (Chen et al. 2010, França et al. 2012).

Aguilar-Rivera et al. (2010) mention that in a complete cycle the sugarcane crop requires 1500 mm of water for its development. However, it is found that there were just two years in which the rainfall reached these levels (1991 and 1993). Therefore favoring the increase of DSI, mainly in the last 11 years of the study period in which there was an important increase in the cultivated hectares and in the PET levels (Fig. 4).

The C and N release in the Huasteca Potosina during sugarcane burning represented practically 0.05 % and 1 % of the total losses of these nutrients in Mexico (Ordóñez et al. 2012). In addition, when comparing these values with those of the other countries such as Brazil and Australia it is observed that the carbon and nitrogen release increased substantially (more than 100 % of the Mexican losses; UNFCCC 2014). The direct affectation is reflected in the land use change due to a growing sugarcane surface displacing natural vegetation as was demonstrated by Peralta-Rivero et al. (2014) for the Mexican Huasteca Region and as it is shown in this research, particularly for the Huasteca Potosina (Fig. 3).

Actually, various strategies exist to manage the C and N in soils. In these particular cases it is imperative to increase the amount of these nutrients in the surface with a minimum soil disturbance. As noted by Lal (2006), it should be in a proportional way to the increase of cultivated hectares. Also, there is an important rise in the quantity of N during the agricultural practices as recommended by SAGARPA (2009a, 2009b, 2009c) from 20 kg/ha to 190 kg/ha. Finally, it is necessary to use sustainable methods to take advantage of the sugarcane’s waste potential in the availability of C, N and water in the soils.


When the cane production level has not a consistent growth related with the increase of the harvested area, then is relevant to analyze the climate and soil conditions under which the cane crop take place. The aptitude levels of sugarcane fluctuated from “not suitable” to “medium” during the period 1990-2010 in the Huasteca Potosina region. Consequently, the production per hectare never reached the optimal levels (>80 t/ha). Furthermore, the weather settings played an important role in each of the studied periods of sugarcane burning where the water flux to atmosphere is related with the potential evapotranspiration. In this study, the potential evapotranspiration (PET) had three types of behavior: it grows when the DSI decrease is helping to keep a water balance (ranged from 1990 to 1994), it decreases when the DSI increases, which affects the water soil availability (years 1995, 1996, 1997 and 2009), and it increases when DSI increases too (1998, 1999 and in the period 2000-2008 and 2010) due probably to an increase in cultivated hectares and the rise of precipitation levels. In the latter two cases the combustion process during sugarcane burning was favored due to soil water loss and the increase of the crop surface, respectively. The analysis that represents the soil-atmosphere water flux was the potential rainfall availability (P-PET). Since 1994 and practically until 2006 the water flux was mostly evaporated to the atmosphere, because P-PET was less than their average value of the study period (10.8 × 102 mm). The conditions in which there was an increase in the C and N release were favored due to a decrease in the amount of water in the soil and consecuently favouring the combustion process of sugarcane during the harvest period. In the years 2007, 2008 and 2010 the annual values of P-PET were greater than the average value due to an increase of precipitation levels. These values did not affect in an important way the release of C and N in the cane burning due to the increase of crop hectares and the number of acres burned. The correlation of analyzed parameters allowed identification of important phenomena, such as the -0.1 value obtained between the production per hectare (with a negative trend) and the crop surface (with a positive slope). Due to the relation between cane production and aptitude levels, as mentioned above, it was possible to identify C and N losses, which was justified with the analysis of DSI, PET and PRA in the study period. Likewise, during the period 1990-1997 the DSI and PET had a correlation value close to one (0.8) due to PET increase. This parameter depends strongly on water availability. In contrast, for the same period, the correlation coefficient between the C and N release and P-PET was close to zero (0.28) due to the decrease in precipitation levels (from 1512 mm to 913 mm) during 1993 and 1994, which helped to increase increase nutrient emissions. Finally, it is relevant to mention agricultural methods and land use strategies to increase the harvested sugarcane without affectation to soil health and climatic conditions. One of these methods is the green harvesting in which there is no release of nutrients to the atmosphere, a method that only involves leaf cutting another method is the use of cane species with a minor percentage of waste to be burned. Although, it is important to analyze the economic cost for farmers in order to know if there really will be a real benefit to them. In contrast, to San Luis Potosí, other Mexican states like Jalisco have fewer hectares to grow sugarcane but still produce more cane per hectare than Huasteca Potosina region. Thus, a land use strategy could be the increase of soil nutrients (mainly N and C) instead of the cultivated surface to support the sugarcane grow.


This research was made through the program “Programa Multidisciplinario de Posgrado en Ciencias Ambientales”, financed by Consejo Nacional de Ciencia y Tecnología, Mexico. Special thanks to Cecilia Barrera-Gutiérrez from the Universidad Autónoma de San Luis Potosí, who provided help in the dry season analysis.


Aguilar-Rivera N. (2013). Ficha técnica del cultivo de Caña de Azúcar. Sistema Potosino de Vigilancia Epidemiológica. Technical Report. Veracruz, Mexico, 20 pp. [online].,_FICHA_T%C3%89CNICA.pdf 21/08/2014.

Aguilar-Rivera N., Galindo M.G. and Fortanelli M.J. (2012). Evaluación agroindustrial del cultivo de caña de azúcar (Saccharum officinarum L.) mediante imágenes SPOT 5 HVR en la Huasteca México. Revista de la Facultad de Agronomía, La Plata 111, 64-74.

Aguilar-Rivera N., Galindo M.G., Contreras S.C. and Fortanelli M.J. (2010). Zonificación productiva cañera en Huasteca Potosina, México. Agronomía Trop. 60, 139-154.

Aguilar-Rivera N., Olvera V.A. and Galindo M.G. (2013). Evaluación de aptitud de tierras de cultivo de caña de azúcar en la Huasteca Potosina, México, por técnicas geomáticas. Revista de Geografía Norte Grande 55, 141-156.

Algara-Siller M. and Contreras-Servin C. (2009). Drought and climate change on the Huasteca region, San Luis Potosí, Mexico. In: Environmental Engineering and Management (M. Theophanides and T. Theophanides, Eds.). ATINER, Athens, Greece, pp. 287-301.

Bautista F., Bautista D. and Delgado-Carranza C. (2009). Calibration of the equations of Hargreaves and Thorn-thwaite to estimate the potential evapotranspiration in semi-arid and subhumid tropical climates for regional applications. Atmósfera 22, 331-348.

Chen L-W.A., Verburg P., Shackelford A., Zhu D., Susfalk R., Chow J.C. and Watson J.G. (2010). Moisture effects on carbon and nitrogen emission from burning of wildland biomass. Atmos. Chem. Phys. Discuss. 10, 7985-8007. DOI: 10.5194/acp-10-6617-2010

CLICOM (2014). Datos climáticos diarios del CLICOM del Servicio Meteorológico Nacional con gráficas del Centro de Investigación Científica y de Educación Superior de Ensenada [online]. 15/03/2014.

CNIDTACHP (2014). Arreglos topológicos y densidad de siembra de caña de azúcar en la Huasteca Potosina. Technical Report. Centro Nacional de Investigación y Desarrollo Tecnológico para la Agroindustria Cañera de la Huasteca Potosina, A.C. [online]. 14/07/2015.

CONABIO (2012). Distribución de la población en México por municipio, 1990, escala 1:250000. Comisión Nacional Para el Conocimiento y Uso de la Biodiversidad. México D.F.: Datos estadísticos del 2010, Instituto Nacional de Geografía y Estadística [online]. 12/12/2014.

CONADESUCA (2014). Reporte de avance de producción de caña y azúcar. Ciclo azucarero 2013/14. Reporte No.36. Semana No. 41. Del 5 al 12 de julio del 2014 - Días corrientes de zafra: 254. Final Report. Comité Nacional para el Desarrollo Sustentable de la Caña de Azúcar. Mexico City, 32 pp.

Dai A, Trenberth K.E. and Qian T. (2004). A global dataset of Palmer drought severity index for 1870-2002: Relationship with soil moisture and effects of surface warming. American Meteorological Society. J. Hydrometeorol. 5, 117-1130. DOI: 10.1175/JHM-386.1.

De Carvalho L.C., Evangelista A.W.P., Oliveira K.M.G., Silva B.M., Alves M de C., Júnior A de S. and Mirando W.L. (2013). FAO Penman-Monteith equation for reference Evapotranspiration from missing data. IDESIA 31, 39-47.

FIRCO (2012). Informe de rendición de cuentas 2006-2012. México. Report. Fideicomiso de Riesgo Compartido. Etapa 1. Mexico City [online]. 11/07/2015.

Franca D de A., Longo K.M., Soares N.T.G., Santos J.C., Freitas S.R., Rudorff B.F.T., Cortez E.V., Anselmo E. and Carvalho Jr. A. (2012). Pre-harvest sugarcane burning: Determination of emission factors through laboratory measurements. Atmosphere 3, 164-180. DOI: 10.3390/atmos3010164

Geissen V., Sánchez-Hernández R., Kampichler C., Ramos-Reyes R., Sepúlveda-Lozada A., Ochoa-Goana S., de Jong B.H.J., Huerta-Lwanga E. and Hernández-Daumas S. (2009). Effects of land-use change on some properties of tropical soils - An example from Southeast Mexico. Geoderma 151, 87-97. DOI:10.1016/j.geoderma.2009.03.011

Gutiérrez E.S., Dávila-Vázquez G., Sánchez J. de A. and Torres J.J.D. (2013). Assessment of hydric balance through climatic variables, in the Cazones River Basin, Veracruz, Mexico. Rev. Ambient. água 8, 104-117. DOI: 10.4136/1980-993X

Haynes R.J. (1997). The concept of soil quality and its applicability to sugarcane production. Proc. S. Afr. Sug. Technol. Ass. 71, 9-14.

Hernández C.M.E and Valdez M.G. 2004. Impactos, vulnerabilidad y adaptación En: Cambio climático: una visión desde México (J. Martínez y A. Fernández Bremauntz, Comp.). Secretaría del Medio Ambiente y Recursos Naturales, Mexico City, 315-325 pp.

Hernández-Valencia I. and López-Hernández D. (2002). Pérdida de nutrimentos por la quema de vegetación en una sabana de Trachypogon Rev. Biol. Trop. 50, 3-4.

IPCC (1996). Revised 1996 guidelines for national greenhouse inventories. Vol. 2, Intergovernmental Panel on Climate Change / Organisation for Economic Co-operation and Development / International Energy Agency. United Kingdom Meteorological Office, Bracknell Published for the Intergovernmental Panel on Climate Change. Manual. Bracknell, England, 84 pp. [online]. 7/05/2014.

IPCC (2000). Good practice guidance and uncertainty management in national greenhouse gas inventories. Manual. IPCC national greenhouse gas inventories programme. Intergovernmental Panel on Climate Change Institute for Global Environmental Strategies. Japan, 517 pp. [online]. 7/05/2014.

Lal R. (2006). Carbon management in agricultural soils. Mitigation and adaptation strategies for global change 12, 303-322. DOI: 10.1007/s11027-006-9036-7.

Mazzoli-Rocha F., Carvalho G., Lanzetti M, Valenga S.S., Silva L.F., Saldiva P.H.N., Zin W.A. and Faffe D.S. (2014). Respiratory toxicity of repeated exposure to particles produced by traffic and sugar cane burning. Respiratory physiology and neurobiology 191, 106-113. DOI: 10.1016/j.resp.2013.11.004.

Ordóñez D.J.A. and Hernández T.T. (2006). Obtención de factores de emisión en el sector agrícola para disminuir incertidumbre en el Inventario Nacional de Emisiones de Gases de Efecto Invernadero. Instituto Nacional de Ecología y Cambio Climático. Report. Mexico City, 22 pp.

Ordóñez D.J.A.B., Hernández T.T. , Galicia N.A., Jiménez Álvarez I.F., Ordóñez D.M.J., Torres O.J.F., León G.J.D., Nuñez R.A., Hernández Q.M., González F.A. and Carrera J.A. (2012). Actualización del Inventario Nacional de Gases de Efecto Invernadero 1990-2010, para el sector Agricultura. Report, Mexico City, 294 pp., [online]. 8/06/2014.

Ortega A. and Corvalán E. (2007). Diagnóstico de suelos. Instituto Nacional de Tecnología Agropecuaria Estación Experimental Agropecuaria Salta. Salta, Argentina, 5 pp. [online]. 5/05/2014.

Peralta-Rivero C., Contreras-Servin C. , Galindo-Mendoza M.G., Mas C.J.F. and Algara-Siller M. (2014). Analysis of land use and land cover changes and evaluation of natural generation and potential restoration areas in the Mexican Huasteca Region. Open Journal of Forestry 4, 124-135. DOI: 10.4236/ojf.2014.42018.

Pineda L., Carbajal N. and Medina R.E. (2006). Region-alization and classification of bioclimatic zones in the central-northeastern region of Mexico using principal component analysis (PCA). Atmósfera 20, 133-145.

Quintero-Durán R. (2008). Efectos de la aplicación de elementos menores en caña de azúcar en suelos del valle del río Cauca. Cali, Colombia, pp. 18-25. [online]. no20_2008_p18-26.pdf 10/09/2014. [ Links ]

Rein P., Turner P. and Mathias K. (2011). Good management practices manual for the cane sugar industry (final). PGBI House-International Finance Corporation. South Africa, Africa, 696 pp. [online]. .

Reyes H., Aguilar R.M., Aguirre R.J.R. and Trejo VI. (2005). Cambios en la cubierta vegetal y uso de suelo en el área del proyecto Pujal-Coy, San Luis Potosí, Mexico, 1973-2002. Boletín del Instituto de Geografía 59, 26-42. DOI: .

SAGARPA (2007). Programa Nacional de la Agroindustria de la Caña de Azúcar. 2007-2012. Secretaría de Agricultura, Ganadería, Desarrollo Rural, Pesca y Alimentación.Mexico City, 10 pp. [online]. 28/06/2015.

SAGARPA (2009a). Digitalización del campo cañero en México para alcanzar la agricultura de precisión de la caña de azúcar. Report. Ingenio Alianza Popular, S.A. de C.V Programa Nacional de la Agroindustria de la Caña de Azúcar. Servicio de Información Agroalimentaria y Pesquera. Colegio de Postgraduados. Secretaría de Agricultura, Ganadería, Desarrollo Rural, Pesca y Alimentación. San Luis Potosí, Mexico, 11 pp. [online]. 4/04/2014

SAGARPA (2009b). Digitalización del campo cañero en México para alcanzar la agricultura de precisión de la caña de azúcar. Report. Ingenio Plan de Ayala, S.A. de C.V. Programa Nacional de la Agroindustria de la Caña de Azúcar. Servicio de Información Agroalimentaria y Pesquera. Colegio de Postgraduados. Secretaría de Agricultura, Ganadería, Desarrollo Rural, Pesca y Alimentación. San Luis Potosí, Mexico, 10 pp. [online]. 4/04/2014.

SAGARPA (2009c). Digitalización del campo cañero en México para alcanzar la agricultura de precisión de la caña de azúcar. Report. Ingenio Pla de San Luis, S.A. de C.V. Programa Nacional de la Agroindustria de la Caña de Azúcar. Servicio de Información Agroalimentaria y Pesquera. Colegio de Postgraduados. Secretaría de Agricultura, Ganadería, Desarrollo Rural, Pesca y Alimentación. San Luis Potosí, Mexico. 10 pp. [online]. 4/04/2014.

Sánchez-Hernández R., Ramos-Reyes R., Geissen V., Mendoza-Palacios J. de D., de la Cruz-Lázaro E., Salcedo Pérez E. and Palma-López D. (2011). Contenido de carbonos en suelo con diferentes usos agropecuarios en el trópico mexicano. Terra Latinoamericana 29, 211-219.

Sandhu H.S., Gilbert R.A., Kingston G., Subiros J.F., Morgan K., Rice R.W. and Davis L. (2013). Effects of sugarcane harvest method on microclimate in Florida and Costa Rica. Agr. Forest Meteorol. 177, 101-109. DOI:10.1016/j.agrformet.2013.04.011.

Scheele E.A. and Kruger D. (2006). Global anthropogenic methane and nitrous oxide emissions. The Energy Journal, Multi-Greenhouse Gas Mitigation and Climate Policy Special Issue 27, 33-44.

Scherr J. and Staphit S. (2009). Mitigating climate change. Through food and land use. Worldwatch Report 179. Ecoaagriculture Partners and Worldwatch Institute. Washington, D.C., USA, 48 pp.

SEGOB (2014). Programa Nacional de la Agroindustria de la Caña de Azúcar 2014-2018. Diario Oficial de la Federación. 2 de mayo de 2014. Secretaría de Gobernación. Mexico City [online]. 21/04/2015.

SEMARNAT (2009). Impacto del cambio climático en las tierras y sus características. Report. Secretaría del Medio Ambiente y Recursos Naturales. Mexico City, 104 pp.

SEMARNAT (2013). Informe de la Situación del Medio Ambiente en México. Compendio de Estadísticas Ambientales. Indicadores Clave y de Desempeño Ambiental. Report. Secretaría del Medio Ambiente y Recursos Naturales. Mexico City, 20 pp.

SENASICA (2010). Boletín climático fitosanitario. Boletín No. 1. Climatología fitosanitaria. Sistema Nacional de Vigilancia Epidemiológica Fitosanitaria. Servicio Nacional de Sanidad, Inocuidad, y Calidad Agroalimentaria. San Luis Potosí, Mexico, 3 pp. [online]. 3/09/2015.

Sentíes-Herrera H.E., Gómez-Merino F.C., Valdez-Balero A., Silva-Rojas H.V. and Trejo-Téllez L.I. (2014). The agro-Industrial sugarcane system in Mexico: Current status, challenges and opportunities. J. Agr. Sci. 6, 26-54.

SIACON (2010). Sistema de información agroalimentaria de consulta. Sistema de Información Agroalimentaria y Pesquera (SIAP) [online]. 13/08/2014.

Swenson S. and Wahr J. (2005). Estimating Large-Scale Precipitation Minus Evapotranspiration from GRACE Satellite Gravity Measurements. J. Hydrometeorol. 7, 252-270. DOI: .

Tukimat N.N.A., Harun S. and Shahid S. (2012). Comparison of different methods in estimating potential evapotranspiration at muda irrigation scheme of Malaysia. J. Agr. Rural Dev. Trop. 113, 77-85.

UNFCCC (2014). Greenhouse gas inventory data - Detailed data by Party. United Nations Framework Convention on Climate Change. [online]. 12/06/2014.

Vilaboa A.I. and Barroso L.A. (2013). Contaminación ambiental por quema de caña de azúcar. Un estudio exploratorio en la región central del estado de Veracruz. Congress Memories. Think Green 2013: Crecimiento verde, retos y oportunidades para México. Veracruz, Veracruz, Mexico. 7 al 9 de octubre, 2013 [online]. 17/06/2015.

Ir a la Página