Land Surface Temperature continuous time series from SSM/I passive microwave radiometers, over northern latitudes (Latitudes > 45°N), during summer snow free period (2000-2011), with links to NetCDF files

Remote sensing instruments are key players to map land surface temperature (LST) at large temporal and spatial scales. In this paper, we present how we combine passive microwave and thermal infrared data to estimate LST during summer snow-free periods over northern high latitudes. The methodology is based on the SSM/I-SSMIS 37 GHz measurements at both vertical and horizontal polarizations on a 25 km × 25 km grid size. LST is retrieved from brightness temperatures introducing an empirical linear relationship between emissivities at both polarizations as described in Royer and Poirier (2010). This relationship is calibrated at pixel scale, using cloud-free independent LST data from MODIS instruments. The SSM/I-SSMIS and MODIS data are synchronized by fitting a diurnal cycle model built on skin temperature reanalysis provided by the European Centre for Medium-Range Weather Forecasts (ECMWF). The resulting temperature dataset is provided at 25 km scale and at an hourly time step during the ten-year analysis period (2000-2011). This new product was locally evaluated at five experimental sites of the EU-PAGE21 project against air temperature measurements and meteorological model reanalysis, and compared to the MODIS LST product at both local and circumpolar scale. The results giving a mean RMSE of the order of 2.2 K demonstrate the usefulness of the microwave product, which is unaffected by clouds as opposed to thermal infrared products and offers a better resolution compared to model reanalysis.

Assessing the role of Mediterranean evergreen oaks canopy cover in land surface albedo and temperature using a remote sensing-based approach

Modifications in vegetation cover can have an impact on the climate through changes in biogeochemical and biogeophysical processes. In this paper, the tree canopy cover percentage of a savannah-like ecosystem (montado/dehesa) was estimated at Landsat pixel level for 2011, and the role of different canopy cover percentages on land surface albedo (LSA) and land surface temperature (LST) were analysed. A modelling procedure using a SGB machine-learning algorithm and Landsat 5-TM spectral bands and derived vegetation indices as explanatory variables, showed that the estimation of montado canopy cover was obtained with good agreement (R2 = 78.4%). Overall, montado canopy cover estimations showed that low canopy cover class (MT_1) is the most representative with 50.63% of total montado area. MODIS LSA and LST products were used to investigate the magnitude of differences in mean annual LSA and LST values between contrasting montado canopy cover percentages. As a result, it was found a significant statistical relationship between montado canopy cover percentage and mean annual surface albedo (R2 = 0.866, p < 0.001) and surface temperature (R2 = 0.942, p < 0.001). The comparisons between the four contrasting montado canopy cover classes showed marked differences in LSA (χ2 = 192.17, df = 3, p < 0.001) and LST (χ2 = 318.18, df = 3, p < 0.001). The highest montado canopy cover percentage (MT_4) generally had lower albedo than lowest canopy cover class, presenting a difference of −11.2% in mean annual albedo values. It was also showed that MT_4 and MT_3 are the cooler canopy cover classes, and MT_2 and MT_1 the warmer, where MT_1 class had a difference of 3.42 °C compared with MT_4 class. Overall, this research highlighted the role that potential changes in montado canopy cover may play in local land surface albedo and temperature variations, as an increase in these two biogeophysical parameters may potentially bring about, in the long term, local/regional climatic changes moving towards greater aridity.

Surface temperature measurements near Ny-Ålesund, Svalbard, Norway in 2008 with links to data files

