## versión On-line ISSN 0718-9516

### J. Soil Sci. Plant Nutr. vol.12 no.2 Temuco  2012

#### http://dx.doi.org/10.4067/S0718-95162012000200009

Journal of Soil Science and Plant Nutrition, 2012, 12 (2), 295-302

Simple linear algorithm to estimate the space-time variability of precipitable water in the Araucanía Region, Chile

L. Morales-Salinas1*, J.C. Parra-Aravena2,3, F. Lang-Tasso1, R. Abarca-Del Río4, E. Jorquera-Fontena5

1Departament of Environmental Sciences and Renewable Natural Resources, Faculty of Agronomic Sciences, Universidad de Chile. P.O. Box 1004, La Pintana, Santiago.*Corresponding author: lmorales@renare.uchile.c
2Department of Physical Sciences, Faculty of Engineering, Sciences and Administration, Universidad de La Frontera, Casilla 54-D, Temuco, Chile.
3Center for Optics and Photonics, Universidad de Concepción, Casilla 4016, Concepción, Chile.
4Departament of Geophysics, Faculty of Physical and Mathematical Sciences, Universidad de Concepción, P.O. Box 160-C, Concepción, Chile.
5
Doctorate Programme in Natural Resources Sciences. Universidad de La Frontera. Casilla 54-D, Temuco, Chile.

Abstract

This work proposes a linear algorithm which relates precipitable water (PW) with altitude and distance from the coast, on a regional scale. For precipitable water, the ERA-Interim (ECMWF) database for the period 1989-2010 was used, while a digital elevation model was used for the geographical variables. The results obtained indicate statistical significance between the two variables (p< 0.01), with a confidence level of 99%. The model makes it possible to describe the mean monthly space-time variability of precipitable water, stratifying the study area into clearly differentiated zones.

Keywords: linear algorithm, precipitable water, digital elevation model.

1. Introduction

Atmospheric water vapour is the principal contributor to the greenhouse effect and plays a key role in our understanding of the Earth's climate. Precipitable water (PW) is the amount of vertically integrated water vapour and can be expressed in g cm-2, or as the height of an equivalent column of liquid water in cm. Traditionally, there have been three atmospheric water vapour observing systems on a global scale: radio soundings, general circulation models and satellite observations. PW is an important component of the hydrological cycle and is adopted as an input variable in global climatological studies. Moreover, it has the potential to support hydrological, biospheric and atmospheric modelling efforts, on both local and regional scales, since it is widely used in energy budget and evapotranspiration studies (Hadjimitsis et al., 2011). Furthermore, PW is an essential requirement in atmospheric correction of high spatial resolution satellite data, and it is also necessary for enhancing the precision of land surface temperature estimates obtained from satellite data (Wan, 1999; Morales et al., 2002; Morales and Parra, 2002; Parra et al, 2006; Vera et al, 2010).

