Simulated dynamics of net primary productivity (NPP) for outdoor livestock feeding coefficients driven by climate change scenarios in México

 

A. I. MONTERROSO RIVAS, J. D. GÓMEZ DÍAZ, M. L. TOLEDO MEDRANO, J. A. TINOCO RUEDA
Departamento de Suelos, Universidad Autónoma Chapingo,
Km. 38.5 Carretera México–Texcoco, Chapingo, 56230, Estado de México.
Corresponding autor: A.I. Monterroso Rivas; e–mail: aimrivas@correo.chapingo.mx

C. CONDE ÁLVAREZ and C. GAY GARCÍA
Centro de Ciencias de la Atmósfera, Universidad Nacional Autónoma de México,
Circuito Exterior, Ciudad Universitaria, México, D. F., 04510, México

 

Received May 18, 2009; accepted May 13, 2010

 

RESUMEN

En este trabajo se aplicó el concepto de productividad primaria neta (NPP, por sus siglas en inglés) como una forma de estimar la capacidad que tienen los ecosistemas de producir materia seca (DM, por sus siglas en inglés) que puede estar disponible para que el ganado cubra sus requerimientos de forraje. El método permite simular bajo condiciones de cambio climático el impacto posible que sobre la NPP y la materia seca se podrán observar para el país a un horizonte de tiempo dado. El concepto además fue aplicado sobre los coeficientes de agostadero actuales y bajo escenarios de cambio climático, permitiendo observar los posibles cambios en superficie requerida para la alimentación sostenible del ganado. Así, se encontró que en México se tienen valores de la NPP que van desde 0 y hasta 50 000 kgDM/ha/año. La DM, por su parte, se encuentra en el rango que va de 0 hasta 25 000 kgDM/ha/año. Los coeficientes de agostadero para el escenario actual no cambian considerablemente con relación a los presentados por COTECOCA, el organismo oficial encargado de su determinación. Bajo los escenarios de cambio climático se encontró que los estados de Baja California, Baja California Sur, Coahuila, Colima, Jalisco, Nuevo León, Puebla, Querétaro, San Luis Potosí, Sonora, Tamaulipas, Veracruz y Yucatán podrán ser los más impactados por posibles cambios en las condiciones actuales y verán modificados sus coeficientes de agostadero en el futuro. A partir de lo anterior se discuten las implicaciones que se podrán observar y se analizan las alternativas propuestas por el gobierno federal para el sector ganadero del país.

 

ABSTRACT

In this paper the concept of Net Primary Productivity (NPP) is used as a way to estimate the capacity of the ecosystem to produce dry matter which may be available for livestock to meet the forage requirements. The method allows the simulation of the possible impact on NPP and dry matter (DM), under climate change conditions observable for the country in a given time horizon. The concept was also used for current coefficients of rangeland and under current climate change scenarios, thus allowing to observe possible changes on the required surface for a sustainable cattle nutrition. It was found that México has NPP values ranging from 0 to 50 000 kgDM / ha / year. Dry matter is in the range that goes from 0 to 25 000 kgDM / ha / year. The coefficients of rangeland for the current scenario do not change significantly compared to those presented by COTECOCA, the official body determining such values. The states that present a greater impact due to current conditions under climate change scenarios were Baja California, Baja California Sur, Coahuila, Colima, Jalisco, Nuevo León, Puebla, Querétaro, San Luis Potosí, Sonora, Tamaulipas, Veracruz and Yucatán; this group will also have changes in its coefficients of rangeland in the future. Considering all of the above we discuss the implications that can be observed and we analyze the alternatives proposed by the federal government for the country's livestock sector.

Keywords: dry matter, rangeland, cattle nutrition, scenarios, sustainability.

 

1. Introduction

Between 14 and 22% of global greenhouse gases (CO₂ equivalent) come from livestock activity (FAO, 2006); specifically it is accountable for 9% of the CO₂ emissions (mainly due to deforestation), 37% of methane emissions (due to digestion of ruminants) and 65% of nitrous oxide (from manure). Livestock takes place in almost any climate and area of the world, from intensive production systems to those where the animals graze freely. However, it is estimated that livestock, especially meat production worldwide, will increase by double by 2050 (FAO, 2006), thus increasing the environmental impact of this activity, as well as food demand –such as in grassland areas. In this sense, the Intergovernmental Panel on Climate Change (IPCC, 2001) found that the combination of increased CO₂ concentrations and temperature changes along with precipitation will have significant impacts on grassland areas, with increased production in humid regions but decreases in arid and semiarid regions. This situation was confirmed years later in the Fourth Assessment Report (IPCC, 2007).

