Long-term (45 years) diversified surface meteorological records from Thumba Equatorial Rocket Launching Station (TERLS), India, were collected and analysed to study the long-term changes in the overall climatology, climatology pertained to a particular observational time, mean daily climatology in temperature, inter-annual variability in temperature, interannual variability in surface pressure, and rainfall for the main Indian seasons—South West and North East monsoons and inter-annual mean monthly anomaly structure in temperature. Results on various analyses show strong and vivid features contributed by climate change for this South Peninsular Indian Arabian Sea Coastal Station, and this paper may be a first time venture which discusses climate change imparted perturbations in several meteorological parameters in different time domains, like a specific time, daily, monthly, and interannually over a station. Being a coastal rocket launching station, climatic change information is crucial for long-term planning of its facilities as well as for various rocket range operational demands.

1. Introduction

Climate change refers to a statistically strong and significant variation in either the average state of the climate or in its variability, persisting for an extended period, typically in decades or longer. Short-period oscillations are statistically insignificant in the scenario of long-term climate change context. Climate change may be due to natural internal processes on earth (atmospheric, seismic, volcanic, and oceanic), external forces (variation in solar activity, rotation and revolution of earth), and more recently anthropogenic activities. It is now well established that anthropogenic activities cause intensification of the greenhouse effect and thereby contribute to the climate change. The paper by Oreskes [1] highlighted the listings in the ISI database on climate change from several abstracts published in refereed scientific journals between 1993 and 2003.

In India, the first time document on the concern of greenhouse effect and possible mitigation measures was published in the Science Reporter, CSIR, by Ayyar [2]. Several papers were published in a special session in Current Science (volume 90, 2006) focussing the concerns over the socioeconomic and scientific aspects attributed to climate change on Indian scenario. To develop high-resolution climate change simulations over the Indian region, Providing Regional Climates for Impact Studies (PRECIS) modelling system was applied [3], and it was observed that the warming is monotonously widespread over India with substantial spatial differences in the projected rainfall changes.

Different data sets have been used in order to examine the monotonic change in the level of global temperature by various investigators [47]. From the diversified methodologies adopted, it is observed that the global surface temperature has risen by 0.74°C with an error of ±0.18°C when estimated by a linear trend over the last 100 years (1906–2005) and 0.07°C ± 0.02°C per decade. An executive summary is made available by Solomon et al. [8] on climate change in the fourth assessment report of the Intergovernmental Panel on Climate Change (IPCC). Climate change features are local or regional highly influenced by agricultural practices, urbanisation, heat island effect, changing irrigation, aerosol modifications, desertification, deforestation, and decay processes of organic materials [914]. Long-term changes in maximum temperature, minimum temperature, and their range have been analyzed over two stations in Cyprus by Price et al. [15] who established a reduction in range and thereby a possibility in increased cloud cover and tropospheric aerosols. Recently seasonal variability in mean sea level pressure extremes over the 11 zones of India and adjoining sea area has been studied by Singh et al. [16].

