Influence of landscape spatial heterogeneity on the non-linear estimation of leaf area index from moderate spatial resolution remote sensing data

The monitoring of earth surface dynamic processes requires global observations of the structure and the functioning of vegetation. Moderate resolution sensors (with pixel size ranging from 250 m to 7 km) provide frequent estimates of biophysical variables to characterize vegetation such as the leaf area index (LAI). However, the computation of LAI from moderate resolution remote sensing data induces a scaling bias on the LAI estimate if the moderate resolution pixel is heterogeneous and if the transfer function that relates remote sensing data to LAI is non-linear.This study provides a model to evaluate and correct the scaling bias. The model is built first for a univariate semi-empirical transfer function relating LAI directly to NDVI. The scaling bias is a function of (i) the degree of non-linearity of the transfer function quantified by its second derivative and (ii) the spatial heterogeneity of the moderate resolution pixel quantified by the variogram of the high spatial resolution (20 m) NDVI image. Then, the model is extended to a bivariate transfer function where LAI is related to red and near infrared reflectances. The scaling bias depends on (i) the Hessian matrix of the transfer function and (ii) the variograms and cross variogram of the red and near infrared reflectances.The scaling bias is investigated on several distinct landscapes from the VALERI database. Adjusting for scaling bias is critical on crop sites which are the most heterogeneous sites at the landscape level. Regarding the univariate transfer function, the magnitude of the scaling bias increases rapidly with pixel size until this size is larger than the typical spatial scale of the data. For the bivariate transfer function, it results from the addition of several components that may add up or cancel each other out. It is thus more difficult to analyze.The accuracy of the model to estimate the scaling bias is discussed. It depends mainly on the ability of the variograms and cross variogram to represent the local dispersion variances and covariance within the moderate resolution pixel. The model is generally highly accurate at 1000 m spatial resolution for the univariate transfer function and less accurate for the bivariate transfer function.     

Stereo Image Analysis of Non-Lambertian Surfaces

Stereo image analysis is based on establishing correspondences between a pair of images by determining similarity measures for potentially corresponding image parts. Such similarity criteria are only strictly valid for surfaces with Lambertian (diffuse) reflectance characteristics. Specular reflections are viewpoint dependent and may thus cause large intensity differences at corresponding image points. In the presence of specular reflections, traditional stereo approaches are often unable to establish correspondences at all, or the inferred disparity values tend to be inaccurate, or the established correspondences do not belong to the same physical surface point. The stereo image analysis framework for non-Lambertian surfaces presented in this contribution combines geometric cues with photometric and polarimetric information into an iterative scheme that allows to establish stereo correspondences in accordance with the specular reflectance behaviour and at the same time to determine the surface gradient field based on the known photometric and polarimetric reflectance properties. The described approach yields a dense 3D reconstruction of the surface which is consistent with all observed geometric and photopolarimetric data. Initially, a sparse 3D point cloud of the surface is computed by traditional blockmatching stereo. Subsequently, a dense 3D profile of the surface is determined in the coordinate system of camera 1 based on the shape from photopolarimetric reflectance and depth technique. A synthetic image of the surface is rendered in the coordinate system of camera 2 using the illumination direction and reflectance properties of the surface material. Point correspondences between the rendered image and the observed image of camera 2 are established with the blockmatching technique. This procedure yields an increased number of 3D points of higher accuracy, compared to the initial 3D point cloud. The improved 3D point cloud is used to compute a refined dense 3D surface profile. These steps are iterated until convergence of the 3D reconstruction. An experimental evaluation of our method is provided for areas of several square centimetres of forged and cast iron objects with rough surfaces displaying both diffuse and significant specular reflectance components, where traditional stereo image analysis largely fails. A comparison to independently measured ground truth data reveals that the root-mean-square error of the 3D reconstruction results is typically of the order 30–100 μm at a lateral pixel resolution of 86 μm. For two example surfaces, the number of stereo correspondences established by the specular stereo algorithm is several orders of magnitude higher than the initial number of 3D points. For one example surface, the number of stereo correspondences decreases by a factor of about two, but the 3D point cloud obtained with the specular stereo method is less noisy, contains a negligible number of outliers, and shows significantly more surface detail than the initial 3D point cloud. For poorly known reflectance parameters we observe a graceful degradation of the accuracy of 3D reconstruction.     

