Modeling Indian Ocean circulation: Bay of Bengal fresh water plume and the Arabian Sea Mini Warm Pool

The Indian subcontinent divides the north Indian Ocean into two tropical basins, namely the Arabian Sea and the Bay of Bengal. The Arabian Sea has high salinity whereas the salinity of the Bay of Bengal is much lower due to the contrast in freshwater forcing of the two basins. The freshwater received by the Bay in large amounts during the summer monsoon through river discharge is flushed out annually by ocean circulation. After the withdrawal of the summer monsoon, the Ganga – Brahmaputra river plume flows first along the Indian coast and then around Sri Lanka into the Arabian Sea creating a low salinity pool in the southeastern Arabian Sea (SEAS). In the same region, during the pre-monsoon months of February – April, a warm pool, known as the Arabian Sea Mini Warm Pool (ASMWP), which is distinctly warmer than the rest of the Indian Ocean, takes shape. In fact, this is the warmest region in the world oceans during this period. Simulation of the river plume and its movement as well as its implications to thermodynamics has been a challenging problem for models of Indian Ocean. Here we address these issues using an ocean general circulation model – first we show that the model is capable of reproducing fresh plumes in the Bay of Bengal as well as its movement and then we use the model to determine the processes that lead to formation of the ASMWP. Hydrographic observations from the western Bay of Bengal have shown the presence of a fresh plume along the northern part of the Indian coast during summer monsoon. The Indian Ocean model when forced by realistic winds and climatological river discharge reproduces the fresh plume with reasonable accuracy. The fresh plume does not advect along the Indian coast until the end of summer monsoon. The North Bay Monsoon Current, which flows eastward in the northern Bay, separates the low salinity water from the more saline southern parts of the bay and thus plays an important role in the fresh water budget of the Bay of Bengal. The model also reproduces the surge of the fresh-plume along the Indian coast, into the Arabian Sea during northeast monsoon. Mechanisms that lead to the formation of the Arabian Sea Mini Warm Pool are investigated using several numerical experiments. Contrary to the existing theories, we find that salinity effects are not necessary for the formation of the ASMWP. The orographic effects of the Sahyadris (Western Ghats) and resulting reduction in wind speed leads to the formation of the ASMWP. During November – April, the SEAS behave as a low-wind heatdominated regime where the evolution of sea surface temperature is solely determined by atmospheric forcing. In such regions the evolution of surface layer temperature is not dependent on the characteristics of the subsurface ocean such as the barrier layer and temperature inversion.

Solar energy potential assessment using GIS

Renewable energy resources are those having a cycling time less than 100 years and are renewed by the nature and their supply exceeds the rate of consumption. Renewable energy systems use resources that are constantly replaced in nature and are usually less polluting. In order to tap the potential of renewable energy sources, there is a need to assess the availability of resources spatially as well as temporally. Geographic Information Systems (GIS) along with Remote Sensing (RS) helps in mapping on spatial and temporal scales of the resources and demand. The spatial database of resource availability and the demand would help in the regional energy planning. This paper discusses the application of geographical information system (GIS) to map the solar potential in Karnataka state, India. Regions suitable for tapping solar energy are mapped on the basis of global solar radiation data, and this analysis provides a picture of the potential. The study identifies that Coastal parts of Karnataka with the higher global solar radiation is ideally suited for harvesting solar energy. The potential analysis reveals that, maximum global solar radiation is in districts such as Uttara Kannada and Dakshina Kannada. Global solar radiation in Uttara Kannada during summer, monsoon and winter are 6.31, 4.40 and 5.48 kWh/sq.m, respectively. Similarly, Dakshina Kannada has 6.16, 3.89 and 5.21 kWh/sq.m during summer, monsoon and winter.

Tracing the source of bottled water using stable isotope techniques

