Landsat Sensors Research Articles

Extending Airborne Lidar-Derived Estimates of Forest Canopy Cover and Height Over Large Areas Using kNN With Landsat Time Series Data

Airborne light detection and ranging (lidar) data provide an accurate and consistent means to obtain reliable forest canopy cover (CC) and height measurements, which are important in determining forest stand structure, volume, and biomass. Extending CC and height measurements over larger areas by integration with satellite imagery increases the value of airborne lidar data. A typical approach has been to use multiple regression, machine-learning, or regression tree methods to determine relationships between the forest structure variables measured in the lidar data and available single-date or multitemporal Landsat sensor reflectance data, and if the relationships are strong enough, to extend those variables over much larger areas than is typically covered by lidar. Such methods can be difficult to apply because of the complexity of the extending models and algorithms, and the long processing times, which may become prohibitive. One machine-learning approach, which uses the k-nearest neighbor (kNN) algorithm, is presented here in a British Columbia, Canada forest environment. Our goal was to simplify the estimation of lidar-derived forest structure variables with available Landsat time series data and compare the results of the kNN model to the traditional multiple regression results, and to those obtained with more complex and computationally demanding random forest (RF) methods. We develop and test the kNN model with airborne lidar-derived estimates of forest CC and height in 1846 relatively young and mature forests. The best kNN model produced estimates of airborne lidar-derived CC in validation sample ( ${\text n} = 1132$ ) of mature forest with three Landsat time series variables with an ${R}^{2} = 0.74$ and RMSE of approximately 10%. This result is comparable to the results obtained in earlier work using more complex machine-learning approaches (in approximately 1/10 the time). Younger (i.e., recently disturbed) forest CC estimation was less successful in the kNN model because of the high degree of structural variability in these forests. Extending airborne lidar estimates of forest height to larger areas was also possible though, as expected, less successful than CC estimation using the kNN approach.

Assessing the role of climate and resource management on groundwater dependent ecosystem changes in arid environments with the Landsat archive

Abstract Groundwater dependent ecosystems (GDEs) rely on near-surface groundwater. These systems are receiving more attention with rising air temperature, prolonged drought, and where groundwater pumping captures natural groundwater discharge for anthropogenic use. Phreatophyte shrublands, meadows, and riparian areas are GDEs that provide critical habitat for many sensitive species, especially in arid and semi-arid environments. While GDEs are vital for ecosystem services and function, their long-term (i.e. ~ 30 years) spatial and temporal variability is poorly understood with respect to local and regional scale climate, groundwater, and rangeland management. In this work, we compute time series of NDVI derived from sensors of the Landsat TM, ETM +, and OLI lineage for assessing GDEs in a variety of land and water management contexts. Changes in vegetation vigor based on climate, groundwater availability, and land management in arid landscapes are detectable with Landsat. However, the effective quantification of these ecosystem changes can be undermined if changes in spectral bandwidths between different Landsat sensors introduce biases in derived vegetation indices, and if climate, and land and water management histories are not well understood. The objective of this work is to 1) use the Landsat 8 under-fly dataset to quantify differences in spectral reflectance and NDVI between Landsat 7 ETM + and Landsat 8 OLI for a range of vegetation communities in arid and semiarid regions of the southwestern United States, and 2) demonstrate the value of 30-year historical vegetation index and climate datasets for assessing GDEs. Specific study areas were chosen to represent a range of GDEs and environmental conditions important for three scenarios: baseline monitoring of vegetation and climate, riparian restoration, and groundwater level changes. Google's Earth Engine cloud computing and environmental monitoring platform is used to rapidly access and analyze the Landsat archive along with downscaled North American Land Data Assimilation System gridded meteorological data, which are used for both atmospheric correction and correlation analysis. Results from the cross-sensor comparison indicate a benefit from the application of a consistent atmospheric correction method, and that NDVI derived from Landsat 7 and 8 are very similar within the study area. Results from continuous Landsat time series analysis clearly illustrate that there are strong correlations between changes in vegetation vigor, precipitation, evaporative demand, depth to groundwater, and riparian restoration. Trends in summer NDVI associated with riparian restoration and groundwater level changes were found to be statistically significant, and interannual summer NDVI was found to be moderately correlated to interannual water-year precipitation for baseline study sites. Results clearly highlight the complementary relationship between water-year PPT, NDVI, and evaporative demand, and are consistent with regional vegetation index and complementary relationship studies. This work is supporting land and water managers for evaluation of GDEs with respect to climate, groundwater, and resource management.