The ground surface temperature is one of the key parameters that determine the thermal regime of permafrost soils in arctic regions. Due to remoteness of most permafrost areas, monitoring of the land surface temperature (LST) through remote sensing is desirable. However, suitable satellite platforms such as MODIS provide spatial resolutions, that cannot resolve the considerable small-scale heterogeneity of the surface conditions characteristic for many permafrost areas. This study investigates the spatial variability of summer surface temperatures of high-arctic tundra on Svalbard, Norway. A thermal imaging system mounted on a mast facilitates continuous monitoring of approximately 100 x 100 m of tundra with a wide variability of different surface covers and soil moisture conditions over the entire summer season from the snow melt until fall. The net radiation is found to be a control parameter for the differences in surface temperature between wet and dry areas. Under clear-sky conditions in July, the differences in surface temperature between wet and dry areas reach up to 10K. The spatial differences reduce strongly in weekly averages of the surface temperature, which are relevant for the soil temperature evolution of deeper layers. Nevertheless, a considerable variability remains, with maximum differences between wet and dry areas of 3 to 4K. Furthermore, the pattern of snow patches and snow-free areas during snow melt in July causes even greater differences of more than 10K in the weekly averages. Towards the end of the summer season, the differences in surface temperature gradually diminish. Due to the pronounced spatial variability in July, the accumulated degree-day totals of the snow-free period can differ by more than 60% throughout the study area. The terrestrial observations from the thermal imaging system are compared to measurements of the land surface temperature from the MODIS sensor. During periods with frequent clear-sky conditions and thus a high density of satellite data, weekly averages calculated from the thermal imaging system and from MODIS LST agree within less than 2K. Larger deviations occur when prolonged cloudy periods prevent satellite measurements. Futhermore, the employed MODIS L2 LST data set contains a number of strongly biased measurements, which suggest an admixing of cloud top temperatures. We conclude that a reliable gap filling procedure to moderate the impact of prolonged cloudy periods would be of high value for a future LST-based permafrost monitoring scheme. The occurrence of sustained subpixel variability of the summer surface temperature is a complicating factor, whose impact needs to be assessed further in conjunction with other spatially variable parameters such as the snow cover and soil properties.

Downscaled land surface temperatures obtained by a particle smoothing approach, with links to NetCDF files

This dataset contains continuous time series of land surface temperature (LST) at spatial resolution of 300m around the 12 experimental sites of the PAGE21 project (grant agreement number 282700, funded by the EC seventh Framework Program theme FP7-ENV-2011). This dataset was produced from hourly LST time series at 25km scale, retrieved from SSM/I data (André et al., 2015, doi:10.1016/j.rse.2015.01.028) and downscaled to 300m using a dynamic model and a particle smoothing approach. This methodology is based on two main assumptions. First, LST spatial variability is mostly explained by land cover and soil hydric state. Second, LST is unique for a land cover class within the low resolution pixel. Given these hypotheses, this variable can be estimated using a land cover map and a physically based land surface model constrained with observations using a data assimilation process. This methodology described in Mechri et al. (2014, doi:10.1002/2013JD020354) was applied to the ORCHIDEE land surface model (Krinner et al., 2005, doi:10.1029/2003GB002199) to estimate prior values of each land cover class provided by the ESA CCI-Land Cover product (Bontemps et al., 2013) at 300m resolution . The assimilation process (particle smoother) consists in simulating ensemble of LST time series for each land cover class and for a large number of parameter sets. For each parameter set, the resulting temperatures are aggregated considering the grid fraction of each land cover and compared to the coarse observations. Miniminizing the distance between the aggregated model solutions and the observations allow us to select the simulated LST and the corresponding parameter sets which fit the observations most closely. The retained parameter sets are then duplicated and randomly perturbed before simulating the next time window. At the end, the most likely LST of each land cover class are estimated and used to reconstruct LST maps at 300m resolution using ESA CCI-Land Cover. The resulting temperature maps on which ice pixels were masked, are provided at daily time step during the nine-year analysis period (2000-2009).

Sea surface temperature and the subsequent freshwater survival rate of some salmon stocks: a surprising link between the climate of land and sea

Previous consideration of the relationship between climate and the survival rate of Pacific salmon eggs and fry has been confined to effects of large variation in the ambient freshwater environment; e.g., stream discharge, temperature, turbidity. This analysis shows sea surface temperatures during the last year of life of maturing adult salmon are also strongly associated with the subsequent survival rate of salmon eggs and fry is fresh water, presumably through development of the future eggs or sperm. In several stocks of three species of North American salmon, the association between the "marine" climate and egg survival is stronger than, or additive to, any estimated climatic association in fresh water. This apparent and surprising link between fresh water and the distant ocean has some interesting and complex implications for management of future salmon production.

Model-based changes in global annual mean surface temperature change (Delta T_g) and radiative forcing due to land ice albedo changes (Delta R_[LI]) over the last 5 Myr, supplementary material

