Open access

Climate Changes of the Recent Past in the South American Continent: Inferences Based on Analysis of Borehole Temperature Profiles

Written By

Valiya M. Hamza and Fábio P. Vieira

Submitted: 25 November 2010 Published: 12 September 2011

DOI: 10.5772/23363

From the Edited Volume

Climate Change - Geophysical Foundations and Ecological Effects

Edited by Juan Blanco and Houshang Kheradmand

Chapter metrics overview

2,758 Chapter Downloads

View Full Metrics

1. Introduction

A detailed understanding of the nature of past climate changes is important in assessment of the effects of global warming trends identified in meteorological records (Hansen and Lebedeff, 1987). Nevertheless, there are large uncertainties in the reconstruction of the climate history of times prior to the period of instrumental records, there being considerable difficulties in experimental determination of past climate changes. In this context geothermal methods, based on results of temperature logs in boreholes stands out as one of the few methods that allow direct measurement of thermal signals in the subsurface induced by climate changes of the past.

Geothermal methods have been employed during the last few decades in extracting information on climate changes of the recent past for several regions of the northern hemisphere (e.g. Cermak, 1971; Lachenbruch et al, 1982; Beltrami et al, 1992; Bodri and Cermak, 1995; Duchkov and Sokolova, 1998; Harris and Chapman, 2001; Majorowicz and Safanda, 2001; Golovanova et al, 2001; Beltrami and Burlon, 2004). Nevertheless, very few attempts have been made in using geothermal data for examining climate variations in low latitudes of the southern hemisphere. Among the studies in this category are the works carried out in Australia (Cull, 1979, 1980; Torok and Nicholls, 1996; Taniguchi et al, 1999a, 1999b; Appleyard, 2005), Brazil (Hamza, 1991; Hamza,1998; Hamza et al, 1991; Cavalcanti and Hamza, 2001; Cerrone and Hamza, 2003) and South Africa (Tyson et al, 1998; Jones et al, 1999). Much of the work carried out in Brazil remain as publications of limited access, such as internal reports (Hamza et al, 1978; Eston et al, 1982; Hamza et al, 1987), academic theses (Vitorello, 1978; Araújo, 1978; Santos, 1986; Ribeiro, 1988; Del Rey, 1989; Cavalcanti, 2003) and meeting proceedings (Ribeiro, 1991; Souza et al, 1991; Hamza and Cavalcanti, 2001; Cerrone and Hamza, 2003; and Conceição and Hamza, 2006). In the present work we provide a synthesis of these earlier works, with emphasis on progress obtained during the period 2006 – 2011.

As prelude to the discussion of results presented in this work we provide first a brief description of the sources and characteristics of the geothermal data employed and the criteria used for data selection. Details of the methods used for extracting information on past climate are set out in the next section. The results obtained in model simulations of temperature-depth profiles are classified into groups, representative of the major geographic zones. Finally, the climate history of South America, deduced from geothermal data, is compared with results of geothermal climate reconstructions from other continental areas.


2. Characteristics of the data base

According to the recent compilations carried out by the National Observatory (Observatório Nacional – ON/MCT) in Brazil geothermal measurements have been carried out in over 5000 localities in South America (Hamza et al, 2010; Vieira and Hamza, 2010). Most of the earlier data were acquired as parts of basic research projects for heat flow determinations and also as parts of applied research projects for oil exploration and geothermal energy assessments. The focus of data acquisition in the earlier works has been on determining temperature gradients in the deeper parts of the boreholes.

The characteristics of these data sets are variable, depending on the methods used for primary data acquisition. Of these, only the ones acquired using the so-called conventional (CVL) method provide direct information on the vertical distribution of temperatures at shallow depths and hence are potentially suitable for climate related investigations. The conventional method has been employed for geothermal studies in 134 localities, which is slightly more than 10% of the overall data set. It includes mainly temperature logs in bore holes and wells and thermal property measurements on samples representative of local geologic formations. In some cases estimation of radiogenic heat production was also carried out. The details of the experimental techniques employed for temperature and thermal conductivity measurements have been discussed in academic theses (e.g. Vitorello, 1978; Araújo, 1978; Del Rey, 1989) and in publications dealing with heat flow measurements (Hamza et al, 1987; Hamza and Muñoz, 1996; Gomes and Hamza, 2005). A direct evaluation of the quality of data acquired in the earlier works is a difficult task since the experimental techniques used for temperature and thermal conductivity measurements have undergone substantial changes over the last few decades.

The sources of conventional data sets employed in the present work may be considered as falling into five main groups:

Figure 1.

Localities of geothermal measurements in South America (Vieira and Hamza, 2010). Letters in the legend refer to Koppen climate classification: A – Tropical; B – Dry; C – Temperate; D – Continental; E - Polar (Koppen, 1936; Pidwirny, 2006).

2.1. Selection criteria for climate studies

The characteristics of the conventional (CVL) data set were examined carefully to screen out records with indications of possible perturbations arising from non-climatic effects. Also, it was necessary to eliminate those which do not provide fairly reliable determinations of both the steady and the transient components of the subsurface thermal field. In an attempt to guarantee the reliability of the data set the following quality assurance conditions were imposed:

  1. The depth of borehole is sufficiently large that the lower section of the thermal profile allows a reliable determination of the geothermal gradient, free of the effects of recent climate changes. Order of magnitude calculations indicate that surface temperature changes of the last centuries would penetrate to depths of nearly 150meters, in a medium with a thermal diffusivity of 10-6 m2/s. Thus boreholes of at least 200 m deep are necessary for a reliable determination of the local geothermal gradient. The choice of this depth limit is rather arbitrary since the possibility that low amplitude climate signals of earlier periods are present at larger depths of up to several hundreds of meters cannot entirely be ruled out. However it is a reasonable compromise for examining subsurface thermal effects of ground surface temperature (GST) variations of the last few centuries;

  2. The temperature-depth profile is free from the presence of any significant non-linear features in the bottom parts of the borehole, usually indicative of advection heat transfer by fluid movements, either in the surrounding formation or in the borehole itself;

  3. The time elapsed between cessation of drilling and measurements in boreholes is at least an order of magnitude large compared to the duration of drilling, minimizing thereby the influence of eventual thermal perturbations generated during the drilling activity;

  4. The lithologic sequences encountered in the borehole have relatively uniform thermal properties and are of sufficiently large thickness that the gradient changes related to variations in thermal properties does not lead to systematic errors in the procedure employed for extracting the climate related signal; and

  5. The elevation changes at the site and in the vicinity of the borehole are relatively small so that the topographic perturbation of the subsurface temperature field at shallow depths is not significant.