The Indian Summer Monsoon (ISM) precipitation recharges ground water aquifers in a large portion of the Indian subcontinent. Monsoonal precipitation over the Indian region brings moisture from the Arabian Sea and the Bay of Bengal (BoB). A large difference in the salinity of these two reservoirs, owing to the large amount of freshwater discharge from the continental rivers in the case of the BoB and dominating evaporation processes over the Arabian Sea region, allows us to distinguish the isotopic signatures in water originating in these two water bodies. Most bottled water manufacturers exploit the natural resources of groundwater, replenished by the monsoonal precipitation, for bottling purposes. The work presented here relates the isotopic ratios of bottled water to latitude, moisture source and seasonality in precipitation isotope ratios. We investigated the impact of the above factors on the isotopic composition of bottled water. The result shows a strong relationship between isotope ratios in precipitation (obtained from the GNIP data base)/bottled water with latitude. The approach can be used to predict the latitude at which the bottled water was manufactured. The paper provides two alternative approaches to address the site prediction. The limitations of this approach in identifying source locations and the uncertainty in latitude estimations are discussed. Furthermore, the method provided here can also be used as an important forensic tool for exploring the source location of bottled water from other regions. Copyright (C) 2011 John Wiley & Sons, Ltd.

Climatic fluctuations during the LIA and post-LIA in the Kumaun Lesser Himalaya, India: Evidence from a 400 y old stalagmite record

This paper presents the first stable isotope (delta O-18 and delta C-13) data of a similar to 400 years (1590-2006 AD) long annual to decadal-resolution speleothem record collected from the Indian Lesser Himalaya. The data show a variation from -2.7 to -5.9 parts per thousand in delta O-18 and -5.3 to -8.8 parts per thousand in delta C-13. The isotopic analyses indicate that the climate during this period can be divided into two stages: a wet phase during the Little Ice Age (LIA) (1590-1850 AD) and comparatively dry phase during the post-LIA after 1850 AD. However, the record also documents the minor dry events during the LIA and a wet episode after the LIA. Within the age uncertainty, the dry spells during the LIA are linked with the historical drought events in the Indian subcontinent and similar latitudes. The isotopic record is consistent with a number of previous studies in the areas influenced by the Westerlies but appears to be conflicting to the regions, dominated by the Indian Summer Monsoon (ISM). This may be due to the possible changes in the strength of Westerlies in the study area and added by negative anomaly of North Atlantic Oscillation (NAO) during the LIA. (C) 2012 Elsevier Ltd and INQUA. All rights reserved.

A note on new indices for the equatorial Indian Ocean oscillation

It is now well known that there is a strong association of the extremes of the Indian summer monsoon rainfall (ISMR) with the El Nio and southern oscillation (ENSO) and the Equatorial Indian Ocean Oscillation (EQUINOO), later being an east-west oscillation in convection anomaly over the equatorial Indian Ocean. So far, the index used for EQUINOO is EQWIN, which is based on the surface zonal wind over the central equatorial Indian Ocean. Since the most important attribute of EQUINOO is the oscillation in convection/precipitation, we believe that the indices based on convection or precipitation would be more appropriate. Continuous and reliable data on outgoing longwave radiation (OLR), and satellite derived precipitation are now available from 1979 onwards. Hence, in this paper, we introduce new indices for EQUINOO, based on the difference in the anomaly of OLR/precipitation between eastern and western parts of the equatorial Indian Ocean. We show that the strong association of extremes of the Indian summer monsoon with ENSO and EQUINOO is also seen when the new indices are used to represent EQUINOO.

Seasonal variability of rainfall recorded in growth bands of the Giant African Land Snail Lissachatina fulica (Bowdich) from India

Seasonal rainfall patterns in Bangalore, India, have been reconstructed using stable isotopic ratios in the growth bands of Giant African Land Snail shells. The present study was conducted at Bangalore, India which receives rain during the summer months. The oxygen isotopic record in the rainwater samples collected during different months covering the period of the summer monsoon of the year 2008 is compared with the isotopic ratio in the gastropod growth bands deposited simultaneously. The chronology of the shell growth band is independently established assuming the growth rate observed in a chamber experiment maintaining similar relative humidity and temperature conditions. A consistent pattern observed in the isotopic ratio in the gastropod growth bands and rainwater is demonstrated and provides a novel approach for precipitation reconstruction at seasonal and weekly time scales. This approach of using isotopic ratios in the gastropod growth bands for rainfall can serve as a substitute for filling gaps in rainfall data and for cases where no rain records are available. In addition, they can be used to determine the frequencies and magnitudes of dry spells from the past records. (C) 2013 Elsevier B.V. All rights reserved.

Effect of SST Variation on ITCZ in APE Simulations