TIME SERIES ANALYSIS OF REMOTE SENSING OBSERVATIONS FOR CITRUS CROP GROWTH STAGE AND EVAPOTRANSPIRATION ESTIMATION

Satellite based earth observation (EO) platforms have proved capability to spatio-temporally monitor changes on the earth's surface. Long term satellite missions have provided huge repository of optical remote sensing datasets, and United States Geological Survey (USGS) Landsat program is one of the oldest sources of optical EO datasets. This historical and near real time EO archive is a rich source of information to understand the seasonal changes in the horticultural crops. Citrus (Mandarin / Nagpur Orange) is one of the major horticultural crops cultivated in central India. Erratic behaviour of rainfall and dependency on groundwater for irrigation has wide impact on the citrus crop yield. Also, wide variations are reported in temperature and relative humidity causing early fruit onset and increase in crop water requirement. Therefore, there is need to study the crop growth stages and crop evapotranspiration at spatio-temporal scale for managing the scarce resources. In this study, an attempt has been made to understand the citrus crop growth stages using Normalized Difference Time Series (NDVI) time series data obtained from Landsat archives (<a href="http://earthexplorer.usgs.gov/"target="_blank">http://earthexplorer.usgs.gov/</a>). Total 388 Landsat 4, 5, 7 and 8 scenes (from year 1990 to Aug. 2015) for Worldwide Reference System (WRS) 2, path 145 and row 45 were selected to understand seasonal variations in citrus crop growth. Considering Landsat 30 meter spatial resolution to obtain homogeneous pixels with crop cover orchards larger than 2 hectare area was selected. To consider change in wavelength bandwidth (radiometric resolution) with Landsat sensors (i.e. 4, 5, 7 and 8) NDVI has been selected to obtain continuous sensor independent time series. The obtained crop growth stage information has been used to estimate citrus basal crop coefficient information (Kcb). Satellite based Kcb estimates were used with proximal agrometeorological sensing system observed relevant weather parameters for crop ET estimation. The results show that time series EO based crop growth stage estimates provide better information about geographically separated citrus orchards. Attempts are being made to estimate regional variations in citrus crop water requirement for effective irrigation planning. In future high resolution Sentinel 2 observations from European Space Agency (ESA) will be used to fill the time gaps and to get better understanding about citrus crop canopy parameters.

TIME SERIES ANALYSIS OF REMOTE SENSING OBSERVATIONS FOR CITRUS CROP GROWTH STAGE AND EVAPOTRANSPIRATION ESTIMATION

Abstract. Satellite based earth observation (EO) platforms have proved capability to spatio-temporally monitor changes on the earth's surface. Long term satellite missions have provided huge repository of optical remote sensing datasets, and United States Geological Survey (USGS) Landsat program is one of the oldest sources of optical EO datasets. This historical and near real time EO archive is a rich source of information to understand the seasonal changes in the horticultural crops. Citrus (Mandarin / Nagpur Orange) is one of the major horticultural crops cultivated in central India. Erratic behaviour of rainfall and dependency on groundwater for irrigation has wide impact on the citrus crop yield. Also, wide variations are reported in temperature and relative humidity causing early fruit onset and increase in crop water requirement. Therefore, there is need to study the crop growth stages and crop evapotranspiration at spatio-temporal scale for managing the scarce resources. In this study, an attempt has been made to understand the citrus crop growth stages using Normalized Difference Time Series (NDVI) time series data obtained from Landsat archives (http://earthexplorer.usgs.gov/). Total 388 Landsat 4, 5, 7 and 8 scenes (from year 1990 to Aug. 2015) for Worldwide Reference System (WRS) 2, path 145 and row 45 were selected to understand seasonal variations in citrus crop growth. Considering Landsat 30 meter spatial resolution to obtain homogeneous pixels with crop cover orchards larger than 2 hectare area was selected. To consider change in wavelength bandwidth (radiometric resolution) with Landsat sensors (i.e. 4, 5, 7 and 8) NDVI has been selected to obtain continuous sensor independent time series. The obtained crop growth stage information has been used to estimate citrus basal crop coefficient information (Kcb). Satellite based Kcb estimates were used with proximal agrometeorological sensing system observed relevant weather parameters for crop ET estimation. The results show that time series EO based crop growth stage estimates provide better information about geographically separated citrus orchards. Attempts are being made to estimate regional variations in citrus crop water requirement for effective irrigation planning. In future high resolution Sentinel 2 observations from European Space Agency (ESA) will be used to fill the time gaps and to get better understanding about citrus crop canopy parameters.