The sites of these selected boreholes are distributed in the eastern parts of Brazil (in the states of Santa Catarina, Paraná, São Paulo, Minas Gerais, Rio de Janeiro and Bahia), in the Amazon region, central cordillera in Colombia and eastern cordillera in Peru. The majority of the selected temperature logs are from boreholes with depths greater than 200 meters. Some log data for depths less than 200 m were also considered, as these are found to provide complementary information on subsurface temperature fields at shallow depths which may be compared with those encountered in areas where deeper boreholes are situated. Most of the data sets acquired during the decades of 1970 and 1980 have temperature measurements at depth intervals of ten to twenty meters. In more recent logs measurements have been made at intervals of two meters. Such recent logs may be considered as capable of providing more robust estimates of the background temperature gradients. Data from boreholes with depths shallower than 150 meters were not considered in the present work in view of the potential uncertainties in the determination of the local undisturbed gradient and consequent difficulties in extraction of the climate signal. On the other hand, information on climate changes of the recent past available in such logs may be used in obtaining qualitative estimates of GST changes. Data acquired at shallow depths of less than 20 meters were excluded from analyses for climate changes, avoiding thereby eventual perturbing effects of diurnal and seasonal variations in the reconstruction of surface temperature history.

In an earlier study Hamza et al (2007) reported results of geothermal measurements for a number of sites in the Brazilian territory. In the present we have included results of additional studies carried out during the period of 2006 to 2011. The data sets have been classified into groups, designated as subtropical highlands, subtropical humid zones (of the interior and of the coastal areas), tropical Amazon region, semi-arid zones and cordilleran regions in western parts of the continent. The area extents of these geographic sectors are in large part similar to the prevailing climate zones indicated in Figure (1). Typical examples of temperature profiles encountered in boreholes in Brazil, Peru and Colombia are illustrated in Figure (2). In this figure some of the temperature-depth profiles have been shifted laterally to convenient positions along the temperature axis, to avoid overlap and to allow for easy visualization. Consequently, the temperature axis in Figure (2) displays only relative values.

Figure 2.

Vertical distributions of temperatures in selected localities in Brazil, Colombia and Peru. The climate changes of the recent past are often considered responsible for the non-linear features in temperature logs.

The vertical distributions of temperatures in boreholes, some examples of which are illustrated in Figure (2) reveal several remarkable features. For example, the temperature depth profiles at shallow depths are consistently concave towards the depth axis, which is indicative of surface warming events of relatively recent times. The widespread occurrence of such temperature-depth profiles in almost all of the major geographic zones, irrespective of the local geological complexities and changes in soil type, is considered a clear indication that the observed features are generated by surface warming events of large spatial dimensions. Temperature profiles that are convex towards the depth axis, and hence characteristic of cooling events, were not encountered.


3. Methods employed in data analysis

Three different methods have been employed in data analysis: forward models, inversion methods and signal backstripping methods. Brief descriptions of these methods are provided in the following sections.

3.1. Forward models

The basics of forward modeling approach has been discussed extensively in the literature (Birch, 1948; Cermak, 1971; Vasseur et al, 1983; Lachenbruch et al, 1986, 1988). For the case where surface temperature variation can be represented by a power law relation, analytic solutions are readily available. Thus, for a linear (or ramp type) change in surface temperature the relation between the amplitude of the climate signal (ΔT) and the time elapsed (t) at any depth (z) is given by the relation (Carslaw and Jaeger, 1959):


where i2 erfc is the second integral of the complementary error function and κ the thermal diffusivity of the medium. The best fitting ramp function is obtained by inverting the above relation using iterative procedures such as linearized Newton’s method. The iterative procedure for this model, referred to as Ramp Inversion (Chisholm and Chapman, 1992; Golovanova et al, 2001; Roy et al, 2002), allows simultaneous determination of the magnitude of surface temperature change and the period of its occurrence. The main limitation of the forward model approach is that it resolves mainly the first-order features in the GST history. This is a consequence of the implicit assumption that the bottom parts of the log data, employed in determination of background temperature gradients, are free of transient perturbations.

The basic steps in forward model approach include identification and separation of the steady and transient components present in temperature profiles at shallow depths. The steady state component is determined by the flow of heat from the interior of the Earth while the transient component is induced by downward propagation of a climate related thermal signal induced at the surface. Usually a linear fit to the deeper portion of the log data, where the climate perturbation is practically absent, allows determination of the local temperature gradient. However some care is necessary in selecting the depth interval for determination of the gradient. If the temperature gradient is calculated using a small subset of data from the lowermost part of the borehole its standard deviation (σG) is likely to be relatively large, a consequence of the large root mean square (rms) deviation associated with the small number of data points. Progressive inclusion of data from the overlying parts in least square analysis leads to a steady initial decrease in σG, as the estimation of gradient becomes more robust. However, as more data are included from shallower depths (where non-linear features are present) this tendency is reversed and σG increases. In the present work, the depth corresponding to the minimum value of σG is considered as indicative of the top of the unperturbed zone. The background gradient determined for the depth interval below this zone is used for calculating the steady component of the temperature field. Subtracting it from the observed temperatures allows determination of the transient component. A similar procedure has also been employed by Roy et al (2002) in the separation of steady and transient components of temperature profiles in the Indian subcontinent.

3.2. Inverse models

In the inverse problem approach (Tarantola and Valette, 1982) a priori information is explicitly incorporated in constraining the solution. The functional space inversion (FSI) method discussed by Shen and Beck (1991, 1992) and Shen et al (1992) makes use of the non-linear least squares theory in solving the one dimensional heat conduction equation in a layered half space. The algorithm employed finds the model that minimizes the misfit function:


where d and d0 are respectively the calculated and observed temperatures, m and m0 the calculated and a priori model parameters and Cd and Cm the covariance matrices of d0 and m0. The term Cd indicates the uncertainty in the observed temperature-depth data while Cm indicates uncertainty in the a priori model. The selection of appropriate values of a priori standard deviations for the temperature (σd0) and thermal conductivity (σk0) data are important in determining the solutions.

