Air Pollution Modeling from Remotely Sensed Data Using Regression Techniques

There is a need for timely information about changes in the air pollution levels in cities for adopting precautionary measures. Keeping this in view, an attempt has been made to develop a model which will be useful to obtain air quality information directly from remotely sensed data easily and quickly. For this study pixel values, vegetation indices and urbanization index from IRS P6 LISS IV and Landsat ETM+ images were used to develop regression based models with Air Pollution Index (API), which were calculated from in-situ air pollutant information. It was found that among the 12 parameters of IRS, highest correlation exists between pixel values in NIR (Near Infra-Red) band (Pearson correlation ?0.77) and Normalized Difference Vegetation Index (NDVI) (Pearson correlation ?0.68) and both have inverse relationship with API. In case of Landsat, the highest correlation was observed in SWIR (Short Wave Infra-Red) band (Pearson correlation ?0.83) and NIR (Pearson correlation ?0.78). Both single and multivariate regression models were calibrated from best correlated variables from IRS and Landsat. Among all the models, multivariate regression model from Landsat with four most correlated variables gave the most accurate air pollution image. On comparison between the API modeled and API interpolated images, 90.5 % accuracy was obtained.     

Evaluating uncertainty in Landsat-derived postfire recovery metrics due to terrain,soil, and shrub type variations in southern California

ABSTRACT

Temporal trajectories of apparent vegetation abundance based on the multi-decadal Landsat image series provide valuable information on the postfire recovery of chaparral shrublands, which tend to mature within one decade. Signals of change in fractional shrub cover (FSC) extracted from time-sequential Normalized Difference Vegetation Index (NDVI) data can be systematically biased due to spatial variation in shrub type, soil substrate, or illumination differences associated with topography. We evaluate the effects of these variables in Landsat-derived metrics of FSC and postfire recovery, based upon three chaparral sites in southern California which contain shrub community ecotones, complex terrain, and soil variations. Detailed validations of prefire and postfire FSC are based on high spatial resolution ortho-imagery; cross-stratified random sampling is used for variable control. We find that differences in the composition and structure of shrubs (inferred from ortho-imagery) can substantially influence FSC-NDVI relations and impact recovery metrics. Differences in soil type have a moderate effect on the FSC-NDVI relation in one of the study sites, while no substantial effects were observed due to variation of terrain illumination among the study sites. Arithmetic difference recovery metrics – based on NDVI values that were not normalized with unburned control plots – correlate in a moderate but significant manner with a change in FSC (R ² values range 0.47–0.59 at two sites). Similar regression coefficients resulted from using Landsat visible reflectance data alone. The lowest correlations to FSC resulted from Soil-Adjusted Vegetation Index (SAVI) and are attributed to the effects of the soil-adjustment factor in sparsely vegetated areas. The Normalized Burn Ratio and Normalized Burn Ratio 2 showed a moderate correlation to FSC. This study confirms the utility of Landsat NDVI data for postfire recovery evaluation and implies a need for stratified analysis of postfire recovery in some chaparral landscapes.     

Estuarine sediment deposition during wetland restoration: A GIS and remote sensing modeling approach

Restoration is currently underway in the industrial salt flats of San Francisco Bay, California. Remote sensing of suspended sediment concentration and other GIS predictor variables were used to model sediment deposition within recently restored ponds. Suspended sediment concentrations were calibrated to reflectance values from Landsat TM 5 and ASTER satellite image data using three statistical techniques—linear regression, multivariate regression, and Artificial Neural Network (ANN) regression. Multivariate and ANN regressions using ASTER proved to be the most accurate methods, yielding r² values of 0.88 and 0.87, respectively. Predictor variables such as sediment grain size and tidal frequency were used in the Marsh Sedimentation (MARSED) model for predicting deposition rates. MARSED results show a root mean square deviation (RMSD) of 66.8 mm (<1σ) between modeled and field observations. This model was applied to a pond breached in November 2010 and indicated that the pond will reach sediment equilibrium levels after 60 months of tidal inundation.     

