Comparative evaluation of horizontal accuracy of elevations of selected ground control points from ASTER and SRTM DEM with respect to CARTOSAT-1 DEM: a case study of Shahjahanpur district,Uttar Pradesh,India

Digital elevation model (DEM) data of Shuttle Radar Topography Mission (SRTM) are distributed at a horizontal resolution of 90 m (30 m only for US) for the world, Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) DEM data provide 30 m horizontal resolution, while CARTOSAT-1 (IRS-P5) gives 2.6 m horizontal resolution for global coverage. SRTM and ASTER data are available freely but 2.6 m CARTOSAT-1 data are costly. Hence, through this study, we found out a horizontal accuracy for selected ground control points (GCPs) from SRTM and ASTER with respect to CARTOSAT-1 DEM to implement this result (observed from horizontal accuracy) for those areas where the 2.6-m horizontal resolution data are not available. In addition to this, the present study helps in providing a benchmark against which the future DEM products (with horizontal resolution less than CARTOSAT-1) with respect to CARTOSAT-1 DEM can be evaluated. The original SRTM image contained voids that were represented digitally as ?140; such voids were initially filled using the measured values of elevation for obtaining accurate DEM. Horizontal accuracy analysis between SRTM- and ASTER-derived DEMs with respect to CARTOSAT-1 (IRS-P5) DEM allowed a qualitative assessment of the horizontal component of the error, and the appropriable statistical measures were used to estimate their horizontal accuracies. The horizontal accuracy for ASTER and SRTM DEM with respect to CARTOSAT-1 were evaluated using the root mean square error (RMSE) and relative root mean square error (R-RMSE). The results from this study revealed that the average RMSE of 20 selected GCPs was 2.17 for SRTM and 2.817 for ASTER, which are also validated using R-RMSE test which proves that SRTM data have good horizontal accuracy than ASTER with respect to CARTOSAT-1 because the average R-RMSE of 20 GCPs was 3.7 × 10^?4 and 5.3 × 10^?4 for SRTM and ASTER, respectively.     

Assessment of the EUMETSAT LSA-SAF evapotranspiration product for drought monitoring in Europe

Evapotranspiration is a key parameter for water stress assessment as it is directly related to the moisture status of the soil-vegetation system and describes the moisture transfer from the surface to the atmosphere. With the launch of the Meteosat Second Generation geostationary satellites and the setup of the Satellite Application Facilities, it became possible to operationally produce evapotranspiration data with high spatial and temporal evolution over the entire continents of Europe and Africa. In the frame of this study we present an evaluation of the potential of the evapotranspiration (ET) product from the EUMETSAT Satellite Application Facility on Land Surface Analysis (LSA-SAF) for drought assessment and monitoring in Europe.To assess the potential of this product, the LSA-SAF ET was used as input for the ratio of ET to reference evapotranspiration (ET₀), the latter estimated from the ECMWF interim reanalysis. In the analysis two case studies were considered corresponding to the drought episodes of spring/summer 2007 and 2011. For these case studies, the ratio ET/ET₀ was compared with meteorological drought indices (SPI, SPEI and Sc-PDSI for 2007 and SPI for 2011) as well as with the anomalies of the fraction of absorbed photosynthetic active radiation (fAPAR) derived from remote sensing data. The meteorological and remote sensing indicators were taken from the European Drought Observatory (EDO) and the CARPATCLIM climatological atlas.Results show the potential of ET/ET₀ to characterize soil moisture variability, and to give additional information to fAPAR and to precipitation distribution for drought assessment. The main limitations of the proposed ratio for drought characterization are discussed, including options to overcome them. These options include the use of filters to discriminate areas with a low percentage vegetation cover or areas that are not in their growing period and the use of evapotranspiration without water restriction (ET_wwr), obtained as output of the LSA-SAF model instead of ET₀. The ET/ET_wwr ratio was tested by comparing its accumulated values per growing period with the winter wheat yield values per country published by Eurostat. The results point to the potential of using the remote sensing based LSA-SAF evapotranspiration and the ET/ET_wwr ratio for vegetation monitoring at large scale, especially in areas where data is generally lacking.     

Application of S-SEBI model for crop evapotranspiration using Landsat-8 data over parts of North India