Sub-pixel water temperature estimation from thermal-infrared imagery using vectorized lake features

Water skin temperature derived from thermal infrared satellite data are used in a wide variety of studies. Many of these studies would benefit from frequent, high spatial resolution (100 m pixels) thermal imagery but currently, at any given location, such data are only available every few weeks from spaceborne sensors such as ASTER. Lower spatial resolution (1 km pixels) thermal imagery is available multiple times per day at any given location, from several sensors such as MODIS on board both the AQUA and TERRA satellite platforms. In order to fully exploit lower spatial resolution imagery, a sub-pixel unmixing technique has been developed and tested at Quesnel Lake, British Columbia, Canada. This approach produces accurate, frequent high spatial resolution water skin temperature maps by exploiting a priori knowledge of water boundaries derived from vectorized water features. The pixel water-fraction maps are then input to a gradient descent algorithm to solve the mixed pixel ground leaving radiance equation for sub-pixel water temperature. Ground-leaving radiance is estimated from standard temperature and emissivity data products for pure pixels and a simple regression technique to estimate atmospheric effects. In this test case, MODIS 1 km thermal imagery was used along with 1:50,000 water features to create a high-resolution (100 m) water skin temperature map. This map is compared to a concurrent ASTER temperature image and found to be within 1 °C of the ASTER skin temperature 99% of the time. This is a considerable improvement over the 2.55 °C difference between the original MODIS product and ASTER image due to land temperature contamination. The algorithm is simple, effective, and unlocks a largely untapped resource for limnological and hydrological studies.     

Appropriate scale of soil moisture retrieval from high resolution radar imagery for bare and minimally vegetated soils

This research investigates the appropriate scale for watershed averaged and site specific soil moisture retrieval from high resolution radar imagery. The first approach involved filtering backscatter for input to a retrieval model that was compared against field measures of soil moisture. The second approach involved spatially averaging raw and filtered imagery in an image-based statistical technique to determine the best scale for site-specific soil moisture retrieval. Field soil moisture was measured at 1225 m² sites in three watersheds commensurate with 7 m resolution Radarsat image acquisition. Analysis of speckle reducing block median filters indicated that 5 × 5 filter level was the optimum for watershed averaged estimates of soil moisture. However, median filtering alone did not provide acceptable accuracy for soil moisture retrieval on a site-specific basis. Therefore, spatial averaging of unfiltered and median filtered power values was used to generate backscatter estimates with known confidence for soil moisture retrieval. This combined approach of filtering and averaging was demonstrated at watersheds located in Arizona (AZ), Oklahoma (OK) and Georgia (GA). The optimum ground resolution for AZ, OK and GA study areas was 162 m, 310 m, and 1131 m respectively obtained with unfiltered imagery. This statistical approach does not rely on ground verification of soil moisture for validation and only requires a satellite image and average roughness parameters of the site. When applied at other locations, the resulting optimum ground resolution will depend on the spatial distribution of land surface features that affect radar backscatter. This work offers insight into the accuracy of soil moisture retrieval, and an operational approach to determine the optimal spatial resolution for the required application accuracy.     

Effect of remote sensing spatial resolution on interpreting tower-based flux observations

