A comparison of two models with Landsat data for estimating above ground grassland biomass in Inner Mongolia,China

Two models, artificial neural network (ANN) and multiple linear regression (MLR), were developed to estimate typical grassland aboveground dry biomass in Xilingol River Basin, Inner Mongolia, China. The normalized difference vegetation index (NDVI) and topographic variables (elevation, aspect, and slope) were combined with atmospherically corrected reflectance from the Landsat ETM+ reflective bands as the candidate input variables for building both models. Seven variables (NDVI, aspect, and bands 1, 3, 4, 5 and 7) were selected by the ANN model (implemented in Statistica 6.0 neural network module), while six (elevation, NDVI, and bands 1, 3, 5 and 7) were picked to fit the MLR function after a stepwise analysis was executed between the candidate input variables and the above ground dry biomass. Both models achieved reasonable results with RMSEs ranging from 39.88% to 50.08%. The ANN model provided a more accurate estimation (RMSE_r = 39.88% for the training set, and RMSE_r = 42.36% for the testing set) than MLR (RMSE_r = 49.51% for the training, and RMSE_r = 53.20% for the testing). The final above ground dry biomass maps of the research area were produced based on the ANN and MLR models, generating the estimated mean values of 121 and 147 g/m², respectively.     

The importance of spatial autocorrelation,extent and resolution in predicting forest bird occurrence

Concerns about declines in forest biodiversity underscore the need for accurate estimates of the distribution and abundance of organisms at large scales and at resolutions that are fine enough to be appropriate for management. This paper addresses three major objectives: (i) to determine whether the resolution of typical air photo-derived forest inventory is sufficient for the accurate prediction of site occupancy by forest birds. We compared prediction success of habitat models using air photo variables to models with variables derived from finer resolution, ground-sampled vegetation plots. (ii) To test whether incorporating spatial autocorrelation into habitat models via autologistic regression increases prediction success. (iii) To determine whether landscape structure is an important factor in predicting bird distribution in forest-dominated landscapes. Models were tested locally (Greater Fundy Ecosystem [GFE]) using cross-validation, and regionally using an independent data set from an area located ca. 250 km to the northwest (Riley Brook [RB]). We found significant positive spatial autocorrelation in the residuals of at least one habitat model for 76% (16/21) of species examined. In these cases, the logistic regression assumption of spatially independent errors was violated. Logistic models that ignored spatial autocorrelation tended to overestimate habitat effects. Though overall prediction success was higher for autologistic models than logistic models in the GFE, the difference was only significantly improved for one species. Further, the inclusion of spatial covariates did little to improve model performance in the geographically discrete study area. For 62% (13/21) of species examined, landscape variables were significant predictors of forest bird occurrence even after statistically controlling for stand-level variability. However, broad spatial extents explained less variation than local factors. In the GFE, 76% (16/21) of air photo and 81% (17/21) of ground plot models were accurate enough to be of practical utility (AUC > 0.7). When applied to RB, both model types performed effectively for 55% (11/20) of the species examined. We did not detect an overall difference in prediction success between air photo and ground plot models in either study area. We conclude that air photo data are as effective as fine resolution vegetation data for predicting site occupancy for the majority of species in this study. These models will be of use to forest managers who are interested in mapping species distributions under various timber harvest scenarios, and to protected areas planners attempting to optimize reserve function.     

The derivation of species response curves with Gaussian logistic regression is sensitive to sampling intensity and curve characteristics

We investigated quantitatively the sensitivity of plant species response curves to sampling characteristics (number of plots, occurrence and frequency of species), along a simulated pH gradient. We defined 54 theoretical unimodal response curves, issued from combinations of six values for optimum (opt = 3, 4, …, 8), three values for tolerance (tol = 0.5, 1.0, and 1.5, sensu ter Braak and Looman [ter Braak, C.J.F., Looman, C.W.N., 1986. Weighted averaging, logistic regression and the Gaussian response model. Vegetatio 65, 3–11]), and three values for maximum probability of presence (p_max = 0.05, 0.20, and 0.50). For each of these 54 theoretical response curves, we built artificial binary data sets (presence/absence) to test the influence of species occurrence, frequency, or number of available plots. With real data extracted from EcoPlant, a phytoecological database for French forests [Gégout, J.-C., Coudun, Ch., Bailly, G., Jabiol, B., 2005. EcoPlant: a forest sites database linking floristic data with soil characteristics and climatic conditions. J. Veg. Sci. 16, 257–260], we compared the ecological response of 50 plant species to soil pH, based first on a small data set (100 randomly sampled plots), and then based on the whole data set available (3810 plots).     