FEASIBILITY STUDY OF LANDSAT-8 IMAGERY FOR RETRIEVING SEA SURFACE TEMPERATURE (CASE STUDY PERSIAN GULF)

Abstract. Sea surface temperature (SST) is one of the critical parameters in marine meteorology and oceanography. The SST datasets are incorporated as conditions for ocean and atmosphere models. The SST needs to be investigated for various scientific phenomenon such as salinity, potential fishing zone, sea level rise, upwelling, eddies, cyclone predictions. On the other hands, high spatial resolution SST maps can illustrate eddies and sea surface currents. Also, near real time producing of SST map is suitable for weather forecasting and fishery applications. Therefore satellite remote sensing with wide coverage of data acquisition capability can use as real time tools for producing SST dataset. Satellite sensor such as AVHRR, MODIS and SeaWIFS are capable of extracting brightness values at different thermal spectral bands. These brightness temperatures are the sole input for the SST retrieval algorithms. Recently, Landsat-8 successfully launched and accessible with two instruments on-board: (1) the Operational Land Imager (OLI) with nine spectral bands in the visual, near infrared, and the shortwave infrared spectral regions; and (2) the Thermal Infrared Sensor (TIRS) with two spectral bands in the long wavelength infrared. The two TIRS bands were selected to enable the atmospheric correction of the thermal data using a split window algorithm (SWA). The TIRS instrument is one of the major payloads aboard this satellite which can observe the sea surface by using the split-window thermal infrared channels (CH10: 10.6 μm to 11.2 μm; CH11: 11.5 μm to 12.5 μm) at a resolution of 30 m. The TIRS sensors have three main advantages comparing with other previous sensors. First, the TIRS has two thermal bands in the atmospheric window that provide a new SST retrieval opportunity using the widely used split-window (SW) algorithm rather than the single channel method. Second, the spectral filters of TIRS two bands present narrower bandwidth than that of the thermal band on board on previous Landsat sensors. Third, TIRS is one of the best space born and high spatial resolution with 30 m. in this regards, Landsat-8 can use the Split-Window (SW) algorithm for retrieving SST dataset. Although several SWs have been developed to use with other sensors, some adaptations are required in order to implement them for the TIRS spectral bands. Therefore, the objective of this paper is to develop a SW, adapted for use with Landsat-8 TIRS data, along with its accuracy assessment. In this research, that has been done for modelling SST using thermal Landsat 8-imagery of the Persian Gulf. Therefore, by incorporating contemporary in situ data and SST map estimated from other sensors like MODIS, we examine our proposed method with coefficient of determination (R2) and root mean square error (RMSE) on check point to model SST retrieval for Landsat-8 imagery. Extracted results for implementing different SW's clearly shows superiority of utilized method by R2 = 0.95 and RMSE = 0.24.

FEASIBILITY STUDY OF LANDSAT-8 IMAGERY FOR RETRIEVING SEA SURFACE TEMPERATURE (CASE STUDY PERSIAN GULF)

Sea surface temperature (SST) is one of the critical parameters in marine meteorology and oceanography. The SST datasets are incorporated as conditions for ocean and atmosphere models. The SST needs to be investigated for various scientific phenomenon such as salinity, potential fishing zone, sea level rise, upwelling, eddies, cyclone predictions. On the other hands, high spatial resolution SST maps can illustrate eddies and sea surface currents. Also, near real time producing of SST map is suitable for weather forecasting and fishery applications. Therefore satellite remote sensing with wide coverage of data acquisition capability can use as real time tools for producing SST dataset. Satellite sensor such as AVHRR, MODIS and SeaWIFS are capable of extracting brightness values at different thermal spectral bands. These brightness temperatures are the sole input for the SST retrieval algorithms. Recently, Landsat-8 successfully launched and accessible with two instruments on-board: (1) the Operational Land Imager (OLI) with nine spectral bands in the visual, near infrared, and the shortwave infrared spectral regions; and (2) the Thermal Infrared Sensor (TIRS) with two spectral bands in the long wavelength infrared. The two TIRS bands were selected to enable the atmospheric correction of the thermal data using a split window algorithm (SWA). The TIRS instrument is one of the major payloads aboard this satellite which can observe the sea surface by using the split-window thermal infrared channels (CH10: 10.6 μm to 11.2 μm; CH11: 11.5 μm to 12.5 μm) at a resolution of 30 m. The TIRS sensors have three main advantages comparing with other previous sensors. First, the TIRS has two thermal bands in the atmospheric window that provide a new SST retrieval opportunity using the widely used split-window (SW) algorithm rather than the single channel method. Second, the spectral filters of TIRS two bands present narrower bandwidth than that of the thermal band on board on previous Landsat sensors. Third, TIRS is one of the best space born and high spatial resolution with 30 m. in this regards, Landsat-8 can use the Split-Window (SW) algorithm for retrieving SST dataset. Although several SWs have been developed to use with other sensors, some adaptations are required in order to implement them for the TIRS spectral bands. Therefore, the objective of this paper is to develop a SW, adapted for use with Landsat-8 TIRS data, along with its accuracy assessment. In this research, that has been done for modelling SST using thermal Landsat 8-imagery of the Persian Gulf. Therefore, by incorporating contemporary in situ data and SST map estimated from other sensors like MODIS, we examine our proposed method with coefficient of determination (R2) and root mean square error (RMSE) on check point to model SST retrieval for Landsat-8 imagery. Extracted results for implementing different SW's clearly shows superiority of utilized method by R<sup>2</sup> = 0.95 and RMSE = 0.24.