In the present paper, features on surface climate change over Thumba Equatorial Rocket Launching Station (TERLS), a coastal station near the South tip of Peninsular India, are critically examined. TERLS was established as an international rocket launching station dedicated to the United Nations for conducting meteorological and sounding rocket experiments on the geomagnetic equator. Meteorology group, TERLS has one of the heaviest responsibilities in the range to cater (i) weather forecasting and data needs of sounding rocket launch operations (in planning, pad operations and launching), (ii) climatological briefing to working crews including scientists to plan experiments, (iii) Organise various met. Observations, (iv) processing, analyses, storage and retrieval of various met. data (Class I surface observatory, autoweather stations, sonic anemometers, 50 m micro meteorological flux tower facility, meteorological balloon launches, rocket chaff derived wind data) and provide the data for users demand. For the planning of a scientific mission, the variability of climate is very crucial and thereby the importance of the present study on the climate change perspective. Importance of weather factors affecting rocket operations can be found elsewhere in detail [17, 18]. Steady increase in temperature anomaly at different levels in the troposphere over TERLS was investigated and documented in the Space research in India [19]. Joseph and Simon [20] have investigated weakening trend of the South West monsoon current through Peninsular India. Dash et al. [21] reported the decreasing frequency of monsoon depressions over the Indian region, and Krishnakumar et al. [22] have analysed significant decrease in South West monsoon rainfall while increase in post monsoon season over the whole Kerala State (also a meteorological subdivision), India. From the voluminous data bank, the present study analyses daily, monthly, seasonal, and interannual variability on different meteorological parameters based on both eye reading instruments and non-instrumental practices in cloud and visibility observations as per World Meteorological Organization (WMO) practices. Till the inception of long-term uninterrupted data collection through instrument-based cloud and visibility measurements, data collection from such non-instrumental means solely vested on conventional methods in surface meteorological observations adopted worldwide is to be accepted for a reasonable climatology extraction. The present work examines the imprints made by the climatic change on several meteorological parameters in different time domains, like a specific time, daily, monthly, and interannually unlike the climate change associated variations on a single meteorological parameter like temperature or rainfall alone have dealt in most of the climatic change studies.

2. Data and Methodology

A class I surface meteorological observatory was established at TERLS (8°32′′N/76° 52′′E, 6.7 amsl) by India Meteorological Department (IMD) approximately 70 m away from the coastline of the Arabian Sea coast (Figure 1). TERLS is situated hardly 8 km West North West of Thiruvananthapuram, the capital city of the Kerala state of India. Various meteorological data from 0830 IST (Indian Standard Time, IST = GMT + 0530 hr) observations for a period of 45 years (1964–2008), except for pressure where it is hourly data (1973–2008), are used in order to investigate changes encountered in meteorological parameters due to climate change. The data used are obtained from the following measurements at 0830 IST (1964–2008) namely, (i) maximum temperature (past 24 hr reported at 0830 IST), (ii) minimum temperature (past 24 hr reported at 0830 IST), (iii) past 24 hr rainfall (reported at 0830 IST), (iv) dry-bulb temperature at 0830 IST, (v) wet-bulb temperature at 0830 IST, (vi) relative humidity at 0830 IST, (vii) soil temperature at 0830 IST, (viii) non-instrumental eye observations on clouds and visibility at 0830 IST, (ix) wind speed and direction at 0830 IST, and (x) hourly surface pressure (1973–2008). Table 1 provides details on various instrumentation, exposure, maintenance, data span and the lowest and the highest in maximum absolute error values among various monthly statistics from 1964 to 2008.

The above mentioned 45 years of data have partitioned into two climatological slots of 30 years each, namely, 1964 to 1993 and from 1979 to 2008, respectively, called as past data for the former and present data for the later. Maximum absolute error values in the total data set (1962–2008) are the summation of errors [23] of past and present data sets. Individual error is computed as the ratio of standard deviation to the number of observations of a particular data set. Trenberth [24] has established the climate shift as a prominent upward trend in temperature after 1976. In order to include this observation influence, half portion (15 years) period from 1979 to 1993 is made common in the latter half of past data and in the former part of present data. The difference obtained among present data and past data has been represented as Δ in the discussions throughout. To see genuineness in findings attributed to climate change, Δ values are compared with absolute errors so as to see that Δ do not fall within absolute errors of corresponding meteorological parameters.

3. Results and Discussion

3.1. Δ Climatology of the Past 24 hr Observations