Validation comparisons between satellite-based surface energy balance models and tower-based flux measurements over heterogeneous landscapes can be strongly influenced by the spatial resolution of the remote sensing inputs. In this paper, a two-source energy balance model developed to use thermal and visible /near-infrared remotely sensed data is applied to Landsat imagery collected during the 2004 Soil Moisture Experiment (SMEX04) conducted in southern Arizona. Using a two dimensional flux-footprint algorithm, modeled surface fluxes are compared to tower measurements at three locations in the SMEX04 study area: two upland sites, and one riparian site. The effect of pixel resolution on evaluating the performance of the land surface model and interpreting spatial variations of land surface fluxes over these heterogeneous areas is evaluated. Three Landsat scenes were examined, one representing the dry season and the other two representing the relatively wet monsoon season. The model was run at three resolution scales: namely the Landsat visible/near-infrared band resolution (30 m), the Landsat 5 thermal band resolution (120 m), and 960 m, which is nominally the MODIS thermal resolution at near-nadir. Comparisons between modeled and measured fluxes at the three tower sites showed good agreement at the 30 m and 120 m resolutions — pixel scales at which the source area influencing the tower measurement (∼ 100 m) is reasonably resolved. At 960 m, the agreement is relatively poor, especially for the latent heat flux, due to sub-pixel heterogeneity in land surface conditions at scales exceeding the tower footprint. Therefore in this particular landscape, thermal data at 1-km resolution are not useful in assessing the intrinsic accuracy of the land-surface model in comparison with tower fluxes. Furthermore, important spatial patterns in the landscape are lost at this resolution. Currently, there are no definite plans supporting high resolution thermal data with regular global coverage below ∼ 700 m after Landsat 5 and ASTER fail. This will be a serious problem for the application and validation of thermal-based land-surface models over heterogeneous landscapes.     

基于TM影像的典型内陆淡水湿地水体提取研究

水是维系湿地生态系统稳定和健康的决定性因子,利用卫星遥感影像快速、准确地提取湿地水体信息已经成为湿地调查、研究与保护的重要手段。鉴于TM遥感影像具有较高的空间分辨率、波谱分辨率、极为丰富的信息量、较高的定位精度和相对较低的价格,其必然成为近一段时期湿地调查、研究与保护的重要数据源之一。研究基于TM遥感影像,运用多种方法针对典型内陆淡水湿地的水体信息进行了提取实验,通过对实验结果的分析得出：在面积的准确性、提取的准确度以及视觉效果3种指标下,光谱分类法较其它方法效果要好,其次为单波段阈值分析法与植被指数法,较差的是多波段谱间关系法与水体指数法;影响提取效果的主要原因是湿地水体提取不够完全,这是由影像的分辨率及湿地特殊的水文条件所造成的,采用像元分解及多源遥感数据融合技术将成为提高水体提取精度的重要手段。     

Merging complementary remote sensing datasets in the context of snow water equivalent reconstruction

Time series of snow covered area (SCA) estimates from the Moderate Resolution Imaging Spectroradiometer (MODIS) and Landsat Enhanced Thematic Mapper (ETM+) were merged with a spatially explicit snowmelt model to reconstruct snow water equivalent (SWE) in the Rio Grande headwaters (3419 km²). A linear optimization scheme was used to derive SCA estimates that preserve the statistical moments of the higher spatial resolution (i.e. 30 m) ETM+ data and resolve the superior temporal signal (i.e. ∼ daily) of the MODIS data. It was found that merging the two SCA products led to an 8% decrease and an 18% increase in the basinwide SWE in 2001 and 2002, respectively, compared to the SWE estimated from ETM+ only. Relative to SWE simulations using only ETM+ data, the hybrid SCA estimates reduced the mean absolute SWE error by 17 and 84% in 2001 and 2002, respectively; errors were determined using intensive snow survey data and two separate methods of scaling snow survey field measurements of SWE to the 1-km model pixel resolution. SWE bias for both years was reduced by 49% and skewness was reduced from − 0.78 to 0.49. These results indicate that the hybrid SWE was closer to being an unbiased estimate of the measured SWE and errors were distributed more normally. The accuracy of the SCA estimates is likely dependent on the vegetation fraction.     

Extending the MODIS 1 km ocean colour atmospheric correction to the MODIS 500 m bands and 500 m chlorophyll-a estimation towards coastal and estuarine monitoring