It is still an open question how equilibrium warming in response to increasing radiative forcing - the specific equilibrium climate sensitivity S - depends on background climate. We here present palaeodata-based evidence on the state dependency of S, by using CO2 proxy data together with a 3-D ice-sheet-model-based reconstruction of land ice albedo over the last 5 million years (Myr). We find that the land ice albedo forcing depends non-linearly on the background climate, while any non-linearity of CO2 radiative forcing depends on the CO2 data set used. This non-linearity has not, so far, been accounted for in similar approaches due to previously more simplistic approximations, in which land ice albedo radiative forcing was a linear function of sea level change. The latitudinal dependency of ice-sheet area changes is important for the non-linearity between land ice albedo and sea level. In our set-up, in which the radiative forcing of CO2 and of the land ice albedo (LI) is combined, we find a state dependence in the calculated specific equilibrium climate sensitivity, S[CO2,LI], for most of the Pleistocene (last 2.1 Myr). During Pleistocene intermediate glaciated climates and interglacial periods, S[CO2,LI] is on average ~ 45 % larger than during Pleistocene full glacial conditions. In the Pliocene part of our analysis (2.6-5 Myr BP) the CO2 data uncertainties prevent a well-supported calculation for S[CO2,LI], but our analysis suggests that during times without a large land ice area in the Northern Hemisphere (e.g. before 2.82 Myr BP), the specific equilibrium climate sensitivity, S[CO2,LI], was smaller than during interglacials of the Pleistocene. We thus find support for a previously proposed state change in the climate system with the widespread appearance of northern hemispheric ice sheets. This study points for the first time to a so far overlooked non-linearity in the land ice albedo radiative forcing, which is important for similar palaeodata-based approaches to calculate climate sensitivity. However, the implications of this study for a suggested warming under CO2 doubling are not yet entirely clear since the details of necessary corrections for other slow feedbacks are not fully known and the uncertainties that exist in the ice-sheet simulations and global temperature reconstructions are large.

On the measurement of the surface energy budget over a land surface during the summer monsoon

The measurement of surface energy balance over a land surface in an open area in Bangalore is reported. Measurements of all variables needed to calculate the surface energy balance on time scales longer than a week are made. Components of radiative fluxes are measured while sensible and latent heat fluxes are based on the bulk method using measurements made at two levels on a micrometeorological tower of 10 m height. The bulk flux formulation is verified by comparing its fluxes with direct fluxes using sonic anemometer data sampled at 10 Hz. Soil temperature is measured at 4 depths. Data have been continuously collected for over 6 months covering pre-monsoon and monsoon periods during the year 2006. The study first addresses the issue of getting the fluxes accurately. It is shown that water vapour measurements are the most crucial. A bias of 0.25% in relative humidity, which is well above the normal accuracy assumed the manufacturers but achievable in the field using a combination of laboratory calibration and field intercomparisons, results in about 20 W m(-2) change in the latent heat flux on the seasonal time scale. When seen on the seasonal time scale, the net longwave radiation is the largest energy loss term at the experimental site. The seasonal variation in the energy sink term is small compared to that in the energy source term.

Role of predictors in downscaling surface temperature to river basin in India for IPCC SRES scenarios using support vector machine

In this paper, downscaling models are developed using a support vector machine (SVM) for obtaining projections of monthly mean maximum and minimum temperatures (T-max and T-min) to river-basin scale. The effectiveness of the model is demonstrated through application to downscale the predictands for the catchment of the Malaprabha reservoir in India, which is considered to be a climatically sensitive region. The probable predictor variables are extracted from (1) the National Centers for Environmental Prediction (NCEP) reanalysis dataset for the period 1978-2000, and (2) the simulations from the third-generation Canadian Coupled Global Climate Model (CGCM3) for emission scenarios A1B, A2, B1 and COMMIT for the period 1978-2100. The predictor variables are classified into three groups, namely A, B and C. Large-scale atmospheric variables Such as air temperature, zonal and meridional wind velocities at 925 nib which are often used for downscaling temperature are considered as predictors in Group A. Surface flux variables such as latent heat (LH), sensible heat, shortwave radiation and longwave radiation fluxes, which control temperature of the Earth's surface are tried as plausible predictors in Group B. Group C comprises of all the predictor variables in both the Groups A and B. The scatter plots and cross-correlations are used for verifying the reliability of the simulation of the predictor variables by the CGCM3 and to Study the predictor-predictand relationships. The impact of trend in predictor variables on downscaled temperature was studied. The predictor, air temperature at 925 mb showed an increasing trend, while the rest of the predictors showed no trend. The performance of the SVM models that are developed, one for each combination of predictor group, predictand, calibration period and location-based stratification (land, land and ocean) of climate variables, was evaluated. In general, the models which use predictor variables pertaining to land surface improved the performance of SVM models for downscaling T-max and T-min