Maximum temperature, minimum temperature, and 24 hr accumulated rainfall are the only 3 meteorological parameters registered at 0830 IST having effect of the past 24 hr variability behaviour. Climatology for various parameters has been derived for past data and present data and presented in Tables 2, 3, 4, and 5 in which months that constitute the different seasons are mentioned. Derived Δ values are also incorporated in the above tables. Comparison between past data and present data maximum temperature climatology shows marked increase in present data climatology in relation to the past data counterpart, and thereby, positive values alone are observed in Δ. Δ values vary from +0.01°C in the North East monsoon (NE monsoon) month November to +0.43°C in the pre- South West monsoon (SW monsoon) month May. The NE monsoon season comprises October and November, whereas the four months from June to September designate SW monsoon season, which are the main two Indian seasons. In minimum temperature, all Δ values are positive with maximum value in the winter month February with a value +0.53°C. The observations among maximum and minimum temperature climatology derived in past data and present data along with Δ values can infer that the minimum daily temperatures have generally increased at a larger rate than the maximum daily temperatures. Karl et al. [25] studied the possible cause of increased urbanisation effects and Rebetez, and Beniston [26] justified the result of increased cloud cover in the large increase in minimum temperature and small increase in maximum temperatures.

From each individual daily maximum and minimum temperature, mean daily temperatures ((maximum temperature + minimum temperature)/2) are computed and further derive the past data and present data climatologies of mean temperatures. Only positive Δ values are observed throughout the year, which clearly indicate that the temperature over this tropical coastal station is increasing. Temperature climatologies have shown the tropical monsoonal characteristic of two maxima and two minima in monthly variation.

The diurnal temperature range is calculated from individual maximum and minimum temperature value and calculated past data and present data temperature range climatologies. Range climatologies are featured by annual variation. The prominent feature of negative Δ values in range is seen, which shows that the difference between maximum temperature and minimum temperature in a day is decreasing. Also, the larger rate of increase in minimum temperatures than maximum temperatures results in these negative Δ values in the diurnal temperature range.

The rainfall Δ shows decrease in rainfall in present data for the summer months April and May. Well-marked negative Δ values of about 21 mm are observed in the SW monsoon months (June–August). An increase in rainfall activity is detected in the SW monsoon month September and in the NE monsoon months October and November (debted to monsoonal climate) with maximum value of +36.50 mm in October. The station weather is dominated by convective mesoscale systems during the NE monsoon. Probably this investigation poses the doubt about whether such meso-scale weather phenomena like thunderstorm formations are on increasing trend during NE Monsoon and thereby an increase in Δ.

3.2. Climatology of a Particular Observation at 0830 IST

Analyses are carried out on different parameters collected at 0830 IST. Δ in dry bulb temperature, wet bulb temperature, and dew point temperature obtained from the two climatologies show only positive values. It ranges from +0.60°C in July to +1.01°C in January for dry bulb temperature. Variations in Δ among the SW monsoon months are less in general for all the above three variables, even though a maximum value in Δ is observed in dry bulb temperature. The less variation of wet bulb and dew point temperatures may be as a result of moisture laden atmosphere during the period which imparts corresponding variations in the moisture parameters like wet bulb temperature and dew point temperature.

Station relative humidity values at 0830 IST show more than 70% throughout the year with maximum in SW monsoon of the order 89%. Low negative values in Δ relative humidity are detected as the air becomes drier in present data compared to past. As the air is getting dried, both dew point and the wet bulb temperatures are supposed to show negative Δ values, but the feature is not detected. This may be inferred as the dry bulb temperature is the key parameter in elevating Δ values of both dew point and wet bulb temperatures, and the meagre relative humidity decrease in present data than past data is insignificant to impart Δ changes in dew point and wet bulb temperatures. Climatologies of RH also have shown two maxima and two minima, whereas the secondary minima around the NE monsoon is shallow in character.

Significant positive Δ values are observed among the two climatologies in soil temperature with a low value of +0.99°C to a high of +2.35°C in July and January, respectively. Perhaps the Δ increase observed in the winter months January and February may contribute reduction in radiational cooling which leads to a reduction in fog formation probability over the station. Soil temperature climatologies also have shown tropical monsoonal type with two maxima and minima.