Grasslands, and generally any other plant, are sensitive to the presence of CO₂ and climate change: impacts on stability and resilience of plant communities can be observed (Mitchell and Cillag, 2001). It has also been confirmed that changes in forage quality and feeding behavior will occur (Easter ling et al., 2007). Furthermore it was found that species not palatable for cattle breeding may spread, which may have effects of less nutritional value of these pastures for livestock (DEFRA, 2000).

It has also been identified that climate variability may directly affect livestock, particularly under El Niño and the Southern Oscillation (ENSO) where changes in vegetation's behavior and water availability take place (Gitay et al., 2001). In semiarid regions there is a risk of severe degeneration of vegetation, which may lead to changes in land degradation, and along with less rain, loss of grazing lands and pastures (Zheng et al., 2002).

Plant output and vegetal species composition in grasslands is highly correlated to the amount of precipitation (Knapp and Smith, 2001). It is expected for the amount of rain in different regions in México to decrease under climate change conditions (Conde et al., 2008). In addition to the latter, it has also been stated that the most vulnerable vegetation types in the country are grasslands, scrubs and oak forests, which have the highest rates of change projected for 2050. Moreover, between 53 and 62% of vegetal communities will be exposed to climate conditions different to the current ones (INE–SEMARNAT, 2006). Out of 198 million hectares of national surface, approximately 16% is for agricultural use, 23% are forests and rainforests and the remaining 61% is rangeland area (PND, 2007), which implies that the national livestock surface is just over one hundred million hectares (COTECOCA, 2002). In addition to the threat of climate change on vegetation, deterioration due to land use caused by the introduction and expansion of livestock must be considered as well (INE–SEMARNAT, 2006). Therefore, the following objectives were proposed for this study 1) evaluate vegetation's capacity to remain as a food source for livestock, 2) evaluate the net primary productivity and rates of rangeland under present and climate change conditions as a measure to estimate the vulnerability of the sector, and 3) promote awareness on natural capacity and the impact on livestock that can be observed due to climate change.

 

2. Methods

2.1 Data

The geographical distribution of monthly and annual temperature and precipitation as climate variables for the base scenario in the period from 1950–2000 was taken from Conde et al. (2008). Soil information and maps (INEGI, 1998) as well as dominant vegetation distribution in the country (SEMARNAT–UNAM, 2001) was also used for this study. Areas influenced by climate were drawn from a cartographic map overlapping of annual isohyets and isotherms, as suggested by Gómez et al. (2008). More than 2000 areas influenced by climate were drawn in all the country, each with their respective average values of temperature and precipitation per month and year. Soil information and vegetation distribution was considered as a base reference to estimate humidity balances.

2.2 Estimation of evaporation and evapotranspiration

The balance of humidity and evapotranspiration were obtained using the method of Thornthwaite (1957) amended version III (Monterroso and Gómez, 2003; Estrada–Berg et al., 2008). The first step in calculating the balance is to define soil's capacity to retain humidity, considering soil type and vegetal cover (CAP). The following step is to calculate potential evapotranspiration (PET), which is the maximum amount of water that can be lost by a continuous layer of vegetation covering the entire area, when an unlimited amount of water is supplied to the soil (Ortiz, 1987). The PET estimate depends on the monthly mean temperature and annual heat index, according to the following:

Where i_j is the heat index for the month j and tm_j is the mean temperature of month j in °C. The annual heat index is obtained according to:

The monthly and annual heat index is dimensionless. In getting the value of (2) the daily value of evaporation (EV) in millimeters can be estimated according to formula (3):

For monthly temperature values higher than 26.5 °C, the EV value must be obtained from tables. From (3) the value of potential evaporation (PET) can be calculated by multiplying the evaporation by a latitude correction factor (cf), taken from the tables.

Next step to estimate the hydric balance is to consider the soil's humidity storage capacity (CAP) as well as the potential loss of accumulated water (APWL), which depends on the difference between precipitation and potential evapotranspiration:

Dunne and Leopold (1978) define storage as the amount of water that a given unit of soil is capable to retain and that is available for plants, depending on the texture and depth of the soil. The humidity deficit (DEF) occurs when precipitation is less than potential evapotranspiration and can be defined as the difference between the amount of water that can potentially be lost in the atmosphere and the amount of water actually lost (Estrada–Berg et al., 2008):

The last step consists in obtaining the real evapotranspiration (ETR), water surplus and runoff, as presented below:

Finally, the water balance estimated three indexes: humidity, aridity and precipitation. An approach to a climate classification using the tables and indexes concluded in eight groups of climates: perhumid, humid, sub–humid humid, sub–humid dry, semi–arid humid, semi–arid, semi–arid dry and arid.

2.3 Net primary productivity (NPP)

Plants, like any other organism, must use energy in biomass production, growth and reproduction. This energy is spent in the respiration process, which in turn releases energy. The energy remaining after respiration and stored as organic matter is the net primary productivity (NPP) (Smith and Smith, 2001). There are two groups of methods to calculate NPP; the indirect and the direct. Direct methods are more accurate since they use a sampling of multiple crops from the vegetation directly from the field and then it is analyzed in the laboratory to obtain the net primary productivity. However, when working in large areas, direct methods are very costly and very time–consuming. To solve this, there are indirect methods or empirical methods, where equations are applied according to climate and biological variables in order to estimate NPP.

Within this group, one of the most frequently used equations is the one proposed by Rosenzweig (1968) and it has been tested in México (Serrano, 1987, Ritter et al., 1999, Alba et al., 2003). Rosenzweig designed a logarithmic equation where the vegetation's real evapotranspiration (RET) can be used to predict NPP. The principle is that RET simultaneously expresses water availability and sun energy, essential and indispensable factors in photosynthesis. The equation as stated and used in this study is expressed as follows

In this equation, NPP refers to the net primary productivity (gMS/m²/year); RET is the real evapotranspiration in mm / year and DM is dry matter. To apply the formula, RET estimate of the hydric balance previously described was used.

2.4 Dry matter and herbaceous estimation

Considering plants use NPP to produce leaves, thicken stems, branches, roots and to form flowers and fruits, a percentage of the NPP corresponding to the amount of herbaceous production for each of the vegetation's communities was defined. We analyzed the types of vegetation in the National Forest Inventory (SEMARNAT–UNAM, 2001) and a percentage of herbaceous production was assigned to each vegetation community according to what is presented in Table I. Grassland formations are those with higher incidence of herbs (up to 80% of the NPP becomes grass) and on the contrary, forests provide less herbaceous production due to their structure and composition.

A use factor was later assigned, taking into consideration that animals do not eat all the grass produced in each zone, and we are additionally including species that are not eaten by animals. The use factor was set at 0.5, a normal value for ecosystems (Novillo et al., 2006). Thus, the amount of food available for animal consumption and cattle and cynegetic load is really below to the total NPP, resulting in smaller initial values for dry usable matter.

2.5 Outdoor livestock feeding coefficients

Rangeland is defined as a land with a capacity to produce forage for livestock and wild animals. The rangeland coefficient (RC) is the animal unit– area ratio adequate to maintain a farm production economically and productively permanent, without damaging the natural resource. In other words it is the necessary surface to maintain one animal unit, at its maximum and permanent production in a given area, compatible with the preservation of natural resources. In México, the Technical Advisory Committee for Determining Rangeland Coefficients (COTECOCA) has been the institution responsible for estimating the national RC and has proposed to estimate the RC according to the following function:

where RC refers to rangeland coefficient (AU / ha / year), DM is dry matter and AU is one animal unit. It is estimated that an AU eats as much as 3% of its body weight. COTECOCA (1976) states that the daily consumption of one animal unit weighing approximately 450 kg, is 13.5 kg in average, resulting in an annual consumption of 4927.5 kg. It is important to point out that one animal unit stands for an adult cow along with its calf of less than seven months. Based on (12), consumption can be standardized at 4.925 kg / DM per animal unit per year; thus, to estimate it, it is necessary to only obtain the dry matter production value in each of the different ecosystems of the country, taken from the previous chapter

2.6 Climate change scenarios

