Effects of woody elements on simulated canopy reflectance: Implications for forest chlorophyll content retrieval

An important bio-indicator of actual plant health status, the foliar content of chlorophyll a and b (Cab), can be estimated using imaging spectroscopy. For forest canopies, however, the relationship between the spectral response and leaf chemistry is confounded by factors such as background (e.g. understory), canopy structure, and the presence of non-photosynthetic vegetation (NPV, e.g. woody elements)—particularly the appreciable amounts of standing and fallen dead wood found in older forests. We present a sensitivity analysis for the estimation of chlorophyll content in woody coniferous canopies using radiative transfer modeling, and use the modeled top-of-canopy reflectance data to analyze the contribution of woody elements, leaf area index (LAI), and crown cover (CC) to the retrieval of foliar Cab content. The radiative transfer model used comprises two linked submodels: one at leaf level (PROSPECT) and one at canopy level (FLIGHT). This generated bidirectional reflectance data according to the band settings of the Compact High Resolution Imaging Spectrometer (CHRIS) from which chlorophyll indices were calculated. Most of the chlorophyll indices outperformed single wavelengths in predicting Cab content at canopy level, with best results obtained by the Maccioni index ([R₇₈₀ − R₇₁₀] / [R₇₈₀ − R₆₈₀]). We demonstrate the performance of this index with respect to structural information on three distinct coniferous forest types (young, early mature and old-growth stands). The modeling results suggest that the spectral variation due to variation in canopy chlorophyll content is best captured for stands with medium dense canopies. However, the strength of the up-scaled Cab signal weakens with increasing crown NPV scattering elements, especially when crown cover exceeds 30%. LAI exerts the least perturbations. We conclude that the spectral influence of woody elements is an important variable that should be considered in radiative transfer approaches when retrieving foliar pigment estimates in heterogeneous stands, particularly if the stands are partly defoliated or long-lived.     

A LiDAR-derived canopy density model for tree stem and crown mapping in Australian forests

The retrieval of tree and forest structural attributes from Light Detection and Ranging (LiDAR) data has focused largely on utilising canopy height models, but these have proved only partially useful for mapping and attributing stems in complex, multi-layered forests. As a complementary approach, this paper presents a new index, termed the Height-Scaled Crown Openness Index (HSCOI), which provides a quantitative measure of the relative penetration of LiDAR pulses into the canopy. The HSCOI was developed from small footprint discrete return LiDAR data acquired over mixed species woodlands and open forests near Injune, Queensland, Australia, and allowed individual trees to be located (including those in the sub-canopy) and attributed with height using relationships (r² = 0.81, RMSE = 1.85 m, n = 115; 4 outliers removed) established with field data. A threshold contour of the HSCOI surface that encompassed ∼ 90% of LiDAR vegetation returns also facilitated mapping of forest areas, delineation of tree crowns and clusters, and estimation of canopy cover. At a stand level, tree density compared well with field measurements (r² = 0.82, RMSE = 133 stems ha^− 1, n = 30), with the most consistent results observed for stem densities ≤ 700 stems ha^− 1. By combining information extracted from both the HSCOI and the canopy height model, predominant stem height (r² = 0.91, RMSE = 0.77 m, n = 30), crown cover (r² = 0.78, RMSE = 9.25%, n = 30), and Foliage & Branch Projective Cover (FBPC; r² = 0.89, RMSE = 5.49%, n = 30) were estimated to levels sufficient for inventory of woodland and open forest structural types. When the approach was applied to forests in north east Victoria, stem density and crown cover were reliably estimated for forests with a structure similar to those observed in Queensland, but less so for forests of greater height and canopy closure.     

Towards estimation of canopy foliar biomass with spectral reflectance measurements