National and regional obligations to control and maintain water quality have led to an increase in coastal and estuarine monitoring. A potentially valuable tool is high temporal and spatial resolution satellite ocean colour data. NASA's MODIS-Terra and -Aqua can capture data at both 250 m and 500 m spatial resolutions and the existence of two sensors provides the possibility for multiple daily passes over a scene. However, no robust atmospheric correction method currently exists for these data, rendering them unusable for quantitative monitoring applications. Therefore, this paper presents an automatic and dynamic atmospheric correction approach allowing the determination of ocean colour. The algorithm is based around the standard MODIS 1 km atmospheric correction, includes cloud masking and is applicable to all of the visible 500 m bands. Comparison of the 500 m ocean colour data with the standard 1 km data shows good agreement and these results are further supported by in situ data comparisons. In addition, a novel method to produce 500 m chlorophyll-a estimates is presented. Comparisons of the 500 m estimates with the standard MODIS OC3M algorithm and to in situ data from a near-coast validation site are given.     

Multi-temporal monitoring of wetland water levels in the Florida Everglades using interferometric synthetic aperture radar (InSAR)

Interferometric synthetic aperture radar (InSAR) techniques can successfully detect phase variations related to the water level changes in wetlands and produce spatially detailed high-resolution maps of water level changes. Despite the vast details, the usefulness of the wetland InSAR observations is rather limited, because hydrologists and water resources managers need information on absolute water level values and not on relative water level changes. We present an InSAR technique called Small Temporal Baseline Subset (STBAS) for monitoring absolute water level time series using radar interferograms acquired successively over wetlands. The method uses stage (water level) observation for calibrating the relative InSAR observations and tying them to the stage's vertical datum. We tested the STBAS technique with two-year long Radarsat-1 data acquired during 2006-2008 over the Water Conservation Area 1 (WCA1) in the Everglades wetlands, south Florida (USA). The InSAR-derived water level data were calibrated using 13 stage stations located in the study area to generate 28 successive high spatial resolution maps (50 m pixel resolution) of absolute water levels. We evaluate the quality of the STBAS technique using a root mean square error (RMSE) criterion of the difference between InSAR observations and stage measurements. The average RMSE is 6.6 cm, which provides an uncertainty estimation of the STBAS technique to monitor absolute water levels. About half of the uncertainties are attributed to the accuracy of the InSAR technique to detect relative water levels. The other half reflects uncertainties derived from tying the relative levels to the stage stations' datum.     

Utility of an image-based canopy reflectance modeling tool for remote estimation of LAI and leaf chlorophyll content at the field scale

This paper presents a physically-based approach for estimating critical variables describing land surface vegetation canopies, relying on remotely sensed data that can be acquired from operational satellite sensors. The REGularized canopy reFLECtance (REGFLEC) modeling tool couples leaf optics (PROSPECT), canopy reflectance (ACRM), and atmospheric radiative transfer (6SV1) model components, facilitating the direct use of at-sensor radiances in green, red and near-infrared wavelengths for the inverse retrieval of leaf chlorophyll content (C_ab) and total one-sided leaf area per unit ground area (LAI). The inversion of the canopy reflectance model is constrained by assuming limited variability of leaf structure, vegetation clumping, and leaf inclination angle within a given crop field and by exploiting the added radiometric information content of pixels belonging to the same field. A look-up-table with a suite of pre-computed spectral reflectance relationships, each a function of canopy characteristics, soil background effects and external conditions, is accessed for fast pixel-wise biophysical parameter retrievals. Using 1 m resolution aircraft and 10 m resolution SPOT-5 imagery, REGFLEC effectuated robust biophysical parameter retrievals for a corn field characterized by a wide range in leaf chlorophyll levels and intermixed green and senescent leaf material. Validation against in-situ observations yielded relative root-mean-square deviations (RMSD) on the order of 10% for the 1 m resolution LAI (RMSD = 0.25) and C_ab (RMSD = 4.4 μg cm^− 2) estimates, due in part to an efficient correction for background influences. LAI and C_ab retrieval accuracies at the SPOT 10 m resolution were characterized by relative RMSDs of 13% (0.3) and 17% (7.1 μg cm^− 2), respectively, and the overall intra-field pattern in LAI and C_ab was well established at this resolution. The developed method has utility in agricultural fields characterized by widely varying distributions of model variables and holds promise as a valuable operational tool for precision crop management. Work is currently in progress to extend REGFLEC to regional scales.     