A new method for generating a clear-sky Landsat composite for cropland from cloud-contaminated Landsat-7 and Landsat-8 images

A new method was developed in this study for producing a clear-sky Landsat composite for cropland from cloud-contaminated Landsat images acquired in a short time period. It used Thiel–Sen regression to normalize all Landsat scenes to a MODIS image to make all Landsat images radiometrically consistent and comparable. Pixel selection criteria combining the modified maximum vegetation index and the modified minimum visible reflectance selection methods were designed to enhance the pixel selection of land/water over cloud/shadow in the image compositing. The advantages of the method include (1) avoiding complicated atmospheric corrections but with reliable surface reflectance results, (2) being insensitive to errors induced by image co-registration uncertainties between Landsat and MODIS images, (3) avoiding the lack of samples for the regression analysis using the full Landsat scenes (rather than overlay regions), and (4) enhancing cloud/shadow detection. The composite image has MODIS-like surface reflectance, thus making MODIS algorithms applicable for retrieving biophysical parameters. The method was automatically implemented on a set of 13 cloud-contaminated (>39%) Landsat-7 (Scan-Line Corrector-Off) and Landsat-8 scenes acquired during peak growing season in a crop region of Manitoba, Canada. The result was a 95.8% cloud-free image. The method can also substantially increase the usage of cloud-contaminated Landsat data.     

Spatially Explicit Estimation of Leaf Area Index Using EO-1 Hyperion and Landsat ETM+ Data: Implications of Spectral Bandwidth and Shortwave Infrared Data on Prediction Accuracy in a Tropical Montane Environment

This study evaluated the utility of narrowband (EO-1 Hyperion) and broadband (Landsat ETM+) remote sensing data for the estimation of leaf area index (LAI) in a tropical environment in Sulawesi, Indonesia. LAI was inferred from canopy gap fraction measurements taken in natural tropical forest and cocoa plantations. Single and multiple spectral bands and spectral indices were used as predictor variables in reduced major axis (RMA) and ordinary least squares (OLS) regression models. The predictive power of most regression models was notably higher when employing narrowband data instead of broadband data. Highly significant relationships between LAI and spectral reflectance were observed near the red-edge region and in most shortwave infrared (SWIR) bands. In contrast to most near-infrared (NIR) narrow bands, the correlation between SWIR reflectance and LAI was not confounded when including both vegetation types and did not suffer from saturation. The results demonstrate that leaf area index of a challenging tropical environment can be estimated with satisfactory accuracy from hyperspectral remote sensing data.     

Estimation of Crown Closure and Tree Density Using Landsat TM Satellite Images in Mixed Forest Stands

The objective of this study was to investigate the relationship between crown closure and tree density in mixed forest stands using Landsat Thematic Mapper (TM) reflectance values (TM 1- TM 5 and TM 7) and six vegetation indices (SR, DVI, SAVI, NDVI, TVI and NLI). In this study, multiple regression analysis was used to estimate the relationships between the crown closure and tree density (number of tree stems per hectare) using reflectance values and vegetation indices (VIs). The results demonstrated that the model that used SR and DVI had the best performances in terms of crown closure (R²?=?0.674) and the model that used the DVI and SAVI had the best performances in terms of tree density (R²?=?0.702). The regression model that used TM 1, TM 3 together with TM 4 showed the performances of the crown closure (R²?=?0.610) and the regression model that used TM 1 showed the performances of the tree density (0.613). Results obtained from this research show that vegetation indices (VIs) were a better predictor of crown closure and tree density than other TM bands.     

Comparison of the Landsat Surface Reflectance Climate Data Record (CDR) and manually atmospherically corrected data in a semi-arid European study area