Canopy foliar biomass, defined as the product of leaf dry matter content and leaf area index, is an important measurement for global biogeochemical cycles. This study explores the potential for retrieving foliar biomass in green canopies using a spectral index, the Normalized Dry Matter Index (NDMI). This narrow-band index is based on absorption at the C-H bond stretch overtone and is correlated with leaf dry matter content in fresh green leaves. PROSPECT and SAIL model simulations suggest that the NDMI at the canopy scale is able to minimize the effects of leaf thickness and leaf water content and to maximize sensitivity to variation in canopy foliar biomass. The simulation outputs were analyzed with an ANOVA, and 87% of the variation in the NDMI is explained by leaf dry matter content. The NDMI was linearly related to foliar biomass (g cm^− 2) from model simulations (R² = 0.97). The NDMI calculated from spectral reflectances for one to four stacked leaves was also correlated with total leaf biomass (R² = 0.59). These results suggest that it may be possible to determine foliar biomass from airborne and satellite-borne imaging spectrometers, such as NASA's HyspIRI mission.     

Continental-scale net radiation and evapotranspiration estimated using MODIS satellite observations

Evapotranspiration (ET) is a major pathway for water loss from many ecosystems, and its seasonal variation affects soil moisture and net ecosystem CO₂ exchange. We developed an algorithm to estimate ET using a semi-empirical Priestley-Taylor (PT) approach, which can be applied at a range of spatial scales. We estimated regional net radiation (R_net) at monthly time scales using MODerate resolution Imaging Spectroradiometer (MODIS) albedo and land surface temperature. Good agreement was found between satellite-based estimates of monthly R_net and field-measured R_net, with a RMSE of less than 30 W m^− 2. An adjustable PT coefficient was parameterized as a function of leaf area index and soil moisture based on observations from 27 AmeriFlux eddy covariance sites. The biome specific optimization using tower-based observations performed well, with a RMSE of 17 W m^− 2 and a correlation of 0.90 for predicted monthly latent heat. We implemented the approach within the hydrology module of the CASA biogeochemical model, and used it to estimate ET at a 1 km spatial resolution for the conterminous United States (CONUS). The RMSE of modeled ET was reduced to 21.1 mm mon^− 1, compared to 27.1 mm mon^− 1 in the original CASA model. The monthly ET rates averaged over the Mississippi River basin were similar to those derived using GRACE satellite measurements and river discharge data. ET varied substantially over the CONUS, with annual mean values of 110 ± 76 mm yr^− 1 in deserts, 391 ± 176 mm yr^− 1 in savannas and grasslands, and 840 ± 234 mm yr^− 1 in broadleaf forests. The PT coefficient was the main driver for the spatial variation of ET in arid areas, whereas R_net controlled ET when mean annual precipitation was higher than approximately 400 mm yr^− 1.     

NDVI responses to the forest canopy and floor from spring to summer observed by airborne spectrometer in eastern Siberia

A popular method of satellite-based monitoring of the photosynthetic potential of vegetation is to calculate the normalised difference vegetation index (NDVI) from measurements of the red (RED) and near-infrared (NIR) bands. Enormous amounts of vegetation information have been obtained over continental to global areas based on NDVI derived from NOAA-AVHRR, Terra/Aqua-MODIS, and SPOT-VEGETATION satellite observations. In eastern Siberia, where sparse boreal forests are dominant, the lack of landscape-scale canopy-reflectance observations impedes interpretation of how NDVI seasonality is controlled by the forest canopy and floor status. We discuss the NDVI of the canopy and floor separately based on airborne spectral reflectance measurements and simultaneous airborne land surface images acquired around Yakutsk, Siberia, using a hedgehopping aircraft from spring to summer 2000. The aerial land surface images (4402 scenes) were visually classified into four types according to the forest condition: no-green canopy and snow floor (Type 1), green canopy and snow floor (Type 2), no-green canopy and no-snow floor (Type 3), and green canopy and no-snow floor (Type 4). The spectral reflectance from 350 to 1200 nm was then calculated for these four types. Type 1 had almost no difference in reflectance between the RED and NIR bands, and the resultant NDVI was slightly negative (− 0.03). Although Type 2 showed a significant difference between the two bands because of canopy greenness, the resultant NDVI was rather low (0.17) because of high reflection from the snow cover on the floor. In Type 3, the significant difference between the two bands was mainly caused by the greenness of the floor, and the NDVI was relatively high (0.45). The NDVI for Type 4 was the highest (0.75) among the four types. The contributions of reflectance from the forest canopy and floor to the total reflectance were tested with a forest radiative transfer model. The reflectance difference between NIR and RED bands (NIR − RED) of Type 4 (15.6%) was approximately double the differences of Type 2 (7.0%) and of Type 3 (7.9%), suggesting half-and-half contributions of forest canopy greenness and floor greenness to the total greenness. The result also suggested that the satellite-derived NDVI in the larch forest around Yakutsk reaches 85% of the maximum NDVI owing to the forest floor greenness, and only the other 15% of the increase in NDVI should be attributed to the canopy foliation. These results quantitatively reveal that the NDVI depends considerably on forest floor greenness and snow cover in addition to canopy greenness in the case of relatively sparse forest in Siberia.     