Estimation and monitoring of crop evapotranspiration (ET_c) or consumptive water use over large-area holds the key to irrigation management plans and regional drought preparedness. The objective of this study was to estimate ET_c by applying the simplified-surface energy balance index (S-SEBI) model to Landsat-8 data for the 2014–2015 period in parts of North India. An average ET_c was estimated 2.72 and 2.47 in mm day^?1 with 0.22, 0.18 standard deviation and 0.11, 0.07 standard error for Kharif and Rabi crops, respectively. On validation part, a close relationship was observed between S-SEBI derived and scintillometer observed evaporative fraction with 0.85 correlation coefficient and 0.86 agreement index. The statistical analysis also endorses the results accuracy and reliability with 0.026 and 0.602, relative root-mean square errors and model efficiency for wheat crop, respectively. The study showed that normalized difference vegetation index and LST are closely related and serve as a proxy for qualitative representation of ET_c.     

Assessing the performance of a large-scale irrigation system by estimations of actual evapotranspiration obtained by Landsat satellite images resampled with cubic convolution

Remote sensing techniques allow monitoring the Earth surface and acquiring worthwhile information that can be used efficiently in agro-hydrological systems. Satellite images associated to computational models represent reliable resources to estimate actual evapotranspiration fluxes, ET_a, based on surface energy balance. The knowledge of ET_a and its spatial distribution is crucial for a broad range of applications at different scales, from fields to large irrigation districts. In single plots and/or in irrigation districts, linking water volumes delivered to the plots with the estimations of remote sensed ET_a can have a great potential to develop new cost-effective indicators of irrigation performance, as well as to increase water use efficiency. With the aim to assess the irrigation system performance and the opportunities to save irrigation water resources at the “SAT Llano Verde” district in Albacete, Castilla-La Mancha (Spain), the Surface Energy Balance Algorithm for Land (SEBAL) was applied on cloud-free Landsat 5 Thematic Mapper (TM) images, processed by cubic convolution resampling method, for three irrigation seasons (May to September 2006, 2007 and 2008). The model allowed quantifying instantaneous, daily, monthly and seasonal ET_a over the irrigation district. The comparison between monthly irrigation volumes distributed by each hydrant and the corresponding spatially averaged ET_a, obtained by assuming an overall efficiency of irrigation network equal to 85%, allowed the assessment of the irrigation system performance for the area served by each hydrant, as well as for the whole irrigation district. It was observed that in all the investigated years, irrigation volumes applied monthly by farmers resulted generally higher than the corresponding evapotranspiration fluxes retrieved by SEBAL, with the exception of May, in which abundant rainfall occurred. When considering the entire irrigation seasons, it was demonstrated that a considerable amount of water could have been saved in the district, respectively equal to 26.2, 28.0 and 16.4% of the total water consumption evaluated in the three years.     

Evaluating five remote sensing based single-source surface energy balance models for estimating daily evapotranspiration in a humid subtropical climate

In the last two decades, a number of single-source surface energy balance (SEB) models have been proposed for mapping evapotranspiration (ET); however, there is no clear guidance on which models are preferable under different conditions. In this paper, we tested five models-Surface Energy Balance Algorithm for Land (SEBAL), Mapping ET at high Resolution with Internalized Calibration (METRIC), Simplified Surface Energy Balance Index (S-SEBI), Surface Energy Balance System (SEBS), and operational Simplified Surface Energy Balance (SSEBop)—to identify the single-source SEB models most appropriate for use in the humid southeastern United States. ET predictions from these models were compared with measured ET at four sites (marsh, grass, and citrus surfaces) for 149 cloud-free Landsat image acquisition days between 2000 and 2010. The overall model evaluation statistics showed that SEBS generally outperformed the other models in terms of estimating daily ET from different land covers (e.g., the root mean squared error (RMSE) was 0.74 mm day⁻¹). SSEBop was consistently the worst performing model and overestimated ET at all sites (RMSE = 1.67 mm day⁻¹), while the other models typically fell in between SSEBop and SEBS. However, for short grass conditions, SEBAL, METRIC, and S-SEBI appear to work much better than SEBS. Overall, our study suggests that SEBS may be the best SEB model in humid regions, although it may require modifications to work better over short vegetation.     

Estimation of Partial Pressure of Carbon Dioxide and Air-Sea Fluxes in Hooghly Estuary Based on In Situ and Satellite Observations

