Modeling urban heat islands in heterogeneous land surface and its correlation with impervious surface area by using night-time ASTER satellite data in highly urbanizing city,Delhi-India

The present study is an assessment and identification of urban heat island (UHI) in the environment of one of the fastest urbanizing city of India, Delhi Metropolis, employing satellite image of ASTER and Landsat 7 ETM+ in the thermal infrared region 3–14 μm. Temporal (2001 and 2005) ASTER datasets were used to analyze the spatial structure of the thermal urban environment subsequently urban heat island (UHI) in relation to the urban surface characteristics and land use/land cover (LULC). The study involves derivation of parameters governing the surface heat fluxes, constructing statistics of ASTER thermal infrared images along with validation through intensive in situ measurements. The average images reveal spatial and temporal variations of land surface temperature (LST) of night-time and distinct microclimatic patterns. Central Business District (CBD) of Delhi, (Connaught Place, a high density built up area), and commercial/industrial areas display heat islands condition with a temperature greater than 4 °C compared to the suburbs. The small increase in surface temperature at city level is mainly attributed to cumulative impact of human activities, changes in LULC pattern and vegetation density. In this study the methodology takes into account spatially-relative surface temperatures and impervious surface fraction value to measure surface UHI intensity between the urban land cover and rural surroundings. Both the spatial and temporal variation in surface temperature associated with impervious surface area (ISA) has been evaluated to assess the effect of urbanization on the local climate.     

Mesospheric OH temperatures: Simultaneous ground-based and SABER OH measurements over Millstone Hill

We present rotational temperature measurements of the mesospheric OH emission layer using a meridional imaging spectrograph at Millstone Hill (42.6°N, 72.5°W). The system is equipped with a state-of-the-art bare-CCD detector and can yield simultaneous quasi-meridional images of the mesospheric OH and O₂ intensity and temperature fields at 87 and 94 km altitude during the course of each night. A cross-validation study of the rotational OH temperature measurements obtained on 61 nights during the autumnal months of 2005–2007 was undertaken with near-simultaneous kinetic temperature measurements made by the SABER instrument aboard the NASA TIMED satellite during overpasses of Millstone Hill. Excellent agreement was obtained between the two datasets with the small differences being attributable to differences in the spatial and temporal averaging inherent between the two datasets.     

Change in refractivity of the atmosphere and large variation in TEC associated with some earthquakes,observed from GPS receiver

The present study reports the analysis of GPS based TEC for our station Surat (21.16°N, 72.78°E) located at the northern crest of equatorial anomaly region in India at times close to some earthquake events (M ? 5) during the year 2009 in India and its neighbouring regions. The TEC data used in the study are obtained from GPS Ionospheric Scintillation and TEC Monitoring (GISTM) system. The TEC data has been analysed corresponding to 11 earthquakes in low solar activity period and quiet geomagnetic condition. We found that, out of 11 cases of earthquakes (M > 5) there were seven cases in which enhancement in TEC occurred on earthquake day and in other four cases there was depletion in TEC on earthquake day. The variation in refractivity prior to earthquake was significant for the cases in which the epicentre lied within a distance of 600 km from the receiving station. By looking into the features on temporal enhancement and depletion of TEC a prediction was made 3–2 days prior to an earthquake (on 28 October 2009 in Bhuj – India). The paper includes a brief discussion on the method of potentially identifying an impending earthquake from ionospheric data.     

Equatorial Ionospheric Anomaly (EIA) and comparison with IRI model during descending phase of solar activity (2005–2009)