Response of the Sundarbans coastline to sea level rise and decreased sediment flow: A remote sensing assessment

The Sundarbans is the world's largest remaining single block of mangrove forest, covering approximately 1 million ha (~ 10,000 km²) of the Ganges-Brahmaputra delta along the coastal areas of India and Bangladesh. Sea level rise and alteration of water flows of the Himalayan headwaters are among the major disturbances threatening these coastal areas. But very few studies exist on the dynamics or current status of the Sundarbans coastline. We used Landsat images spanning from 1973 to 2010, and an algorithm that we developed, to consistently estimate the spatiotemporal dynamics of erosion and accretion for four different time intervals and the whole study period. Our results show that the direction and extent of erosion and accretion rates varied throughout the different periods. Erosion was the highest in the 1973-1979 interval, with 23.2 km² year⁻¹ of land loss. However, that rate substantially declined in the following periods, reaching a rate of 7-10 km² year⁻¹. Accretion showed a rate of 10 km² year⁻¹ between 1973 and 1989, but substantially declined to ~ 4 km² year⁻¹ between 1989 and 2010. Accretion rate has declined in the recent years but erosion rate has remained relatively high. As a result the delta front has undergone a net erosion of ~ 170 km² of coastal land in the 37 years of our study period. These numbers are significantly higher than the previously reported rates and magnitudes of erosion in this area. The methods and maps developed in this study may be helpful in management planning of this vulnerable coastline.     

Influence of woody elements of a Norway spruce canopy on nadir reflectance simulated by the DART model at very high spatial resolution

A detailed sensitivity analysis investigating the effect of woody elements introduced into the Discrete Anisotropic Radiative Transfer (DART) model on the nadir bidirectional reflectance factor (BRF) for a simulated Norway spruce canopy was performed at a very high spatial resolution (modelling resolution 0.2 m, output pixel size 0.4 m). We used such a high resolution to be able to parameterize DART in an appropriate way and subsequently to gain detailed understanding of the influence of woody elements contributing to the radiative transfer within heterogeneous canopies. Three scenarios were studied by modelling the Norway spruce canopy as being composed of i) leaves, ii) leaves, trunks and first order branches, and finally iii) leaves, trunks, first order branches and small woody twigs simulated using mixed cells (i.e. cells approximated as composition of leaves and/or twigs turbid medium, and large woody constituents). The simulation of each scenario was performed for 10 different canopy closures (CC = 50-95%, in steps of 5%), 25 leaf area index (LAI = 3.0-15.0 m² m^− 2, in steps of 0.5 m² m^− 2), and in four spectral bands (centred at 559, 671, 727, and 783 nm, with a FWHM of 10 nm). The influence of woody elements was evaluated separately for both, sunlit and shaded parts of the simulated forest canopy, respectively. The DART results were verified by quantifying the simulated nadir BRF of each scenario with measured Airborne Imaging Spectroradiometer (AISA) Eagle data (pixel size of 0.4 m). These imaging spectrometer data were acquired over the same Norway spruce stand that was used to parameterise the DART model.The Norway spruce canopy modelled using the DART model consisted of foliage as well as foliage including robust woody constituents (i.e. trunks and branches). All results showed similar nadir BRF for the simulated wavelengths. The incorporation of small woody parts in DART caused the canopy reflectance to decrease about 4% in the near-infrared (NIR), 2% in the red edge (RE) and less than 1% in the green band. The canopy BRF of the red band increased by about 2%. Subsequently, the sensitivity on accounting for woody elements for two spectral vegetation indices, the normalized difference vegetation index (NDVI) and the angular vegetation index (AVI), was evaluated. Finally, we conclude on the importance of including woody elements in radiative transfer based approaches and discuss the applicability of the vegetation indices as well as the physically based inversion approaches to retrieve the forest canopy LAI at very high spatial resolution.     