Structural changes in lake functioning induced from nutrient loading and climate variability

Climate variability is increasingly recognized as an important regulatory factor, capable of influencing the structural properties of aquatic ecosystems. Lakes appear to be particularly sensitive to the ecological impacts of climate variability, and several long time series have shown a close coupling between climate, lake thermal properties and individual organism physiology, population abundance, community structure, and food web dynamics. Thus, understanding the complex interplay among meteorological forcing, hydrological variability, and ecosystem functioning is essential for improving the credibility of model-based water resources/fisheries management. Our objective herein is to examine the relative importance of the ecological mechanisms underlying plankton seasonal variability in Lake Washington, Washington State (USA), over a 35-year period (1964–1998). Our analysis is founded upon an intermediate complexity plankton model that is used to reproduce the limiting nutrient (phosphate)–phytoplankton–zooplankton–detritus (particulate phosphorus) dynamics in the lake. Model parameterization is based on a Bayesian calibration scheme that offers insights into the degree of information the data contain about model inputs and allows obtaining predictions along with uncertainty bounds for modeled output variables. The model accurately reproduces the key seasonal planktonic patterns in Lake Washington and provides realistic estimates of predictive uncertainty for water quality variables of environmental management interest. A principal component analysis of the annual estimates of the underlying ecological processes highlighted the significant role of the phosphorus recycling stemming from the zooplankton excretion on the planktonic food web variability. We also identified a moderately significant signature of the local climatic conditions (air temperature) on phytoplankton growth (r = 0.41), herbivorous grazing (r = 0.38), and detritus mineralization (r = 0.39). Our study seeks linkages with the conceptual food web model proposed by Hampton et al. [Hampton, S.E., Scheuerell, M.D., Schindler, D.E., 2006b. Coalescence in the Lake Washington story: interaction strengths in a planktonic food web. Limnol. Oceanogr. 51, 2042–2051.] to emphasize the “bottom-up” control of the Lake Washington plankton phenology. The posterior predictive distributions of the plankton model are also used to assess the exceedance frequency and confidence of compliance with total phosphorus (15 μg L⁻¹) and chlorophyll a (4 μg L⁻¹) threshold levels during the summer-stratified period in Lake Washington. Finally, we conclude by underscoring the importance of explicitly acknowledging the uncertainty in ecological forecasts to the management of freshwater ecosystems under a changing global environment.     

Modeling Species' Distributions to Improve Conservation in Semiurban Landscapes: Koala Case Study

Abstract:  Models of species' distributions are commonly used to inform landscape and conservation planning. In urban and semiurban landscapes, the distributions of species are determined by a combination of natural habitat and anthropogenic impacts. Understanding the spatial influence of these two processes is crucial for making spatially explicit decisions about conservation actions. We present a logistic regression model for the distribution of koalas (  Phascolarctos cinereus ) in a semiurban landscape in eastern Australia that explicitly separates the effect of natural habitat quality and anthropogenic impacts on koala distributions. We achieved this by comparing the predicted distributions from the model with what the predicted distributions would have been if anthropogenic variables were at their mean values. Similar approaches have relied on making predictions assuming anthropogenic variables are zero, which will be unreliable if the training data set does not include anthropogenic variables close to zero. Our approach is novel because it can be applied to landscapes where anthropogenic variables are never close to zero. Our model showed that, averaged across the study area, natural habitat was the main determinant of koala presence. At a local scale, however, anthropogenic impacts could be more important, with consequent implications for conservation planning. We demonstrated that this modeling approach, combined with the visual presentation of predictions as a map, provides important information for making decisions on how different conservation actions should be spatially allocated. This method is particularly useful for areas where wildlife and human populations exist in close proximity.     

Incorporating low-resolution historic species location data decreases performance of distribution models