This study contributes to the quality assessment of atmospherically corrected Landsat surface reflectance data that are routinely generated by the Landsat Ecosystem Disturbance Adaptive Processing System (LEDAPS). This dataset, named Landsat Surface Reflectance Climate Data Record (Landsat CDR), is available at global scale and offers unprecedented opportunities to land monitoring and management services that require atmospherically corrected Earth observation (EO) data. Our assessment is based on the comparison of the Landsat CDR data against a set of Landsat and DEIMOS-1 images processed to a high degree of accuracy using an industry-standard atmospheric correction algorithm (ATCOR-2). The software package has been used for many years and its correction procedures can be considered consolidated and well-established. The dataset of Landsat and DEIMOS-1 images was acquired over a semi-arid agricultural area located in Lower Austria and was independently corrected by using a manual fine-tuning of ATCOR-2 parameters to reach the highest possible accuracy. Results show a very good correspondence of the surface reflectance in each of the six reflective spectral channels as well as for the NDVI (Normalized Difference Vegetation Index). An additional comparison against a NDVI time series from MODIS revealed also a good correspondence. Coefficients of determination (R²) between the two multi-year and multi-seasonal Landsat/DEIMOS datasets range between 0.91 (blue band) and 0.98 (nIR, SWIR-1 and SWIR-2). The results obtained for our semi-arid test site in Austria confirm previous findings and suggest that automatic atmospheric procedures, such as the one implemented by LEDAPS are accurate enough to be used in land monitoring services that require consistent multi-temporal surface reflectance data.     

Integrating age in the detection and mapping of incongruous patches in coffee (Coffea arabica) plantations using multi-temporal Landsat 8 NDVI anomalies

The development of cost-effective, reliable and easy to implement crop condition monitoring methods is urgently required for perennial tree crops such as coffee (Coffea arabica), as they are grown over large areas and represent long term and higher levels of investment. These monitoring methods are useful in identifying farm areas that experience poor crop growth, pest infestation, diseases outbreaks and/or to monitor response to management interventions. This study compares field level coffee mean NDVI and LSWI anomalies and age-adjusted coffee mean NDVI and LSWI anomalies in identifying and mapping incongruous patches across perennial coffee plantations. To achieve this objective, we first derived deviation of coffee pixels from the global coffee mean NDVI and LSWI values of nine sequential Landsat 8 OLI image scenes. We then evaluated the influence of coffee age class (young, mature and old) on Landsat-scale NDVI and LSWI values using a one-way ANOVA and since results showed significant differences, we adjusted NDVI and LSWI anomalies for age-class. We then used the cumulative inverse distribution function (α ≤ 0.05) to identify fields and within field areas with excessive deviation of NDVI and LSWI from the global and the age-expected mean for each of the Landsat 8 OLI scene dates spanning three seasons. Results from accuracy assessment indicated that it was possible to separate incongruous and healthy patches using these anomalies and that using NDVI performed better than using LSWI for both global and age-adjusted mean anomalies. Using the age-adjusted anomalies performed better in separating incongruous and healthy patches than using the global mean for both NDVI (Overall accuracy = 80.9% and 68.1% respectively) and for LSWI (Overall accuracy = 68.1% and 48.9% respectively). When applied to other Landsat 8 OLI scenes, the results showed that the proportions of coffee fields that were modelled incongruent decreased with time for the young age category and while it increased for the mature and old age classes with time. We concluded that the method could be useful for the identification of anomalous patches using Landsat scale time series data to monitor large coffee plantations and provide an indication of areas requiring particular field attention.     

Temporal replication of the national land cover database using active machine learning

Regional and national level land cover datasets, such as the National Land Cover Database (NLCD) in the United States, have become an important resource in physical and social science research. Updates to the NLCD have been conducted every 5 years since 2001; however, the procedure for producing a new release is labor-intensive and time-consuming, taking 3 or 4 years to complete. Furthermore, in most countries very few, if any, such releases exist, and thus there is high demand for efficient production of land cover data at different points in time. In this paper, an active machine learning framework for temporal updating (or backcasting) of land cover data is proposed and tested for three study sites covered by the NLCD. The approach employs a maximum entropy classifier to extract information from one Landsat image using the NLCD, and then replicate the classification on a Landsat image for the same geographic extent from a different point in time to create land cover data of similar quality. Results show that this framework can effectively replicate the land cover database in the temporal domain with similar levels of overall and within class agreement when compared against high resolution reference land cover datasets. These results demonstrate that the land cover information encapsulated in the NLCD can effectively be extracted using solely Landsat imagery for replication purposes. The algorithm is fully automated and scalable for applications at landscape and regional scales for multiple points in time.     

