Spectral Mixture Analysis as a Unified Framework for the Remote Sensing of Evapotranspiration

• Main Authors: Daniel Sousa, Christopher Small
• Format: Article
• Language: English
• Published: MDPI AG 2018-12-01
• Series: Remote Sensing
• Subjects: spectral mixture analysis; evapotranspiration; surface energy balance
• Online Access: https://www.mdpi.com/2072-4292/10/12/1961

This study illustrates a unified, physically-based framework for mapping landscape parameters of evapotranspiration (<i>ET</i>) using spectral mixture analysis (SMA). The framework integrates two widely used approaches by relating radiometric surface temperature to subpixel fractions of substrate (<i>S</i>), vegetation (<i>V</i>), and dark (<i>D</i>) spectral endmembers (EMs). Spatial and temporal variations in these spectral endmember fractions reflect process-driven variations in soil moisture, vegetation phenology, and illumination. Using all available Landsat 8 scenes from the peak growing season in the agriculturally diverse Sacramento Valley of northern California, we characterize the spatiotemporal relationships between each of the <i>S</i>, <i>V</i>, <i>D</i> land cover fractions and apparent brightness temperature (<i>T</i>) using bivariate distributions in the <i>ET</i> parameter spaces. The dark fraction scales inversely with shortwave broadband albedo (&#961; &lt; &#8722;0.98), and show a multilinear relationship to <i>T</i>. Substrate fraction estimates show a consistent (&#961; &#8776; 0.7 to 0.9) linear relationship to <i>T</i>. The vegetation fraction showed the expected triangular relationship to <i>T</i>. However, the bivariate distribution of <i>V</i> and <i>T</i> shows more distinct clustering than the distributions of Normalized Difference Vegetation Index (<i>NDVI</i>)-based proxies and <i>T</i>. Following the Triangle Method, the <i>V</i> fraction is used with <i>T</i> to compute the spatial maps of the <i>ET</i> fraction (<i>EF</i>; the ratio of the actual total <i>ET</i> to the net radiation) and moisture availability (<i>Mo</i>; the ratio of the actual soil surface evaporation to potential <i>ET</i> at the soil surface). <i>EF</i> and <i>Mo</i> estimates derived from the <i>V</i> fraction distinguish among rice growth stages, and between rice and non-rice agriculture, more clearly than those derived from transformed <i>NDVI</i> proxies. Met station-based reference <i>ET</i> &amp; soil temperatures also track vegetation fraction-based estimates of <i>EF</i> &amp; <i>Mo</i> more closely than do <i>NDVI</i>-based estimates of <i>EF</i> &amp; <i>Mo</i>. The proposed approach using <i>S</i>, <i>V</i>, <i>D</i> land cover fractions in conjunction with T (SVD+T) provides a physically-based conceptual framework that unifies two widely-used approaches by simultaneously mapping the effects of albedo and vegetation abundance on the surface temperature field. The additional information provided by the third (Substrate) fraction suggests a potential avenue for <i>ET</i> model improvement by providing an explicit observational constraint on the exposed soil fraction and its moisture-modulated brightness. The structures of the <i>T</i>, <i>EF</i> &amp; <i>Mo</i> vs SVD feature spaces are complementary and that can be interpreted in the context of physical variables that scale linearly and that can be represented directly in process models. Using the structure of the feature spaces to represent the spatiotemporal trajectory of crop phenology is possible in agricultural settings, because variations in the timing of planting and irrigation result in continuous trajectories in the physical parameter spaces that are represented by the feature spaces. The linear scaling properties of the SMA fraction estimates from meter to kilometer scales also facilitate the vicarious validation of <i>ET</i> estimates using multiple resolutions of imagery.