Developing robust species distribution models is important as model outputs are increasingly being incorporated into conservation policy and management decisions. A largely overlooked component of model assessment and refinement is whether to include historic species occurrence data in distribution models to increase the data sample size. Data of different temporal provenance often differ in spatial accuracy and precision. We test the effect of inclusion of historic coarse-resolution occurrence data on distribution model outputs for 187 species of birds in Australian tropical savannas. Models using only recent (after 1990), fine-resolution data had significantly higher model performance scores measured with area under the receiver operating characteristic curve (AUC) than models incorporating both fine- and coarse-resolution data. The drop in AUC score is positively correlated with the total area predicted to be suitable for the species (R² = 0.163-0.187, depending on the environmental predictors in the model), as coarser data generally leads to greater predicted areas. The remaining unexplained variation is likely to be due to the covariate errors resulting from resolution mismatch between species records and environmental predictors. We conclude that decisions regarding data use in species distribution models must be conscious of the variation in predictions that mixed-scale datasets might cause.     

Modeling freshwater fish distributions using multiscale landscape data: A case study of six narrow range endemics

Species distribution models (SDMs) have become integral tools in scientific research and conservation planning. Despite progress in the assessment of various statistical models for use in SDMs, little has been done in way of evaluating appropriate ecological models. In this paper, we evaluate the multiscale filter framework as a suitable theoretical model for predicting freshwater fish distributions in the upper Green River system (Ohio River drainage), USA. The spatial distributions of six fishes with contrasting biogeographies were modeled using boosted regression trees and multiscale landscape data. Species biogeography did not appear to affect predictive performance and all models performed well statistically with receiver operating characteristic area under the curve (AUC) ranging from 0.87 to 0.98. Predictive maps show accurate estimations of ranges for five of six species based on historical collections. The relative influence of each type of environmental feature and spatial scale varied markedly with between species. A hierarchical effect was detected for narrowly distributed species. These species were highly influenced by soil composition at larger spatial scales and land use/land cover (LULC) patterns at more proximal scales. Conversely, LULC pattern was the most influential feature for widely distributed at all spatial scales. Using multiscale data capable of capturing hierarchical landscape influences allowed production of accurate predictive models and provided further insight into factors controlling freshwater fish distributions.     

Comparison and ranking of different modelling techniques for prediction of site index in Mediterranean mountain forests

Forestry science has a long tradition of studying the relationship between stand productivity and abiotic and biotic site characteristics, such as climate, topography, soil and vegetation. Many of the early site quality modelling studies related site index to environmental variables using basic statistical methods such as linear regression. Because most ecological variables show a typical non-linear course and a non-constant variance distribution, a large fraction of the variation remained unexplained by these linear models. More recently, the development of more advanced non-parametric and machine learning methods provided opportunities to overcome these limitations. Nevertheless, these methods also have drawbacks. Due to their increasing complexity they are not only more difficult to implement and interpret, but also more vulnerable to overfitting. Especially in a context of regionalisation, this may prove to be problematic. Although many non-parametric and machine learning methods are increasingly used in applications related to forest site quality assessment, their predictive performance has only been assessed for a limited number of methods and ecosystems.In this study, five different modelling techniques are compared and evaluated, i.e. multiple linear regression (MLR), classification and regression trees (CART), boosted regression trees (BRT), generalized additive models (GAM), and artificial neural networks (ANN). Each method is used to model site index of homogeneous stands of three important tree species of the Taurus Mountains (Turkey): Pinus brutia, Pinus nigra and Cedrus libani. Site index is related to soil, vegetation and topographical variables, which are available for 167 sample plots covering all important environmental gradients in the research area. The five techniques are compared in a multi-criteria decision analysis in which different model performance measures, ecological interpretability and user-friendliness are considered as criteria.When combining these criteria, in most cases GAM is found to outperform all other techniques for modelling site index for the three species. BRT is a good alternative in case the ecological interpretability of the technique is of higher importance. When user-friendliness is more important MLR and CART are the preferred alternatives. Despite its good predictive performance, ANN is penalized for its complex, non-transparent models and big training effort.     

Development of TRIPLEX-Management model for simulating the response of forest growth to pre-commercial thinning

In order to simulate forest growth response to pre-commercial thinning (PCT), TRIPLEX1.0 - a process-based model designed to predict forest growth as well as carbon (C) and nitrogen (N) dynamics - was modified and improved to also simulate managed forest ecosystem thinning practices. A three-parameter Weibull distribution model was integrated to simulate thinning treatments within the newly developed TRIPLEX-Management model. The thinning intensity component within the model allows users to simulate thinning treatments by applying basal area, stand density and volume to quantify thinning intensity. Natural mortality decreased following thinning due to an increase in growing space for residual stems. Predicted litterfall pools also increased after thinning events took place. The TRIPLEX-Management model was tested against published observational data for Jack Pine (Pinus banksiana Lamb.) stands subjected to PCT in Northwestern Ontario, Canada. The coefficients of determination (R²) between the predicted and observed variables including stand density, mean DBH (diameter at breast height), the quadratic mean DBH, total volume and merchantable volume as well as belowground, aboveground, and total biomass ranged from 0.50 to 0.88 (n = 20, P < 0.001) with the exception of mean tree height (R² = 0.25, n = 20, P < 0.05). Overall, the Willmott index of agreement between predicted and observed variables ranged from 0.97 to 1.00. Results show that the TRIPLEX-Management model is generally capable of simulating growth response to PCT for Jack Pine stands.     