An empirical model is developed and used with remotely sensed predictors: sea surface temperature (SST) and chlorophyll-a concentration (Chl-a), to compute surface water partial pressure of carbon dioxide (pCO_2w) and air-sea fluxes of CO₂ in the Hooghly estuary and its adjacent coastal oceans. In situ observations used here were based on measurements carried out in this region during winter and summer periods in 2008. The estimated pCO_2w compares well with the in situ observations at root mean square error ±18 μatm. In winter, estimated pCO_2w ranges between 320 and 500 μatm with large values (>400 μatm) on the south-western and south-eastern flanks of the coastal domain and lower values (340–375 μatm) on the main-channel. In summer, it remained spatially uniform at 450 μatm. Extrapolation of the results over the study region based on the Moderate Imaging Specroradiometer (MODIS) measured SST and Chl-a suggests that the region is a strong source of atmospheric CO₂ during the summer with net release of 0.095 Tg C year^?1 (equivalent to mean flux of 90 molC m^?2 year^?1) and is a weak source during the winter with net release of 0.006 Tg C yr^?1 (0.5 molC m^?2 year^?1) from the geographical extent of 6000 Km² area.     

Pasture enclosures increase soil carbon dioxide flux rate in Semiarid Rangeland,Kenya

Background

Pasture enclosures play an important role in rehabilitating the degraded soils and vegetation, and may also influence the emission of key greenhouse gasses (GHGs) from the soil. However, no study in East Africa and in Kenya has conducted direct measurements of GHG fluxes following the restoration of degraded communal grazing lands through the establishment of pasture enclosures. A field experiment was conducted in northwestern Kenya to measure the emission of CO₂, CH₄ and N₂O from soil under two pasture restoration systems; grazing dominated enclosure (GDE) and contractual grazing enclosure (CGE), and in the adjacent open grazing rangeland (OGR) as control. Herbaceous vegetation cover, biomass production, and surface (0–10 cm) soil organic carbon (SOC) were also assessed to determine their relationship with the GHG flux rate.

Results

Vegetation cover was higher enclosure systems and ranged from 20.7% in OGR to 40.2% in GDE while aboveground biomass increased from 72.0 kg DM ha^?1 in OGR to 483.1 and 560.4 kg DM ha^?1 in CGE and GDE respectively. The SOC concentration in GDE and CGE increased by an average of 27% relative to OGR and ranged between 4.4 g kg^?1 and 6.6 g kg^?1. The mean emission rates across the grazing systems were 18.6 μg N m^?2 h^?1, 50.1 μg C m^?2 h^?1 and 199.7 mg C m^?2 h^?1 for N₂O, CH₄, and CO₂, respectively. Soil CO₂ emission was considerably higher in GDE and CGE systems than in OGR (P?<?0.001). However, non-significantly higher CH₄ and N₂O emissions were observed in GDE and CGE compared to OGR (P?=?0.33 and 0.53 for CH₄ and N₂O, respectively). Soil moisture exhibited a significant positive relationship with CO₂, CH₄, and N₂O, implying that it is the key factor influencing the flux rate of GHGs in the area.

Conclusions

The results demonstrated that the establishment of enclosures in tropical rangelands is a valuable intervention for improving pasture production and restoration of surface soil properties. However, a long-term study is required to evaluate the patterns in annual CO₂, N₂O, CH₄ fluxes from soils and determine the ecosystem carbon balance across the pastoral landscape.

     

A comparison of geographic datasets and field measurements to model soil carbon using random forests and stepwise regressions (British Columbia,Canada)

We used geographic datasets and field measurements to examine the mechanisms that affect soil carbon (SC) storage for 65 grazed and non-grazed pastures in southern interior grasslands of British Columbia, Canada. Stepwise linear regression (SR) modeling was compared with random forest (RF) modeling. Models produced with SR performed better than those produced using RF models (r² = 0.56–0.77 AIC = 0.16–0.30 for SR models; r² = 0.38–0.53 and AIC = 0.18–0.30 for RF models). The factors most significant when predicting SC were elevation, precipitation, and the normalized difference vegetation index (NDVI). NDVI was evaluated at two scales using: (1) the MOD 13Q1 (250 m/16-day resolution) NDVI data product from the moderate resolution imaging spectro-radiometer (MODIS) (NDVI_MODIS), and (2) a handheld multispectral radiometer (MSR, 1 m resolution) (NDVI_MSR) in order to understand the potential for increasing model accuracy by increasing the spatial resolution of the gridded geographic datasets. When NDVI_MSR data were used to predict SC, the percentage of the variance explained by the model was greater than for models that relied on NDVI_MODIS data (r² = 0.68 for SC for non-grazed systems, modeled with SR based on NDVI_MODIS data; r² = 0.77 for SC for non-grazed systems, modeled with SR based on NDVI_MSR data). The outcomes of this study provide the groundwork for effective monitoring of SC using geographic datasets to enable a carbon offset program for the ranching industry.     