Validating MODIS-derived land surface evapotranspiration with in situ measurements at two AmeriFlux sites in a semiarid region

Reducing uncertainties in the estimation of land surface evapotranspiration (ET) from remote-sensing data is essential to better understand earth-atmosphere interactions. This paper demonstrates the applicability of temperature-vegetation index triangle (T-s-VI) method in estimating regional ET and evaporative fraction (EF, defined as the ratio of latent heat flux to surface available energy) from MODIS/Terra and MODIS/Aqua products in a semiarid region. We have compared the satellite-based estimates of ET and EF with eddy covariance measurements made over 4 years at two semiarid grassland sites: Audubon Ranch (AR) and Kendall Grassland (KG). The lack of closure in the eddy covariance measured surface energy components is shown to be more serious at MODIS/Aqua overpass time than that at MODIS/Terra overpass time for both AR and KG sites. The T-s-VI-derived EF could reproduce in situ EF reasonably well with BIAS and root-mean-square difference (RMSD) of less than 0.07 and 0.13, respectively. Surface net radiation has been shown to be systematically overestimated by as large as about 60 W/m(2). Satisfactory validation results of the T-s-VI-derived sensible and latent heat fluxes have been obtained with RMSD within 54 W/m(2). The simplicity and yet easy use of the T-s-VI triangle method show a great potential in estimating regional ET with highly acceptable accuracy that is of critical significance in better understanding water and energy budgets on the Earth. Nevertheless, more validation work should be carried out over various climatic regions and under other different land use/land cover conditions in the future.

Explicit representation of subgrid heterogeneity in a GCM land-surface scheme

Essery, R L H, Best, M J, Betts, R A, Cox, P M & Taylor, C M, Explicit representation of subgrid heterogeneity in a GCM land-surface scheme. Journal of Hydrometeorology 4, pp 530-543 (2003).

A Bayesian approach to estimate sensible and latent heat over vegetated land surface

Sensible and latent heat fluxes are often calculated from bulk transfer equations combined with the energy balance. For spatial estimates of these fluxes, a combination of remotely sensed and standard meteorological data from weather stations is used. The success of this approach depends on the accuracy of the input data and on the accuracy of two variables in particular: aerodynamic and surface conductance. This paper presents a Bayesian approach to improve estimates of sensible and latent heat fluxes by using a priori estimates of aerodynamic and surface conductance alongside remote measurements of surface temperature. The method is validated for time series of half-hourly measurements in a fully grown maize field, a vineyard and a forest. It is shown that the Bayesian approach yields more accurate estimates of sensible and latent heat flux than traditional methods.

Climate impacts of recent multidecadal changes in Atlantic Ocean Sea Surface Temperature: A multimodel comparison

During the twentieth century sea surface temperatures in the Atlantic Ocean exhibited prominent multidecadal variations. The source of such variations has yet to be rigorously established—but the question of their impact on climate can be investigated. Here we report on a set of multimodel experiments to examine the impact of patterns of warming in the North Atlantic, and cooling in the South Atlantic, derived from observations, that is characteristic of the positive phase of the Atlantic Multidecadal Oscillation (AMO). The experiments were carried out with six atmospheric General Circulation Models (including two versions of one model), and a major goal was to assess the extent to which key climate impacts are consistent between the different models. The major climate impacts are found over North and South America, with the strongest impacts over land found over the United States and northern parts of South America. These responses appear to be driven by a combination of an off-equatorial Gill response to diabatic heating over the Caribbean due to increased rainfall within the region and a Northward shift in the Inter Tropical Convergence Zone (ITCZ) due to the anomalous cross-equatorial SST gradient. The majority of the models show warmer US land temperatures and reduced Mean Sea Level Pressure during summer (JJA) in response to a warmer North Atlantic and a cooler South Atlantic, in line with observations. However the majority of models show no significant impact on US rainfall during summer. Over northern South America, all models show reduced rainfall in southern hemisphere winter (JJA), whilst in Summer (DJF) there is a generally an increase in rainfall. However, there is a large spread amongst the models in the magnitude of the rainfall anomalies over land. Away from the Americas, there are no consistent significant modelled responses. In particular there are no significant changes in the North Atlantic Oscillation (NAO) over the North Atlantic and Europe in Winter (DJF). Additionally, the observed Sahel drying signal in African rainfall is not seen in the modelled responses. Suggesting that, in contrast to some studies, the Atlantic Multidecadal Oscillation was not the primary driver of recent reductions in Sahel rainfall.