Application of a Random Forest algorithm to predict spatial distribution of the potential yield of Ruditapes philippinarum in the Venice lagoon, Italy

We present a modelling framework that combines machine learning techniques and Geographic Information Systems to support the management of an important aquaculture species, Manila clam (Ruditapes philippinarum). We use the Venice lagoon (Italy), the first site in Europe for the production of R. philippinarum, to illustrate the potential of this modelling approach. To investigate the relationship between the yield of R. philippinarum and a set of environmental factors, we used a Random Forest (RF) algorithm. The RF model was tuned with a large data set (n = 1698) and validated by an independent data set (n = 841). Overall, the model provided good predictions of site-specific yields and the analysis of marginal effect of predictors showed substantial agreement among the modelled responses and available ecological knowledge for R. philippinarum. The most influent environmental factors for yield estimation were percentage of sand in the sediment, salinity, and water depth. Our results agree with findings from other North Adriatic lagoons. The application of the fitted RF model to continuous maps of all the environmental variables allowed estimates of the potential yield for the whole basin. Such a spatial representation enabled site-specific estimates of yield in different farming areas within the lagoon. We present a possible management application of our model by estimating the potential yield under the current farming distribution and comparing it to a proposed re-organization of the farming areas. Our analysis suggests a reduction of total yield is likely to result from the proposed re-organization.     

The MONICA model: Testing predictability for crop growth, soil moisture and nitrogen dynamics

A fundamentally revised version of the HERMES agro-ecosystem model, released under the name of MONICA, was calibrated and tested to predict crop growth, soil moisture and nitrogen dynamics for various experimental crop rotations across Germany, including major cereals, sugar beet and maize. The calibration procedure also included crops grown experimentally under elevated atmospheric CO₂ concentration. The calibrated MONICA simulations yielded a median normalised mean absolute error (nMAE) of 0.20 across all observed target variables (n = 42) and a median Willmott's Index of Agreement (d) of 0.91 (median modelling efficiency (ME): 0.75). Although the crop biomass, habitus and soil moisture variables were all within an acceptable range, the model often underperformed for variables related to nitrogen. Uncalibrated MONICA simulations yielded a median nMAE of 0.27 across all observed target variables (n = 85) and a median d of 0.76 (median ME: 0.30), also showing predominantly acceptable results for the crop biomass, habitus and soil moisture variables. Based on the convincing performance of the model under uncalibrated conditions, MONICA can be regarded as a suitable simulation model for use in regional applications. Furthermore, its ability to reproduce the observed crop growth results in free-air carbon enrichment experiments makes it suited to predict agro-ecosystem behaviour under expected future climate conditions.     

An integrated assessment method of urban drainage system: A case study in Shenzhen City, China

In recent years, the urban drainage system in China is facing the dual pressure of renovation and construction. This requires that the integrated assessment for the planning and operation of the urban drainage system is obligatory. To evaluate the urban drainage system, an integrated assessment methodology based on the analytic hierarchy process (AHP), integrated simulation, and fuzzy assessment is established. This method is a multi-criteria decision adding app roach to the assessment of the urban drainage system comprehensively. Through the integration of the Storm Water Management Model (SWMM), a simple wastewater treatment plant model, and a surface water quality model, an integrated modelling system for the urban drainage system is developed and applied as a key tool for assessment. Using the established method, a case study in Shenzhen City has been implemented to evaluate and compare two urban drainage system reno vation plans, the distributed plan and the centralized plan. Because of the particularity of this case study, the established method is not applied entirely. Considering the water environ mental impact, ecological impact, technological feasibility, and economic cost, the integrated performance of the distri buted plan is better. As shown in this case study, the proposed method is found to be both effective and practical.     

Obtaining Environmental Favourability Functions from Logistic Regression