The main advantage of FSI formulation is that it does not predetermine the steady state temperature profile. Instead, both the steady state and transient profiles are estimated simultaneously. In addition, it allows consideration of the vertical distribution of thermo-physical properties and their uncertainties as model parameters, allowing thereby determination of a more detailed GST history where it is possible to identify second order features. Also, FSI formulation includes as model parameters all variables that govern the conductive thermal regime (background heat flow density, the GST history as a function of time and thermal properties as functions of depth).

As prelude to the presentation of the results obtained by the functional space inversion method we provide brief descriptions of the steps taken in data processing and analysis. These include specifying the depth intervals and time periods used for inversion, setting a priori standard deviations of temperature and thermal conductivity data sets and measures taken to minimize the undesirable consequences of the null hypothesis in the inversion scheme.

The estimates of the depth at which the thermal regime is supposedly untouched by the GST variations and the time limit beyond which GST variations cannot be recovered from the given temperature log data, were chosen in accordance with the depth extent of the available temperature log data. In particular, the depth estimate is set to be greater than the deepest data point because the calculated data are projected (interpolated) from the finite element solution. As for the time limit there is no harm in setting a value compatible with or greater than the depth extent of the borehole (Shen and Beck, 1992). On the other hand, use of a shorter-than-necessary time span would end up in “telescoping” the GST history. Unless there are independent evidences indicating a rapid return to unperturbed conditions it seems prudent to assume that this return take place gradually. Use of shorter time spans leads to slight reductions in the magnitude and duration of the cooling events of the earlier periods. As pointed out by Shen and Beck (1991; 1992) the results of GST history, determined by functional space inversion, is sensitive to a priori standard deviations of thermal conductivity. The preferred values for the standard deviations are based on considerations of the trade-off between consistency of the solution and data resolution. In the present work we have used 50mK for standard deviation of the temperature data (σd0) and 1W/m/K for standard deviation of the thermal conductivity (σk0). These values are of the same order of magnitude as those adopted by Safanda and Rajver (2001) and Golovanova et al (2001).

FSI inversion makes use of an a priori null hypothesis in obtaining robust estimates of the prior steady state. However, this built-in feature can potentially lead to undesirable results when inversions are attempted for determining GST history of periods not correlated with the subsurface temperature data. Thus, in the absence of suitable temperature data for shallow depths the inversion scheme generates artificial values for the late part of GST history. It is clear that acquisition of reliable temperature data for shallow depths is important in determining GST history of the last few decades. On the other hand, the results also show that occurrence of artificial cooling trends for the decades prior to 1970 are possible only in cases where temperature measurements are restricted to depths greater than 100m.

3.3. Signal back stripping approach

Both the conventional and the Bayesian inversion methods have inherent difficulties in identification of individual thermal signals originating from climate variations that are episodic. In such cases methods based on signal back stripping approach are more convenient. Following the standard practice we also assume that the residual temperature profile represents is a superposition of individual perturbations. The characteristics of such perturbations vary, mainly as a result of differences in the magnitudes and time periods of individual GST episodes. The essence of the procedure adopted in the present work can be understood by considering magnitudes of temperature perturbations at two conveniently selected depths z1 and z2 in the residual profile:


Dividing (3a) by (3b) and designating the ratio dT1/dT2 by δ we have:


Note that equation (4) does not depend on the magnitude of the temperature perturbation but only on the selected value of the time period. Iterative methods may now be employed for determining the appropriate value of δ for the depth interval. The magnitude of this perturbation can be obtained from the relation:


In applying this procedure it is necessary to start with results of the lowermost section of the residual temperature profile. If the borehole is sufficiently deep it is fairly reasonable to assume that the residual temperatures of the deeper parts retain the effects of only the earliest perturbation. The magnitude and period of this earliest perturbation can be determined through the use of equations (3) and (5). It also opens up the possibility of removing the effects of this particular perturbation from the original residual temperature profile. The result is a back stripped profile free of the effects of the earliest perturbation. The procedure is repeated successively for identifying and removing the perturbations arising from later climate episodes.


4. Estimates of surface temperature changes

4.1. Forward model results

A summary of the results obtained in fitting forward models to the observational data discussed in this work, is presented in Table (1). It includes magnitudes of the GST change (ΔT) and their duration (t) as well as the values of the root mean square (rms) misfit between the model and the observational data. For reasons of brevity, we present here only the results for the ramp function model. Climate changes inferred on the basis of this model are grouped together for the five major geographic zones: subtropical highlands (elevations >400 meters) subtropical lowlands (elevations <400m), humid regions, semi arid zones, tropical rain forest and cordilleran regions.

The data in table (1) reveal differences in GST between the major geographic zones. For example, the highland areas seem to be characterized by changes in GST in the range of 2 to 3.8oC, the mean value being 3.3oC. The lowlands of the interior also seem to have nearly similar GST changes. On the other hand, the humid zones have relatively low GST values, in the range of 0.4 to 1.8oC, with a mean of 1.1oC. Low GST values are also found in the semiarid zones in the northeast, in the interval of 1 to 2oC. Also, the times of duration of climate changes are different for the geographic zones. Thus, the subtropical highlands, subtropical lowlands and humid zones are characterized by climate episodes with durations varying from 40 to 120yrs, the mean value being 80 yrs. The semi arid zones on the other hand have climate variations with relatively larger duration, the mean value being 150yrs.

The vertical distributions of transient components derived from forward model results (which are considered as indicative of temperature signals resulting from climate changes) are illustrated in Figure (3). Note that the magnitudes of the transient components decrease rapidly with depth. At depths greater than 150 meters these values fall below the experimental detection limits of temperature changes in boreholes. The period of GST changes fall in the range of 130 to 260 years. Hence the beginning of the climate change seems to have taken place during the period of approximately 1750 to 1850. A closer examination of the results illustrated in Figure (3) reveals some marked differences in the GST values within individual geographic zones. The primary reason for the occurrence of such intra-zonal variations is unknown at the moment, but it is likely that they are related to microclimatic histories of individual sites.