Using LiDAR and quickbird data to model plant production and quantify uncertainties associated with wetland detection and land cover generalizations

Spatiotemporal data from satellite remote sensing and surface meteorology networks have made it possible to continuously monitor global plant production, and to identify global trends associated with land cover/use and climate change. Gross primary production (GPP) and net primary production (NPP) are routinely derived from the Moderate Resolution Imaging Spectroradiometer (MODIS) onboard satellites Terra and Aqua, and estimates generally agree with independent measurements at validation sites across the globe. However, the accuracy of GPP and NPP estimates in some regions may be limited by the quality of model input variables and heterogeneity at fine spatial scales. We developed new methods for deriving model inputs (i.e., land cover, leaf area, and photosynthetically active radiation absorbed by plant canopies) from airborne laser altimetry (LiDAR) and Quickbird multispectral data at resolutions ranging from about 30 m to 1 km. In addition, LiDAR-derived biomass was used as a means for computing carbon-use efficiency. Spatial variables were used with temporal data from ground-based monitoring stations to compute a six-year GPP and NPP time series for a 3600 ha study site in the Great Lakes region of North America. Model results compared favorably with independent observations from a 400 m flux tower and a process-based ecosystem model (BIOME-BGC), but only after removing vapor pressure deficit as a constraint on photosynthesis from the MODIS global algorithm. Fine-resolution inputs captured more of the spatial variability, but estimates were similar to coarse-resolution data when integrated across the entire landscape. Failure to account for wetlands had little impact on landscape-scale estimates, because vegetation structure, composition, and conversion efficiencies were similar to upland plant communities. Plant productivity estimates were noticeably improved using LiDAR-derived variables, while uncertainties associated with land cover generalizations and wetlands in this largely forested landscape were considered less important.     

Effects of soil moisture and water depth on ERS SAR backscatter measurements from an Alaskan wetland complex

We conducted a preliminary investigation of the response of ERS C-band SAR backscatter to variations in soil moisture and surface inundation in wetlands of interior Alaska. Data were collected from 5 wetlands over a three-week period in 2007. Results showed a positive correlation between backscatter and soil moisture in sites dominated by herbaceous vegetation cover (r = 0.74, p < 0.04). ERS SAR backscatter was negatively correlated to water depth in all open (non-forested) wetlands when water table levels were more than 6 cm above the wetland surface (r = − 0.82, p < 0.001). There was no relationship between backscatter and soil moisture in the forested (black spruce-dominated) wetland site. Our preliminary results show that ERS SAR data can be used to monitor variations in hydrologic conditions in high northern latitude wetlands (including peatlands), particularly sites with sparse tree cover.     

Impact of DEM accuracy and resolution on topographic indices

Topography is an important land-surface characteristic that affects most aspects of the water balance in a catchment, including the generation of surface and sub-surface runoff; the flow paths followed by water as it moves down and through hillslopes and the rate of water movement. All of the spatially explicit fully distributed hydraulic and hydrological models use topography (represented by the DEM of the area modelled) to derive bathymetry. DEM is also used to derive some other key information critical in fully distributed hydraulic and hydrological models.With high-resolution DEMs such as LiDAR (Light Detection and Ranging) becoming more readily available and also with the advancements in computing facilities which can handle these large data sets, there is a need to quantify the impact of using different resolution DEMs (e.g. 1 m against 10 m or 25 m) on hydrologically important variables and the loss of accuracy and reliability of the results as we move from high resolution to coarser resolution.The results from statistical analysis carried out to compare field survey elevations with the LiDAR DEM-derived elevations, show that there are small differences between the two data sets but LiDAR DEM is a reasonably good representation of the actual ground surface compared to other commonly used DEMs derived from contour maps.The results from the analysis clearly show that the accuracy and resolution of the input DEM have serious implications on the values of the hydrologically important spatial indices derived from the DEM. The result also indicates that the loss of details by re-sampling the higher resolution DEM to coarser resolution are much less compared to the details captured in the commonly available coarse resolution DEM derived from contour maps. Topographic indices based on contour derived DEMs should be used with caution and where available, the higher resolution DEM should be used instead of the coarse resolution one.     