The ionospheric variability at equatorial and low latitude region is known to be extreme as compared to mid latitude region. In this study the ionospheric total electron content (TEC), is derived by analyzing dual frequency Global Positioning System (GPS) data recorded at two stations separated by 325 km near the Indian equatorial anomaly region, Varanasi (Geog latitude 25°, 16^/ N, longitude 82°, 59^/ E, Geomagnetic latitude 16°, 08^/ N) and Kanpur (Geog latitude 26°, 18^/ N, longitude 80°, 12^/ E, Geomagnetic latitude 17°, 18^/ N). Specifically, we studied monthly, seasonal and annual variations as well as solar and geomagnetic effects on the equatorial ionospheric anomaly (EIA) during the descending phase of solar activity from 2005 to 2009. It is found that the maximum TEC (EIA) near equatorial anomaly crest yield their maximum values during the equinox months and their minimum values during the summer. Using monthly averaged peak magnitude of TEC, a clear semi-annual variation is seen with two maxima occurring in both spring and autumn. Results also showed the presence of winter anomaly or seasonal anomaly in the EIA crest throughout the period 2005–2009 only except during the deep solar minimum year 2007–2008. The correlation analysis indicate that the variation of EIA crest is more affected by solar activity compared to geomagnetic activity with maximum dependence on the solar EUV flux, which is attributed to direct link of EUV flux on the formation of ionosphere and main agent of the ionization. The statistical mean occurrence of EIA crest in TEC during the year from 2005 to 2009 is found to around 12:54 LT hour and at 21.12° N geographic latitude. The crest of EIA shifts towards lower latitudes and the rate of shift of the crest latitude during this period is found to be 0.87° N/per year. The comparison between IRI models with observation during this period has been made and comparison is poor with increasing solar activity with maximum difference during the year 2005.     

A morphological study of GPS-TEC data at Agra and their comparison with the IRI model

The temporal and seasonal variations of Total Electron Content (TEC) are studied at Agra (Geographic Lat. 27.17°N, Long. 78.89°E, Dip: 41.4°), India, which is in the equatorial anomaly region, for a period of 12 months from 01 January to 31 December, 2007 using a Global Positioning System (GPS) receiver. The mean TEC values show a minimum at 0500 h LT (LT = UT + 5.5 h) and a peak value at about 1400 h LT. The lowest TEC values are observed in winter whereas largest values are observed in equinox and summer. Anomalous variations are found during the period of magnetic disturbances. These results are compared with the TEC derived from IRI-2007 using three different options of topside electron density, NeQuick, IRI01-corr, and IRI-2001. A good agreement is found between the TEC obtained at Agra and those derived from IRI models.     

Modeling medium-scale TEC structures,observed by Belgian GPS receivers network

GALOCAD project “Development of a Galileo Local Component for the nowcasting and forecasting of atmospheric disturbances affecting the integrity of high precision Galileo applications” aims to perform a detailed study on ionospheric small- and medium-scale structures and to assess the influence of these structures on the reliability of Galileo precise positioning applications. GPS-derived TEC (total electron content) is obtained from the Belgium Dense Network (BDN), consisting of 67 permanent GPS stations. An empirical 3-D model is developed for studying these ionospheric structures. The model, named LLT model, described temporal variations of TEC in latitude/longitude frame (46°, 52°)N and (−1°, 11°)E. The spatial variations of TEC are modeled by Tchebishev base functions, while the temporal variations are described by a trigonometric basis. To fit the model to the data, the observed area is divided into bins with (1° × 1°) geographic scale and 6 min on time axis. LLT model is made flexible, with varying number of coefficients along each axis. This allows different degree of smoothing, which is the key element of the present approach. Model runs with higher number of coefficients, capturing in details medium-scale TEC structures are subtracted from results obtained with smaller number of coefficients; the latter represent the background ionosphere. The residual structures are localized and followed as they travel across the observed area. In this way, the size, velocity, and direction of the irregular structures are obtained.     

Response of the EIA ionosphere to the 7-8 May 2005 geomagnetic storm

In this paper, response of low latitude ionosphere to a moderate geomagnetic storm of 7–8 May 2005 (SSC: 1920 UT on 7 May with Sym-H minimum, ∼−112 nT around 1600 UT on 8 May) has been investigated using the GPS measurements from a near EIA crest region, Rajkot (Geog. 22.29°N, 70.74°E, Geomag.14°), India. We found a decrease in total electron content (TEC) in 12 h after the onset of the storm, an increase during and after 6 h of Sym-H deep minimum with a decrease below its usual-day level on the second day during the recovery phase of the storm. On 8 May, an increase of TEC is observed after sunset and during post-midnight hours (maximum up to 170%) with the formation of ionospheric plasma bubbles followed by a nearly simultaneous onset of scintillations at L-band frequencies following the time of rapid decrease in Sym-H index (−30 nT/h around 1300 UT).     