Towards decadal soil salinity mapping using Landsat time series data

Salinization is one of the major soil problems around the world. However, decadal variation in soil salinization has not yet been extensively reported. This study exploited thirty years (1985–2015) of Landsat sensor data, including Landsat-4/5 TM (Thematic Mapper), Landsat-7 ETM+ (Enhanced Thematic Mapper Plus) and Landsat-8 OLI (Operational Land Imager), for monitoring soil salinity of the Yellow River Delta, China. The data were initially corrected for atmospheric effects, and then matched the spectral bands of EO-1 (Earth Observing One) ALI (Advanced Land Imager). Subsequently, soil salinity maps were derived with a previously developed PLSR (Partial Least Square Regression) model. On intra-annual scale, the retrievals showed that soil salinity increased in February, stabilized in March, and decreased in April. On inter-annual scale, soil salinity decreased within 1985–2000 (−0.74gkg−1/10a, p<0.001), and increased within 2000–2015 (0.79gkg−1/10a, p<0.001). Our study presents a new perspective for use of multiple Landsat data in soil salinity retrieval, and further the understanding of soil salinization development over the Yellow River Delta.

Including land cover change in analysis of greenness trends using all available Landsat 5, 7, and 8 images: A case study from Guangzhou, China (2000–2014)

Abstract Remote sensing has proven a useful way of evaluating long-term trends in vegetation “greenness” through the use of vegetation indices like Normalized Differences Vegetation Index (NDVI) and Enhanced Vegetation Index (EVI). In particular, analyses of greenness trends have been performed for large areas (continents, for example) in an attempt to understand vegetation response to climate. These studies have been most often used coarse resolution sensors like Moderate Resolution Image Spectroradiometer (MODIS) and Advanced Very High Resolution Radiometer (AVHRR). However, trends in greenness are also important at more local scales, particularly in and around cities as vegetation offers a variety of valuable ecosystem services ranging from minimizing air pollution to mitigating urban heat island effects. To explore the ability to monitor greenness trends in and around cities, this paper presents a new way for analyzing greenness trends based on all available Landsat 5, 7, and 8 images and applies it to Guangzhou, China. This method is capable of including the effects of land cover change in the evaluation of greenness trends by separating the effects of abrupt and gradual changes, and providing information on the timing of greenness trends. An assessment of the consistency of surface reflectance from Landsat 8 with past Landsat sensors indicates biases in the visible bands of Landsat 8, especially the blue band. Landsat 8 NDVI values were found to have a larger bias than the EVI values; therefore, EVI was used in the analysis of greenness trends for Guangzhou. In spite of massive amounts of development in Guangzhou from 2000 to 2014, greenness was found to increase, mostly as a result of gradual change. Comparison of the greening magnitudes estimated from the approach presented here and a Simple Linear Trend (SLT) method indicated large differences for certain time intervals as the SLT method does not include consideration for abrupt land cover changes. Overall, this analysis demonstrates the importance of considering land cover change when analyzing trends in greenness from satellite time series in areas where land cover change is common.

Mangrove response to environmental change in Australia's Gulf of Carpentaria.