The general atmospheric circulation models (GCM) used in the study were ECHAM–5, HADGEM–1 and GFDL–CM–2.0 which will be named for simplicity as Echam, Hadley and GFDL; the socio–economic scenarios tested were A2 and B2 (Nakicenovic et al., 2000) for 2030 and 2050 time horizons. The outputs of these models have a spatial resolution of 2.5° × 2.5° and the change rates were taken from Conde et al. (2008). The reasons for change in each model were applied on the base maps of climate and water balance previously described and hence the new cartography for assessed climate change scenarios could be obtained

 

3. Results

3.1 Net primary productivity (NPP)

The NPP represents the net flow of carbon from the atmosphere into green plants per unit of area and time. In México, this parameter is in the range that goes from 0 to 50 000 kilograms of dry matter per hectare per year (kgDM / ha / year). The variation is largely related to the ecological conditions of each region, the amount of energy received, soil vegetation cover and the amount of precipitation. Figure 1 shows the distribution of this parameter in México for the base scenario and the general circulation models for B2 scenario and time horizon to 2050.

3.2 Dry matter estimation

Figure 2 shows the estimated range of DM for the country and the area occupied. About 38% of the land area has a DM production ranging from 0–200 kgDM / ha / year, mainly located in the northern region and in small areas in the central part of the country. It is followed by the range of 500–1000 kgDM / ha / year with approximately 15% of the area. In 76% of the national area it is produced less than 1000 kgDM / ha / year, indicating that three quarters of the country's natural pasture availability is limited. Areas that are located in this range correspond to the dry and semi–arid regions of the country, representing about 50% of national surface (Monterroso and Gómez, 2003), and also to temperate forest areas of closed canopy that limits the production of herbaceous plants. The greater production of DM from herbaceous is generated in grasslands, and due to the ecological character of each site, the range can go from 5000 to 50 000 kgDM / ha / year. The plant communities of tropical and temperate forests have a reduced herbaceous layer and the production of DM in these areas has no higher values, since the canopy is closed and allowed no sunlight, the possibility that the herbaceous layer develops optimally is limited. Thus, the range in which lie these forests, in terms of production of herbaceous plants, goes from 300 to 5000 kgDM / ha / year. The bushes and vegetation of arid and semi–arid areas provide a production of DM from grass in ranges that can go from 0 to 1500 kgDM / ha / year, depending on prevailing rainfall, temperature and solar radiation received, and the type of soil.

3.3 Outdoor livestock feeding coefficients

It was found that the rangeland coefficient (RC) that dominates the country and the base scenario is greater than 20 ha per animal unit, covering about 61% of the total area (Table II). These areas are associated with drier zones and soil as limiting factor for the development of herbaceous plants, that correspond in part to the very arid climate zones estimated by Monterroso and Gómez (2003), plus semiarid areas with dry climate and also soil as limiting factor. On the contrary, the RC of less than 5 hectares covers 15.5% of the national surface and is distributed in small strips, associated mainly to grassland areas in the south of the country where the climate is hot and humid.

The RC bigger than 100 ha per animal unit, although it is distributed throughout most part of the country, is located mainly in the north of the country, being Baja California Sur, Coahuila and Sonora the states with larger areas, associated mainly to arid and very dry areas where the grass production is very limited. The production of DM by herbaceous plants in this range is less than 90 kgDM / ha / year. Greater detail of these results is presented in Annex 2.

3.4 Climate change scenarios

The rates of change obtained from each of the general circulation models differ in the variations of temperature and precipitation throughout the year for the different regions of the country. All models predict an increase in temperature, although different for each of the regions. In general, the GFDL presents overall smaller increments but generalized throughout the country and every month. The biggest increase is in Echam, and Hadley model presents intermediate values. This relationship remains across the two time horizons and in both scenarios, differing in each. The increase is higher for 2050 and the A2 scenario. The estimates for precipitation differ between models, while in general the GFDL rates of change have increases in most months, specially in the summer when most of the rain happens, this provision applies to both scenarios and time horizons covered. Instead, Hadley and Echam models presented a trend of lower rainfall, which is accentuated in the A2 scenario for 2050. These conditions have a direct effect on the analysis of NPP, DM and outdoor livestock feeding coefficients. Figure 3 shows the percentage of changes of NPP compared to the baseline scenario.