An introduction of mountain-based GPS radio occultation experiments in China

The mountain-based GPS radio occultation is a novel approach to lower atmospheric profiling. The experiments of the mountain-based GPS radio occultation were conducted on the top of Mt. Yaogu (29.38°N, 113.68°E, ∼1240 m) on December 17, 2003, and on the top of Mt. Jiugong (29.39°N, 114.65°E, ∼1550 m) on July 24, 2004. Based on these observation data, the scientific data processing software has been developed and is used to retrieve successfully the atmospheric refractivity profiles. The validation experiment was performed on the top of Mt. Wuling (40.60°N, 117.48°E, ∼2118 m) during August 1–29, 2005. Collocated automatic weather station and the radiosondes nearby were operated simultaneously for the comparison campaign. Results show that the radio occultation technique obtained about 40 profiles every day with the receiver antenna pointing to the south. Comparisons show that the refractivity measured by occultation agree well with those by the radiosondes, but not well with those by the automatic weather station due to their much different geographic locations of measurements. Results of these experiments suggest that the mountain-based GPS radio occultation is an economic reliable novel technique monitoring temporal and spatial variations of local lower atmospheric environments.     

The solar eclipse and its associated ionospheric TEC variations over Indian stations on January 15, 2010

An annular solar eclipse occurred over the Indian subcontinent during the afternoon hours of January 15, 2010. This event was unique in the sense that solar activity was minimum and the eclipse period coincides with the peak ionization time at the Indian equatorial and low latitudes. The number of GPS receivers situated along the path of solar eclipse were used to investigate the response of total electron content (TEC) under the influence of this solar eclipse. These GPS receivers are part of the Indian Satellite Based Augmentation System (SBAS) named as ‘GAGAN’ (GPS Aided Geo Augmented Navigation) program. The eight GPS stations located over the wide range of longitudes allows us to differentiate between the various factors induced due to solar eclipse over the equatorial and low latitude ionosphere. The effect of the eclipse was detected in diurnal variations of TEC at all the stations along the eclipse path. The solar eclipse has altered the ionospheric behavior along its path by inducing atmospheric gravity waves, localized counter-electrojet and attenuation of solar radiation intensity. These three factors primarily control the production, loss and transport of plasma over the equatorial and low latitudes. The localized counter-electrojet had inhibited the equatorial ionization anomaly (EIA) in the longitude belt of 72°E–85°E. Thus, there was a negative deviation of the order of 20–40% at the equatorial anomaly stations lying in this ‘inhibited EIA region’. The negative deviation of only 10–20% is observed for the stations lying outside the ‘inhibited EIA region’. The pre-eclipse effect in the form of early morning enhancement of TEC associated with atmospheric gravity waves was also observed during this solar eclipse. More clear and distinctive spatial and temporal variations of TEC were detected along the individual satellite passes. It is also observed that TEC starts responding to the eclipse after 30 min from start of eclipse and the delay of the maximum TEC deviation from normal trend with respect to the maximum phase of the eclipse was close to one hour in the solar eclipse path.     

Study of TEC changes during geomagnetic storms occurred near the crest of the equatorial ionospheric ionization anomaly in the Indian sector

GPS satellites data obtained at Bhopal (23.16° N, 77.36° E, geomagnetic latitude 14.23° N) India were analyzed to study the TEC changes during several geomagnetic storms (−300 nT < Dst < −50 nT) occurred in 2005–2007. We had segregated the storms according to the Dst value, i.e. moderate storms (−100 nT < Dst ? −50 nT), strong storms (−150 nT < Dst < −100 nT), and severe storms (Dst less than −150 nT). Total of 21 geomagnetic storms (10 moderate, 9 strong, 2 severe) are considered for the present study. Deviation in vertical total electron content (VTEC) during the main phase of the storm was found to be associated with the prompt penetration of electric field originated due to the under-shielding and over-shielding conditions for almost all geomagnetic storms discussed in this paper. For most of the storms VTEC shows the positive percentage deviation during the main phase while it shows positive as well as the negative deviation during the recovery phase of the storms. The −80% deviation in VTEC was found for geomagnetic storm occurred on July 17, 2005 and the negative trend continued for recovery phase of the storm. This was mainly due to the thermospheric composition changes by Joule heating effect at auroral latitudes that generate electric field disturbance at low latitudes. Traveling ionospheric disturbances (TIDs) were responsible for the formation of wave like nature in VTEC for the storms occurred on May 15, 2005, whereas it was not observed for storm occurred on August 24, 2005.     