Across their range, mangroves are responding to coastal environmental change. However, separating the influence of human activities from natural events and processes (including that associated with climatic fluctuation) is often difficult. In the Gulf of Carpentaria, northern Australia (Leichhardt, Nicholson, Mornington Inlet, and Flinders River catchments), changes in mangroves are assumed to be the result of natural drivers as human impacts are minimal. By comparing classifications from time series of Landsat sensor data for the period 1987–2014, mangroves were observed to have extended seawards by up to 1.9 km (perpendicular to the coastline), with inland intrusion occurring along many of the rivers and rivulets in the tidal reaches. Seaward expansion was particularly evident near the mouth of the Leichhardt River, and was associated with peaks in river discharge with LiDAR data indicating distinct structural zones developing following each large rainfall and discharge event. However, along the Gulf coast, and particularly within the Mornington Inlet catchment, the expansion was more gradual and linked to inundation and regular sediment supply through freshwater input. Landward expansion along the Mornington Inlet catchment was attributed to the combined effects of sea level rise and prolonged periods of tidal and freshwater inundation on coastal lowlands. The study concluded that increased amounts of rainfall and associated flooding and sea level rise were responsible for recent seaward and landward extension of mangroves in this region.

Perspectives on monitoring gradual change across the continuity of Landsat sensors using time-series data

Abstract There are many types of changes occurring over the Earth's landscapes that can be detected and monitored using Landsat data. Here we focus on monitoring “within-state,” gradual changes in vegetation in contrast with traditional monitoring of “abrupt” land-cover conversions. Gradual changes result from a variety of processes, such as vegetation growth and succession, damage from insects and disease, responses to shifts in climate, and other factors. Despite the prevalence of gradual changes across the landscape, they are largely ignored by the remote sensing community. Gradual changes are best characterized and monitored using time-series analysis, and with the successful launch of Landsat 8 we now have appreciable data continuity that extends the Landsat legacy across the previous 43 years. In this study, we conducted three related analyses: (1) comparison of spectral values acquired by Landsats 7 and 8, separated by eight days, to ensure compatibility for time-series evaluation; (2) tracking of multitemporal signatures for different change processes across Landsat 5, 7, and 8 sensors using anniversary-date imagery; and (3) tracking the same type of processes using all available acquisitions. In this investigation, we found that data representing natural vegetation from Landsats 5, 7, and 8 were comparable and did not indicate a need for major modification prior to use for long-term monitoring. Analyses using anniversary-date imagery can be very effective for assessing long term patterns and trends occurring across the landscape, and are especially good for providing insights regarding trends related to long-term and continuous trends of growth or decline. We found that use of all available data provided a much more comprehensive level of understanding of the trends occurring, providing information about rate, duration, and intra- and inter-annual variability that could not be readily gleaned from the anniversary date analyses. We observed that using all available clear Landsat 5–8 observations with the new Continuous Change Detection and Classification (CCDC) algorithm was very effective for illuminating vegetation trends. There are a number of potential challenges for assessing gradual changes, including atmospheric impacts, algorithm development and visualization of the changes. One of the biggest challenges for studying gradual change will be the lack of appropriate data for validating results and products.

Early spring post-fire snow albedo dynamics in high latitude boreal forests using Landsat-8 OLI data.

Taking advantage of the improved radiometric resolution of Landsat-8 OLI which, unlike previous Landsat sensors, does not saturate over snow, the progress of fire recovery progress at the landscape scale (< 100m) is examined. High quality Landsat-8 albedo retrievals can now capture the true reflective and layered character of snow cover over a full range of land surface conditions and vegetation densities. This new capability particularly improves the assessment of post-fire vegetation dynamics across low- to high- burn severity gradients in Arctic and boreal regions in the early spring, when the albedos during recovery show the greatest variation. We use 30 m resolution Landsat-8 surface reflectances with concurrent coarser resolution (500m) MODIS high quality full inversion surface Bidirectional Reflectance Distribution Functions (BRDF) products to produce higher resolution values of surface albedo. The high resolution full expression shortwave blue sky albedo product performs well with an overall RMSE of 0.0267 between tower and satellite measures under both snow-free and snow-covered conditions. While the importance of post-fire albedo recovery can be discerned from the MODIS albedo product at regional and global scales, our study addresses the particular importance of early spring post-fire albedo recovery at the landscape scale by considering the significant spatial heterogeneity of burn severity, and the impact of snow on the early spring albedo of various vegetation recovery types. We found that variations in early spring albedo within a single MODIS gridded pixel can be larger than 0.6. Since the frequency and severity of wildfires in Arctic and boreal systems is expected to increase in the coming decades, the dynamics of albedo in response to these rapid surface changes will increasingly impact the energy balance and contribute to other climate processes and physical feedback mechanisms. Surface radiation products derived from Landsat-8 data will thus play an important role in characterizing the carbon cycle and ecosystem processes of high latitude systems.

Optimizing the number of training areas for modeling above-ground biomass with ALS and multispectral remote sensing in subtropical Nepal