Also, the depth distribution of the transient components indicates that the magnitude of warming event is relatively smaller (in the range of 1.4 to 2.20C) for the semi-arid zone in north eastern parts of Brazil. The close agreement between the results for the localities in Table (1) is considered as indication that local changes in vegetation cover and soil types have only a minor influence on the surface thermal budget in semi-arid zones. The GST change in this region seems to have had its beginning during the time period of 1850 to 1900, significantly earlier than the corresponding periods for other geographic zones.

Climate ZoneLocalityCoordinatesΔT
Climate variationrms
DurationYear of Onset
Subtropical HighlandsÁ. Lindóia220 29’/460 38’3.210518774.4
Amparo22° 43’/46° 46’2.811518525.5
Araras22° 21’/47° 22’3.44019426.4
Cosmópolis230 43’/470 12’3.87519075.5
Itapira220 28’/460 43’3.840194210.2
Rafard23° 00’/47° 31’3.86019226.8
Jacarei230 18’/450 57’2.410518803.5
São Paulo23° 34’/46° 44’2.110519057.2
Serra Negra22° 36’/46° 32’2.010518773.9
P. Caldas210 55’/460 25’2.050192613.7
Teresópolis22° 26’/ 42° 57’2.05519455.5
L. Muller28º 40´/49º 30´2.680189510.0
Subtropical LowlandsSeropédica22º 46´/43º 39´3.8145186023.1
Miracema22º 01'/41º 06'3.63019697.3
Campos21º 46'/41º 17'3.02019806.7
Itapemirim19° 06’/410 04’3.665191014.1
S. Sebastião23° 48’/45° 25’2.711018807.8
Subtropical HumidItu230 15’/470 19’1.26019223.2
Jundiaí230 10’/460 52’1.25019327.2
Papanduva260 23’/500 08’1.290188515.0
Cach. Sul30º 00´/52º 55´1.870190512.0
Maricá22° 54’/42° 45’0.41519852.9
Rio Bonito21° 25’/42° 12’1.84019607.7
Semi-AridArraial12° 32’/42° 50’1.48518908.4
Caraiba09° 28’/39° 50’1.411018659.3
Jacobina11° 11’/40° 31’1.915018288.7
Poço Fora090 41’/390 51’1.810518704.9
Tropical Rain ForestManaus02º 57´/60º 01´1.340196015.0
Belém01º 27´/48º 27´2.25519539.4
Salinópolis00º 38´/47º 20´2.16519438.6
Dom Eliseu04º 17´/47º 34´1.96019489.2
Cordilleran RegionsPen 74281º 10´/04º 17´4.020017797.4
Lobitos81º 16´/04º 27´2.520017798.3
Lomitos70º 39´/17º 16´2.020017797.2
Bogotá74º 05´/04º 38´2.815018609.5

Table 1.

Results of Ramp Inversions of GST changes, for selected localities in South America. ΔT is the magnitude of ramp change, and rms the root mean square misfit.

Figure 3.

Magnitudes of climate changes and its vertical distributions, deduced from forward model fits to borehole temperature profiles at selected localities in South America.

A careful examination of figure (3) reveals that some of the transient temperature profiles do have small negative values for depth intervals corresponding to the lower parts of the transient sections. The negative residuals are usually considered as arising from cooling events prior to recent warming episodes. However, the magnitudes of such events appear as subdued features, because of the implicit assumption in forward model approach that the bottom parts of the borehole are free of transient perturbations. Results of numerical simulations indicate that occurrence of negative residuals is quite sensitive to subtle changes in the value adopted for the background temperature gradient. This is the main limitation of the forward models, which resolves mainly for the first-order features in the GST history. We conclude by noting that residual temperature profiles similar to those obtained in the present work were also reported by Golovanova et al (2001) for the Urals region and Roy et al (2002) for the Indian Peninsula. A complementary analysis of this problem is provided in section (4.2) below, where we discuss the residual temperature profiles derived using the inversion method.

4.2. Results of functional space inversion

The method of Functional space inversion (FSI) was employed in analysis of temperature-depth profiles for 14 localities distributed over three main geographic zones of Brazil. The criteria used in the selection of profiles included availability of both thermal property data of subsurface layers and supplementary information on the history of changes in the vegetation cover. In discussing the results it is important to point out that the FSI method provides more detailed information on the GST history than that provided by the forward model approach. However, in comparing the GST histories of several localities it is convenient to work with deviations from the site specific mean rather than the absolute value. In the present work, GST deviations are calculated by subtracting the model results from the a posteriori estimate of the site specific mean.

The results of FSI method, illustrated in figure (4), reveal several characteristic features in the GST history of the study area. Foremost among these are the indications that surface temperatures have increased by as much as 1 to 4oC, over the last century. This observation is in reasonable agreement with the results of the ramp function model discussed in the previous section. However, in all localities the warming events seem to be preceded by cooling episodes occurring over the time period of approximately 1700 to 1900. The amplitudes of the cooling events are much less, falling generally in the range of 0.5 to 1oC. For time periods prior to the 17th century the resolving power of the inversion method is poor, a consequence of the limitations in sensitivity and precision of sensors used for temperature measurements in boreholes. Hence FSI model calculations for periods prior to 1700 may not necessarily be representative of true climate history.

Even though the data set is poor the results of Figure (4) seem to indicate that the warming trends are less pronounced in semi-arid regions relative to the tropical humid regions. The results also reveal a small time shift in the occurrence of climate warming in highlands regions compared with that for semiarid regions. For example, the warming event in the highlands region had its beginning during the time period of 1850 to 1900 while that for semi-arid regions seems to have had its beginning during the period of 1670 to 1860. Also the durations of both warming and cooling episodes appear to be relatively smaller in highland regions when compared with those for the semi-arid regions. Such differences have important implications for understanding evolutionary trends in climate history of eastern Brazil.

The vertical distributions of transient components derived from FSI inversions are similar to those found for the forward model results, illustrated in Figure (3). The magnitudes of the transient components decrease rapidly with depth. However, the residual temperatures become significantly negative at depth intervals corresponding to the lower parts of the transient sections. Such negative residuals are indicative of the occurrence of cooling events prior to recent warming episodes. We conclude that the vertical distributions of residual temperatures in FSI method are better representations of subsurface transient thermal regimes than those derived from the forward model approach.