Logistic regression is a statistical tool widely used for predicting species’ potential distributions starting from presence/absence data and a set of independent variables. However, logistic regression equations compute probability values based not only on the values of the predictor variables but also on the relative proportion of presences and absences in the dataset, which does not adequately describe the environmental favourability for or against species presence. A few strategies have been used to circumvent this, but they usually imply an alteration of the original data or the discarding of potentially valuable information. We propose a way to obtain from logistic regression an environmental favourability function whose results are not affected by an uneven proportion of presences and absences. We tested the method on the distribution of virtual species in an imaginary territory. The favourability models yielded similar values regardless of the variation in the presence/absence ratio. We also illustrate with the example of the Pyrenean desman’s (Galemys pyrenaicus) distribution in Spain. The favourability model yielded more realistic potential distribution maps than the logistic regression model. Favourability values can be regarded as the degree of membership of the fuzzy set of sites whose environmental conditions are favourable to the species, which enables applying the rules of fuzzy logic to distribution modelling. They also allow for direct comparisons between models for species with different presence/absence ratios in the study area. This makes them more useful to estimate the conservation value of areas, to design ecological corridors, or to select appropriate areas for species reintroductions. Received: June 2005 / Revised: July 2005     

Incorporating spatial autocorrelation into cellular automata model: An application to the dynamics of Chinese tamarisk (Tamarix chinensis Lour.)

Spatial autocorrelation (SAC) is frequently encountered in most spatial data in ecology. Cellular automata (CA) models have been widely used to simulate complex spatial phenomena. However, little has been done to examine the impact of incorporating SAC into CA models. Using image-derived maps of Chinese tamarisk (Tamarix chinensis Lour.), CA models based on ordinary logistic regression (OLCA model) and autologistic regression (ALCA model) were developed to simulate landscape dynamics of T. chinensis. In this study, significant positive SAC was detected in residuals of ordinary logistic models, whereas non-significant SAC was found in autologistic models. All autologistic models obtained lower Akaike's information criterion corrected for small sample size (AICc) values than the best ordinary logistic models. Although the performance of ALCA models only satisfied the minimum requirement, ALCA models showed considerable improvement upon OLCA models. Our results suggested that the incorporation of the autocovariate term not only accounted for SAC in model residuals but also provided more accurate estimates of regression coefficients. The study also found that the neglect of SAC might affect the statistical inference on underlying mechanisms driving landscape changes and obtain false ecological conclusions and management recommendations. The ALCA model is statistically sound when coping with spatially structured data, and the adoption of the ALCA model in future landscape transition simulations may provide more precise probability maps on landscape transition, better model performance and more reasonable mechanisms that are responsible for landscape changes.     

Statistical inference on patch-specific survival and movement rates from marked animals

Empirical estimates of patch-specific survival and movement rates are needed to parametrize spatially explicit population models, and for inference on the effects of habitat quality and fragmentation on populations. Data from radio-marked animals, in which both the fates and habitat locations of animals are known over time, can be used in conjunction with continuous-time proportional hazards models to obtain inferences on survival rates. Discrete-time conditional logistic models may provide inference on both survival and movement rates. We use Monte Carlo simulation to investigate accuracy of estimates of survival from both approaches, and movement rates from conditional logistic regression, for two habitats. Bias was low (relative bias < 0.04) and interval coverage accurate (close to the nominal 0.95) for estimates of habitat effect on survival based on proportional hazards. Bias was high ( relative bias 0.60) and interval coverage poor ( = 0.26 vs. nominal 0.95) for estimates of habitat effect based on conditional logistic regression; bias was especially influenced by heterogeneity in survival and the shape of the hazard function, whereas both bias and coverage were affected by ‘memory’ effects in movement patterns. Bias estimates of movement rate was low ( relative bias < 0.05), but interval coverage was poor ( = 0.48–0.80), possibly as a result of poor performance of a Taylor series estimate of variance. An example is provided from a radio-telemetry study of 47 wintering American woodcock (Scolopax minor), illustrating practical difficulties in field studies to parametrize these models. We also discuss extensions of continuous-time models to explicitly include a movement process, and further examine tradeoffs between continuous and discrete models.     

River network structure shapes interannual feedbacks between adult sea lamprey migration and larval habitation