Wind speed (both scalar and resultant) is decreased in present climatologies. Mathematically resultant wind speed shall be determined as the square root of the sum of the squares of zonal and meridional wind components. The encountered Δ values that rise to a maximum value of the order 0.5 m/s are due to decrease in wind speeds in the present data climatology in comparison with past data climatology. Wind steadiness factor values at 0830 IST are computed as the ratio of mean resultant wind speed to mean scalar wind (average value of scalar wind speeds) expressed in %, and the two slot climatologies are made. Wind steadiness factor is a general meteorological technique in order to derive the quantitative estimate of steadiness of wind direction [27]. The feature of losing controllability of blowing winds from a specific direction (computed from wind steadiness) as the demand of a particular season might be justified as a result of heat island effect of the very fast developing urbanisation hardly 8 km East South East of the observatory and the increasing vehicle populations and thereby drastic modifications in at 0830 IST observations which belong to land breeze prevalence over the station in the non-SW monsoonal months (up to a maximum of 20% decrease). Throughout the day, strong seabreeze is in prevalence in the SW monsoon months. Even in the SW monsoon which is dominantly associated with unidirectional wind flow, slight deteriorations in the steadiness of the wind directions have observed. Both the above observations brought out the phenomenal changes in the steadiness of offshore winds (land breeze) at 0830 IST during the non-SW monsoonal months at greater extend and onshore wind (sea breeze) at 0830 IST during the SW monsoon steadiness in smaller grades. The feature of very prominent decrease in direction steadiness in preset data to the past is clearly detected. As negative Δ dominance exists in the non-SW monsoon months where change over of land-sea breezes are so pronounced for the station, this observation may be further studied in order to check any possible modifications in the meso-scale meteorological phenomena like sea/land breeze, thunderstorms due to associated thermal and moisture environment changes.

Visibility and cloud climatologies are also studied. As these observational procedures are non-instrumental, the data is too subjective and quality observation is solely vested on the experience of the weather observer. The observation shows, better visibility codes in the present data climatology. A visibility code of 96 and above signifies that a well-defined object (like a mobile phone tower) at least 1 km far can be distinctly differentiable. Among various cloud climatology, dominant positive Δ values are existed in high clouds through out the year which may be pointed out as one of the reason for subdued albedo effect which enhances greenhouse effect, and hereby surface temperature escalation. This may be the result of number of high cloud types (Ci and Cc) are less compared to the varieties of low and middle cloud types, which enables easy estimation for the observer while making counts on high clouds. The increase in high cloud climatology values in present data may be the result of increase in air traffic resulted contrails over the region which results an increase in cirrus cloud formations due to tropical dynamics unlike microphysical processes in mid-latitudes as in Minnis et al. [28], Zerefos et al. [29]. Eleftheratos et al. [30] have studied the variability of cirrus clouds over the tropics governed by the activity of El Nino Southern Oscillation, the Quasi Biennial Oscillation and solar activity. Possible radiation budget alterations, modifications at the time of occurrence of minimum and maximum temperatures are to be critically examined as a part of future investigations. Prominent two maxima (secondary maximum is very shallower compared to the primary) and two minima were seen in middle, high and total cloud climatology.

3.3. Mean Daily Climatology in Temperature

Mean of daily maximum, minimum, mean and range in temperature between daily maximum and minimum temperatures are computed and depicted for the 365 days in a year (Figures 2, 3, 4 and 5) along with their standard deviations. The increase in mean daily values in present data temperatures are quite obvious and standard deviation values are fallen within 2.7°C for all the four mean daily climatologies. A difference statistics for 365 days between present data to past data is made and is shown in Figure 6. Out of 365 days, more than 80% mean daily maximum, minimum and mean temperatures are distributed as greater than zero occurrences in the difference statistics. An elevated frequency of occurrence of more than 65% under less than zero class is a concrete feature which implies that the range among maximum and minimum temperature is alarmingly diminishing.

3.4. Interannual Variability