Remote estimation of water quality parameters of Himalayan lake (Kashmir) using Landsat 8 OLI imagery

The overarching aim of this study was to derive simple and accurate algorithms for the retrieval of water quality parameters for the Wular Lake using Landsat 8 OLI satellite data. The water quality parameters include pH, COD, DO, alkalinity, hardness, chloride, TDS, total suspended solids (TSS), turbidity, electric conductivity and phosphate. Regression analysis was performed using atmospherically corrected true reflectance values of original OLI bands, images after applying enhancement techniques (NDVI, principal components) and the values of the water quality parameters at different sample locations to obtain the empirical relationship. Most of the parameters were well correlated with single OLI bands with R² greater than 0.5, whereas phosphate showed a good correlation with NDVI image. The parameters like pH and DO showed a good relation with the principal component I and IV, respectively. The high concentration of pH, COD, turbidity and TSS and low concentration of DO infers the anthropogenic impact on lake.     

Formulation of Time Series Vegetation Index from Indian Geostationary Satellite and Comparison with Global Product

To study impact of climate change on vegetation time series vegetation index has a vital role to know the behaviour of vegetation dynamics over a time period. INSAT 3A CCD (Charged Couple Device) is the only geostationary sensor to acquire regular coverage of Asia continent at 1 km × 1 km spatial resolution with high temporal frequency (half-an-hour). A formulation of surface reflectances in red, near infrared (NIR), short wave infrared (SWIR) and NDVI from INSAT 3A CCD has been defined and integrated in the operational chain. The atmospheric correction of at-sensor reflectances using SMAC (Simple Model for Atmospheric Correction) model improved the NDVI by 5–40% and also increased its dynamic range. The temporal dynamics of 16-day NDVI composite at 0500 GMT for a growing year (June 2008–March 2009) showed matching profiles with reference to global products (MODIS TERRA) over known land targets. The root mean square deviation (RMSD) between the two was 0.14 with correlation coefficient (r) 0.84 from 200 paired datasets. This inter-sensor cross-correlation would help in NDVI calibration to add continuity in long term NDVI database for climate change studies.     

Local-scale agricultural drought monitoring with satellite-based multi-sensor time-series

ABSTRACT

Globally, drought constitutes a serious threat to food and water security. The complexity and multivariate nature of drought challenges its assessment, especially at local scales. The study aimed to assess spatiotemporal patterns of crop condition and drought impact at the spatial scale of field management units with a combined use of time-series from optical (Landsat, MODIS, Sentinel-2) and Synthetic Aperture Radar (SAR) (Sentinel 1) data. Several indicators were derived such as Normalized Difference Vegetation Index (NDVI), Normalized Difference Moisture Index (NDMI), Land Surface Temperature (LST), Tasseled cap indices and Sentinel-1 based backscattering intensity and relative surface moisture. We used logistic regression to evaluate the drought-induced variability of remotely sensed parameters estimated for different phases of crop growth. The parameters with the highest prediction rate were further used to estimate thresholds for drought/non-drought classification. The models were evaluated using the area under the receiver operating characteristic curve and validated with in-situ data. The results revealed that not all remotely sensed variables respond in the same manner to drought conditions. Growing season maximum NDVI and NDMI (70–75%) and SAR derived metrics (60%) reflect specifically the impact of agricultural drought. These metrics also depict stress affected areas with a larger spatial extent. LST was a useful indicator of crop condition especially for maize and sunflower with prediction rates of 86% and 71%, respectively. The developed approach can be further used to assess crop condition and to support decision-making in areas which are more susceptible and vulnerable to drought.     