Comparison of MODIS gross primary production estimates for forests across the U.S.A. with those generated by a simple process model, 3-PGS

Since 2000, NASA's Moderate Resolution Imaging Spectro-radiometer (MODIS) has provided 1 × 1 km estimates of 8-day gross primary production (GPP). The MODIS algorithm computes GPP as a simple function of absorbed photosynthetically active radiation and a regionally assigned light-use conversion efficiency (LUE) that is reduced if temperature or atmospheric vapor pressure deficits are suboptimal. We compared MODIS-derived GPP estimates for forested areas across the United States of America (U.S.A.) with those generated by the 3-PGS (Physiological Principles Predicting Growth using Satellite data) model, the latter of which considers spatial variation in available soil water storage capacity (ASWC) and nitrogen content. We expected seasonal and annual MODIS GPP values to be in close agreement with those derived from the 3-PGS model in regions with adequate precipitation, soil water storage, and moderately fertile soils. 3-PGS was initially run with STATSGO-derived soils information provided by the Oak Ridge National Laboratory. The analysis was expanded to include sensitivity analyses with ASWC set at 50, 100, 300, and 400 mm to identify areas within nine major ecoregions where drought might prove to be a major limitation on GPP. The majority of forests across the U.S.A. were relatively insensitive to large variations in ASW storage. In areas where ASWC was assumed < 200 mm and average annual rainfall was < 100 mm yr^− 1, GPP was predicted to be reduced by > 60%. There was generally good agreement (within 20%) between MODIS and 3-PGS estimates of forest GPP across the U.S.A. GPP predicted by the MODIS model was generally higher in ecoregions with substantial drought and with relatively low soil fertility. The latter, which influences LUE, was more than twice as important as soil drought.     

Application of hyperspectral vegetation indices to detect variations in high leaf area index temperate shrub thicket canopies

Accurate measurement of leaf area index (LAI), an important characteristic of plant canopies directly linked to primary production, is essential for monitoring changes in ecosystem C stocks and other ecosystem level fluxes. Direct measurement of LAI is labor intensive, impractical at large scales and does not capture seasonal or annual variations in canopy biomass. The need to monitor canopy related fluxes across landscapes makes remote sensing an attractive technique for estimating LAI. Many vegetation indices, such as Normalized Difference Vegetation Index (NDVI), tend to saturate at LAI levels > 4 although tropical and temperate forested ecosystems often exceed that threshold. Using two monospecific shrub thickets as model systems, we evaluated the potential of a variety of algorithms specifically developed to improve accuracy of LAI estimates in canopies where LAI exceeds saturation levels for other indices. We also tested the potential of indices developed to detect variations in canopy chlorophyll to estimate LAI because of the direct relationship between total canopy chlorophyll content and LAI. Indices were evaluated based on data from direct (litterfall) and indirect measurements (LAI-2000) of LAI. Relationships between results of direct and indirect ground-sampling techniques were also evaluated. For these two canopies, the indices that showed the highest potential to accurately differentiate LAI values > 4 were derivative indices based on red-edge spectral reflectance. Algorithms intended to improve accuracy at high LAI values in agricultural systems were insensitive when LAI exceeded 4 and offered little or no improvement over NDVI. Furthermore, indirect ground-sampling techniques often used to evaluate the potential of vegetation indices also saturate when LAI exceeds 4. Comparisons between hyperspectral vegetation indices and a saturated LAI value from indirect measurement may overestimate accuracy and sensitivity of some vegetation indices in high LAI communities. We recommend verification of indirect measurements of LAI with direct destructive sampling or litterfall collection, particularly in canopies with high LAI.     

Impact of footprint diameter and off-nadir pointing on the precision of canopy height estimates from spaceborne lidar