Evaluation of ASTER and MODIS land surface temperature and emissivity products using long-term surface longwave radiation observations at SURFRAD sites

Land surface temperature (LST) and emissivity are key parameters in estimating the land surface radiation budget, a major controlling factor of global climate and environmental change. In this study, Terra Advanced Spaceborne Thermal Emission Reflection Radiometer (ASTER) and Aqua MODerate resolution Imaging Spectroradiometer (MODIS) Collection 5 LST and emissivity products are evaluated using long-term ground-based longwave radiation observations collected at six Surface Radiation Budget Network (SURFRAD) sites from 2000 to 2007. LSTs at a spatial resolution of 90 m from 197 ASTER images during 2000-2007 are directly compared to ground observations at the six SURFRAD sites. For nighttime data, ASTER LST has an average bias of 0.1 °C and the average bias is 0.3 °C during daytime. Aqua MODIS LST at 1 km resolution during nighttime retrieved from a split-window algorithm is evaluated from 2002 to 2007. MODIS LST has an average bias of − 0.2 °C. LST heterogeneity (defined as the Standard Deviation, STD, of ASTER LSTs in 1 × 1 km² region, 11 × 11 pixel in total) and instrument calibration error of pyrgeometer are key factors impacting the ASTER and MODIS LST evaluation using ground-based radiation measurements. The heterogeneity of nighttime ASTER LST is 1.2 °C, which accounts for 71% of the STD of the comparison, while the heterogeneity of the daytime LST is 2.4 °C, which accounts for 60% of the STD. Collection 5 broadband emissivity is 0.01 larger than that of MODIS Collection 4 products and ASTER emissivity. It is essential to filter out the abnormal low values of ASTER daily emissivity data in summer time before its application.     

Validation of AMSR-E soil moisture using L-band airborne radiometer data from National Airborne Field Experiment 2006