The effect of meridional variation of sea surface temperature (SST) on tropical atmospheric circulation is analyzed using Aqua-planet Experiment (APE) simulations. The meridional SST gradient around the narrow SST peak in CONTROL simulation favours a strong and single equatorial Intertropical Convergence Zone (ITCZ, defined by the maximum of zonally averaged total precipitation) in all APE models. In contrast, flat equatorial SST peak (FLAT simulation) favours split/double ITCZs flanking the SST maximum, in the majority of the APE models. Although there is reasonable agreement for SST sensitivity of ITCZ among the APE models in CONTROL, there exists disparity among them in FLAT case. Similarly, while the total and convective precipitation responses are consistent among the models, the large-scale precipitation response shows considerable inter-model variations in FLAT case. The APE intercomparison indicates that the occurrence and positioning of the ITCZ are primarily related to boundary layer moisture convergence as a response to the meridional variation of SST. Furthermore, the meridional gradient of tropospheric temperature is found to be an important factor that can influence the positioning of ITCZ. FLAT SST distribution is found to be similar to the observed distribution over the Indian region during summer season. Models that yield double ITCZs in this case simulate an easterly jet over the equatorial region (similar to 15 degrees equatorward of the ITCZ). This is analogous to the Tropical Easterly Jet (TEJ), which is a unique feature observed over the Indian region during summer monsoon season, with its core at 12 degrees N, equatorward of the seasonal convergence zone centered along 25 degrees N. In these models, positive meridional temperature gradient and the associated easterly shear in the atmosphere strengthened by moisture convergence penetrate up to the upper troposphere, with which TEJ is in thermal wind balance.

Trends in aerosol optical depth over Indian region: Potential causes and impact indicators

The first regional synthesis of long-term (back to similar to 25 years at some stations) primary data (from direct measurement) on aerosol optical depth from the ARFINET (network of aerosol observatories established under the Aerosol Radiative Forcing over India (ARFI) project of Indian Space Research Organization over Indian subcontinent) have revealed a statistically significant increasing trend with a significant seasonal variability. Examining the current values of turbidity coefficients with those reported similar to 50 years ago reveals the phenomenal nature of the increase in aerosol loading. Seasonally, the rate of increase is consistently high during the dry months (December to March) over the entire region whereas the trends are rather inconsistent and weak during the premonsoon (April to May) and summer monsoon period (June to September). The trends in the spectral variation of aerosol optical depth (AOD) reveal the significance of anthropogenic activities on the increasing trend in AOD. Examining these with climate variables such as seasonal and regional rainfall, it is seen that the dry season depicts a decreasing trend in the total number of rainy days over the Indian region. The insignificant trend in AOD observed over the Indo-Gangetic Plain, a regional hot spot of aerosols, during the premonsoon and summer monsoon season is mainly attributed to the competing effects of dust transport and wet removal of aerosols by the monsoon rain. Contributions of different aerosol chemical species to the total dust, simulated using Goddard Chemistry Aerosol Radiation and Transport model over the ARFINET stations, showed an increasing trend for all the anthropogenic components and a decreasing trend for dust, consistent with the inference deduced from trend in Angstrom exponent.

Impact of diurnal forcing on intraseasonal sea surface temperature oscillations in the Bay of Bengal

The diurnal cycle is an important mode of sea surface temperature (SST) variability in tropical oceans, influencing air-sea interaction and climate variability. Upper ocean mixing mechanisms are significant at diurnal time scales controlling the intraseasonal variability (ISV) of SST. Sensitivity experiments using an Ocean General Circulation Model (OGCM) for the summer monsoon of the year 2007 show that incorporation of diurnal cycle in the model atmospheric forcings improves the SST simulation at both intraseasonal and shorter time scales in the Bay of Bengal (BoB). The increase in SST-ISV amplitudes with diurnal forcing is approximate to 0.05 degrees C in the southern bay while it is approximate to 0.02 degrees C in the northern bay. Increased intraseasonal warming with diurnal forcing results from the increase in mixed layer heat gain from insolation, due to shoaling of the daytime mixed layer. Amplified intraseasonal cooling is dominantly controlled by the strengthening of subsurface processes owing to the nocturnal deepening of mixed layer. In the southern bay, intraseasonal variability is mainly determined by the diurnal cycle in insolation, while in the northern bay, diurnal cycle in insolation and winds have comparable contributions. Temperature inversions (TI) develop in the northern bay in the absence of diurnal variability in wind stress. In the northern bay, SST-ISV amplification is not as large as that in the southern bay due to the weaker diurnal variability of mixed layer depth (MLD) limited by salinity stratification. Diurnal variability of model MLD is not sufficient to create large modifications in mixed layer heat budget and SST-ISV in the northern bay.