The changes are more pronounced in the A2 scenario and for 2050, although for each scenario (A2 and B2) patterns of change show great similarity with variations only in the time horizons.

In the A2 scenario for 2030 and 2050 the pattern of variation in the area of the various classes of NPP is very similar. An increase in the bigger class that goes from 40 001–50 000 kgDM / ha / year is observed. This is associated with warm humid areas where no water deficit is expected because they are located in the Gulf of México slope with rain all year. With the increasing of temperature and a proper moisture balance, these zones may enhance their productivity. For all three models an increase from 1.5 to 2.0% is estimated.

In the B2 scenario the impact on productive capacity is lower than in the A2 and follows a similar pattern of involvement in different kinds of DM production, with some variations between classes.

These impacts are transferred to the DM available for grazing, and considering that the NPP is adjusted by a coefficient associated with the type of vegetation, and because it is expected that this remain, the change will be similar for each class, each model scenarios and the time periods considered. This applies to the outdoor livestock feeding coefficients and, for what was discussed in the section, to the variation that they will have.

Figure 4 shows the rangeland coefficients for the base scenario and the socio–economic B2 scenario for 2050 model.

 

4. Discussion and conclusions

Natural ecosystems in México grow in different environmental conditions that characterize them and that have an influence on the production potential in them. The geographic location of each biome provides distinction, especially in primary productivity. In this sense, the method used to obtain the NPP allowed to correlate the climate and environmental variables for future NPP simulation under climate change scenarios. By using this method we could confirm that the estimates that only consider temperature overestimate the NPP modeling (Ritter et al., 1999). Including hydric balance in modeling allowed a better NPP estimate.

Estimating NPP involves considering different environmental and ecological aspects that at a national level may be overestimated. To apply a preliminary hydric balance helped to better estimate NPP and reduce spatial errors implied in the formula's application. It is possible that Rosenzweiz's formula suggests high levels of productivity, so it is recommended to do more studies in order to corroborate its validity and especially its real meaning for more regional scales. In this study it was decided to use the formula given its simplicity of implementation at a national level, although it has not been proven in national studies. It allowed getting more detailed information, with a significant contribution to this work.

In order to obtain more realistic values in estimating the rangeland coefficients, we estimated available dry matter for cattle from NPP values, using herbaceous production factors for plant communities as well as a use factor. This allowed reducing the possible errors in the overestimation of NPP, since herbaceous production factors (Table I) assigned by experts and the use factor (50% of herbaceous production) are reported as normal values for the country's ecosystems.

In this sense, different potential for herbaceous production for cattle consumption was found in the country, where over half of the national area (60%) had low dry matter production that can be used by livestock, mainly located in areas where rainfall is less than 800 millimeters per year. Moreover, the region where there is greater availability of herbaceous plants is found in the states of Tabasco, Veracruz and Chiapas, among some others, however, these areas represent just less than 5% of the national area.

The lowest dry matter production, i.e. for the range from 0 to 200 kg of dry matter per hectare per year, represent about 60% of the national area and are related to restrictive conditions, such as low rainfall and low soil fertility. If the range is widen –for example from 0 to 300 kgDM / year / ha– the national area increases to 65%. For the case of a range from 0 to 500 kgDM / year / ha, the surface increases to 75%. This suggests that in three quarters of the country it can be fed less than one animal unit per hectare per year naturally, living from extensive grazing, i.e. without shed conditions and in a sustainable manner.

Pastures are where better opportunities to produce food for livestock can take place. The different types of grass can be distinguished by the availability of humidity and temperature, being the humid part of the country the one with the largest natural potential. The tropical region presents the greatest potential for livestock, as it has better conditions for plant growth, mainly linked to high solar radiation and high rainfall.

This affects the definition of rangeland coefficients. The table in Annex 1 compares the values presented by COTECOCA (weighted based on their maximum and minimum estimates) with the results obtained in this study. It is noticed that for eight states the values are practically similar (±10%); for 18 states the required area per animal unit is greater, while in six states the estimates of this study are less than the ones from COTECOCA. It is worth noting the differences found in the states of Chiapas, Colima, Hidalgo, Oaxaca, San Luis Potosi, Sinaloa and Veracruz: information reported by COTECOCA does not exceed 10 hectares required per animal unit, whereas in this study we obtained up to 28ha. Such is the case of Chiapas, where the increase in the area required is of about 1481% (going from 1.8 to 28 ha required per animal unit).