Morphological studies on ionospheric VHF scintillations over an Indian low latitude station during a solar cycle period (2001–2010)

The amplitude scintillations data recorded at 244 MHz from the geostationary satellite, FLEETSAT (73°E) at a low latitude station, Waltair (17.7°N, 83.3°E) during the ten year period of high to low solar activity from 2001 to 2010 is considered to study the occurrence characteristics of the VHF scintillations. A close association between the intense scintillations on VHF signals during pre-midnight hours, associated with range type of spread-F on ionograms and a relatively weak and slow fading scintillations during post-midnight hours associated with frequency type of spread-F is observed during the relatively high sunspot years from 2001 to 2004, whereas during the low sunspot years from 2005 to 2010 the scintillation activity as well as spread-F activity are found to be minimum. During both the high and low sunspot years, it is observed that the maximum scintillation activity occurs during equinoctial months followed by winter with the minimum occurrence during summer months. The annual mean percentage occurrence of scintillations is found to be clearly associated with the variations in the annual mean sunspot number. The nocturnal variations in the occurrence of scintillations show the onset of scintillation activity starts from 19:00 h LT with maximum of occurrence around 21:00 h LT. A clear semiannual variation in the occurrence of scintillations is observed during pre-midnight hours with two peaks in equinoctial months of March/April and September/October. The number of scintillation patches observed is found to be more during pre-midnight hours compared to those during post-midnight hours. The most probable scintillation patch duration lies around 30 min. Further, it is also found that the number of scintillation patches with durations of 60 min and more decreases with the increase in the patch duration. It is also observed in general that the scintillation activity is inhibited during geomagnetic disturbed days.     

Diagnostics of equatorial and low latitude ionosphere by TEC mapping over Brazil

The total electron content (TEC) in the equatorial and low-latitude ionosphere over Brazil was monitored in two dimensions by using 2011 data from the ground-based global navigation satellite system (GNSS) receiver network operated by the Brazilian Institute for Geography and Statistics. It was possible to monitor the spatial and temporal variations in TEC over Brazil continuously during both day and night with a temporal interval of 10 min and a spatial resolution of about 400 km. The daytime equatorial ionization anomaly (EIA) and post-sunset plasma enhancement (PS-EIA) were monitored over an area corresponding to a longitudinal extension of 4000 km in South America. Considerable day-to-day variation was observed in EIA and PS-EIA. A large latitudinal and longitudinal gradient of TEC indicated a significant ionospheric range error in application of the GNSS positioning system. Large-scale plasma bubbles after sunset were also mapped over a wide range. Depletions with longitudinally separated by more than 800 km were observed. They were extended by more than 2000 km along the magnetic field lines and drifted eastward. It is expected that 2-dimensional TEC mapping can serve as a useful tool for diagnosing ionospheric weather, such as temporal and spatial variation in the equatorial plasma trough and crest, and particularly for monitoring the dynamics of plasma bubbles.     

Spectral and temporal characteristics of aerosol optical depth over a wet tropical location in North East India

A network of multi wavelength solar radiometer (MWR) stations has been in operation since the 1980s in India for measurement of aerosol optical depth (AOD). This network was augmented recently with the addition of a large number of stations located across the length and breath of India covering a variety of climate regimes. The spectral and temporal variations of aerosol optical depths observed over Dibrugarh located in the North East of India (27.3°N, 94.5°E) are investigated by analyzing the data obtained from a MWR during October 2001–September 2003 using the Langley technique. AOD varies with time of the day, month of the year and season. From January to April and October to December, aerosol optical depth decreases with wavelength whereas during May–September aerosol optical depth has been found to be nearly independent of wavelengths. AOD is higher during pre-monsoon season (March–May) and lower in the monsoon (June–September) season at about all wavelengths. The temporal variation of AOD over Dibrugarh have also been compared with those reported from selected locations in India.     