Another conspicuous feature of the results by the FSI method in Figure (4) is the presence of brief cooling episodes for the recent decades, since 1970. There is a possibility that short period cooling episodes are spurious, a consequence of the null hypothesis employed in the computational process of the FSI method (Shen and Beck, 1992). The potential undesirable effects of the null hypothesis may be minimized by limiting the GST history to periods compatible with reliable subsurface temperature data. In the present work temperature data are available for depths less than 50 meters in all of the data sets employed in inversion schemes. Consequently the cooling trends observed in our results, for the decades prior to 1970, cannot be attributed to the spurious effects of the null hypothesis. Examples of cooling trends similar to those found in the present work can also be seen in the results reported for several localities in Europe and North America (Bodri and Cermak, 1995; Rajver et al, 1998; Stulc et al, 1998; Golovanova et al, 2001).

Figure 4.

History of ground surface temperature variations associated with climate changes of the recent past, deduced on the basis of functional space inversion models fits to borehole temperature profiles, for selected localities in South America (Adapted with modifications from Hamza et al, 2007).

The uncertainties in the estimates of GST in the inversion method can somewhat be improved by carrying out simultaneous inversion of temperature profiles of several sites in the same geographical province. In the present work we have carried out simultaneous inversions only for repeat measurements at Seropédica, in the state of Rio de Janeiro. In view of notable differences in the local soil conditions and vegetation cover, no attempt has been made for carrying out simultaneous inversions of temperature profiles from different locations discussed in the present work. Also, most of these sites are separated by large distances and fall within areas with distinctly different microclimate conditions and geographic characteristics.

A summary of the results of inverse modelling is presented in Table (2). The summary includes maximum and minimum values of ground surface temperatures and their respective times of occurrence. Also given in this table are the difference between the maximum and minimum values of GST, the time elapsed between the maximum and minimum and aposteriori estimates of undisturbed GST. The values of the differences in magnitudes are in the range of 1.3 to 30C, which is slightly lower than the range indicated by the results of the forward model approach. On the other hand, the values of the differences in the time periods are in the range of 62 to 147 years, which is comparable to the ranges indicated by the results of the forward model approach.

LocalityMaximumMinimumDifferenceT0 (°C)Log Year
C. Itapemirim23.7197120.018833.78821.21976
Poço de Fora30.4197628.218222.115428.71976

Table 2.

Results of functional space inversions for selected data sets. Maximum and minimum values of GST variations and their respective years of occurrences are given in columns 2 - 5. Also given are model results for GST deviations and the respective periods. T0 is the aposteriori steady state temperature.

4.3. Results of signal back-stripping approach

Techniques of signal back-stripping were used in analysis of temperature log data from 15 localities in south east Brazil. In the present work we limit the discussion to the results obtained for the locality of Seropédica (Rio de Janeiro). The vertical distributions of temperatures obtained during the various stages of the back stripping process are illustrated in the set of upper and lower panels of Figure (5).

The segment on the left side of the upper panel in this figure refers to results of the first stage of the back stripping process. Here the blue curve indicates the initial reduced temperatures and the red curve the first signal extracted by the back stripping method. The extracted signal points to an episode of climate change that took place 115 years back in time and had a magnitude of 2.90C. Similarly, the segments on the right side of the upper panel and that on the left side of the lower panel of this figure illustrate, respectively, the results of the second and third stages of the back stripping process. According to the results obtained the signal extracted in the second stage points to an episode of climate change that took place 93 years back in time and had a magnitude of 0.20C. On the other hand, the signal extracted in the third stage points to an episode of climate change that took place 21 years back in time and had a magnitude of 0.80C. The segment on the right side of the lower panel illustrates the residual temperatures, after extracting the above mentioned three climate signals. The back stripping process is terminated at this third stage since the magnitudes of the residual temperatures are less than 0.050C, below the sensitivity limit of experimental system used in data acquisition. The overall history of climate change derived from the back stripping method is illustrated in Figure (6).

Figure 5.

Results of signal back stripping method employed in disentangling the climate history of Seropédica (Rio de Janeiro).

Figure 6.

Climate history of Seropedica (Rio de Janeiro) derived from results of the signal back stripping method.


5. Discussion

5.1. Climate changes in South America

The results obtained in the present work have contributed to substantial improvements in data base for climate change in South America. The area covered in the data base includes several geographic zones in Brazil, Colombia and Peru. Classical inverse models were employed in the analysis of temperature logs from over 30 localities and, in addition, Bayesian inverse modelling was carried out for data from 20 selected sites. The model results have allowed determination of the magnitude as well as the duration of ground surface temperature changes in the major geographic zones of South America. The map of recent climate changes derived from geothermal data is presented in Figure (7).

According to the results obtained the magnitude of GST changes are in the range of 2 to 3.5oC but have had their beginning during the early decades of the 20th century. Nearly similar trends are seen in temperature-depth profiles of bore holes in tropical as well as subtropical zones of the interior and coastal areas. The data from semi arid zones also indicate occurrence of surface warming events but the magnitudes are in the range of 1.4 to 2.2oC while the duration of the warming event is larger, extending back into the last decades of the 19th century. The magnitudes of GST variations are relatively large in localities which have undergone recent changes in vegetation cover. Also there are indications that GST changes are practically insignificant in areas of tropical rain forest.

5.2. Global climate changes

The improved data base on GST changes in the South American continent has contributed to a better understanding of climate changes in the southern hemisphere. In this context, it is convenient to examine also its implications for global variations. It is clear that improved assessments of global changes require an integrated analysis of both the geothermal data sets as well as the meteorological data sets (such as HadCRUT3) reported for oceanic and continental regions (Brohan et al, 2006; Jones and Moberg, 2003; Rayner et al, 2003). Major obstacles in such a venture springs from the wide disparities in data density, and it is

Figure 7.

Map of recent changes in ground surface temperatures derived from geothermal data sets, and supplemented with ground data for selected localities.

convenient to adopt procedures that minimize problems arising from the non-homogeneous distribution. In the present case, the surface area of the globe is divided into a regular grid system composed of 5 x 5 degree grid cells and average values of surface temperature variation in the grid elements calculated. Experimental data are available for almost all of the grid elements covering the continental regions and a significant part of the grid system for the oceanic regions. Following the common practice employed in deriving maps of global representations interpolated values were used for grid elements without data.