Inter annual variability has critically examined for annual mean of minimum, maximum and mean temperatures, annual rainfall in the SW and NE monsoon over the station and for the Kerala State and average annual surface pressure in the SW and NE monsoons for the station. Each inter annual variability data has undergone an 11 year moving average with end point extrapolation till the raw data availability to bring out long term trend and the trend line also has shown in each inter annual variability figures similar to Krishnakumar et al. [22] for the interannual variability of rainfall over Kerala. Interpolated values are obtained for each year from the trend and such a data set obtained, has been subjected to statistical significant test. interpolated data has not taken for a particular year, if the raw data point is missing for the corresponding year. This method has enabled no reduction in the degrees of freedom (df). To confirm the statistical significance on various 11 year moving averaged data and also to assess the strength of sample correlation coefficient ( ), linear regression t-test were carried out and the results are presented (Table 6). Generally strong and significant linear trend variations are observed. Even though linear trends of maximum and minimum temperatures have shown strong and significant statistical trends, the derived variable-temperature range has resulted a weak and insignificant test statistic. This may be statistically inferred as a variable derived out of two variables which favours strong and significant test statistic separately need not show the same test static of strong and significant.

3.4.1. Interannual Variability in Temperature

To assess the interannual variability with temperature, annual mean minimum, mean maximum and mean of mean temperatures are found out for each year from 1964 to 2008. Also to ascertain deviation from the normal, each yearly mean value is deducted from the mean of the entire values available for present data and thereby extracted the anomaly in temperature. Mean and standard deviation of daily mean, maximum, minimum and daily temperature range are computed for a year and depicted in Figures 7, 8, 9 and 10. The fitted linear trend line slope (dy/dx) as well as correlation coefficient ( ) values are annotated in graphs. It is seen that from both mean daily temperature as well as from the maximum temperature an abrupt prominent shift or a discernible upward trend is occurred after 1976 may be like the “climate shift” as stated by Trenberth [24]. The 45 year rate in mean daily, maximum and minimum are obtained as +0.585°C, +0.540°C and +0.630°C respectively or with decadal rates of (+ )°C, (+ )°C and (+ )°C accordingly. In the obtained value for the global mean surface temperature as a whole (+ )°C, surface temperatures over land have risen at about double the ocean rate after 1979 (more than (+ )°C per decade versus (+ )°C per decade) as per the Trenberth et al. [31]. From 1950 to 2004, the global land trend in diurnal temperature range was extracted as negative (−0.07°C per decade) although tropical areas are still underrepresented [31]. Kumar et al. [32] showed trends of increase in diurnal temperature range, significant increase in annual maximum temperature, and practically trendless annual minimum temperature over India in the last century. Our inter annual variability results coincide mostly with former studied data period which supports the decrease of diurnal temperature range, and the same is very vivid from Figure 6 discussed earlier. It infers that for this tropical coastal station the obtained decadal rate of (+ )°C, which is in agreement with global ocean rate, may be due to the station being coastal in origin. The error value reduction by ±0.03°C for TERLS can be due to more stabilized single station data. Also rate of increase in trend for minimum temperature is high in comparison to maximum temperature. A decreasing trend in maximum temperature is noticed after 2005 with decreasing standard deviation values. This observation might be correct as a result of shifting the observatory since 1964 to hardly 30 m towards north in 2005 were smooth and free sea-breeze flow which enables abrupt reach of maximum temperature in the diurnal course around the time of sea breeze onset. The feature of attaining the maximum temperature of the day around the time of onset of sea breeze is known [33]. As the station is of coastal in origin, the occurrence of daily maximum temperature is very much linked to the setting and strong establishment of sea breeze front towards the land areas.

Mean monthly anomaly structures in minimum, maximum, and mean temperatures are presented in time section contour diagrams (Figures 11, 12 and 13). All contour depictions show well-marked temperature increase during the course of years and are so pronounced after 1976.

3.4.2. Interannual Variability in SW and NE Monsoon Rainfall and Surface Pressure