Lower-mesospheric inversion layers over brazilian equatorial region using TIMED/SABER temperature profiles

Lower-mesospheric inversion layers (MILs) were studied using the temperature profiles observed by TIMED/SABER over Cariri (7.5°S, 36.5°W), Brazil, in 2005. A total 175 MILs were identified with the maximum occurrence in April and October and the minimum in January and July. The lower MIL is located in a height region from 70 to 90 km, with the peak at around 83 ± 4 km with the temperature of 205 ± 5 K, and the thickness of 4–10 km. The results show large amplitudes of MILs during equinoxes and minimum in solstices, with a clear semiannual variation. A general feature of lower MIL in monthly mean profile was observed twice a year, one from February to May, and the other from August to October with a downward shift of the top level. These results suggest that formation and long persistence of MIL is an important factor to investigate propagation of atmospheric gravity waves in the mesosphere-lower thermosphere (MLT) region.     

CME–flare association during the 23rd solar cycle

The relation between coronal mass ejections (CMEs) and solar flares are statistically studied. More than 10,000 CME events observed by SOHO/LASCO during the period 1996–2005 have been analyzed. The soft X-ray flux measurements provided by the Geostationary Operational Environmental Satellite (GOES), recorded more than 20,000 flares in the same time period. The data is filtered under certain temporal and spatial conditions to select the CME–flare associated events. The results show that CME–flare associated events are triggered with a lift-off time within the range 0.4–1.0 h. We list a set of 41 CME–flare associated events satisfying the temporal and spatial conditions. The listed events show a good correlation between the CME energy and the X-ray flux of the CME–flare associated events with correlation coefficient of 0.76.     

TEC variability near northern EIA crest and comparison with IRI model

Monthly median values of hourly total electron content (TEC) is obtained with GPS at a station near northern anomaly crest, Rajkot (geog. 22.29°N, 70.74°E; geomag. 14.21°N, 144.9°E) to study the variability of low latitude ionospheric behavior during low solar activity period (April 2005 to March 2006). The TEC exhibit characteristic features like day-to-day variability, semiannual anomaly and noon bite out. The observed TEC is compared with latest International Reference Ionosphere (IRI) – 2007 model using options of topside electron density, NeQuick, IRI01-corr and IRI-2001 by using both URSI and CCIR coefficients. A good agreement of observed and predicted TEC is found during the daytime with underestimation at other times. The predicted TEC by NeQuick and IRI01-corr is closer to the observed TEC during the daytime whereas during nighttime and morning hours, IRI-2001 shows lesser discrepancy in all seasons by both URSI and CCIR coefficients.     

Simultaneous bi-station measurements of mesospheric waves from Indian low latitudes

Coordinated mesospheric night-airglow measurements have been carried out from two stations, Gadanki (13.5° N, 79.2° E) and Allahabad (25.5° N, 81.9° E), India during April 2009 to study the common gravity wave features. With two nights of coordinated measurements we find some of the wave periodicities to be similar at the two locations. Simultaneous OH and O₂ intensity measurements over Gadanki reveal these features to be upward propagating gravity waves while the coordinated OH intensity measurements of similar waves from Allahabad show the large spatial extent of these waves.     

The asymmetrical features in electron density during extreme solar minimum