While analyzing the causes for these differences we concluded that they are due to both the cartography used in each study, and the observed changes in land use.

Rangeland coefficients defined by COTECA were for the years 1972–1986, therefore the used cartography fulfilled standards set in those years; however, 20 years have gone by since its publication. In this study, we used cartography from 2001, and the classification and nomenclature of land use coverage are both more detailed and different from those mentioned above. On the other hand we have the changes in land use in the country's ecosystems. The cartography used here was sensitive to such changes occurred since then; they were definitely considered in their spatial representation. The states with greater differences in rangeland coefficients mentioned above –Chiapas for example– agree in having areas covered with forests and rainforests within their boundaries, which are ecosystems that have had major changes in land use since the 80s.

The most favorable rangeland coefficients are located in the Gulf of México, Yucatán Peninsula and South Pacific (particularly in Chiapas), while the less optimal for livestock development are located in the Baja California Peninsula, to the northwest and northeast of the country. This should be regarded as a natural capacity for animal support, because under shed conditions and with greater financial resources, it is possible to increase the coefficients and improve natural capacity.

In applying possible future climate change scenarios, the driest regions of the country increase their vulnerability, and categorically the three models agree on temperature rises, primarily resulting in the reduction of humidity in the soil and affecting the development of herbaceous stratum in ecosystems, thus affecting the capacity of animal support as well. This region and its livestock productivity must be addressed under very specific policies that focus on reducing vulnerability, investing financial resources to compensate for what nature will no longer provide in a natural way. The states that are more vulnerable to these changes are Sonora, Chihuahua and Coahuila.

In the Yucatán Peninsula and specifically in the state of Yucatán, about 13% of its surface has an RC of more than 20 ha / AU / year in the base scenario; change scenarios suggest an increase of 5% up to 18%. In Quintana Roo State, the RC in the range 20–50 ha / AU / year for the base scenario covers 29% and models agree in the possibility of a rise, reaching 55%. In the states of Guerrero, Hidalgo, México, Michoacán, Nayarit and Oaxaca is noted that under change scenarios it is very likely that current conditions are maintained or even slightly improved.

At this point it is important to note that in México the rangeland coefficients determined by COTECOCA are not dynamic and are calculated for natural conditions, i.e. without considering the disturbance caused by misuse or improvements in the conditions of the sites evaluated locally. They are legal and current valid to determine the size of the small cattle property. However, for current management purposes, values here reported are only demonstrative and we point out that they must be adjusted to local and existing vegetation, climate, soil and animal species conditions using the resources of the sites being considered, among other factors. Therefore, climate change scenarios presented here should be taken with caution and its usefulness is for demonstration purposes only in case future research lines are suggested as well as for current schemes in planning and adapting foreseen changes.

In this sense the Special Program on Climate Change (SPCC, 2008) promoted by the federal executive government and that emerged from the National Strategy for Climate Change (NSCC, 2007) shows a strong inclination to focus efforts on reducing emissions through the development of mitigation projects. Although we recognize the efforts described above, we believe that vulnerability must be reduced and to develop better adaptation capacities in the livestock sector in the country.

Specifically in the Program for livestock, two objectives are proposed: reduce livestock vulnerability towards climate change impact and deepen the knowledge on the impacts on and vulnerability of livestock sector towards variability and climate change. As for the second objective we recognize the lack of studies on the subject in the country, so we can only propose research lines, while for the first objective we propose specific guidelines that can be improved.

In this sense, we applaud the federal government's willingness to address climate change impacts in this sector. However, the proposals that focus on reducing vulnerability are perceived as short, since it is necessary to include a biophysical and social dimensions for its evaluation, as well as to determine the adaptation capacity of the sector (Carter, 2007), and to know why and to what it is vulnerable to. It is expected that the results presented here assist in the formulation of public policies, pointing out a starting point for possible impacts on livestock and helping to identify why it would be in some degree of vulnerability. Different climate change models suggest major changes implemented in the number of animals that can be supported in natural ecosystems. The next step, as set out in the SPCC will be to approach producers (especially with low–income) and develop detailed studies that allow to reduce vulnerability and enhance adaptation capacity of the sector.

 

References