基于Landsat TM的城市热岛效应与地表特征参数稳健关系模型

 首先,利用Landsat TM热红外影像结合地面气象观测资料反演地面温度,揭示了济南市夏季城市热岛效应| 然后,基于稳 健的LTS与最小二乘回归（LS）分析探讨了城乡地面热辐射与地表特征参数的线性变化趋势,认为植被指数（NDVI、SAVI和TCG）、 湿度指数（NDMI和TCW）以及近红外反照率与地表温度的变化趋势相反,亮度指数（NDBI和TCB）和可见光反照率与地表温度的变化 趋势一致,而短光波段反照率与地表温度不存在明显相关趋势。研究结果表明,NDMI能很好地解释地表温度变化,且最为稳健; 其次是NDVI、SAVI、TCG和NDBI,它们对地表温度的解释程度高且稳健性较强; 可见光反照率虽能较好解释地表温度,但其稳健性 较差; 近红外反照率、TCW和TCB对地表温度的解释程度和稳健性相对较低。     

Estimation of green grass/herb biomass from airborne hyperspectral imagery using spectral indices and partial least squares regression

The main objective was to determine whether partial least squares (PLS) regression improves grass/herb biomass estimation when compared with hyperspectral indices, that is normalised difference vegetation index (NDVI) and red-edge position (REP). To achieve this objective, fresh green grass/herb biomass and airborne images (HyMap) were collected in the Majella National Park, Italy in the summer of 2005. The predictive performances of hyperspectral indices and PLS regression models were then determined and compared using calibration (n = 30) and test (n = 12) data sets. The regression model derived from NDVI computed from bands at 740 and 771 nm produced a lower standard error of prediction (SEP = 264 g m⁻²) on the test data compared with the standard NDVI involving bands at 665 and 801 nm (SEP = 331 g m⁻²), but comparable results with REPs determined by various methods (SEP = 261 to 295 g m⁻²). PLS regression models based on original, derivative and continuum-removed spectra produced lower prediction errors (SEP = 149 to 256 g m⁻²) compared with NDVI and REP models. The lowest prediction error (SEP = 149 g m⁻², 19% of mean) was obtained with PLS regression involving continuum-removed bands. In conclusion, PLS regression based on airborne hyperspectral imagery provides a better alternative to univariate regression involving hyperspectral indices for grass/herb biomass estimation in the Majella National Park.     

Normalised difference spectral indices and urban land cover as indicators of land surface temperature (LST)

Land cover changes associated with urbanisation modify microclimate, leading to urban heat islands, whereby cities are warmer than the surrounding countryside. Understanding the factors causing this phenomenon could help urban areas adapt to climate change and improve living conditions of inhabitants. In this study, land surface temperatures (LST) of Aarhus, a city in the high latitudes, are estimated from the reflectance of a thermal band (TIRS1; Band 10; 10.60–11.19 μm) of Landsat 8 on five dates in the summer months (one in 2015, and four in 2018). Spectral indices, modelled on the normalised difference vegetation index (NDVI), using all combinations of the first seven bands of Landsat 8 are calculated and their relationships with LST, analysed. Land cover characteristics, in terms of the percentages of tree cover, building cover and overall vegetation cover are estimated from airborne LiDAR data, building footprints and 4-band aerial imagery, respectively. The correlations between LST, the spectral indices and land cover are estimated.The difference in mean temperature between the rural and urban parts of Aarhus is up to 3.96 °C, while the difference between the warmer and colder zones (based on the mean and SD of LST) is up to 13.26 °C. The spectral index using the near infrared band (NIR; Band 5; 0.85-0.88 μm) and a short-wave infrared band (SWIR2; Band 7; 2.11–2.29 μm) has the strongest correlations (r: 0.62 to 0.89) with LST for the whole study area. This index is the inverse of normalised burn ratio (NBR), which has been used for mapping burnt areas. Spectral indices using different combinations of the infrared bands have stronger correlations with LST than the more widely used vegetation indices such as NDVI. The percentage of tree cover has a higher negative correlation (Pearson’s r: -0.68 to -0.75) with LST than overall vegetation cover (r: -0.45 to -0.63). Tree cover and building cover (r: 0.53 to 0.71) together explain up to 68 % of the variation in LST. Modification of tree and building cover may therefore have the potential to regulate urban LST.     

融合多源遥感数据的高分辨率城市植被覆盖度估算