Techniques of spherical harmonic representation (Hamza et al, 2008) were employed in analysis of global variations. The harmonic representation of surface temperature variations (q) may be represented as:


where is the longitude = 90 – ψ, is the colatitude, P´nm(cos ) is the associated Legendre function that is fully normalized and Anm and Bnm the coefficients of the harmonic expansion. The expression for evaluation of P´nm is:


where Pnm is the associated Legendre function given by:


In equation (9) Int(nm/2)is the largest integer that is lower than(nm)/2.

Full normalization of associated Legendre functions (Pnm) requires that the following equations be satisfied:


The coefficients Anm and Bnm are evaluated by fitting the harmonic expansion to the set of experimental data, which are the values of surface temperature changes (q) and their respective geographic coordinates ( and ).

The results obtained on the basis of harmonic representation of the mixed data set are illustrated in the global map of surface temperature changes, in Figure (8). According to features discernible in this figure the thermal signals of climate change in continental areas of North America, Europe, Asia, West Africa and Eastern South America have magnitudes in excess of one degree centigrade, while those in most of the oceanic regions is less than one degree centigrade. On the other hand most of the oceanic regions are characterized by surface temperature variations of less than 0.80C.


6. Conclusions

Much of the information on surface temperature changes of the recent past in South America are derived from historical records of local meteorological stations. Some of the

Figure 8.

Spherical harmonic representation of surface temperature variations associated with climate changes of the recent past. Note the systematic differences in magnitudes of climate change between continental and oceanic areas.

main problems in the use of such records include non-uniformity in the accuracy and precision of sensors employed in data acquisition, gaps in the time series of observations and influence of urban warming effects (most of the meteorological observatories are located in or close to towns or cities). The results reported in the present work, based on geothermal methods, are in general free of such problems.

An important result emerging from studies based on geothermal methods is that the climate was relatively cooler during the 17th and 18th centuries. The climate histories, deduced from geothermal data, are found to be consistent with the results of available meteorological records in Brazil, Colombia and Peru. Comparative studies indicate that the magnitudes and duration of recent climate changes in South America are similar to those found in other continental areas such as North America, Asia and Europe.

Analysis of surface temperature changes on global scale reveal that climate warming effects in the northern hemisphere are more pronounced because of the relatively large areas of continental blocks in these regions. On the other hand, the climate warming effects are relatively subdued in the southern hemisphere, where relatively larger proportions of oceans dominate. It is possible that systematic difference in magnitudes of global warming between continental and oceanic regions is a consequence of large scale mixing of near surface waters with those of deep ocean circulation systems. This mixing process is considered to be responsible for pronounced mitigation of the effects of global warming in the southern hemisphere.