Alba López M. P., M. González Espinoza, N. Ramírez Marcial and M. A. Castillo Santiago, 2003. Determinantes de la distribución de Pinus spp., en la Altiplanicie Central de Chiapas, México. Bol. Soc. Bot. Méx. 73, 7–15.         [ Links ]

Carter T. R., R. N. Jones, X. Lu, S. Bhadwal, C. Conde, L. O. Mearns, B. C. O'Neill, M. D. A. Rounsevell and M. B. Zurek, 2007. New Assessment methods and the characterization of future conditions. Climate Change 2007: Impacts, Adaptation and Vulnerability. Contribution of Working Group II to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change (M.L. Parry, O.F.Canziani, J.P. Palutikof, P.J. van der Linden and C.E. Hanson, Eds.), Cambridge University Press, Cambridge, UK, 133–171.         [ Links ]

Conde C., B. Martínez, O. Sánchez, F. Estrada, A. Fernández, J. Zavala and C. Gay, 2008. Escenarios de Cambio Climático (2030 y 2050) para México y Centro América. Temperatura y Precipitación. Available at: http://www.atmosfera.unam.mx/cclimat/index.php?option=com_content & view=article & id=61 & Itemid=74.         [ Links ]

COTECOCA, 1976. Coeficientes de Agostadero Nacionales. Comisión Técnico Consultiva de Coeficientes de Agostadero (COTECOCA), Secretaría de Agricultura y Ganadería, con base en: Cotecoca, SARH, Memorias de Coeficientes de Agostadero, años 1972–1986, México.         [ Links ]

COTECOCA, 2002. Comisión Técnico Consultiva de Coeficientes de Agostadero, SAGARPA, con base en: COTECOCA, SARH. 1972–1981. Revisado en 2002.         [ Links ]

DEFRA, 2000. Impact of climate change on grasslands and livestock. In: Climate Change and Agriculture in the United Kingdom. Department for Environment, Food and Rural Affairs, London, pp. 43–56.         [ Links ]

Dunne T. and L. Leopold, 1978. Water in environmental planning. W. H. Freeman. USA, 818 pp.         [ Links ]

Easterling W. E., P. K. Aggarwal, P. Batima, K. M. Brander, L. Erda, S. M. Howden, A. Kirilenko, J. Morton, J.–F. Soussana, J. Schmidhuber and F. N. Tubiello, 2007. Food, fibre and forest products. Climate Change 2007: Impacts, Adaptation and Vulnerability. Contribution of Working Group II to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change (M.L. Parry, O.F. Canziani, J.P. Palutikof, P.J. van der Linden and C.E. Hanson, Eds.), Cambridge University Press, Cambridge, UK, 273–313.         [ Links ]

ENACC, 2007. Estrategia Nacional de Cambio Climático. Comisión Intersecretarial de Cambio Climático. Gobierno de la República. México. 157 pp.         [ Links ]

Estrada–Berg, J. B., J. D. Gómez, J. W. Estrada–Berg, A. I. Monterroso and J. A. Tinoco. 2008. El cambio climático en la República Mexicana, aplicaciones metodológicas en el manejo de suelos. Universidad Autónoma Chapingo. México. 143 pp.         [ Links ]

FAO, 2006. Livestock's long shadow, environmental issues and options. Food and Agricultural Organization of the United Nations. Rome. 375 pp.         [ Links ]

Gitay H., S. Brown,W. E. Easterling, B. Jallow, J. Antle, M. Apps, R. Beamish, C. Cerri and Coauthors, 2001. Ecosystems and their services. Climate Change 2001: Impacts, Adaptation and Vulnerability to Climate Change. Contribution of Working Group II to the Third Assessment Report of the Intergovernmental Panel on Climate Change (J.J. McCarthy, O.F. Canziani, N.A. Leary, D.J. Dokken and K.S. White, Eds.), Cambridge University Press, Cambridge, 236–342.         [ Links ]

Gómez J. D., J. Etchevers, A. Monterroso, C. Gay, J. Campo and M. Martínez, 2008. Spatial estimation of mean temperature and precipitation in areas of scarce meteorological information. Atmósfera 21, 35–56.         [ Links ]

INE–SEMARNAT, 2006. México, Tercera Comunicación Nacional ante la Convención Marco de las Naciones Unidas sobre Cambio Climático. Instituto Nacional de Ecología and Secretaría de Medio Ambiente y Recursos Naturales. México. 211 pp.         [ Links ]