Remote sensing-based inventories of above-ground forest biomass (AGB) require a set of training plots representative of the area to be studied, the collection of which is the most expensive part of the analysis. These are time-consuming and costly because the large variety in forest conditions requires more plots to adequately capture this variability. A field campaign in general is challenging and is hampered by the complex topographic conditions, limited accessibility, steep mountainous terrains which increase labor efforts and costs. In addition it is also depend on the ratio between size of study area and number of training plots. In this study, we evaluate the number of training areas (sample size) required to estimate AGB for an area in the southern part of Nepal using airborne laser scanning (ALS), RapidEye and Landsat data. Three experiments were conducted: (i) AGB model performance, based on all the field training plots; (ii) reduction of the sample size, based on the ALS metrics and the AGB distribution; and (iii) prediction of the optimal number of training plots, based on the correlation between the remote sensing and field data. The AGB model was fitted using the sparse Bayesian method. AGB model performance was validated using an independent validation dataset. The effect of the strategies for reducing the sample size was readily apparent for the ALS-based AGB prediction, but the RapidEye and Landsat sensor data failed to capture any such effect. The results indicate that adequate coverage of the variability in tree height and density was an important condition for selecting the training plots. In addition, the ALS-based AGB prediction required the smallest number of training plots and was also quite stable with a small number of field plots.

Spatial prediction of soil surface texture in a semiarid region using random forest and multiple linear regressions

Soil texture is an essential and extremely variable physical property that strongly influences many other soil properties that are highly relevant for agricultural production, e.g., fertility and water retention capacity. In plain areas, terrain properties derived from a digital elevation model are not effective for digital soil mapping, and the variation in the properties of such areas remains a challenge. In this regard, remote sensing can facilitate the mapping of soil properties. The purpose of this study was to evaluate the efficiency of using of data obtained from the Thematic Mapper (TM) sensor of Landsat 5 for digital soil mapping in a semi-arid region, based on multiple linear regression (MLR) and a random forest model (RFM). To this end, 399 samples of the soil surface layer (0–20cm) were used to predict the sand, silt and clay contents, using the bands 1, 2, 3, 4, 5 and 7, the Normalized Difference Vegetation Index (NDVI), the grain size index (GSI), and the relationships between bands 3 and 2, bands 3 and 7, and bands 5 and 7 (clay index) of the Landsat 5 TM sensor as covariates. Among these covariates, only band 1 (b1), the relationship between bands 5 and 7 (b5/b7) for sand, silt and clay, and band 4 (b4) for silt were not significantly correlated according to Pearson's correlation analysis. The validation of the models showed that the properties were best estimated using the RFM, which explained 63% and 56% of the spatial variability of sand and clay, respectively, whereas the MLR explained 30% of the spatial variation of silt. An analysis of the relevance of the variables predicted by the RFM showed that the covariates b3/b7, b5, NDVI and b2 explained most of the variability of the considered properties. The RFM proved to be more advantageous than the MLR with respect to insensitivity to overfitting and the presence of noise in the data. In addition, the RFM produced more realistic distribution maps of the soil properties than did the MLR, taking into account that the estimated values of the soil attributes were in the same range as the calibration data, while the MLR model estimates were out of the range of the calibration data.

A time series analysis of urbanization induced land use and land cover change and its impact on land surface temperature with Landsat imagery

Abstract The relationship between land use and land cover (LULC) patterns and thermal characteristic has long been studied for examining the impact of urbanization on urban thermal environment. Previous studies were inclined to the use of one or a few images without fully considering the temporal domain of the reflective and thermal infrared data on board Landsat sensors. This paper took the Atlanta metropolitan area as a case study to illustrate LULC change and its impact on land surface temperature (LST) variations. The Landsat L1T (Standard terrain correction) images from TM/ETM + from 1984 to 2011 (507 in total) were downloaded through the USGS online portal and consistently calibrated to surface reflectance and brightness temperature (BT). The cloud-, cloud shadow-, and snow-contaminated pixels were excluded in the analysis according to the metadata, and a further screening procedure based on the RIRLS (Robust Iteratively Reweighted Least Squares) technique was implemented. The time series LSTs (TSLSTs) was derived using the single channel algorithm because it required only the parameters of water vapor and land surface emissivity and had a reported error close to 1 K. The LULC classification and change detection was accomplished by using the Continuous Change Detection and Classification (CCDC) algorithm. The TSLSTs were further decomposed into the seasonal and trend components by an additive model. Results showed that the overall LULC classification and change detection accuracies were 89% and 92%, respectively. Urban growth was mainly observed in Fulton County and Gwinnett County. High-intensity urban land had the largest mean LST value of 294.9 K and yearly amplitude of 17.4 K. A comparison of the trend component between urban and non-urban land covers showed a difference of 1.8 K per decade. Next, temporal thermal signatures were created to characterize and quantify the impact of urban LULC changes. Decomposition analysis showed that the conversion of evergreen forest to medium-intensity urban land generated the largest difference in annual LST variation (5.7 K) and the largest trend difference (0.0004 K/day).