Modeling non-stationarity in intensity, duration and frequency of extreme rainfall over India

Significant changes are reported in extreme rainfall characteristics over India in recent studies though there are disagreements on the spatial uniformity and causes of trends. Based on recent theoretical advancements in the Extreme Value Theory (EVT), we analyze changes in extreme rainfall characteristics over India using a high-resolution daily gridded (1 degrees latitude x 1 degrees longitude) dataset. Intensity, duration and frequency of excess rain over a high threshold in the summer monsoon season are modeled by non-stationary distributions whose parameters vary with physical covariates like the El-Nino Southern Oscillation index (ENSO-index) which is an indicator of large-scale natural variability, global average temperature which is an indicator of human-induced global warming and local mean temperatures which possibly indicate more localized changes. Each non-stationary model considers one physical covariate and the best chosen statistical model at each rainfall grid gives the most significant physical driver for each extreme rainfall characteristic at that grid. Intensity, duration and frequency of extreme rainfall exhibit non-stationarity due to different drivers and no spatially uniform pattern is observed in the changes in them across the country. At most of the locations, duration of extreme rainfall spells is found to be stationary, while non-stationary associations between intensity and frequency and local changes in temperature are detected at a large number of locations. This study presents the first application of nonstationary statistical modeling of intensity, duration and frequency of extreme rainfall over India. The developed models are further used for rainfall frequency analysis to show changes in the 100-year extreme rainfall event. Our findings indicate the varying nature of each extreme rainfall characteristic and their drivers and emphasize the necessity of a comprehensive framework to assess resulting risks of precipitation induced flooding. (C) 2014 Elsevier B.V. All rights reserved.

Investigation of orographic features influencing spatial distribution of rainfall over the Western Ghats of India using satellite data

Using remotely sensed Tropical Rainfall Measuring Mission (TRMM) 3B42 rainfall and topographic data from the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) Digital Elevation Model (DEM), the impact of oroghraphical aspects such as topography, spatial variability of elevation and altitude of apexes are examined to investigate capacious summer monsoon rainfall over the Western Ghats (WG) of India. TRMM 3B42 v7 rainfall data is validated with Indian Meteorological Department (IMD) gridded rainfall data at 0.5 degrees resolution over the WG. The analysis of spatial pattern of monsoon rainfall with orography of the WG ascertains that the grade of orographic precipitation depends mainly on topography of the mountain barrier followed by steepness of windward side slope and altitude of the mountain. Longer and broader, i.e. cascaded topography, elevated summits and gradually increasing slopes impel the enhancement in precipitation. Comparing topography of various states of the WG, it has been observed that windward side of Karnataka receives intense rainfall in the WG during summer monsoon. It has been observed that the rainfall is enhanced before the peak of the mountain and confined up to the height about 800m over the WG. In addition to this, the spatial distribution of heavy and very heavy rainfall events in the last 14 years has also been explored. Heavy and very heavy rain events on this hilly terrain are categorized with a threshold of precipitation (R) in the range 150>R>120mmday(-1) and exceeding 150mmday(-1) using probability distribution of TRMM 3B42 v7 rainfall. The areas which are prone to heavy precipitation are identified. The study would help policy makers to manage the hazard scenario and, to improve weather predictions on mountainous terrain of the WG.

Meridional gradients in aerosol vertical distribution over Indian Mainland: Observations and model simulations