A spaceborne lidar mission could serve multiple scientific purposes including remote sensing of ecosystem structure, carbon storage, terrestrial topography and ice sheet monitoring. The measurement requirements of these different goals will require compromises in sensor design. Footprint diameters that would be larger than optimal for vegetation studies have been proposed. Some spaceborne lidar mission designs include the possibility that a lidar sensor would share a platform with another sensor, which might require off-nadir pointing at angles of up to 16°. To resolve multiple mission goals and sensor requirements, detailed knowledge of the sensitivity of sensor performance to these aspects of mission design is required.This research used a radiative transfer model to investigate the sensitivity of forest height estimates to footprint diameter, off-nadir pointing and their interaction over a range of forest canopy properties. An individual-based forest model was used to simulate stands of mixed conifer forest in the Tahoe National Forest (Northern California, USA) and stands of deciduous forests in the Bartlett Experimental Forest (New Hampshire, USA). Waveforms were simulated for stands generated by a forest succession model using footprint diameters of 20 m to 70 m. Off-nadir angles of 0 to 16° were considered for a 25 m diameter footprint diameter.Footprint diameters in the range of 25 m to 30 m were optimal for estimates of maximum forest height (R² of 0.95 and RMSE of 3 m). As expected, the contribution of vegetation height to the vertical extent of the waveform decreased with larger footprints, while the contribution of terrain slope increased. Precision of estimates decreased with an increasing off-nadir pointing angle, but off-nadir pointing had less impact on height estimates in deciduous forests than in coniferous forests. When pointing off-nadir, the decrease in precision was dependent on local incidence angle (the angle between the off-nadir beam and a line normal to the terrain surface) which is dependent on the off-nadir pointing angle, terrain slope, and the difference between the laser pointing azimuth and terrain aspect; the effect was larger when the sensor was aligned with the terrain azimuth but when aspect and azimuth are opposed, there was virtually no effect on R² or RMSE. A second effect of off-nadir pointing is that the laser beam will intersect individual crowns and the canopy as a whole from a different angle which had a distinct effect on the precision of lidar estimates of height, decreasing R² and increasing RMSE, although the effect was most pronounced for coniferous crowns.     

Direct electrochemistry of cytochrome P450 6A1 in mimic bio-membrane and its application for pesticides sensing

The direct electrochemistry of house fly cytochrome P4506A1 (CYP6A1) confined in dioctadecyl dimethyl ammonium bromide (DDAB) film was achieved. The immobilized CYP6A1 displayed a pair of redox peaks with a formal potential of −0.36 mV in pH 7.0 O₂-free phosphate buffers at scan rate of 1 V s⁻¹ and the direct electron transfer of CYP6A1 was characterized by voltammetry. The CYP6A1 in the DDAB film retained its bioactivity and could catalyze the reduction of dissolved oxygen. Upon addition of its substrate aldrin or heptachlor to the air-saturated solution, the reduction peak current of dissolved oxygen increased, which indicates the catalytic behavior of CYP6A1 to its substrates. By amperometry a calibration linear range was obtained to be 9.08 × 10⁻⁶-4.54 × 10⁻⁵ mol L⁻¹ with a sensitivity of 80 μA mM⁻¹ for aldrin or 8.91 × 10⁻⁶-4.46 × 10⁻⁵ mol L⁻¹ with a sensitivity of 66 μA mM⁻¹ for heptachlor. The apparent Michaelis-Menten constant for the electrocatalytic activity of CYP6A1 was found to be 7.468 × 10⁻⁵ mol L⁻¹ for aldrin and 4.316 × 10⁻⁵ mol L⁻¹ for heptachlor. The bioelectrocatalytic products were analysed using gas chromatography (GC) and electron ionization-mass spectrometry (EI-MS). The results confirmed that epoxidation was the main pathways of CYP6A1-mediated organochlorine pesticides oxidation.     

Estimation of net primary productivity by integrating remote sensing data with an ecosystem model