Characterization of Landsat-7 to Landsat-8 reflective wavelength and normalized difference vegetation index continuity.

At over 40 years, the Landsat satellites provide the longest temporal record of space-based land surface observations, and the successful 2013 launch of the Landsat-8 is continuing this legacy. Ideally, the Landsat data record should be consistent over the Landsat sensor series. The Landsat-8 Operational Land Imager (OLI) has improved calibration, signal to noise characteristics, higher 12-bit radiometric resolution, and spectrally narrower wavebands than the previous Landsat-7 Enhanced Thematic Mapper (ETM+). Reflective wavelength differences between the two Landsat sensors depend also on the surface reflectance and atmospheric state which are difficult to model comprehensively. The orbit and sensing geometries of the Landsat-8 OLI and Landsat-7 ETM+ provide swath edge overlapping paths sensed only one day apart. The overlap regions are sensed in alternating backscatter and forward scattering orientations so Landsat bi-directional reflectance effects are evident but approximately balanced between the two sensors when large amounts of time series data are considered. Taking advantage of this configuration a total of 59 million 30m corresponding sensor observations extracted from 6,317 Landsat-7 ETM+ and Landsat-8 OLI images acquired over three winter and three summer months for all the conterminous United States (CONUS) are compared. Results considering different stages of cloud and saturation filtering, and filtering to reduce one day surface state differences, demonstrate the importance of appropriate per-pixel data screening. Top of atmosphere (TOA) and atmospherically corrected surface reflectance for the spectrally corresponding visible, near infrared and shortwave infrared bands, and derived normalized difference vegetation index (NDVI), are compared and their differences quantified. On average the OLI TOA reflectance is greater than the ETM+ TOA reflectance for all bands, with greatest differences in the near-infrared (NIR) and the shortwave infrared bands due to the quite different spectral response functions between the sensors. The atmospheric correction reduces the mean difference in the NIR and shortwave infrared but increases the mean difference in the visible bands. Regardless of whether TOA or surface reflectance are used to generate NDVI, on average, for vegetated soil and vegetation surfaces (0 ≤ NDVI ≤ 1), the OLI NDVI is greater than the ETM+ NDVI. Statistical functions to transform between the comparable sensor bands and sensor NDVI values are presented so that the user community may apply them in their own research to improve temporal continuity between the Landsat-7 ETM+ and Landsat-8 OLI sensor data. The transformation functions were developed using ordinary least squares (OLS) regression and were fit quite reliably (r 2 values >0.7 for the reflectance data and >0.9 for the NDVI data, p-values <0.0001).

Digital mapping of corrosion risk in coastal urban areas using remote sensing and structural condition assessment: case study in cyprus

AbstractAtmospheric corrosion is one of the main factors leading to performance deterioration of reinforced concrete buildings; and, hence, periodic structural condition monitoring is required to assess and repair the adverse effects of corrosion. However, this can become a cumbersome and expensive task to undertake for large populations of buildings, scattered in large urban areas. To optimize the use of available resources, appropriate tools are required for the assessment of corrosion risk of reinforced concrete construction. This paper proposes a framework for the production of digital corrosion risk maps for urban areas; Cyprus was used as a case study. This framework explored multi-temporal satellite remote sensing data from the Landsat sensors as well as corrosion risk factors derived from the results of a recently completed research project, entitled “STEELCOR”. This framework was used to develop two corrosion risk scenarios within Geographical Information Systems, and to produce corrosion risk maps for three coastal cities of Cyprus. The thematic maps indicated that, for slight corrosion damage, the distance of reinforced concrete buildings from the coast was more influential than the building age. While, for significant corrosion damage, the maps indicated that the age of RC buildings was more influential than the distance from the coast.

A New Geostatistical Solution to Remote Sensing Image Downscaling