AMSR-E has been extensively evaluated under a wide range of ground and climate conditions using in situ and aircraft data, where the latter were primarily used for assessing the T_B calibration accuracy. However, none of the previous work evaluates AMSR-E performance under the conditions of flood irrigation or other forms of standing water. Also, it should be mentioned that global soil moisture retrievals from AMSR-E typically utilize X-band data. Here, C-band based AMSR-E soil moisture estimates are evaluated using 1 km resolution retrievals derived from L-band aircraft data collected during the National Airborne Field Experiment (NAFE'06) field campaign in November 2006. NAFE'06 was conducted in the Murrumbidgee catchment area in southeastern Australia, which offers diverse ground conditions, including extensive areas with dryland, irrigation, and rice fields. The data allowed us to examine the impact of irrigation and standing water on the accuracy of satellite-derived soil moisture estimates from AMSR-E using passive microwave remote sensing. It was expected that in fields with standing water, the satellite estimates would have a lower accuracy as compared to soil moisture values over the rest of the domain. Results showed sensitivity of the AMSR-E to changes in soil moisture caused by both precipitation and irrigation, as well as good spatial (average R = 0.92 and RMSD = 0.049 m³/m³) and temporal (R = 0.94 and RMSD = 0.04 m³/m³) agreement between the satellite and aircraft soil moisture retrievals; however, under the NAFE'06 ground conditions, the satellite retrievals consistently overestimated the soil moisture conditions compared to the aircraft.     

Evaluating hyperspectral imaging of wetland vegetation as a tool for detecting estuarine nutrient enrichment

Nutrient enrichment and eutrophication are major concerns in many estuarine and wetland ecosystems, and the need is urgent for fast, efficient, and synoptic ways to detect and monitor nutrients in wetlands and other coastal systems across multiple spatial and temporal scales. We integrated three approaches in a multi-disciplinary evaluation of the potential for using hyperspectral imaging as a tool to assess nutrient enrichment and vegetation responses in tidal wetlands. For hyperspectral imaging to be an effective tool, spectral signatures must vary in ways correlated with water nutrient content either directly, or indirectly via such proxies as vegetation responses to elevated nitrogen. Working in Elkhorn Slough, central California, where intensive farming practices generate considerable runoff of fertilizers and pesticides, we looked first for long- and short-term trends among temporally ephemeral point data for nutrients and other water quality characters collected monthly at 18 water sampling stations since 1988. Second, we assessed responses of the dominant wetland plant, Salicornia virginica (common pickleweed) to two fertilizer regimes in 0.25 m² experimental plots, and measured changes in tissue composition (C, H, N), biomass, and spectral responses at leaf and at canopy scales. Third, we used HyMap hyperspectral imagery (126 bands; 15–19 nm spectral resolution; 2.5 m spatial resolution) for a synoptic assessment of the entire wetland ecosystem of Elkhorn Slough. We mapped monospecific Salicornia patches (～ 56–500 m²) on the ground adjacent to the 18 regular water sampling sites, and then located these patches in the hyperspectral imagery to correlate long-term responses of larger patches to water nutrient regimes. These were used as standards for correlating plant canopy spectral responses with nitrogen variation described by the water sampling program. There were consistent positive relationships between nitrogen levels and plant responses in both the field experiment and the landscape analyses. Two spectral indices, the Photochemical Reflectance Index (PRI) and Derivative Chlorophyll Index (DCI), were correlated significantly with water nutrients. We conclude that hyperspectral imagery can be used to detect nutrient enrichment across three spatial and at least two temporal scales, and suggest that more quantitative information could be extracted with further research and a greater understanding of physiological and physical mechanisms linking water chemistry, plant properties and spectral imaging characteristics.     

Modelling of urban sensible heat flux at multiple spatial scales: A demonstration using airborne hyperspectral imagery of Shanghai and a temperature–emissivity separation approach

Urbanization related alterations to the surface energy balance impact urban warming (‘heat islands’), the growth of the boundary layer, and many other biophysical processes. Traditionally, in situ heat flux measures have been used to quantify such processes, but these typically represent only a small local-scale area within the heterogeneous urban environment. For this reason, remote sensing approaches are very attractive for elucidating more spatially representative information. Here we use hyperspectral imagery from a new airborne sensor, the Operative Modular Imaging Spectrometer (OMIS), along with a survey map and meteorological data, to derive the land cover information and surface parameters required to map spatial variations in turbulent sensible heat flux (Q_H). The results from two spatially-explicit flux retrieval methods which use contrasting approaches and, to a large degree, different input data are compared for a central urban area of Shanghai, China: (1) the Local-scale Urban Meteorological Parameterization Scheme (LUMPS) and (2) an Aerodynamic Resistance Method (ARM). Sensible heat fluxes are determined at the full 6 m spatial resolution of the OMIS sensor, and at lower resolutions via pixel aggregation and spatial averaging. At the 6 m spatial resolution, the sensible heat flux of rooftop dominated pixels exceeds that of roads, water and vegetated areas, with values peaking at ～ 350 W m^? 2, whilst the storage heat flux is greatest for road dominated pixels (peaking at around 420 W m^? 2). We investigate the use of both OMIS-derived land surface temperatures made using a Temperature–Emissivity Separation (TES) approach, and land surface temperatures estimated from air temperature measures. Sensible heat flux differences from the two approaches over the entire 2 × 2 km study area are less than 30 W m^? 2, suggesting that methods employing either strategy maybe practica1 when operated using low spatial resolution (e.g. 1 km) data. Due to the differing methodologies, direct comparisons between results obtained with the LUMPS and ARM methods are most sensibly made at reduced spatial scales. At 30 m spatial resolution, both approaches produce similar results, with the smallest difference being less than 15 W m^? 2 in mean Q_H averaged over the entire study area. This is encouraging given the differing architecture and data requirements of the LUMPS and ARM methods. Furthermore, in terms of mean study Q_H, the results obtained by averaging the original 6 m spatial resolution LUMPS-derived Q_H values to 30 and 90 m spatial resolution are within ～ 5 W m^? 2 of those derived from averaging the original surface parameter maps prior to input into LUMPS, suggesting that that use of much lower spatial resolution spaceborne imagery data, for example from Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) is likely to be a practical solution for heat flux determination in urban areas.     