This paper describes a method for integrating leaf area index (LAI) derived from remote sensing data with an ecosystem model for accurate estimation of net primary productivity (NPP). The ecosystem model used in this study was Sim-CYCLE, with which LAI retrieved from the data acquired by MODIS sensor (MODIS-LAI) was integrated. Global annual NPP was estimated as 59.6 Gt C year⁻¹ by MOD-Sim-CYCLE (Sim-CYCLE after integration of MODIS-LAI), whereas it was 62.7 Gt C year⁻¹ in case of Sim-CYCLE for the year 2001. Both models predicted highest NPP around the equator with another smaller peak occurring around 60°N. These two regions represented the tropical and boreal forests biomes, respectively. The NPP estimated by MOD-Sim-CYCLE exceeded the NPP estimated by Sim-CYCLE in these two regions. Other than the tropical and boreal forests biomes, NPP values estimated by the MOD-Sim-CYCLE were typically lower than Sim-CYCLE across the latitudes. Validations of results in Australia and USA showed that MOD-Sim-CYCLE estimated NPP more accurately than Sim-CYCLE. Our results demonstrate the utility of combining satellite-observation with an ecosystem process model to achieve improved accuracy in estimates and monitoring global net primary productivity.     

The m-pancycle-connectivity of a WK-Recursive network

In this paper, we study the m-pancycle-connectivity of a WK-Recursive network. We show that a WK-Recursive network with amplitude W and level L is strictly (5 × 2^L−1 − 2)-pancycle-connected for W ? 3. That is, each pair of vertices in a WK-recursive network with amplitude greater than or equal to 3 resides in a common cycle of every length ranging from 5 × 2^L−1 − 2 to N, where N is the size of the interconnection network; and the value 5 × 2^L−1 − 2 reaches the lower bound of the problem.     

Predicting temperate conifer forest successional stage distributions with multitemporal Landsat Thematic Mapper imagery

Forest succession is a fundamental ecological process which can impact the functioning of many terrestrial processes, such as water and nutrient cycling and carbon sequestration. Therefore, knowing the distribution of forest successional stages over a landscape facilitates a greater understanding of terrestrial ecosystems. One way of characterizing forest succession over the landscape is to use satellite imagery to map forest successional stages continuously over a region. In this study we use a forest succession model (ZELIG) and a canopy reflectance model (GORT) to produce spectral trajectories of forest succession from young to old-growth stages, and compared the simulated trajectories with those constructed from Landsat Thematic Mapper (TM) imagery to understand the potential of mapping forest successional stages with remote sensing. The simulated successional trajectories captured the major characteristics of observed regional mean succession trajectory with Landsat TM imagery for Tasseled Cap indices based on age information from the Pacific Northwest Forest Inventory and Analysis Integrated Database produced by Pacific Northwest Research Station, USDA Forest Service. Though the successional trajectories are highly nonlinear in the early years of succession, a linear model fits well the regional mean successional trajectories for brightness and greenness due to significant cross-site variations that masked the nonlinearity over a regional scale (R² = 0.8951 for regional mean brightness with age; R² = 0.9348 for regional mean greenness with age). Regression analysis found that Tasseled Cap brightness and greenness are much better predictors of forest successional stages than wetness index based on the data analyzed in this study. The spectral history based on multitemporal Landsat imagery can be used to effectively identify mature and old-growth stands whose ages do not match with remote sensing signals due to change occurred during the time between ground data collection and image acquisition. Multitemporal Landsat imagery also improves prediction of forest successional stages. However, a linear model on a stand basis has a limited predictive power of forest stand successional stages (adjusted R² = 0.5435 using the Tasseled Cap indices from all four images used in this study) due to significant variations in remote sensing signals for stands at the same successional stage. Therefore, accurate prediction of forest successional stage using remote sensing imagery at stand scale requires accounting for site-specific factors influence remotely sensed signals in the future.     

Leaf chlorophyll content retrieval from airborne hyperspectral remote sensing imagery