Multi-year observations from the network of ground-based observatories (ARFINET), established under the project `Aerosol Radiative Forcing over India' (ARFI) of Indian Space Research Organization and space-borne lidar `Cloud Aerosol Lidar with Orthogonal Polarization' (CALIOP) along with simulations from the chemical transport model `Goddard Chemistry Aerosol Radiation and Transport' (GOCART), are used to characterize the vertical distribution of atmospheric aerosols over the Indian landmass and its spatial structure. While the vertical distribution of aerosol extinction showed higher values close to the surface followed by a gradual decrease at increasing altitudes, a strong meridional increase is observed in the vertical spread of aerosols across the Indian region in all seasons. It emerges that the strong thermal convections cause deepening of the atmospheric boundary layer, which although reduces the aerosol concentration at lower altitudes, enhances the concentration at higher elevations by pumping up more aerosols from below and also helping the lofted particles to reach higher levels in the atmosphere. Aerosol depolarization ratios derived from CALIPSO as well as the GOCART simulations indicate the dominance of mineral dust aerosols during spring and summer and anthropogenic aerosols in winter. During summer monsoon, though heavy rainfall associated with the Indian monsoon removes large amounts of aerosols, the prevailing southwesterly winds advect more marine aerosols over to landmass (from the adjoining oceans) leading to increase in aerosol loading at lower altitudes than in spring. During spring and summer months, aerosol loading is found to be significant, even at altitudes as high as 4 km, and this is proposed to have significant impacts on the regional climate systems such as Indian monsoon. (C) 2015 Elsevier Ltd. All rights reserved.

Investigations of aerosol black carbon from a semi-urban site in the Indo-Gangetic Plain region

Long-term (2009-2012) data from ground-based measurements of aerosol black carbon (BC) from a semi-urban site, Pantnagar (29.0 degrees N, 79.5 degrees E, 231 m amsl), in the Indo-Gangetic Plain (IGP) near the Himalayan foothills are analyzed to study the regional characterization. Large variations are seen in BC at both diurnal and seasonal scales, associated with the mesoscale and synoptic meteorological processes, and local/regional anthropogenic activities. BC diurnal variations show two peaks (morning and evening) arising from the combined effects of the atmospheric boundary layer (ABL) dynamics and local emissions. The diurnal amplitudes as well as the rates of diurnal evolution are the highest in winter season, followed by autumn, and the lowest in summer-monsoon. BC exhibits nearly an inverse relation with mixing layer depth in all seasons; being strongest in winter (R-2 = 0.89) and weakest (R-2 = 0.33) in monsoon (July-August). Unlike BC, co-located aerosol optical depths (AOD) and aerosol absorption are highest in spring over IGP, probably due to the presence of higher abundances of aerosols (including dust) above the ABL (in the free troposphere). AOD (500 nm) showed annual peak (>0.6) in May-June, dominated by coarse mode, while fine mode aerosols dominated in late autumn and early winter. Aerosols profiles from CALIPSO show highest values close to the surface in winter/autumn, similar to the feature seen in surface BC, whereas at altitudes > 2 km, the extinction is maximum in spring/summer. WRF-Chem model is used to simulate BC temporal variations and then compared with observed BC. The model captures most of the important features of the diurnal and seasonal variations but significantly underestimated the observed BC levels, suggesting improvements in diurnal and seasonal varying BC emissions apart from the boundary layer processes. (C) 2015 Elsevier Ltd. All rights reserved.

The impact of latent heating on the location and strength of the tropical easterly jet

The tropical easterly jet (TEJ) is a prominent atmospheric circulation feature observed during the Asian summer monsoon. It is generally assumed that sensible heating over the Tibetan Plateau directly influences the location of the TEJ. However, other studies have suggested the importance of latent heating in determining the jet location. In this paper, the relative importance of latent heating on the maintenance of the TEJ is explored through simulations with a general circulation model. The simulation of the TEJ by the Community Atmosphere Model, version 3.1 is discussed in detail. These simulations showed that the location of the TEJ is well correlated with the location of the precipitation. Significant zonal shifts in the location of the precipitation resulted in similar shifts in the zonal location of the TEJ. These zonal shifts had minimal effect on the large-scale structure of the jet. Further, provided that precipitation patterns were relatively unchanged, orography did not directly impact the location of the TEJ. These changes were robust even with changes in the cumulus parameterization. This suggests the potential important role of latent heating in determining the location and structure of the TEJ. These results were used to explain the significant differences in the zonal location of the TEJ in the years 1988 and 2002. To understand the contribution of the latitudinal location of latent heating on the strength of the TEJ, aqua-planet simulations were carried out. It has been shown that for similar amounts of net latent heating, the jet is stronger when heating is in the higher tropical latitudes. This may partly explain the reason for the jet to be very strong during the JJA monsoon season.