Interannual variability in rainfall and surface pressure particularly for SW monsoon and NE monsoon are studied. A decrease in trend in SW monsoon and an increase in trend in NE monsoons are the observed features over the station with decadal variability of (− ) mm and (+ ) mm, respectively (Figures 14 and 15). Analysis by considering initial rainfall values (at the beginning year on the fitted straight line) of South West monsoon, North East monsoon, corresponding slope of them as 817.822 mm, 296.535 mm, −21.5 mm/decade and 39.3 mm/decade, respectively, shows the likelihood exceed of North East monsoon rainfall over South West monsoon over this region after the year 2050 if the current trend continues. Similar analyses are carried out on the slot 1965 to 2007 for the Kerala meteorological subdivision from the longest instrumental rainfall series of the Indian regions (1813–2006) ([34]; http://www.tropmet.res.in/). The actual rainfall amount over Kerala during 2007 is obtained from the quarterly journal the MAUSAM (July 2008 and October 2008 issues). Decadal variations show a decreasing trend in SW monsoon by (− ) mm and increasing trend in NE monsoon by (+ ) mm (Figures 16 and 17).

Interannual variability in surface pressure was derived and depicted for SW monsoon (Figure 18) and NE monsoon (Figure 19). For this, every hourly data available since 1973 was subjected to mean monthly value computation. Number of observations available ranges from 1100 to 2245 in SW monsoon and 500 to 1495 in NE monsoon, respectively, and provided along with the standard deviation values in the above figures. Since atmospheric pressure is a highly stabilized property over a region for a particular season, the decadal trend variability observed is (+ ) mb and (– ) mb for SW monsoon and NE monsoon, respectively. The positive (increase) and negative (decrease) interannual trends observed in SW monsoon and NE monsoon, respectively, might have contributed possible implications in the development of associated synoptic scale meteorological systems (position and orientation of troughs/ridges and low pressures) and thereby corresponding decrease and increase in SW monsoon and NE monsoon rainfall trends.

4. Conclusion

Long-term surface meteorological records from TERLS were collected and analysed to study the long-term changes in overall climatology, climatology pertained to a particular time of observation, mean daily climatology in temperature, interannual variability in temperature, interannual variability in SW and NE monsoon rainfall, and mean monthly anomaly structure in temperature. Results on various analyses show strong and vivid features contributed by climate change. The salient features observed contributed by climate change are as follows:(1) All the temperatures show increase in their behaviour with slight decrease in RH between the present climatology and past climatology.(2) Steadiness of wind decreases especially in non-SW monsoon months, which gives an indirect signature of possible modifications in meso-scale systems.(3)High cloud climatology increases and thereby probable modifications in radiation budget, time of occurrence of maximum and minimum temperatures.(4)Daily statistics for 365 days show more than 80% exceedance in present climatology in comparison to past in minimum, maximum, and mean temperature.(5)Interannual variability analysis in mean temperature for this tropical coastal station shows +1.3°C/100 years which is same as that of global ocean rate.(6)Observed interannual rainfall variability in SW monsoon and North East monsoon shows decrease in trend by –21.45 mm/decade and increase in trend by +39.3 mm/decade, respectively. Prognosis on this observation points out the likelihood exceedance of NE monsoon rainfall over the South West monsoon after the year 2050 for the station. But it is prejudicial to infer that the rainfall increase in NE monsoon can compensate the forecasted decrease in SW monsoon, since NE monsoon is mostly of heavy rainfall in nature as a result of meso-scale convective systems. Such rainfall runoff is very fast than the slow water infiltration feature in association with SW monsoon rainfall, which is essential for irrigation and agricultural needs.(7)Decadal variability in surface pressure has shown +0.06 mb and −0.10 mb for SW monsoon and NE monsoon, respectively.

The imprints due to climate change on various meteorological parameters over this coastal station proved that the region also faces vulnerable situation due to the phenomena, and therefore, necessary engineering methods/socioeconomic steps are to be undertaken.


The availability of this unique data for a continuous period of 45 years is the outcome of concerted and whole hearted efforts put in by the staffs of TERLS/India Meteorological Department. It is a pleasure to acknowledge with esteemed gratitude the team that worked since 1963 at the TERLS, Thumba, whose efforts have yielded the voluminous unique data bank which is available for investigations. The authors are extremely grateful to Dr. K Krishnamoorthy, Director, Space Physics Laboratory, VSSC, and the reviewers for their valuable comments that helped to improve the paper to the final form.