Hyperspectral remote sensing has great potential for accurate retrieval of forest biochemical parameters. In this paper, a hyperspectral remote sensing algorithm is developed to retrieve total leaf chlorophyll content for both open spruce and closed forests, and tested for open forest canopies. Ten black spruce (Picea mariana (Mill.)) stands near Sudbury, Ontario, Canada, were selected as study sites, where extensive field and laboratory measurements were carried out to collect forest structural parameters, needle and forest background optical properties, and needle biophysical parameters and biochemical contents chlorophyll a and b. Airborne hyperspectral remote sensing imagery was acquired, within one week of ground measurements, by the Compact Airborne Spectrographic Imager (CASI) in a hyperspectral mode, with 72 bands and half bandwidth 4.25-4.36 nm in the visible and near-infrared region and a 2 m spatial resolution. The geometrical-optical model 4-Scale and the modified leaf optical model PROSPECT were combined to estimate leaf chlorophyll content from the CASI imagery. Forest canopy reflectance was first estimated with the measured leaf reflectance and transmittance spectra, forest background reflectance, CASI acquisition parameters, and a set of stand parameters as inputs to 4-Scale. The estimated canopy reflectance agrees well with the CASI measured reflectance in the chlorophyll absorption sensitive regions, with discrepancies of 0.06%-1.07% and 0.36%-1.63%, respectively, in the average reflectances of the red and red-edge region. A look-up-table approach was developed to provide the probabilities of viewing the sunlit foliage and background, and to determine a spectral multiple scattering factor as functions of leaf area index, view zenith angle, and solar zenith angle. With the look-up tables, the 4-Scale model was inverted to estimate leaf reflectance spectra from hyperspectral remote sensing imagery. Good agreements were obtained between the inverted and measured leaf reflectance spectra across the visible and near-infrared region, with R² = 0.89 to R² = 0.97 and discrepancies of 0.02%-3.63% and 0.24%-7.88% in the average red and red-edge reflectances, respectively. Leaf chlorophyll content was estimated from the retrieved leaf reflectance spectra using the modified PROSPECT inversion model, with R² = 0.47, RMSE = 4.34 μg/cm², and jackknifed RMSE of 5.69 μg/cm² for needle chlorophyll content ranging from 24.9 μg/cm² to 37.6 μg/cm². The estimates were also assessed at leaf and canopy scales using chlorophyll spectral indices TCARI/OSAVI and MTCI. An empirical relationship of simple ratio derived from the CASI imagery to the ground-measured leaf area index was developed (R² = 0.88) to map leaf area index. Canopy chlorophyll content per unit ground surface area was then estimated, based on the spatial distributions of leaf chlorophyll content per unit leaf area and the leaf area index.     

Remote estimation of chl-a concentration in turbid productive waters — Return to a simple two-band NIR-red model?

Today the water quality of many inland and coastal waters is compromised by cultural eutrophication in consequence of increased human agricultural and industrial activities. Remote sensing is widely applied to monitor the trophic state of these waters. This study investigates the performance of near infrared-red models for the remote estimation of chlorophyll-a concentrations in turbid productive waters and evaluates several near infrared-red models developed within the last 34 years. Three models were calibrated for a dataset with chlorophyll-a concentrations from 0 to 100 mg m⁻³ and validated for independent and statistically different datasets with chlorophyll-a concentrations from 0 to 100 mg m⁻³ and 0 to 25 mg m⁻³ for the spectral bands of the MEdium Resolution Imaging Spectrometer (MERIS) and MODerate resolution Imaging Spectroradiometer (MODIS). The MERIS two-band model estimated chlorophyll-a concentrations slightly more accurately than the more complex models, with mean absolute errors of 2.3 mg m⁻³ for chlorophyll-a concentrations from 0 to 100 mg m⁻³ and 1.2 mg m⁻³ for chlorophyll-a concentrations from 0 to 25 mg m⁻³. Comparable results from several near infrared-red models with different levels of complexity, calibrated for inland and coastal waters around the world, indicate a high potential for the development of a simple universally applicable near infrared-red algorithm.     

The effects of aggregated land cover data on estimating NPP in northern Wisconsin

Ecosystem models are routinely used to estimate net primary production (NPP) from the stand to global scales. Complex ecosystem models, implemented at small scales (< 10 km²), are impractical at global scales and, therefore, require simplifying logic based on key ecological first principles and model drivers derived from remotely sensed data. There is a need for an improved understanding of the factors that influence the variability of NPP model estimates at different scales so we can improve the accuracy of NPP estimates at the global scale. The objective of this study was to examine the effects of using leaf area index (LAI) and three different aggregated land cover classification products-two factors derived from remotely sensed data and strongly affect NPP estimates-in a light use efficiency (LUE) model to estimate NPP in a heterogeneous temperate forest landscape in northern Wisconsin, USA. Three separate land cover classifications were derived from three different remote sensors with spatial resolutions of 15, 30, and 1000 m. Average modeled net primary production (NPP) ranged from 402 gC m^− 2 year^− 1 (15 m data) to 431 gC m^− 2 year^− 1 (1000 m data), for a maximum difference of 7%. Almost 50% of the difference was attributed each to LAI estimates and land cover classifications between the fine and coarse scale NPP estimate. Results from this study suggest that ecosystem models that use biome-level land cover classifications with associated LUE coefficients may be used to model NPP in heterogeneous land cover areas dominated by cover types with similar NPP. However, more research is needed to examine scaling errors in other heterogeneous areas and NPP errors associated with deriving LAI estimates.     