One way of obtaining the PW is through radio soundings or database compilations, such as those extracted from Moderate Resolution Imaging Spectroradiometer (MODIS) and ERA-Interim, an ensemble of meteorological data based on the Reanalysis carried out by the European Centre for Medium-Range Weather Forecasts. ERA-Interim (Dee et al, 2011) considers information from 1 January 1989 to the present, and can be accessed on line (http://www.ecmwf.int/research/era/do/get/era-interim). The information is presented in the form of both daily and monthly average values, and covers the whole planet with a grid resolution of 1.5 degrees of latitude and longitude. For each node of the grid, profiles are available for temperature, pressure, specific humidity and other variables.

The objective of the present communication is to propose a simple linear algorithm for estimating and describing the monthly mean space-time variability of precipitable water, on a regional scale, based on the geographical characteristics of the observation point.

2. Materials and methods

2.1 Study area

The study area is the Araucanía Region, Chile (Figure 1), with a total area of 31,840.3 km2. The study area presents two well differentiated climatic characteristics. The northern part of the Region has a warm temperate climate with a short dry season (less than 4 months), while the part lying to the south of latitude -38° is associated with a wet temperate climate, with Mediterranean influence. The Region presents distinct geographical units: coastal plain, coastal mountain range, central depression, Andean pre-cordillera and Andes Mountains, with average altitudes varying between 200 and 2,000 m.a.s.l.

2.2 Precipitable water

For this study, the average monthly precipitable water values were obtained from data supplied by ERA-Interim for the period 1989-2010. Considering that the grid resolution is around 166.5 x 166.5 km, the data were interpolated to bring them down to a resolution of 20 x 20 km, using kriging interpolation by blocks with a linear variogram. The interpolated information was used to apply downscaling (Fuentes et al., 2000; Wetterhall et al., 2007) to obtain the precipitable water vapour at a resolution of 1 km x 1 km, to coincide with the spatial resolution of one pixel of the NOAA satellite image. The following linear model is proposed for the parameterization of precipitable water as a function of geographical characteristics such as altitude and distance from the coast:

A, B and C are monthly mean coefficients depending on latitude and longitude; ALT is the altitude (m), DIL represents the distance from the coast (m) and g represents the error associated with estimation in g cm-2. Equation 1 is founded because the amount of water vapour in the atmosphere is affected by temperature variations and depends on variations in altitude (ALT), and geographical conditions, such as distance to the Pacific Ocean (DIL) (AGU, 2002; Parwati et al, 2007). The application of the model assumes the use of a digital elevation model, with ALT as a variable. In this study we used GTOPO30 (http://edc.usgs.gov/ products/elevation/gtopo30/gtopo30.html), geo-referenced to the datum WGS84 (Figure 1).

3. Results

Table 1 shows, by month: regression coefficients A, B and C; the coefficient of determination (r2); the statistical significance (p); and the standard error g. It may be seen that the proposal of the algorithm to describe the relation between precipitable water, altitude and distance from the coast is significant. In fact, the statistical significance, on average p< 0.01, indicates that the relation between the variables considered is statistically significant, with a confidence level of 99%.

The value of r2 shows that the model used explains, on average, 83 % of the variability in precipitable water vapour, representing a strong relationship between the variables.

The time variability of coefficients A, B and C of equation (1) for the Araucanía Region of Chile are shown in Figure 2.

In order to obtain categories with uniform values for PW, a non-hierarchical cluster analysis was done using the minimum distance criterion (Perez, 2004). This technique allowed pixels with similar precipitable water vapour values to be considered in a single cluster, while pixels with different values were located in other clusters. Figure 3 shows the result of the cluster analysis carried out on the monthly mean values; it may be clearly observed that the precipitable water is stratified into different zones. These correspond to the coastal area, characterised by the Nahuelbuta coastal range; the central depression; the Andean pre-cordillera; and the Andes Mountains.

The result of the cluster analysis in the study area shows that precipitable water vapour follows space-time trends defined at a mean monthly level, which are very similar to the Koppen climate classification system. This fact is based on the observed dependence on geographical position, altitude and, in the case of this study area, distance from the coast. The same descriptors are used to model the spatial distribution of the climate in the study area, with satisfactory results (Díaz et al, 2010). In the study area, water vapour is transported primarily from the Pacific Ocean (Morales et al, 2009).

Table 2 shows the statistics describing each cluster in monthly values, with columns for the average (AVG), the Range (Wmax - Wmin ) (RAN), and the coefficient of variation in percentage (CV) of the precipitable water for each cluster in the time period considered.

Figure 4 shows the behaviour of the precipitable water for each cluster in the study area. In general terms, all the curves follow a similar comportment over time, with maximum values concentrated between January and March. This coincides with maximum solar radiation at median latitudes in the southern hemisphere, and therefore with a considerable increase in evaporation from the Pacific Ocean.

In addition, it can be seen that the cluster with the greatest precipitable water is cluster 1, the coastal zone, while the clusters with intermediate values are the central depression and the pre-cordillera. Cluster 2 equates to the zone with the lowest precipitable water content, corresponding to the Andes range in the Araucanía Region.

4. Conclusions

A linear algorithm relating precipitable water to altitude and distance from the coast makes it possible to estimate and describe the space-time variability of precipitable water. This study, done on a regional scale, produced a statistical significance between the two variables of p< 0.01, indicating that the relation is statistically significant, with confidence level 99%.

A non-hierarchical cluster analysis, applied to the values for precipitable water, allowed the Araucanía Region to be stratified into different uniform zones, with the highest values of W found in the coastal zone, intermediate values in the central depression and pre-cordillera, and the lowest in the Andes range.

The proposed model offers an operational solution for obtaining PW information as a monthly mean when no daily information from radio soundings or satellite platforms is available. This provides an entry variable on which to base atmospheric corrections to thermal images, e.g. those supplied by the NOAA and TERRA platforms. This allows for applications in agricultural sciences when there is a need to estimate particular biophysical variables, such as thermal inertia and evapotranspiration, which appear to be closely linked to soil surface temperature.

Acknowledgements

J.C. Parra and L. Morales are grateful for financial support from project DI11-0033 of La Frontera University. J.C. Parra acknowledges Grant CONICYT PFB08-024.

References

American Geophysical Union (AGU). 2002. Water vapor in the climate system. Florida Ave., N.W., Washington DC 20009.         [ Links ]

Dee, D.P., S.M, Uppala., A.J, Simmons., P, Berrisford., P, Poli., S, Kobayashi., U, Andrae., M.A, Balmaseda., G, Balsamo., P, Bauer., P. Bechtold., A.C.M, Beljaar., L, van de Berg., J, Bidlot., N, Bormann., C, Delsol., R, Dragani., M, Fuentes., A.J, Geer., L, Haimberger., S.B. Healy., H, Hersbach., E.V, Holm., L, Isaksen., P, Kallberg., M, Kohler., M, Matricardi., A P, McNally., B.M, Monge-Sanz., J.J, Morcrette., B.K, Park., C, Peubey., P,de Rosnay., C, Tavolato., J.N, Thépaut., F,Vitart. 2011. The Era interim reanalysis: configuration and performance of the data assimilation system. Q.J.R. Meteorol. Soc. 137, 553 - 597.         [ Links ]

Díaz M, Diego., Luis, Morales S., Giorgio, Castellaro G., Fernando, Neira R. 2010. Topoclimatic modeling of thermopluviometric variables for the Bío-Bío and La Araucanía Regions, Chile. Chilean Journal of Agricultural Research. 70, 4, 604-615.         [ Links ]

Fuentes, U., D. Heimann. 2000. An improved statistical-dynamical downscaling scheme and its application to the Alpine precipitation climatology. Theor. Appl. Climatol. 65, 119-135.         [ Links ]

Hadjimitsis, D., Z. Mitraka., I. Gazani2., A, Retalis., N. Chrysoulakis., S. Michaelides. 2011. Estimation of spatio-temporal distribution of precipitable water using MODIS and AVHRR data: a case study for Cyprus. Adv. Geosci. 30, 23-29.         [ Links ]

Morales, L., J.C. Parra., J.A. Sobrino. 2002. Estimating land surface temperature in South America from NOAA-AVHRR images and reanalysis data. p. 698-707. Proc. First International Symposium on Recent Advances in Quantitative Remote Sensing. 16-20 September. Universidad de Valencia, Valencia, España.         [ Links ]

Morales L., J. C, Parra. 2002. Estimación de la radiancia emitida por la atmósfera en el infrarrojo térmico usando datos del reanálisis. Revista Mexicana de Física, Vol. 48 S3, pp. 119-121.         [ Links ]

Morales L., C.Mattar, L. Da-Silva and R. Abarca. 2009. Water vapour trends at several tropospheric levels over South America between 1973 and 2003. Il Nuovo Cimento, vol. 124 B, N° 8, pp 869-883.         [ Links ]

Parra, J.C., J.A, Sobrino., L, Morales., G, Castellano., J. Uribe, 2006. Application of a Split-Window Algorithm to estimate land surface temperature from AVHRR-NOAA data. Chilean J. Agric. Res., Vol. 66, Nro. 4, p. 385-392.         [ Links ]

Parwati, S., Totok, S., Hasnaeni. 2007. Based on MO-DIS precipitable water for supporting spatial information over Java Island. Remote Sensing and earth Sciences, 4, 33-45.         [ Links ]

Pérez, C. 2004. Técnicas de análisis multivariante de datos, (Ed. Pearson S.A. Madrid, España, 2004).         [ Links ]

Wan, Z. 1999. MODIS Land-Surface Temperature Algorithm Theoretical Basis Document, Institute for Computational Earth System Science University of California, Santa Barbara. Santa Barbara, CA, 93106-3060.         [ Links ]

Wetterhall, F., Halldin, S., Xu, C.Y. 2007. Seasonality properties of four statistical-downscaling methods in central Sweden. Theor. Appl. Climatol. 2007. 87,1-4, pp 123-137.         [ Links ]

Vera, L., J.C, Parra., L, Morales., C, Mattar., E. Jorquera. 2010. Comparative analysis of Split-Window algorithms for estimating soil temperature. R.C. Suelo Nutr. Veg. 10, (1) 35-39.         [ Links ]

Todo el contenido de esta revista, excepto dónde está identificado, está bajo una Licencia Creative Commons