INEGI, 1998. Mapa de suelos dominantes de la República Mexicana escala 1:250000. Instituto Nacional de Estadística y Geografía. Aguascalientes, México. (122 charts).         [ Links ]

IPCC, 2001. Climate Change 2001: Impacts, Adaptation and Vulnerability. Contribution of Working Group II to the Third Assessment Report of the Intergovernmental Panel on Climate Change (J.J. McCarthy, O.F. Canziani, N.A. Leary, D.J. Dokken and K.S. White, Eds.), Cambridge University Press, Cambridge, UK. 236–342.         [ Links ]

IPCC, 2007. Climate Change 2007: Impacts, Adaptation and Vulnerability. Contribution of Working Group II to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change (M.L. Parry, O.F. Canziani, J.P. Palutikof, P.J. van der Linden and C.E. Hanson, Eds.), Cambridge University Press, Cambridge, UK, 273–313.         [ Links ]

Knapp A. K. and M. D. Smith, 2001. Variation among biomes in temporal dynamics of aboveground primary production. Science 291, 481–484.         [ Links ]

Leith H. and E. Box, 1972. Evapotranspiration and primary production, C.W. Thornthwaite Memorial Model. Papers on Selected Topics in Climatology 2, 37–44.         [ Links ]

Mitchell S. W. and F. Csillag, 2001. Assessing the stability and uncertainty of predicted vegetation growth under climatic variability: northern mixed grass prairie. Ecol. Model. 139, 101–121.         [ Links ]

Monterroso A. I. and J. D. Gómez, 2003. Escenarios climatológicos de la República Mexicana ante el cambio climático. Tomo 1. Comisión Nacional de Zonas Áridas and Universidad Autónoma Chapingo. Chapingo, Méx. 170 pp.         [ Links ]

Novillo Camacho C., J. Beneitez and R. Romero Calcerrada, 2006. Modelo de estimación de la productividad primaria neta para determinar actividades forestales sostenibles a escala regional. En: El acceso a la información espacial y las nuevas tecnologías geográficas, (M.T. Camacho Olmedo, J.A. Cañete Pérez, J.J. Lara Valle, Eds.) Actas del XII Congreso de tecnologías de la información geográfica. Universidad de Granada, España. pp. 339–351.         [ Links ]

Ortíz Solorio C. A., 1987. Elementos de agrometeorología cuantitativa con aplicaciones en la República Mexicana. Departamento de Suelos. Universidad Autónoma Chapingo. Chapingo, México. 327 pp.         [ Links ]

PECC, 2008. Programa Especial de Cambio Climático 2008–2012. Comisión Intersecretarial de Cambio Climático. Gobierno de la República. México. 213 pp.         [ Links ]

PND, 2007. Plan Nacional de Desarrollo 2007–2012. Gobierno de la República. México. 321 pp.         [ Links ]

Ritter Ortiz W., S. A. Guzmán Ruíz and Estrada A. Betancourt, 1999. Un método gráfico para la regionalización agroecológica de recursos naturales. Investigaciones Geográficas 38, 57–68.         [ Links ]

Rosenzweig M. L., 1968. Net primary productivity of terrestrial environments: predictions from climatological data. American Naturalist 102, 67–84.         [ Links ]

SEMARNAT–UNAM, 2001. Mapa de uso de suelo y vegetación 2001. Secretaría de Medio Ambiente y Recursos Naturales and Instituto de Geografía, Universidad Nacional Autónoma de México. México.         [ Links ]

Serrano V., 1987. Las comunidades de roedores desertícolas del Bolsón de Mapimí, Durango. Acta Zool. Méx. (ns) 20, 1–22.         [ Links ]

Smith R. L. and T. M. Smith, 2001. Ecology and field biology. Addison Wesley. 771 pp.         [ Links ]

Thornthwaite C. W. and J. R. Matter, 1957. Instructions and tables for computing potential evapotranspiration and the water balance. Climatology 10, 185–311.         [ Links ]

Zheng Y. Q., G. Yu, Y. F. Qian, M. Miao, X. Zeng and H. Liu, 2002. Simulations of regional climatic effects of vegetation change in China. Q. J. Roy. Meteor. Soc., 128, 2089–2114.         [ Links ]