Application of Remotely-Sensed Data to Estimate a Water Budget for Data-Scarce Endorheic Basins: A Case Study of Lake Urmia basin,Iran

Estimating the water budgets of large basins is a challenge because of the lack of data and information. It becomes more complicated in endorheic basins that consist of separate land and water phases. The application of remotely-sensed data is one solution in this regard. The present study addresses this issue and develops a modeling framework to evaluate a water budget based on remotely-sensed data for endorheic basins. To explore the methodology, Lake Urmia basin was selected as a case study. The lake water level has declined steeply since 1995 and stakeholders have agreed to allocate 3100 MCM of water per year to the lake. This makes it necessary to monitor river inflow into the lake to fulfill the agreement. Gauging stations have been employed around the lake, but they could not account for shortages such as water uptake below the stations. To do this, separate water budgets for the water body and the land were required. More specifically, it was necessary to estimate actual evapotranspiration (ET _a ) from freshwater (E _f ) and saltwater (E _s ) estimated using the SEBAL model. Different methods were applied to estimate soil moisture, groundwater exploitation, and surface-groundwater inflow into the lake. A comparison of the observed and estimated amounts showed good agreement. For instance, the coefficient of determination for the observed/reported and estimated ET _a and E _f were 0.83 and 0.84, respectively. The average annual inflow was estimated to be 2.2 BCM/year for 2002–2008 using the RS model, which is about 84 % of the total inflow from the last recording stations before the lake and shows influence of water exploitation after these stations. Future study should focus on increasing temporal and spatial resolution of the method     

Spatio-temporal patterns of energy exchange and evapotranspiration during an intense drought for drylands in Brazil

The Caatinga biome, located in the northeastern region of Brazil, is the most populated dryland region on the planet and extremely vulnerable to land degradation due to climatological and anthropogenic factors. Energy partitioning substantially influences the local climate and affects the water cycle, which is of utmost importance for the economy and livelihood of the region. Recently, eddy covariance (EC) towers were installed in the area; thus, the scientific community can thoroughly assess the water and energy fluxes over this unique biome. While EC towers have a high degree of accuracy, they only measure energy fluxes over a small land footprint. Given the biome spatial heterogeneity, the use of EC-based techniques has the limitation of not comprehensively representing water and energy fluxes profiles over the entire region. Incorporating remote sensing (RS) data into the landscape analysis is a feasible solution to overcome this issue, given that satellite data can capture the phenomena represented by the EC measurements across large spatial scales. Our research studied the capability of the Surface Energy Balance Algorithm for Land (SEBAL) and MOD16-ET products to represent the EC measurements regarding energy and mass exchange, with an ultimate objective of applying the best approach to assess these fluxes regionally. We applied the SEBAL model using only remote sensing data from the Moderate Resolution Imaging Spectroradiometer (MODIS) sensor. The MOD16-ET model uses a different approach but is also based on MODIS data. Our analysis was based on three years (2014–2016) of data, which was limited to the availability of the EC tower data. We found that for the EC-based measurements, energy balance closure (EBC) achieved an average of 0.84, which is considerably high for the region. This is possibly due to the EC tower being installed on a preserved Caatinga plot, with reduced heterogeneity and higher plant density. When analyzing RS-based products to represent ET profiles in the region, we found that the SEBAL model accurately represented water fluxes during the wet season but not the dry season, whereas the MOD16-ET showed a better agreement with EC-based water fluxes throughout all the seasons. SEBAL inaccuracy in drylands is partially due to the narrow range between the cold and hot pixels in an image, as the algorithm relies on this range for input parameters, especially in the dry season. Therefore, we concluded that MOD16-ET is capable of better-representing water fluxes in the Caatinga region. We analyzed the fluxes regionally and quantified annual ET for the three years. These results are especially relevant for local policymakers on dealing with water and landscape issues in a region where the livelihood and well-being of the population is inextricably bound to water availability.     

Retrieving water surface temperature from archive LANDSAT thermal infrared data: Application of the mono-channel atmospheric correction algorithm over two freshwater reservoirs