Dependence of Indian monsoon rainfall on moisture fluxes across the Arabian Sea and the impact of coupled model sea surface temperature biases

The Arabian Sea is an important moisture source for Indian monsoon rainfall. The skill of climate models in simulating the monsoon and its variability varies widely, while Arabian Sea cold sea surface temperature (SST) biases are common in coupled models and may therefore influence the monsoon and its sensitivity to climate change. We examine the relationship between monsoon rainfall, moisture fluxes and Arabian Sea SST in observations and climate model simulations. Observational analysis shows strong monsoons depend on moisture fluxes across the Arabian Sea, however detecting consistent signals with contemporaneous summer SST anomalies is complicated in the observed system by air/sea coupling and large-scale induced variability such as the El Niño-Southern Oscillation feeding back onto the monsoon through development of the Somali Jet. Comparison of HadGEM3 coupled and atmosphere-only configurations suggests coupled model cold SST biases significantly reduce monsoon rainfall. Idealised atmosphere-only experiments show that the weakened monsoon can be mainly attributed to systematic Arabian Sea cold SST biases during summer and their impact on the monsoon-moisture relationship. The impact of large cold SST biases on atmospheric moisture content over the Arabian Sea, and also the subsequent reduced latent heat release over India, dominates over any enhancement in the land-sea temperature gradient and results in changes to the mean state. We hypothesize that a cold base state will result in underestimation of the impact of larger projected Arabian Sea SST changes in future climate, suggesting that Arabian Sea biases should be a clear target for model development.

The relationship between land–ccean surface temperature contrast and radiative forcing

Previous work has demonstrated that observed and modeled climates show a near-time-invariant ratio of mean land to mean ocean surface temperature change under transient and equilibrium global warming. This study confirms this in a range of atmospheric models coupled to perturbed sea surface temperatures (SSTs), slab (thermodynamics only) oceans, and a fully coupled ocean. Away from equilibrium, it is found that the atmospheric processes that maintain the ratio cause a land-to-ocean heat transport anomaly that can be approximated using a two-box energy balance model. When climate is forced by increasing atmospheric CO2 concentration, the heat transport anomaly moves heat from land to ocean, constraining the land to warm in step with the ocean surface, despite the small heat capacity of the land. The heat transport anomaly is strongly related to the top-of-atmosphere radiative flux imbalance, and hence it tends to a small value as equilibrium is approached. In contrast, when climate is forced by prescribing changes in SSTs, the heat transport anomaly replaces ‘‘missing’’ radiative forcing over land by moving heat from ocean to land, warming the land surface. The heat transport anomaly remains substantial in steady state. These results are consistent with earlier studies that found that both land and ocean surface temperature changes may be approximated as local responses to global mean radiative forcing. The modeled heat transport anomaly has large impacts on surface heat fluxes but small impacts on precipitation, circulation, and cloud radiative forcing compared with the impacts of surface temperature change. No substantial nonlinearities are found in these atmospheric variables when the effects of forcing and surface temperature change are added.

Incorporation of water vapour transfer in the JULES Land Surface Model: implications for key soil variables and land surface fluxes

This study focuses on the mechanisms underlying water and heat transfer in upper soil layers, and their effects on soil physical prognostic variables and the individual components of the energy balance. The skill of the JULES (Joint UK Land Environment Simulator) land surface model (LSM) to simulate key soil variables, such as soil moisture content and surface temperature, and fluxes such as evaporation, is investigated. The Richards equation for soil water transfer, as used in most LSMs, was updated by incorporating isothermal and thermal water vapour transfer. The model was tested for three sites representative of semi-arid and temperate arid climates: the Jornada site (New Mexico, USA), Griffith site (Australia) and Audubon site (Arizona, USA). Water vapour flux was found to contribute significantly to the water and heat transfer in the upper soil layers. This was mainly due to isothermal vapour diffusion; thermal vapour flux also played a role at the Jornada site just after rainfall events. Inclusion of water vapour flux had an effect on the diurnal evolution of evaporation, soil moisture content and surface temperature. The incorporation of additional processes, such as water vapour flux among others, into LSMs may improve the coupling between the upper soil layers and the atmosphere, which in turn could increase the reliability of weather and climate predictions.