Quantitative mapping of arid alluvial fan surfaces using field spectrometer and hyperspectral remote sensing

Mapping and dating of arid and semi-arid alluvial fans are of great importance in many Quaternary studies. Yet the most common mapping method of these features is based on visual, qualitative interpretation of air-photos. In this study we examine the feasibility of mapping arid alluvial surfaces by using airborne hyperspectral reflective remote sensing methodology. This technique was tested on Late Pleistocene to Holocene alluvial fan surfaces located in the hyperarid southern Arava valley, Israel. Results of spectral field measurements showed that the surface reflectance is controlled by two main surficial processes, which are used as relative age criteria: the degree of desert pavement development (gravel coverage %) controls the absorption feature depths, while the rock coating development influences significantly the overall reflectance of the surface, but its effect on the absorption feature depths is limited. We show that as the percent of the surface covered by gravels increases, the absorption feature depth of the common gravels, in this case carbonate at 2.33 μm, increases as well; whereas the absorption features depth of the fine particle in-between the gravels, decrease (hydroxyl and ferric absorption features at 2.21 μm, and 0.87 μm, respectively), as the fines are removed from the surface. Using these correlations we were able to map the surface gravel coverage (%) on the entire alluvial fan, by calculating the gravel coverage (%) in each pixel of the hyperspectral image. The prediction of gravel coverage (%) is with accuracy of ± 15% (e.g. gravel coverage of 50% can be predicted to be 35% to 65%). Using extensive accuracy assessment data, we show that the spectral based mapping maintained high accuracy degree (R² = 0.57 to 0.83). The quantitative methodology developed in this study for mapping alluvial surfaces can be adapted for other surfaces and piedmonts throughout the arid regions of the world.     

Tree cover and height estimation in the Fennoscandian tundra-taiga transition zone using multiangular MISR data

The tundra-taiga transition zone stretches around the northern hemisphere separating boreal forest to the south from treeless tundra to the north. Tree cover and height are important variables to characterize this vegetation transition. Accurate continuous fields of tree cover and height would enable the delineation of the forest extent according to different criterion and provide useful data for change detection of this climatically sensitive ecotone. This study examined if multiangular remote sensing data has potential to improve the accuracy of the tree cover and height estimates in relation to nadir-view data. The satellite data consisted of Multi-angle Imaging SpectroRadiometer (MISR) data at 275 m and 1.1 km resolutions. The study area was located in the Fennoscandian tundra-taiga transition zone, in northernmost Finland. The continuous fields of tree cover and height were estimated using neural networks, which were trained and assessed by high-resolution biotope inventory data. The spectral-angular data together produced lower estimation errors than single band nadir, multispectral nadir or single band multiangular data alone. RMSE of the tree cover estimates reduced from 7.8% (relative RMSE 67.4%) to 6.5% (56.1%) at 275 m resolution, and from 5.4% (49.2%) to 4.1% (36.9%) at 1.1 km resolution, when multispectral nadir data were used together with multiangular data. RMSE of the tree height estimates reduced from 2.3 m (44.3%) to 2.0 m (37.6%) and from 1.8 m (35.4%) to 1.3 m (25.4%), respectively. The largest estimation errors occurred in mires and in areas of dense shrub cover, but the use of multiangular data also reduced estimation errors in these areas. The results suggest that directional information has potential to improve the tree cover and height estimates, and hence the accuracy of the land cover change detection in the tundra-taiga transition zone.