Evapotranspiration on western U.S. rivers estimated using the Enhanced Vegetation Index from MODIS and data from eddy covariance and Bowen ratio flux towers

We combined remote sensing and in-situ measurements to estimate evapotranspiration (ET) from riparian vegetation over large reaches of western U.S. rivers and ET by individual plant types. ET measured from nine flux towers (eddy covariance and Bowen ratio) established in plant communities dominated by five major plant types on the Middle Rio Grande, Upper San Pedro River, and Lower Colorado River was strongly correlated with Enhanced Vegetation Index (EVI) values from the Moderate Resolution Imaging Spectrometer (MODIS) sensor on the NASA Terra satellite. The inclusion of maximum daily air temperatures (T_a) measured at the tower sites further improved this relationship. Sixteen-day composite values of EVI and T_a were combined to predict ET across species and tower sites (r² = 0.74); the regression equation was used to scale ET for 2000-2004 over large river reaches with T_a from meteorological stations. Measured and estimated ET values for these river segments were moderate when compared to historical, and often indirect, estimates and ranged from 851-874 mm yr^− 1. ET of individual plant communities ranged more widely. Cottonwood (Populus spp.) and willow (Salix spp.) stands generally had the highest annual ET rates (1100-1300 mm yr^− 1), while mesquite (Prosopis velutina) (400-1100 mm yr^− 1) and saltcedar (Tamarix ramosissima) (300-1300 mm yr^− 1) were intermediate, and giant sacaton (Sporobolus wrightii) (500-800 mm yr^− 1) and arrowweed (Pluchea sericea) (300-700 mm yr^− 1) were the lowest. ET rates estimated from the flux towers and by remote sensing in this study were much lower than values estimated for riparian water budgets using crop coefficient methods for the Middle Rio Grande and Lower Colorado River.     

Red tide detection and tracing using MODIS fluorescence data: A regional example in SW Florida coastal waters

Near real-time data from the MODIS satellite sensor was used to detect and trace a harmful algal bloom (HAB), or red tide, in SW Florida coastal waters from October to December 2004. MODIS fluorescence line height (FLH in W m^− 2 μm^− 1 sr^− 1) data showed the highest correlation with near-concurrent in situ chlorophyll-a concentration (Chl in mg m^− 3). For Chl ranging between 0.4 to 4 mg m^− 3 the ratio between MODIS FLH and in situ Chl is about 0.1 W m^− 2 μm^− 1 sr^− 1 per mg m^− 3 chlorophyll (Chl = 1.255 (FLH × 10)^0.86, r = 0.92, n = 77). In contrast, the band-ratio chlorophyll product of either MODIS or SeaWiFS in this complex coastal environment provided false information. Errors in the satellite Chl data can be both negative and positive (3-15 times higher than in situ Chl) and these data are often inconsistent either spatially or temporally, due to interferences of other water constituents. The red tide that formed from November to December 2004 off SW Florida was revealed by MODIS FLH imagery, and was confirmed by field sampling to contain medium (10⁴ to 10⁵ cells L^− 1) to high (> 10⁵ cells L^− 1) concentrations of the toxic dinoflagellate Karenia brevis. The FLH imagery also showed that the bloom started in mid-October south of Charlotte Harbor, and that it developed and moved to the south and southwest in the subsequent weeks. Despite some artifacts in the data and uncertainty caused by factors such as unknown fluorescence efficiency, our results show that the MODIS FLH data provide an unprecedented tool for research and managers to study and monitor algal blooms in coastal environments.