The availability of the panchromatic (PAN) band in remote sensing images gives birth to so-called image fusion techniques for increasing the spatial resolution of images to that of the PAN band. The spatial resolution of such spatially sharpened images, such as for the MODIS and Landsat sensors, however, may not be sufficient to provide the required detailed land-cover/land-use information. This paper proposes an area-to-point regression kriging (ATPRK)-based geostatistical solution to increase the spatial resolution of remote sensing images beyond that of any input images, including the PAN band. The new approach is a two-stage approach, including covariate downscaling and ATPRK-based image fusion. The new approach treats the PAN band as the covariate and takes advantages of its textural information. It explicitly accounts for the size of support, spatial correlation, and the point spread function of the sensor and has the characteristic of perfect coherence with the original coarse data. Moreover, the new downscaling approach can be extended readily by incorporating other ancillary information. The proposed approach was examined using both Landsat and MODIS images. The results show that it can produce more accurate sharpened images than four benchmark approaches.

A Generic Algorithm to Estimate LAI, FAPAR and FCOVER Variables from SPOT4_HRVIR and Landsat Sensors: Evaluation of the Consistency and Comparison with Ground Measurements

The leaf area index (LAI) and the fraction of photosynthetically active radiation absorbed by green vegetation (FAPAR) are essential climatic variables in surface process models. FCOVER is also important to separate vegetation and soil for energy balance processes. Currently, several LAI, FAPAR and FCOVER satellite products are derived moderate to coarse spatial resolution. The launch of Sentinel-2 in 2015 will provide data at decametric resolution with a high revisit frequency to allow quantifying the canopy functioning at the local to regional scales. The aim of this study is thus to evaluate the performances of a neural network based algorithm to derive LAI, FAPAR and FCOVER products at decametric spatial resolution and high temporal sampling. The algorithm is generic, i.e., it is applied without any knowledge of the landcover. A time series of high spatial resolution SPOT4_HRVIR (16 scenes) and Landsat 8 (18 scenes) images acquired in 2013 over the France southwestern site were used to generate the LAI, FAPAR and FCOVER products. For each sensor and each biophysical variable, a neural network was first trained over PROSPECT+SAIL radiative transfer model simulations of top of canopy reflectance data for green, red, near-infra red and short wave infra-red bands. Our results show a good spatial and temporal consistency between the variables derived from both sensors: almost half the pixels show an absolute difference between SPOT and LANDSAT estimates of lower that 0.5 unit for LAI, and 0.05 unit for FAPAR and FCOVER. Finally, downward-looking digital hemispherical cameras were completed over the main land cover types to validate the accuracy of the products. Results show that the derived products are strongly correlated with the field measurements (R2 &gt; 0.79), corresponding to a RMSE = 0.49 for LAI, RMSE = 0.10 (RMSE = 0.12) for black-sky (white sky) FAPAR and RMSE = 0.15 for FCOVER. It is concluded that the proposed generic algorithm provides a good basis to monitor the seasonal variation of the vegetation biophysical variables for important crops at decametric resolution.

Multi-Image and Multi-Sensor Change Detection for Long-Term Monitoring of Arid Environments With Landsat Series

An automated procedure has been proposed to monitor by multispectral satellite imagery the cultivation expansion between 1987 and 2013 in the arid environment of the Fayyum Oasis (Egypt), which is subject to land reclamation. A change detection procedure was applied to the four years investigated (1987, 1998, 2003 and 2013). This long-term analysis is based on images from the Landsat series, adopting a classification strategy relying on vegetation index computations. In particular: (a) the consequences of the radiometric differences of three Landsat sensors on the vegetation index values were analyzed using data simulated by a hyperspectral Hyperion image; (b) the problems resulting from harvesting cycles were minimized using five images per year, after a preliminary analysis on the effects deriving from the number of processed images; (c) an accuracy assessment was carried out on the 2003 and 2013 maps using high resolution images for a portion of the investigated area, with an estimated overall accuracy of 91% for the change detection. The method is implemented in a batch procedure and can be applied to other similar environmental contexts, supporting analyses for sustainable development and exploitation of soil and water resources.

Comparing the MLC and JavaNNS Approaches in Classifying Multi-Temporal LANDSAT Satellite Imagery over an Ephemeral River Area

This paper analyzes two pixel-based classification approaches to support the analysis of land cover transformations based on multitemporal LANDSAT sensor data covering a time space of about 24 years. The research activity presented in this paper was carried out using Lama San Giorgio (Bari, Italy) catchment area as a study case, being this area prone to flooding as proved by its geological and hydrological characteristics and by the significant number of floods occurred in the past. Land cover classes were defined in accordance with on the CN method with the aim of characterizing land use based on attitude to generate runoff. Two different classifiers, i.e. Maximum Likelihood Classifier (MLC) and Java Neural Network Simulator (JavaNNS) models, were compared. The Artificial Neural Networks (ANN) approach was found to be the most reliable and efficient when lacking ground reference data and a priori knowledge on input data distribution.