The variations of plasma density in topside ionosphere during 23rd/24th solar cycle minimum attract more attentions in recently years. In this analysis, we use the data of electron density (Ne) from DEMETER (Detection of Electromagnetic Emissions Transmitted from Earthquake Regions) satellite at the altitude of 660–710 km to investigate the solstitial and equinoctial asymmetry under geomagnetic coordinate system at LT (local time) 1030 and 2230 during 2005–2010, especially in solar minimum years of 2008–2009. The results reveal that ΔNe (December–June) is always positive over Southern Hemisphere and negative over northern part whatever at LT 1030 or 2230, only at 0–10°N the winter anomaly occurs with ΔNe (December–June) > 0, and its amplitude becomes smaller with the declining of solar flux from 2005 to 2009. The ΔNe between September and March is completely negative during 2005–2008, but in 2009, it turns to be positive at latitudes of 20°S–40°N at LT 1030 and 10°S–20°N at LT 2230. Furthermore, the solstitial and equinoctial asymmetry index (AI) are calculated and studied respectively, which all depends on local time, latitude and longitude. The notable differences occur at higher latitudes in solar minimum year of 2009 with those in 2005–2008. The equinoctial AI at LT 2230 is quite consistent with the variational trend of solar flux with the lowest absolute AI occurring in 2009, the extreme solar minimum, but the solstitial AI exhibits abnormal enhancement during 2008 and 2009 with bigger AI than those in 2005–2007. Compared with the neutral compositions at 500 km altitude, it illustrates that [O/N₂] and [O] play some roles in daytime and nighttime asymmetry of Ne at topside ionosphere.     

A comparison of different regression models for downscaling Landsat and MODIS land surface temperature images over heterogeneous landscape

Remotely sensed high spatial resolution thermal images are required for various applications in natural resource management. At present, availability of high spatial resolution (<200 m) thermal images are limited. The temporal resolution of such images is also low. Whereas, coarser spatial resolution (∼1000 m) thermal images with high revisiting capability (∼1 day) are freely available. To bridge this gap, present study attempts to downscale coarser spatial resolution thermal image to finer spatial resolution using relationships between land surface temperature (LST) and vegetation indices over a heterogeneous landscape of India. Five regression based models namely (i) Disaggregation of Radiometric Temperature (DisTrad), (ii) Temperature Sharpening (TsHARP), (iii) TsHARP with local variant, (iv) Least median square regression downscaling (LMS_DS) and (v) Pace regression downscaling (PR_DS) are applied to downscale LST of Landsat Thematic Mapper (TM) and Terra MODIS (Moderate Resolution Imaging Spectroradiometer) images. All the five models are first evaluated on Landsat image aggregated to 960 m resolution and downscaled to 480 m and 240 m resolution. The downscale accuracy is achieved using LMS_DS and PR_DS models at 240 m resolution at 0.61 °C and 0.75 °C respectively. MODIS data downscaled from 1000 m to 250 m spatial resolution results root mean square error (RMSE) of 1.43 °C and 1.62 °C for LMS_DS and PR_DS models, respectively. The LMS_DS model is less sensitive to outliers in heterogeneous landscape and provides higher accuracy when compared to other models. Downscaling model is found to be suitable for agricultural and vegetated landscapes up to a spatial resolution of 250 m but not applicable to water bodies, dry river bed sand sandy open areas.     

Measurement of mesopause temperature from hydroxyl nightglow at Kolhapur (16.8°N, 74.2°E), India

The paper reports the nightglow observations of hydroxyl (8–3), (7–2) and (6–2) Meinel band carried out at a low latitude station Kolhapur (16.8°N, 74.2°E, dip latitude 10.6°N), India during November 2002 to May 2005 with the objective of investigating mesopause dynamics based on derived OH rotational temperature. Overall, 132 nights of quality data were collected using filter-tilting photometer and an all-sky scanning photometer. The mean mesopause temperature observed at Kolhapur is 195 ± 11, 196 ± 9 and 195 ± 7 K from OH (8–3), (7–2) and (6–2) band emissions, respectively, using transition probabilities given by Langhoff et al. [Langhoff, S.R., Werner, H.J., Rosmus, P. Theoretical transition probabilities for the OH Meinel system. Journal of Molecular Spectroscopy 118, 507–529, 1986]. Small wave-like variations (periodicities ∼ few hours) existing over long period variations in derived temperatures are also present. A steady decrease of emission intensities from evening to dawn hours has been observed in approximately 59% of nights. No significant change of nightly mean temperatures has been noted. Furthermore, about 62% of observed nightly mean temperatures lie within one error bar of MSISE-90 model predictions.