This work was carried out as part of a research project with funding from Conselho Nacional de Desenvolvimento Científico – CNPq (Project No. 301865/2008-6; Produtividade de Pesquisa - PQ). We thank the topic editor and reviewers for critical comments and suggestions. Dr. Andres Papa of the Geophysics Department (ON-MCT, Rio de Janeiro) provided institutional support. The source code for functional space inversion program was provided by Dr. Paul Shen, University of Western Ontario, Canada


  1. 1. AlfaroC.AlvaradoI.QuinteroW.HamzaV. M.VargasC.BriceñoL. A.2009Preliminary map of geothermal gradients of Colombia., XII Congreso Colombiano de Geologia, 711September, Paípa- Boyacá, Columbia.
  2. 2. AppleyardS. J.2005Late Holocene temperature record from southwestern Australia: evidence of global warming from deep boreholes,Australian J. of Earth Sciences, 52161166
  3. 3. AraújoR. L. C.1978Heat flow measurements in the Alkaline intrusive complex of Poços de Caldas (in Portuguese). Unpublished M.Sc. Thesis, University of São Paulo, São Paulo, Brazil.
  4. 4. BeltramiH.BourlonE.2004Ground warming patterns in the Northern Hemisphere during the past five centuries, Earth Plnet. Sci. Lett., 227169177
  5. 5. BeltramiH.JessopA. M.MareschalJ. C.1992Ground temperature histories in eastern and central Canada from Geothermal measurements: Evidence of climatic changePalaeogeogr. Palaeoclimatol. Pallaeoecol., 98167184
  6. 6. BirchF.1948The effects of Pleistocene climatic variations upon geothermal gradientsAm. J. Sci., 246729760
  7. 7. BodriL.CermakV.1995Climate changes of the last millennium inferred from borehole temperatures: results from the Czech Republic- Part I, Global and Planetary Change, 11111125
  8. 8. BrohanP.KennedyJ. J.HarrisI.TettS. F. B.JonesP. D.2006Uncertainty estimates in regional and global observed temperature changes: a new dataset from 1850. J. Geophysical Research 111, D12106, doi:10.1029/2005JD006548.
  9. 9. CarslawH. S.JaegerJ. C.1959Conduction of Heat in Solidspp, Oxford University Press, New York.
  10. 10. CavalcantiA. S. B.2003Palaeoclimate variations in Brazil based on the geothermal method (in Portuguese), Unpublished M.Sc. Thesis, Observatório Nacional, Rio de Janeiro, Brazil.
  11. 11. CavalcantiA. S. B.HamzaV. M.2001Climate changes of the recent past in southern parts of Brazil (Extended abstract), 7th International Congress of the Brazilian Geophysical Society, Salvador (BA), Brazil, 14
  12. 12. CermakV.1971Underground temperature and inferred climatic temperature of the past millennium, Palaeogrogr., Palaeoclimatol., Palaeoecol., 10119
  13. 13. CerroneB. N.HamzaV. M.2003Climate changes of the recent past in the state of Rio de Janeiro, based on the geothermal method (Extended abstract in Portuguese), 8th International Congress of the Brazilian Geophysical Society- SBGf, Rio de Janeiro, Brazil, 14
  14. 14. ChisholmT. J.ChapmanD. S.1992Climate change inferred from analysis of borehole temperatures: an example from western UtahJ. Geophys. Res., 971415514175
  15. 15. ConceiçãoE.HamzaV. M.2006Subsurface thermal signals of the greenhouse effect in the state of São Paulo (Extended Abstract in Portuguese), 2nd International Symposium of the Brazilian Geophysical Society- SBGf, Natal (RN), Brazil, 16
  16. 16. CullJ. P.1979Climatic corrections to Australian heat flow data. BMR Journal of Australian Geology and Geophysics, 4303307
  17. 17. CullJ. P.1980Geothermal records of climate change in New South Wales, Search, 11(6), 201-203.
  18. 18. Del ReyA. C.1989Hydrogeothermal studies of the regions of Águas de Lindoia, Amparo e Socorro- Northeastern parts of the state of São Paulo (in Portuguese), Unpublished M.Sc. Thesis, University of São Paulo, São Paulo, Brazil.
  19. 19. Del reyA. C.HamzaV. M.1989Terrestrial heat flow variations in the northeastern part of the state of São Paulo: A case for transport of geothermal heat by inter-fracture fluid flows. In: A.E. Beck, G. Garven and L. Stegena (Eds.), Hydrogeological regimes and their subsurface thermal effects. Geophysical Monograph, 47, IUGG 2American Geophysical Union, Washington, 137-148.
  20. 20. DuchkovA. D.SokolovaL. S.1998Investigation of temperature changes at the Earth’s surface via borehole geothermometry, In Problems of climatic reconstruction and environment of the Holocene and Pleistocene in Siberia (in Russian), 1151157Nauk, Novosibirsk, Russia.
  21. 21. EstonS. M.HamzaV. M.BeckerE. A.FurumotoS.1982Geothermal research in exploration of hydrocarbons in the Paraná basin (in Portuguese), Internal Report, 18106Instituto de Pesquisas Tecnológicas, São Paulo, Brazil.
  22. 22. GolovanovaI. V.HarrisR. N.SelezniovaG. V.StulcP.2001Evidence of climate warming in the southern Urals region derived from borehole temperatures and meteorological data, Global and Planetary Change, 29167188
  23. 23. GomesA. J. L.2003Assessment of Geothermal Resources of the state of Rio de Janeiro (in Portuguese), Unpublished M.Sc. Thesis, National Observatory, Rio de Janeiro (Brazil). 138pp.
  24. 24. GomesA. J. L.HamzaV. M.2005Geothermal gradients and heat flow in the state of Rio de Janeiro, Brazilian J. of Geophysics, 234325348
  25. 25. HamzaV. M.1982Terrestrial heat flow in the alkaline intrusive complex of Poços de Caldas, Brazil. Tectonophysics, 834562
  26. 26. HamzaV. M.1998A proposal for continuous recording of subsurface temperatures at the sites of Geomagnetic field Observatories. Rev. Geofísica, 48183198
  27. 27. HamzaV. M.1991Recent climate changes in the southern hemisphere: the geothermal evidence (In Portuguese), Proceedings 2nd Congress of the Brazilian Geophysical Society, Salvador (BA), 971973
  28. 28. HamzaV. M.2006Thermal Signals of Global Warming: Evaluation using geothermal methods (in Portuguese), IV Symposium of the Amzon Region, 1214September, Belém (PA), Brazil.
  29. 29. HamzaV. M.2007Climate changes in the eastern segment of the Amazon region: Inferences based on geothermal methods (in Portuguese), V Symposium of the Amzon Region, 30 October to 1 November, Belém (PA), Brazil.
  30. 30. HamzaV. M.AlfaroC.AlvaradoI.QuinteroW.LópezJ.MonsalveM. L.PulgarínB.MadridC.2009Surface Temperature Variations of the Last Two Centuries in Bogotá, Columbia., XII Congreso Colombiano de Geologia, 711September, Paípa- Boyacá, Columbia.
  31. 31. HamzaV. M.CardosoR. R.VieiraF. P.GuimaraesS. N.2010Geothermal Maps of Brazil (in Portuguese), National Atlas of Brazil- 2010, IBGE.
  32. 32. HamzaV. M.CavalcantiA. S. B.2001Thermal imprints of changes in vegetation cover and climate patterns in borehole temperature profiles in Brazil (Abstract), IAGA-IASPEI Assembly, Hanoi, Vietnam, Book of Abstracts, 413414
  33. 33. HamzaV. M.CavalcantiA. S. B.BenyosefL. C.2007Surface thermal perturbations of the recent past at low latitudes-Inferences based on borehole temperature data from Eastern Brazil,Climates of the Past, 3113
  34. 34. HamzaV. M.CerroneB. N.GomesA. J. L.CardosoR. A.2003A Geothermal Climate Change Observatory in the Southern Hemisphere (Abstract), Proceedings of XXIII General Assembly, IUGG, Sapporo, Japan, v. A, 197
  35. 35. HamzaV. M.EstonS. M.AraujoR. L. C.VitorelloI.UssamiN.1978Brazilian Geothermal Data Collection- Series-1, Publication IPT 11091316
  36. 36. HamzaV. M.FrangipaniA.BeckerE. A.1987Geothermal Maps of Brazil (In Portuguese) Internal Report 25305Instituto de Pesquisas Tecnologicas, São Paulo, Brazil.
  37. 37. HamzaV. M.MuñozM.1996Heat flow map of South AmericaGeothermics599 EOF621 EOF
  38. 38. HamzaV. M.RibeiroF. B.BeckerE. A.1991Recent Climatic Changes in the Southern Hemisphere. XX General Assembly, IUGG, Vienna (Austria).
  39. 39. HamzaV. M.SilvaDias. F. J. S.GomesA. J. L.TercerosZ. G. D.2005Numerical and functional representations of regional heat flow in South AmericaPhysics of the Earth and Planetary Interiors152223256
  40. 40. HansenJ. E.LebedeffS.1987Global trends of measured surface air temperatureJ. Geophys. Res., 92, 13345-13372, doi:10.1029/JD092iD1113345
  41. 41. HarrisR. N.ChapmanD. S.2001Mid-Latitude (30-60N) climatic warming inferred by combining borehole temperatures with surface air temperatures, Geophys. Res. Lett., 28747750
  42. 42. HuangS.PollackH. N.1998Global Borehole Temperature Database for Climate Reconstruction. IGBP PAGES/World Data Center-A for Paleoclimatology Data Contribution Series #1998044NOAA/NGDC Paleoclimatology Program, Boulder CO, USA.
  43. 43. JonesP. D.MobergA.2003Hemispheric and Large-Scale Surface Air Temperature Variations: An Extensive Revision and an Update to 2001J. Clim., 1620623
  44. 44. JonesW. Q. W.TysonP. D.CooperG. P. R.1999Modeling climate change in South Africa from perturbed borehole temperature profiles, Quaternary International, 57/58185192
  45. 45. KoppenW.1936Das geographisca System der Klimate, in: Handbuch der Klimatologie, edited by: Koppen, W. and Geiger, G., 1. C. Gebr, Borntraeger, 144
  46. 46. LachenbruchA. H.ClaudouhosT. T.SaltusR. W.1988Permafrost temperature and the changing climate, Proc. Fifth International Conference on Permafrost, Trondheim, Norway, 19
  47. 47. LachenbruchA. H.SassJ. H.MarshallB. V.1986Changing Climate: Geothermal Evidence from Permafrost in the Alaskan Artic, Science, 234689696
  48. 48. LachenbruchA. H.SassJ. H.MarshallB. V.MosesT. H.Jr 1982Permafrost, heat flow and geothermal regime at Prudhoe Bay, Alaska, J. Geophys. Res., 8793019316
  49. 49. MajorowiczJ. A.SafandaJ.2001Composite surface temperature history from simultaneous inversion of borehole temperatures in western Canadian plainsGlobal and Planetary Change29231239
  50. 50. PidwirnyM. (2006). "Climate Classification and Climatic Regions of the World". Fundamentals of Physical Geography, 2nd Edition.
  51. 51. RajverD.SafandaJ.ShenP. Y.1998The climate record inverted from borehole temperatures in SloveniaTectonophysics291263276
  52. 52. RaynerN. A.ParkerD. E.HortonE. B.FollandC. K.AlexanderL. V.RowellD. P.KentE. C.KaplanA.2003Globally complete analyses of sea surface temperature, sea ice and night marine air temperature, 18712000J. Geophysical Research 108, 4407, doi:10.1029/2002JD002670.
  53. 53. RibeiroF.1988Study of the problem of determination of equilibrium temperatures in boreholes (in Portuguese), Unpublished Ph.D. Thesis, University of São Paulo, São Paulo, Brazil, 187pp.
  54. 54. RibeiroF.1991Reconstruction of the climate change from geothermal profiles: Advances in Theory (in Portuguese), Proceedings 2nd Congress of the Brazilian Geophysical Society, Salvador (BA), 959964
  55. 55. RoyS.HarrisR. N.RaoR. U. M.ChapmanD. S.2002Climate change in India inferred from geothermal observationsJ. Geophys. Res., 107, B75116
  56. 56. SafandaJ.RajverD.2001Signature of last ice age in the present subsurface temperatures in the Czech Republic and Slovenia,Global and Planetary Change29241257
  57. 57. SantosJ.1986Heat flow density measurements in the northern parts of the Paraná basin (in Portuguese), Unpublished M.Sc. Thesis, University of São Paulo, São Paulo, Brazil, 124pp.
  58. 58. SantosJ.HamzaV. M.ShenP. Y.1986A method for measurement of terrestrial heat flow density in water wells, Brazilian J. Geophysics, 44553
  59. 59. ShenP. Y.BeckA. E.1991Least square inversion of borehole temperature measurement in functional space.J. Geophys. Res., 961996519979
  60. 60. ShenP. Y.BeckA. E.1992Palaeoclimate and heat flow density inferred from temperature data in the Superior Province of the Canadian Shield. Global and Planetary Change, 98143165
  61. 61. ShenP. Y.WangK.BeltramiH.MareschalJ.C.1992A comparative study of inverse methods for estimating climatic history from borehole temperature dataPalaeogeographyPalaeoclimatology, Palaeoecology, 98113127
  62. 62. SouzaJ. R. S.AraújoR. L. C.MakinoM.1991Transient variations in shallow thermal profiles associated with microclimatic changes in the Amazon (in Portuguese), Proceedings 2nd Congress of the Brazilian Geophysical Society, Salvador (BA), 965970
  63. 63. StulcP.GolovanovaI. V.SelezniovaG. V.1998Climate change record in the Earth- Example of borehole data analysis in the Urals region, Russia, Phys. Chem. Earth, 2311091114
  64. 64. TaniguchiM.WilliamsonD. R.PeckA. J.1999aDisturbances of temperature-depth profiles due to surface climate change and subsurface water flow: 2. An effect of step increase in surface temperature caused by forest clearing in southwestern Australia, Water Resources Research3515191529
  65. 65. TaniguchiM.WilliamsonD. R.PeckA. J.1999bEstimations of surface temperature and subsurface heat flux following forest removal in the south-west of Western AustraliaHydrological Processes1222052216
  66. 66. TarantolaA.ValetteB.1982Inverse problems = quest for information. J. Geophys., 50159170
  67. 67. TorokS.NichollsN.1996An historical temperature record for Australia. Aust. Met. Mag. 45251260
  68. 68. TysonP. D.MasonS. J.JonesM. Q.CooperG. R.1998Global warming and geothermal profiles: the surface rock-temperature response in South AfricaGeophys. Res. Lett., 2527112713
  69. 69. VasseurG.BernardPh.van de MentebrouchJ.KastY.JolivetJ.1983Holocene palaeotemperatures deduced from geothermal measurements, Paleogeography, Paleoclimatology, Palaeontology, 43237259
  70. 70. VieiraF. P.HamzaV. M.2010Global Heat Loss: New Estimates Using Digital Geophysical Maps and GIS Techniques. IV Symposium of Brazilian Geophysical Society, Brasília, 14 a 17 de novembro de 2010.
  71. 71. VitorelloI.1978Heat flow and radiogenic heat production in Brazil, with implications for thermal evolution of continents.Ph. D. dissertation, Univ. of Michigan, Ann Arbor, U.S.A.
  72. 72. VitorelloI.HamzaV. M.PollackH. N.AraújoR. L. C.1978Geothermal Investigations in Brazil, Rev. BrasGeof., 8(2), 71-89.

Written By

Valiya M. Hamza and Fábio P. Vieira

Submitted: 25 November 2010 Published: 12 September 2011