Adult sea lampreys (Petromyzon marinus) migrating upstream to spawn follow a pheromone released by instream larvae. The size (i.e. flow) of a tributary dilutes the concentration of this pheromone, such that the downstream propagation pattern of larval pheromone must be influenced by patterns in the relative sizes and numbers of confluent tributaries. We developed an individual-based model to explicitly test the resulting hypothesis that river network structure influences the migration decisions of adult lampreys following the larval pheromone, and in turn the distribution of larvae. First, we initialized the model using randomly generated river networks, and found a strong positive relationship between network diameter and larval aggregation. Larvae aggregated over time, and the degree and rate of this aggregation depended on network diameter. Second, we initialized the model using a river network based on the Muskegon River, Michigan, and compared model-generated larval distribution to available field survey data. We found a significant correlation between model-generated larval abundance and field-measured larval densities (r₂ = 0.54; p < 0.0001). We also found an inverse relationship between subwatershed area and the degree to which path-dependent effects influenced larval abundance in that subwatershed. Our results overall suggest that larval distribution across a watershed results from a system of context-dependent interannual feedbacks shaped by network structure and the past migratory and spawning behavior of adults.     

A general process-based mass-balance model for phosphorus/eutrophication as a tool to estimate historical reference values for key bioindicators,as exemplified using data for the Gulf of Riga

This work describes how a general, process-based mass-balance model (CoastMab) for phosphorus for coastal areas may be used as a tool to estimate realistic values of “natural” or preindustrial reference levels of key bioindicators in coastal science, including the Secchi depth, a standard measure of water clarity, the chlorophyll-a concentration, an operational measure of phytoplankton biomass and the concentration of cyanobacteria, a measure of the concentration of harmful algae. The CoastMab-model is an ecosystem model giving monthly predictions to achieve seasonal variations of basin-wide properties. The selected case-study area, the Gulf of Riga, is sensitive to nutrient loading because of its shallowness and low openness towards the Baltic Proper. The morphometry of any coastal area, as given by the size and form parameters, influences all internal processes, such as sedimentation, resuspension, diffusion in water and from sediments to water, biouptake and retention in biota, stratification, mixing and outflow. There has been no mass-balance modeling for nitrogen (N) in this work because empirical data (from the HELCOM database) clearly indicate that the monthly primary production in the Gulf of Riga is regulated by phosphorus (P) – the mean monthly total-N to total-P ratios are well over 7.2 (the Redfield-ratio) and generally higher than 15 for the data used in this study (from 1992 to 2005). At present anthropogenic loads, the average modeled monthly values for Secchi depth, chlorophyll (Chl), cyanobacteria (CB) and total-P (TP) are 3.2 m, 3.8 μg/l, 78 μg/l and 31.3 μg/l, respectively. If 50% of all anthropogenic sources to the Gulf of Riga via rivers, point sources and diffuse sources were to be removed, these values would be 3.6 m, 3.4 μg Chl/l, 63 μg CB/l and 29.1 μg TP/l. If 60% of the anthropogenic phosphorus fluxes to the Baltic Proper were to be omitted and as well as 75% of all direct anthropogenic sources to the Gulf of Riga, the values would be 4.6 m, 2.7 μg Chl/l, 45 μg CB/l and 25.4 μg TP/l. These values represent the “natural” reference levels and it is not realistic to expect that remedial measures would improve the conditions more than that. Using the CoastWeb-model, similar calculations can be made for any given coastal area and the data necessary for such calculations are discussed in this work.     

Estimation of annual potential evapotranspiration at regional scale based on the effect of moisture on soil respiration

Potential evapotranspiration (PET) is an important component of water cycle. For traditional models derived from the principle of aerodynamics and the surface energy balance, its calculation always includes many parameters, such as net radiation, water vapor pressure, air temperature and wind speed. We found that it can be acquired in an easier way in specific regions. In this study, a new PET model (PETP model) derived from two empirical models of soil respiration was evaluated using the Penman-Monteith equation as a standard method. The results indicate that the PETP model estimation concur with the Penman-Monteith equation in sites where annual precipitation ranges from 717.71 mm to 1727.37 mm (R² = 0.68, p = 0.0002), but show large discrepancies in all sites (R² = 0.07, p = 0.1280). Then we applied our PETP model at the global scale to the regions with precipitation higher than 700 mm using 2.5° CMAP data to obtain the annual PET for 2006. As expected, the spatial pattern is satisfactory overall, with the highest PET values distributed in the lower latitudes or coastal regions, and with an average of 1292.60 ± 540.15 mm year⁻¹. This PETP model provides a convenient approach to estimate PET at regional scales.