Water surface temperature is a key element in characterizing the thermodynamics of waterbodies, and for irregularly-shaped inland reservoirs, LANDSAT thermal infrared images are the best alternative yet for the retrieval of this parameter. However, images must be corrected mainly for atmospheric effects in order to be fully exploitable. The objective of this study is to validate the mono-channel correction algorithm for single-band thermal infrared LANDSAT data as put forward by Jiménez-Muñoz et al. (2009). Two freshwater reservoirs in continental France were selected as study sites, and best use was made of all accessible image and field data. Results obtained are satisfactory and in accordance with the literature: r² values are above 0.90 and root-mean-square error values are comprised between 1 and 2 °C. Moreover, paired Wilcoxon signed rank tests showed a highly significant difference between field and uncorrected image data, a very highly significant difference between uncorrected and corrected image data, and no significant difference between field and corrected image data. The mono-channel algorithm is hence recommended for correcting archive LANDSAT single-band thermal infrared data for inland waterbody monitoring and study.     

基于蒸渗仪的蒸散量时间尺度扩展方法对比

蒸散量（ET,Evapotranspiration）时间尺度扩展方法模拟效果受多种因素影响。本文利用蒸渗仪实测的冬小麦生育期ET数据,对比分析基于蒸发比、作物系数和冠层阻力的不同ET时间尺度扩展方法的估算效果。结果表明,对小时到日尺度的ET时间扩展而言,利用上午时段（10：00-11：00）的ET实测数据,基于冠层阻力的日ET尺度扩展值估算效果要优于基于蒸发比和作物系数的效果,而利用下午时段（14：00-15：00）的ET实测数据,基于作物系数的日ET尺度扩展值估算效果要优于基于蒸发比和冠层阻力的效果,对于从典型日到生育期尺度的ET时间扩展,基于作物系数的生育期ET尺度扩展值估算效果要好于基于蒸发比和冠层阻力的效果。应根据当地的气候气象条件、作物种植状况、时间尺度扩展类型、土壤水分、ET观测设施种类等因素,择优适宜的ET时间尺度扩展方法。     

Detection of Changes in Glacier Mass Balance Using Satellite and Meteorological Data in Tirungkhad Basin Located in Western Himalaya

Himalayan glaciers and their mass balance are poorly sampled through direct mass balance measurements. Thus, based on Landsat datasets of ETM⁺ (2000), ETM⁺ (2006) and TM (2011), mass balance studies of 32 glaciers was carried out using accumulation area ratio (AAR) method in the Tirungkhad river basin, a tributary of Satluj River, located in western Himalayan region. The overall specific mass balance was negative varying from ?27 cm (2000) to ?41 cm (2011). Out of 32 glaciers, 27 glaciers (81.2 %) showed negative mean mass balance and 5 glaciers (18.7 %) showed positive mean mass balance. Mean of specific mass balance for the year 2000, 2006 and 2011 was found to be ?48 cm, ?55 cm and ?0.61 cm respectively, in case of glaciers with negative mass balance while in case of glaciers with positive mass balance, it was 0.67 cm (2000), 0.56 cm (2006) and 0.47 cm (2011). The investigations suggested a loss of ?0.034 km³ of glacial ice for 2000, 0.036 km³ for 2006 and 0.038 km³ for 2011 respectively. The negative mass balance of the glaciers since 2000 correlates well with the increasing trend of annual mean temperature and decreasing trend of precipitation (snow water equivalent (SWE) and rainfall). Based on Mann Kendall test the temperature and SWE trends were significant at 95 % confidence level, however, the rainfall trend was insignificant.     

Quantification of LST and CO2 levels using Landsat-8 thermal bands on urban environment

Land surface temperature (LST) is an important aspect in global to regional change studies, for control of climate change and balancing of high temperature. Urbanization is one of the influencing factors increasing land surface and atmospheric temperature, by the emission of greenhouse gases (e.g. CO₂, NO and methane). In the present study, LST was derived from Landsat-8 of multitemporal data sets to analyse the spatial structure of the urban thermal environment in relation to the urban surface characteristics and land use–land cover (LULC). LST is influenced by the greenhouse gases i.e. CO₂ plays an important role in increasing the earth’s surface temperature. In order to provide the evidence of influence of CO₂ on LST, the relationship between LST, air temperature and CO₂ was analysed. Landsat-8 satellite has two thermal bands, 10 and 11. These bands were used to accurately to calculate the temperature over the study area. Results showed that the strength of correlation between ground monitoring data and satellite data was high. Based on correlation values of each month April (R² = 0.994), May (R² = 0.297) and June (R² = 0.934), observed results show that band 10 was significantly correlating with air temperature. Relationship between LST and CO₂ levels were obtained from linear regression analysis. band 11 was correlating significantly with CO₂ values in each of the months April (R² = 0.217), May (R² = 0.914) and June, (R² = 0.934), because band 11 is closer to the 15-micron band of CO₂. From the results, it was observed that band 10 can be used for calculating air temperature and band 11 can be used for estimation of greenhouse gases.     