Earthquakes of the Nepal Himalaya : towards a physical model of the seismic cycle

Home to hundreds of millions of souls and land of excessiveness, the Himalaya is also the locus of a unique seismicity whose scope and peculiarities still remain to this day somewhat mysterious. Having claimed the lives of kings, or turned ancient timeworn cities into heaps of rubbles and ruins, earthquakes eerily inhabit Nepalese folk tales with the fatalistic message that nothing lasts forever. From a scientific point of view as much as from a human perspective, solving the mysteries of Himalayan seismicity thus represents a challenge of prime importance. Documenting geodetic strain across the Nepal Himalaya with various GPS and leveling data, we show that unlike other subduction zones that exhibit a heterogeneous and patchy coupling pattern along strike, the last hundred kilometers of the Main Himalayan Thrust fault, or MHT, appear to be uniformly locked, devoid of any of the “creeping barriers” that traditionally ward off the propagation of large events. The approximately 20 mm/yr of reckoned convergence across the Himalaya matching previously established estimates of the secular deformation at the front of the arc, the slip accumulated at depth has to somehow elastically propagate all the way to the surface at some point. And yet, neither large events from the past nor currently recorded microseismicity nearly compensate for the massive moment deficit that quietly builds up under the giant mountains. Along with this large unbalanced moment deficit, the uncommonly homogeneous coupling pattern on the MHT raises the question of whether or not the locked portion of the MHT can rupture all at once in a giant earthquake. Univocally answering this question appears contingent on the still elusive estimate of the magnitude of the largest possible earthquake in the Himalaya, and requires tight constraints on local fault properties. What makes the Himalaya enigmatic also makes it the potential source of an incredible wealth of information, and we exploit some of the oddities of Himalayan seismicity in an effort to improve the understanding of earthquake physics and cipher out the properties of the MHT. Thanks to the Himalaya, the Indo-Gangetic plain is deluged each year under a tremendous amount of water during the annual summer monsoon that collects and bears down on the Indian plate enough to pull it away from the Eurasian plate slightly, temporarily relieving a small portion of the stress mounting on the MHT. As the rainwater evaporates in the dry winter season, the plate rebounds and tension is increased back on the fault. Interestingly, the mild waggle of stress induced by the monsoon rains is about the same size as that from solid-Earth tides which gently tug at the planets solid layers, but whereas changes in earthquake frequency correspond with the annually occurring monsoon, there is no such correlation with Earth tides, which oscillate back-and-forth twice a day. We therefore investigate the general response of the creeping and seismogenic parts of MHT to periodic stresses in order to link these observations to physical parameters. First, the response of the creeping part of the MHT is analyzed with a simple spring-and-slider system bearing rate-strengthening rheology, and we show that at the transition with the locked zone, where the friction becomes near velocity neutral, the response of the slip rate may be amplified at some periods, which values are analytically related to the physical parameters of the problem. Such predictions therefore hold the potential of constraining fault properties on the MHT, but still await observational counterparts to be applied, as nothing indicates that the variations of seismicity rate on the locked part of the MHT are the direct expressions of variations of the slip rate on its creeping part, and no variations of the slip rate have been singled out from the GPS measurements to this day. When shifting to the locked seismogenic part of the MHT, spring-and-slider models with rate-weakening rheology are insufficient to explain the contrasted responses of the seismicity to the periodic loads that tides and monsoon both place on the MHT. Instead, we resort to numerical simulations using the Boundary Integral CYCLes of Earthquakes algorithm and examine the response of a 2D finite fault embedded with a rate-weakening patch to harmonic stress perturbations of various periods. We show that such simulations are able to reproduce results consistent with a gradual amplification of sensitivity as the perturbing period get larger, up to a critical period corresponding to the characteristic time of evolution of the seismicity in response to a step-like perturbation of stress. This increase of sensitivity was not reproduced by simple 1D-spring-slider systems, probably because of the complexity of the nucleation process, reproduced only by 2D-fault models. When the nucleation zone is close to its critical unstable size, its growth becomes highly sensitive to any external perturbations and the timings of produced events may therefore find themselves highly affected. A fully analytical framework has yet to be developed and further work is needed to fully describe the behavior of the fault in terms of physical parameters, which will likely provide the keys to deduce constitutive properties of the MHT from seismological observations.