准确获取城市植被覆盖定量信息对城市生态环境评价,城市规划及可持续城市发展具有重要意义.遥感技术的发展为获取区域及全球植被覆盖信息提供了有效手段,目前基于单传感器、单时相遥感数据的城市植被覆盖度估算方法得到较为广泛的应用.然而,由于城市地表覆盖的复杂性、植被类型的多样性,在一定程度上影响了城市植被覆盖信息提取的精度.为此...     

融合高时空分辨率数据估算植被净初级生产力

植被净初级生产力NPP(Net Primary Production)遥感估算与分析,有赖于高时空分辨率的遥感数据,但目前中高分辨率的遥感数据受卫星回访周期及天气的影响,在中国南方地区难以获取连续时间序列的数据,从而影响了高精度的区域植被净初级生产力的遥感估算。为此,提出一种基于多源遥感数据时空融合技术与CASA模型估算高时空分辨率NPP的方法。首先,利用多源遥感数据,即Landsat8 OLI数据与MODIS13Q1数据,采用遥感数据时空融合方法,获得了时间序列的Landsat8 OLI融合数据;然后,基于Landsat8 OLI时空融合数据,并采用CASA模型,以长株潭城市群核心区为例,进行区域植被NPP的遥感估算。研究结果表明,基于时间序列Landsat融合数据估算的30m分辨率的NPP具有良好的空间细节信息,且估算值与实测值的相关系数达0.825,与实测NPP数据保持了较好的一致性。     

Wetland Monitoring,Serving as an Index of Land Use Change-A Study in Samaspur Wetlands,Uttar Pradesh,India

Wetlands are among the most productive ecosystems in the world and any alterations might lead to changes in their bio-physical, socio-economic and climatic conditions. Wetland dynamics as an index of land use change were studied. Satellite remote sensing was utilized to understand the periodic and seasonal dynamics of Samaspur wetlands using Landsat and RESOURCESAT-1 temporal data. Index-based (i.e., Normalized Difference Water Index (NDWI) and Normalized Difference Vegetation Index (NDVI)) classification resulted in meaningful discrimination of wetland classes. Results indicate (i) effective water spread areas have increased to optimum capacity at 1990 due to the influence of Sharda canal, (ii) expansion of the agricultural area has led to reduction of the wetland buffer area, and (iii) increase in vegetation biomass due to pesticide-fertilizer runoff and sedimentation load. We also reiterate (i) free availability of the Landsat satellite data in public domain facilitating such monitoring studies and (ii) availability and utility of SWIR band information in wetland classification exercise. The study concludes that policy-driven measures have both long and short term impacts on land use and its natural wetland ecosystems; and the characterizing the later serves as indictor of the former and perhaps vice versa.     

Fusing MODIS with Landsat 8 data to downscale weekly normalized difference vegetation index estimates for central Great Basin rangelands,USA

Data fused from distinct but complementary satellite sensors mitigate tradeoffs that researchers make when selecting between spatial and temporal resolutions of remotely sensed data. We integrated data from the Moderate Resolution Imaging Spectroradiometer (MODIS) sensor aboard the Terra satellite and the Operational Land Imager sensor aboard the Landsat 8 satellite into four regression-tree models and applied those data to a mapping application. This application produced downscaled maps that utilize the 30-m spatial resolution of Landsat in conjunction with daily acquisitions of MODIS normalized difference vegetation index (NDVI) that are composited and temporally smoothed. We produced four weekly, atmospherically corrected, and nearly cloud-free, downscaled 30-m synthetic MODIS NDVI predictions (maps) built from these models. Model results were strong with R² values ranging from 0.74 to 0.85. The correlation coefficients (r ≥ 0.89) were strong for all predictions when compared to corresponding original MODIS NDVI data. Downscaled products incorporated into independently developed sagebrush ecosystem models yielded mixed results. The visual quality of the downscaled 30-m synthetic MODIS NDVI predictions were remarkable when compared to the original 250-m MODIS NDVI. These 30-m maps improve knowledge of dynamic rangeland seasonal processes in the central Great Basin, United States, and provide land managers improved resource maps.