Biomass estimation of Sonneratia caseolaris (l.) Engler at a coastal area of Hai Phong city (Vietnam) using ALOS-2 PALSAR imagery and GIS-based multi-layer perceptron neural networks

This study tested the use of machine learning techniques for the estimation of above-ground biomass (AGB) of Sonneratia caseolaris in a coastal area of Hai Phong city, Vietnam. We employed a GIS database and multi-layer perceptron neural networks (MLPNN) to build and verify an AGB model, drawing upon data from a survey of 1508 mangrove trees in 18 sampling plots and ALOS-2 PALSAR imagery. We assessed the model’s performance using root-mean-square error, mean absolute error, coefficient of determination (R²), and leave-one-out cross-validation. We also compared the model’s usability with four machine learning techniques: support vector regression, radial basis function neural networks, Gaussian process, and random forest. The MLPNN model performed well and outperformed the machine learning techniques. The MLPNN model-estimated AGB ranged between 2.78 and 298.95 Mg ha^?1 (average = 55.8 Mg ha^?1); below-ground biomass ranged between 4.06 and 436.47 Mg ha^?1 (average = 81.47 Mg ha^?1), and total carbon stock ranged between 3.22 and 345.65 Mg C ha^?1 (average = 64.52 Mg C ha^?1). We conclude that ALOS-2 PALSAR data can be accurately used with MLPNN models for estimating mangrove forest biomass in tropical areas.     

Accuracy of GPS-derived relative positions as a function of interstation distance and observing-session duration

 Ten days of GPS data from 1998 were processed to determine how the accuracy of a derived three-dimensional relative position vector between GPS antennas depends on the chord distance (denoted L) between these antennas and on the duration of the GPS observing session (denoted T). It was found that the dependence of accuracy on L is negligibly small when (a) using the `final' GPS satellite orbits disseminated by the International GPS Service, (b) fixing integer ambiguities, (c) estimating appropriate neutral-atmosphere-delay parameters, (d) 26 km ≤ L ≤ 300 km, and (e) 4 h ≤T ≤ 24 h. Under these same conditions, the standard error for the relative position in the north–south dimension (denoted S _n and expressed in mm) is adequately approximated by the equation S _n =k _n /T ^0.5 with k _n =9.5 ± 2.1 mm · h^0.5 and T expressed in hours. Similarly, the standard errors for the relative position in the east–west and in the up-down dimensions are adequately approximated by the equations S _e =k _e /T ^0.5 and S _u =k _u /T ^0.5, respectively, with k _e =9.9 ± 3.1 mm · h^0.5 and k _u =36.5 ± 9.1 mm · h^0.5. Received: 5 February 2001 / Accepted: 14 May 2001     

Testing of global pressure/temperature (GPT) model and global mapping function (GMF) in GPS analyses

Several sources of a priori meteorological data have been compared for their effects on geodetic results from GPS precise point positioning (PPP). The new global pressure and temperature model (GPT), available at the IERS Conventions web site, provides pressure values that have been used to compute a priori hydrostatic (dry) zenith path delay z _h estimates. Both the GPT-derived and a simple height-dependent a priori constant z _h performed well for low- and mid-latitude stations. However, due to the actual variations not accounted for by the seasonal GPT model pressure values or the a priori constant z _h, GPS height solution errors can sometimes exceed 10 mm, particularly in Polar Regions or with elevation cutoff angles less than 10 degrees. Such height errors are nearly perfectly correlated with local pressure variations so that for most stations they partly (and for solutions with 5-degree elevation angle cutoff almost fully) compensate for the atmospheric loading displacements. Consequently, unlike PPP solutions utilizing a numerical weather model (NWM) or locally measured pressure data for a priori z _h, the GPT-based PPP height repeatabilities are better for most stations before rather than after correcting for atmospheric loading. At 5 of the 11 studied stations, for which measured local meteorological data were available, the PPP height errors caused by a priori z _h interpolated from gridded Vienna Mapping Function-1 (VMF1) data (from a NWM) were less than 0.5 mm. Height errors due to the global mapping function (GMF) are even larger than those caused by the GPT a priori pressure errors. The GMF height errors are mainly due to the hydrostatic mapping and for the solutions with 10-degree elevation cutoff they are about 50% larger than the GPT a priori errors.