Article Analysis of the Scaling Effects in the Area-Averaged Fraction of Vegetation Cover Retrieved Using an NDVI-Isoline-Based Linear Mixture Model

The spectral unmixing of a linear mixture model (LMM) with Normalized Difference Vegetation Index (NDVI) constraints was performed to estimate the fraction of vegetation cover (FVC) over the earth’s surface in an effort to facilitate long-term surface vegetation monitoring using a set of environmental satellites. Although the integrated use of multiple sensors improves the spatial and temporal quality of the data sets, area-averaged FVC values obtained using an LMM-based algorithm suffer from systematic biases caused by differences in the spatial resolutions of the sensors, known as scaling effects. The objective of this study is to investigate the scaling effects in area-averaged FVC values using analytical approaches by focusing on the monotonic behavior of the scaling effects as a function of the spatial resolution. The analysis was conducted based on a resolution transformation model introduced recently by the authors in the accompanying paper (Obata et al., 2012). The maximum value of the scaling effects present in FVC values was derived analytically and validated numerically. A series of derivations identified the error bounds (inherent uncertainties) of the averaged FVC values caused by the scaling effect. The results indicate a fundamental difference between the NDVI and the retrieved FVC from NDVI, which should be noted for accuracy improvement of long-term observation datasets.


Introduction
Time series analysis of biophysical and climatological parameters is indispensable for studies of the Earth's climate, land, and hydrology systems and the interactions among these systems.Such analysis often relies on the information retrieved from reflectance spectra measured by satellite instruments [1,2].The fraction of vegetation cover (FVC) is one method used to measure the horizontal spread of vegetation and is defined as the ratio of the vegetated area to the entire pixel area in a satellite image [3].The FVC can be estimated from a reflectance spectrum via empirical or physical model-based techniques, such as the vegetation index (VI) model [4], spectral mixture analysis (linear mixture model, LMM) [3,5,6], and radiative transfer models [7].Several types of FVC datasets have been produced using historical and current optical sensor data (Table 1).For example, the European Commission provided a dataset called fCover that was derived from an ensemble of satellites.The program used to compile this dataset is known as the Carbon Cycle and Change in Land Observation Products (CYCLOPS) with SPOT-VEGETATION [8].The technique used in this work is based on neural networks and uses simulated data from radiative transfer models [9,10] to train the neural network.NOAA initiated a program called the Global Vegetation Processing System (GVPS), which has produced a data product known as the global vegetation fraction (GVF).The dataset covers over 30 years of observation term and used the Normalized Difference Vegetation Index (NDVI) [11] to estimate the FVC.Several other programs have been initiated more recently, including the compilation of the dataset called fCover from ENVISAT-MERIS, which uses a neural network with a radiative transfer model with eleven spectral bands [12]; the FVC products determined using Variable Multiple Endmember Spectral Mixture Analysis (VESMA) of data from a geostational satellite; MSG-SEVIRI [13], and the Vegetation Continuous Field (VCF) from MODIS [14] based on a regression tree algorithm.fCover [8] fCover [12] FVC [13] GVF [11] VCF [ Long-term observations of biophysical quantities spanning several decades generally require integrated use of an ensemble of satellite sensors [15].Such efforts have been widely applied in numerous studies involving land remote sensing applications [16][17][18][19][20][21].Discrepancies between the characteristics of the data sets, such as the spectral band configuration or the spatial resolution, however, introduce systematic biases into the parameter retrievals.These discrepancies can decrease the accuracy of the land observations as well as the predictions based on the datasets; thus, the effects of such discrepancies must be investigated thoroughly to identify mechanisms for rectifying the discrepancies.This study contributes to an improved understanding of the effects of spatial resolution discrepancies on the calculated area-averaged FVC values [22].
Conceptual approaches and measurement models related to the scale issues have been introduced by Strahler et al. [23] and Woodstock and Strahler [24] in the context of remote sensing.Comprehensive and intensive reviews on the scale issues and related problems have been provided by several researchers [23,[25][26][27][28][29][30].Although various types of problems have been investigated by numerous studies, monotonic behavior has not been fully analyzed.This work shed light on the monotonic aspects of the FVC retrieval algorithm within a framework of the modifiable area unit problem [31] originated from the field of spatial analysis, which was later introduced to the remote sensing community [25,26,28,29].
Scaling effects have been investigated extensively in the context of biophysical parameter retrievals, such as calculations of the vegetation index [32][33][34][35] leaf area index (LAI) [36][37][38] and FVC [22].Several studies have attempted to derive an analytical expression for the differences between the area-averaged NDVI values at two extreme resolutions [33,35,39], the coarsest and finest resolution, called the lumped and distributed cases [34].These investigations reported that an NDVI value can never reach its maximum or minimum values at an intermediate resolution between the two extreme cases.Our previous studies [40,41] showed that the area-averaged values of the NDVI reach maximum and minimum values at the two extreme resolutions under certain conditions.A proof was provided based on the monotonic behavior of the scaling effects.One question arose in the case of FVC calculations: Do the area-averaged FVC values change monotonically with the spatial resolution, similar to the case of the NDVI?The present study attempts to answer this question.
Spectral unmixing under NDVI constraints (the NDVI-isoline-based LMM) has been used to estimate FVC values.The technique has the advantages of both the VI and LMM approaches [35,42,43].The retrieved FVC values may be degraded by scaling effects [22] because the NDVI (used in the algorithm) suffers from such effects.The sources of the scaling effects in the FVC and their monotonic behavior were previously investigated analytically [44].Although that study analyzed several fundamental properties of the FVC scaling effects, a geometrical interpretation of the monotonic behavior and a thorough investigation of the key parameters that significantly influenced the error bounds were not performed.The present study investigates the FVC changes as a function of the spatial resolution, which corresponds to the number of pixels within a fixed area, as an extension of our previous study [41].The analysis of the FVC scaling effects is based on the resolution transform model proposed in our previous work [41] that is about the scaling effects in the NDVI.In that work we have shown that the NDVI varies monotonically as a function of spatial resolution in a certain resolution sequence.As an extension of the work, we focus on the scale issues in the FVC that uses NDVI in the retrieval algorithm.Although the FVC retrieval algorithm discussed in this study uses NDVI as a parameter, there is no guarantee that the FVC varies monotonically as the NDVI does.This study addresses this issue.
The structure of this work is as follows.After the brief explanation of theoretical background, we clarify the monotonicity of scaling effects in FVC (Section 3).We then explain a geometric relationship between the monotonic trend and endmember spectra by numerical experiments using simulated endmember spectra in Section 4, followed by a validation exercise in Section 5. Derivation of the maximum error bounds of the scaling effects and numerical simulation will be explained in Section 6.The discussion and conclusion sections (Sections 7 and 8, respectively) follow.

NDVI-Isoline-Based LMM
A reflectance spectrum from a target pixel may be represented as a linear sum of the spectra corresponding to a set of representative surface types (endmember spectra) in a LMM.We assume here that a target field consists of vegetation and non-vegetation classes.Red and near-infrared (NIR) bands are considered.Under these assumptions, a modeled spectrum, ρ ρ ρ = (ρ r , ρn ) (where the subscripts r and n indicate red and NIR spectra, respectively), may be written as where ρ ρ ρ v = (ρ v,r , ρv,n ) and ρ ρ ρ s = (ρ s,r , ρs,n ) are vegetation and non-vegetation endmember spectra in the model (the subscripts v and s indicate the vegetation and soil components, that is, non-vegetation land areas, respectively).The weight of a vegetation endmember, ω, represents the FVC in the LMM.A unity constraint is imposed in Equation (1).To retrieve the FVC, in general, ω may be determined by minimizing a distance measure (e.g., the Euclidean norm) between the modeled spectrum and the measured spectrum [3].The condition that determines the FVC then permits the vegetation index to be substituted for the reflectance spectrum.The choice of NDVI for defining the condition is known as the NDVI-isoline-based LMM [42,45].The estimated value of the FVC may be expressed analytically by solving for ω after equating a modeled and a target NDVI where v is the NDVI from the target spectrum (ρ ρ ρ = (ρ r , ρ n )).The variables vv and vs are the NDVI from the vegetation and non-vegetation endmember spectra in the model, respectively.The variable η is the ratio of the L-1 norm of the vegetation endmember to the non-vegetation endmember.
Under the two-endmember assumption, if an endmember spectra assumed in a model rigorously corresponds to the true endmember spectra over a target field, the spatially averaged FVC will be independent of the spatial resolution (the averaged FVC values must be true values) [44]; however, accurately estimating the endmember spectra based on satellite data is difficult due to the surface heterogeneity or the existence of multiple endmember spectra within a pixel [45][46][47][48].Therefore, the FVC values calculated using this algorithm are biased by the differences between the modeled and true endmember spectra.In this case, the FVC retrieval algorithm becomes a nonlinear function of the input spectra, and this nonlinearity is a source of scaling effects.
The analytical form of the biased FVC under the two-endmember assumption (the target field is modeled under the same assumption) can be derived according to the following steps.First, we assume that a measured spectrum, ρ ρ ρ, can be expressed as a combination of true endmember spectra (ρ ρ ρ v = (ρ v,r , ρ v,n ) and ρ ρ ρ s = (ρ s,r , ρ s,n )) with a true FVC value ω given by We then rewrite the FVC estimate given in Equation (2) using Equation (4) as the target spectrum, where

Scaling Effects of the FVC on the Landsat7-ETM+
Scaling effects in the calculation of an area-averaged FVC can be seen in the data from Landsat7-ETM+.In order to clearly show the scaling effects, we conducted the following experiment using the data acquired on 7 July 2001 (Path:109, Row:36), which covers both urban and suburban regions in the Aichi prefecture in Japan.The spectral data of the size of 64 × 64 pixels (extracted from the original scene) was used in this experiment (Figure 1(a)).The spectral data was aggregated using an average-based moving block window into a data set comprising images at multiple resolutions (2 i × 2 i (i = 0, 1, • • • , 5)), as shown in Figure 2. The vegetation endmember spectra assumed in the algorithm was (0.05, 0.35).The effects of the endmember spectra on the scaling effects were demonstrated by assuming either of two non-vegetation spectra, (0.12, 0.14) and (0.18, 0.22), denoted as the "soil endmem.-1"and "soil endmem.-2",respectively.The reflectance spectra of the target area and endmember spectra assumed in the model are plotted in Figure 1(b).The vegetation spectrum is denoted by the green circle, whereas the soil endmem.-1and -2 spectra are denoted by the dark and light brown circles, respectively.
The averaged FVC values (ω j ) were computed as a function of the spatial resolution, which corresponds to the number of pixels within a fixed area (j).The scaling effects in the averaged FVC, calculated as the relative differences between the ωj and ω1 , were computed according to The relative differences obtained from the pair of spectra comprising the vegetation and soil endmem.-1are indicated by the filled circles in Figure 3.The averaged FVC values as a function of the spatial resolution, for the soil endmem.-2,are indicated by the filled squares (Figure 3).Using the pair of vegetation and soil endmem.-1spectra, the averaged FVC increased as a function of the number of pixels (coarser to finer).By contrast, the averaged FVC decreased when the spectral pair comprising the vegetation and soil endmem.-2spectra was assumed in the algorithm.The results clearly indicate that the magnitude and trends (increasing or decreasing) of the scaling effects depend heavily on the endmember spectra assumed in the model.The scaling effects also depend on the endmember spectra reflected from the target area.The mechanism underlying the scaling effects observed in the calculated FVC has not been previously elucidated.The present study aims to address this very issue.

Resolution Transformation Model
The FVC estimate, ω, presented in Equation ( 5), incurs bias errors if ρ ρ ρ v ̸ = ρ ρ ρ v or ρ ρ ρ s ̸ = ρ ρ ρ s .These discrepancies are the source of the scaling effects.To analyze the scaling effects of ω and the monotonic behavior of the scaling effects as a function of the spatial resolution, we introduced the resolution transformation rule described in our previous work [41].Resolution transformations are conducted according to a simple rule for pixel partitioning that interconnects the area-averaged values of biophysical parameters at different spatial resolutions.A variable α in the model represents the fraction of a pixel remaining after the partitioning process.The definition of the model is described in detail in [41].
The reflectance spectrum for a pixel k at the j-th resolution level can be expressed as ρ ρ ρ j,k .Under the two-endmember assumption, the spectrum can also be expressed as the FVC (true value) for pixel k at the j-th resolution, ω j,k , with endmember spectra for the target field (ρ ρ ρ v and ρ ρ ρ s ) given by The estimated FVC for the k-th pixel at the j-th resolution, ωj,k , can be expressed as The area-averaged FVC values estimated using the algorithm for the j-th resolution level become Note that the areas of all pixels are equal in an above equation.A set of appropriate weights should be used for area averaging processes if the pixel size is not identical [41].The variables assumed for the partitioning and unmixing processes are described in Figure 4.In the figure, a target spectrum derived from data at the 1st resolution level (ρ ρ ρ 1,1 ) is expressed as a function of the endmember spectra corresponding to the vegetation and non-vegetation coverage over a target field (ρ ρ ρ v and ρ ρ ρ s , respectively) and the FVC over the target field (ω 1,1 ).The FVC estimate at the 1st resolution level (ω 1,1 = ω1 ) is obtained by unmixing based on the spectrum ρ ρ ρ 1,1 and the endmember spectra assumed in the model (ρ ρ ρ v and ρ ρ ρ s ).The resolution transformation is performed by partitioning the data at the first resolution level into two sub-pixels.The partition ratios across the pixel pair are represented by α and (1 − α).The target spectra in each pixel (ρ ρ ρ 2,k ) are functions of the FVC for the pixel (ω 2,k ) and the endmember spectra (ρ ρ ρ v and ρ ρ ρ s ).FVC estimates for each pixel at the second resolution level (ω 2,k ) can be obtained by unmixing the target spectra ρ ρ ρ 2,k based on the endmember spectra assumed in the model, ρ ρ ρ v and ρ ρ ρ s .The area-averaged FVC for the second resolution level is a linear sum of ω2,1 and ω2,2 with weights α and (1 − α).The averaged FVC estimate at the j-th resolution level can be expressed in terms of the sum of the average FVC value at the (j − 1)-th resolution level and the factor ∆ω, according to By substituting 2 into j in Equation ( 11), ∆ω can be transformed as follows: Recall that the transformation of the spatial resolution can be modeled by repeatedly applying this "unit" resolution transformation from level 1 to 2. ∆ω in Equation ( 12) thus becomes a key to analyzing the monotonicity as a function of the spatial resolution.For instance, if the sign of ∆ω is invariant under the resolution transformation, the average FVC values certainly vary monotonically.

Monotonicity of the Area-Averaged FVC
The variable ∆ω provides a basis for analyzing the monotonicity of the area-averaged FVC.∆ω is a function of (1) the endmember spectra used in the model; (2) the target field (ρ ρ ρ v , ρ ρ ρ s , ρ ρ ρ v , and ρ ρ ρ s ); and (3) the parameters that vary for each partitioning event (during a resolution transform), such as the true value of the FVC in the data at the second resolution level (ω 2,k ) or the partition ratio (α).We therefore focus on the sign of ∆ω and the dependence of the above variables on ∆ω.

Partial Derivative of ∆ω with Respect to ω 2,1
The sign and dependence of the inputs on ∆ω were analyzed by deriving the partial derivative of ∆ω with respect to the FVC for a 1st pixel in the 2nd resolution level, ω 2,1 , using Equation ( 12), where ω1 is independent of ω 2,1 .The derivative terms on the right-hand side of Equation ( 13) can be obtained by taking the derivative of ω2,k (for k = 1, 2) with respect to ω 2,k , where β, γ, and ϕ are defined by The derivative of ω2,2 with respect to ω 2,1 can be expressed as The variable ω 1,1 can be expressed as Therefore, ω 2,2 can be written as Accordingly, Equation ( 16) yields the following result (after some rearrangements), By setting k = 1 in Equation ( 14) and by substituting Equation (19) into Equation ( 13), the partial derivative of ∆ω with respect to ω 2,1 can be obtained as where

Analysis of ∂∆ω/∂ω 2,1
The behavior of ∆ω can be understood from an analysis of the changes in the sign of ∂∆ω/∂ω 2,1 .For example, if ∂∆ω/∂ω 2,1 changes from negative to positive, the sign of ∆ω should be positive, regardless of the parameters that determine the spatial characteristics, such as the true FVC, ω 2,k , or the pixel partition ratio, α.In this case, the sign depends only on the endmember spectra (ρ ρ ρ v , ρ ρ ρ s , ρ ρ ρ v , and ρ ρ ρ s ).First, the factor β in Equation ( 20) can be rearranged to give the following expression Because the four ratios in the parentheses in Equation ( 22) are equal to the ratio vegetation index (RVI), the two parentheses are inevitably both positive.Thus, Equation ( 22) is positive as well.The sign of ∂∆ω/∂ω 2,1 depends only on the second factor in the right-hand side of Equation (20).The factor may be formulated as h(ω 2,1 , ω 2,2 ), If ω 2,1 = ω 2,2 , the function h(ω 2,1 , ω 2,2 ), ∂∆ω/∂ω 2,1 , and ∆ω are equal to zero, then the FVC is scale-invariant.Likewise, if γ = 0, the above function will be identical to 0, hence, a certain relationship is guaranteed between the modeled and true endmember spectra, represented by the following expression, Equation (24) indicates that a vector spanned by the assumed endmember spectra (∆ρ ρ ρ) necessarily becomes parallel to the vector spanned by the true endmember spectra (∆ρ ρ ρ), as illustrated in Figure 5. Empty and filled squares represent the true vegetation and non-vegetation endmember spectra, respectively.Similarly, the circular symbols indicate the vegetation and non-vegetation endmember spectra, respectively, assumed in the algorithm.
For ω 2,1 ̸ = ω 2,2 and γ ̸ = 0, however, the function h(ω 2,1 , ω 2,2 ) is not equal to zero, and the sign of the function depends on all variables.In this case, Equation ( 23) can be expressed as Figure 5. Illustration of the relationship between the vectors describing the true endmember spectra and the assumed endmember spectra (∆ρ ρ ρ and ∆ρ ρ ρ) in the red-NIR reflectance space.Empty and filled squares indicate the vegetation and non-vegetation endmember spectra in the target field, and empty and filled circles indicate the assumed endmember spectra in the LMM.When the two vectors are parallel, scaling effects are not observed in the average FVC.
If 0 ≤ ζ < 1, the last parenthesis in Equation (25) will be negative, hence the sign of ∂∆ω/∂ω 2,1 will vary from negative to positive as follows. If The sign of ∂∆ω/∂ω 2,1 for ζ < 0, however, does not change monotonically with ω 2,1 because the sign of the last factor in Equation ( 25) cannot be unique (dependent on ω 2,1 ).Finally, if ζ = 1, which corresponds to the case of γ = 0, the average FVC will be independent of the spatial resolution.

Behavior of ∆ω
The sign and dependence of the input variables on ∆ω can be summarized as Consequently, the magnitude of the relationship between ω2 and ω1 can be written as Note that if ζ is equal to or greater than zero, the relationship between ω2 and ω1 is independent of the true value of the fractional abundances of the endmember spectra (ω 2,k ) or on the fractional areas by which one pixel is divided into two (α).On one hand, the relationship depends only on the spectra of the endmember components assumed in the algorithm and the true spectra of the target field (ρ ρ ρ v , ρ ρ ρ s , ρ ρ ρ v , and ρ ρ ρ s ).

Monotonicity of the Area-Averaged FVC
As derived in the previous subsection, the sign of the factor ∆ω is invariant during any resolution transform if the endmember spectra satisfy certain conditions (ζ ≥ 0).In other words, the relationship between the FVC values before and after partitioning should be invariant under each step of the resolution transform.The relationship results in monotonic changes in the average FVC for ζ ≥ 0 as a function of the resolution transform sequence.Recall that the resolution transform sequence generated by the repeated application of a simple partitioning rule is referred to as a "resolution class".The findings from the derivation presented in this section are summarized in the following theorems (under the two-endmember assumption).On the one hand, the average FVC values can change non-monotonically for ζ < 0 because the relationship between the average FVC values at any two resolutions will vary.

Geometrical Relationship between ζ and the Endmember Spectra over the Delta Reflectance Space
The variable ζ may be interpreted geometrically as a factor that determines the monotonic behavior within the scope of monotonicity (ζ ≥ 0).By applying Equation ( 26), we can derive equations that comprehensively describe the relation between ζ and the endmember spectra as follows: where det(∆ρ ρ ρ, ∆ρ ρ ρ) represents the determinant of a matrix consisting of the vectors ∆ρ ρ ρ and ∆ρ ρ ρ.The geometrical relationship is described over a delta-reflectance space in which the x-axis represents the difference between the red reflectance of the vegetation and non-vegetation classes (delta-red reflectance).The y-axis in this space represents the delta-NIR reflectance (Figure 6). Figure 6.A geometrical interpretation of the variable ζ, the endmember spectra in the model, ∆ρ ρ ρ (colored dashed arrows), and the target field, ∆ρ ρ ρ (a solid arrow).If ζ > 1, the determinant, det(∆ρ ρ ρ, ∆ρ ρ ρ), becomes positive, indicating that the vector over the endmember spectra to be estimated, ∆ρ ρ ρ, should be described in regions of the left-hand side of the vector corresponding to the difference between the true vectors, ∆ρ ρ ρ (a solid vector), indicated in blue.That is, ∆ρ ρ ρ may be expressed as a vector rotated counterclockwise relative to ∆ρ ρ ρ.Similarly, if 0 ≤ ζ < 1, the determinant is less than unity and ∆ρ ρ ρ may be expressed as a vector rotated clockwise relative to ∆ρ ρ ρ.

Delta-red reflectance
Delta-NIR reflectance 0 If ζ > 1, then the determinant, det(∆ρ ρ ρ, ∆ρ ρ ρ) becomes positive, indicating that the vector expressed over the basis of the endmember spectra to be estimated, ∆ρ ρ ρ, should be described in regions of the left-hand side of the vector corresponding to the difference between the true vectors, ∆ρ ρ ρ (the solid vector in Figure 6), indicated by a blue dashed vector over the delta-red and delta-NIR reflectance space.In other words, ∆ρ ρ ρ may be expressed as a vector rotated counterclockwise from ∆ρ ρ ρ.Similarly, if ζ < 1, the determinant is less than unity, and ∆ρ ρ ρ may be expressed as a vector rotated clockwise relative to ∆ρ ρ ρ, as illustrated by the red dashed line in Figure 6.Note that these interpretations fall within the limited cases that display monotonicity in the averaged FVC (ζ ≥ 0).
Next, we present a set of numerical examples that demonstrate the relationship between the monotonicity (ζ value) and the endmember spectra (including cases in which non-monotonicity (ζ < 0) holds) over the delta-reflectance space.The factor ζ is explicitly a function of ρ ρ ρ v and ρ ρ ρ s and ∆ρ ρ ρ, given by ρ ρ ρ v and ρ ρ ρ s .In this simulation, two vectors, ∆ρ ρ ρ and ∆ρ ρ ρ, were considered (in Equation ( 32)).∆ρ ρ ρ was held constant, and the vector ∆ρ ρ ρ was varied from 0 to 2π in 0.10 radian increments.The origins of the two key vectors ∆ρ ρ ρ and ∆ρ ρ ρ have been shifted relative to the origin of the coordinate system for comparison purpose.The lengths of these vectors were set to be equal to facilitate interpretation.Note that the lengths do not affect monotonicity trend.An example of the results is shown in Figure 7.The black arrows indicate ∆ρ ρ ρ and the gray arrows indicate ρ ρ ρ v and ρ ρ ρ s .The colored arrows indicate ∆ρ ρ ρ for the respective trends: green arrows correspond to the non-decreasing trends of the averaged FVC (0 ≤ ζ < 1).Red arrows fall into the range of non-increasing (ζ > 1).Blue arrows indicate the non-monotonic case (ζ < 0).In Figure 7, ∆ρ ρ ρ, denoted by the black arrow, is located just on the boundary of two regions (indicated by the red and green arrows).Note that the boundaries between the trends (indicated by the arrows of different colors) are determined by the three vectors, ∆ρ ρ ρ, ρ ρ ρ v , and ρ ρ ρ s , which match the representation of Figure 6. Figure 8 shows the vector ∆ρ ρ ρ = (∆ρ r , ∆ρ n ) in a vector chart for all nine cases of the endmember spectra summarized in Table 2.The results show the influence of the change in ∆ρ ρ ρ on the relation.The results show that the delta-reflectance ranges of the decreasing and increasing trends depended significantly on the reflectances of the non-vegetation endmember.For instance, although no obvious differences are observed among the results shown in Figure 8(a,d,g) for the different vegetation endmember spectra assumed, the contribution of the non-decreasing (green arrows) components increased with the red reflectance in the non-vegetation endmember, as shown in Figure 8(a-c).The blue arrows corresponding to the non-monotonic behavior are situated between ρ ρ ρ v and ρ ρ ρ s , indicated by the gray arrows (Figure 8).Non-monotonicity in the averaged FVC is rarely observed because the angle ∆ρ ρ ρ relative to the x-axis often exceeds that of ρ ρ ρ v ; however, if the non-vegetation endmember reflects extremely lower radiances (for example, in coastal and watershed regions) ∆ρ ρ ρ may fall among the blue vectors, as shown in Figure 8(a,d,g).Even if ∆ρ ρ ρ points in the opposite direction (−∆ρ ρ ρ), the ζ value is invariant under Equation ( 26).This property results in a point reflection of ∆ρ ρ ρ with respect to the monotonicity, as shown in Figure 8.

Validation of the Analytical Implications
In this section, we introduce a set of numerical experiments in an effort to validate the results derived here in a practical application.Equation ( 31) implies an important fact, that the FVC values produced under a resolution transformation do not display variations for ζ = 1.This result implies that the endmember spectra assumed in the algorithm provide FVC values that are less sensitive to scaling effects.In the section below, the numerical experiments are explained to provide evidence for a practical demonstration that demonstrates the utility of our analysis.
The major difficulties associated with this validation lie in the difficulty of estimating ζ using actual satellite data.ζ is a function of the "true" endmember spectra for both vegetation and non-vegetation surfaces, and a rigorous estimation/determination of these endmember spectra is generally impossible.Therefore, we circumvent this difficulty by effectively varying the value of ζ according to the following procedure.
The satellite data presented in Figure 1 were used to simulate the FVC values at various spatial resolution levels (for a total of seven resolution cases).In this experiment, we processed the data using different pairs of assumed endmember spectra.The non-vegetation endmember spectra (ρ ρ ρ s ) were varied as shown in Table 3, whereas the spectrum of the vegetation endmember was held fixed at (0.05, 0.35), the spectrum presented in Figure 1.The variations on ρ ρ ρ s were implemented to effectively introduce variations in ζ.Note that the estimation of ζ is impossible because the "true" endmembers are unknown; however, the variations in ρ ρ ρ s most likely shift the value of ζ, as anticipated.Figure 9 shows the scaling effects in the FVC values as a function of the endmember spectra assumed.The figure clearly indicates that (1) the trend in the FVC values as a function of the spatial resolution changes as the value of ζ changes (ζ depends on ρ ρ ρ s in this experiment); and (2) FVC is less sensitive to ζ for some pairs of endmember spectra, which agrees well with the implications of the results derived above.Table 3. Red and NIR reflectances of the assumed non-vegetation endmember spectra in the validation experiment.Red reflectance 0.12 0.13 0.14 0.15 0.16 0.17 0.18 NIR reflectance 0.14 0.16 0.17 0.18 0.19 0.20 0.22 Figure 9. Scaling effects in the FVC, determined for various endmember spectra corresponding to the non-vegetation surface (Table 3) assumed in the algorithm.0.12 0.13 0.14 0.15 0.16 0.17 0.18 6. Maximum Variations Induced by the Scaling Effects Observed in the Averaged FVC

Derivation of the Maximum Variations
The magnitude of the scaling effects depend on the true FVC value within a target area under a fixed pair of assumed endmember spectra.In this section, we focus on the maximum scaling effects in the FVC as a function of the true FVC.The magnitude of the scaling effects is measured as the difference between the two extreme resolutions (the lumped and distributed cases).This difference can be considered to provide the bounds on the errors resulting from the scaling effects when the assumed endmember spectra meet the conditions for monotonicity, as described in Theorem 1 (ζ ≥ 0).
The focus of this discussion is on the differences between the FVC values at the coarsest and finest resolution.Note that the finest resolution corresponds to the case in which all pixels are composed of only one type of surface.In this case, the spectrum of each pixel may be represented by a single endmember spectrum (either vegetation or non-vegetation).Comparing the two extremes, the scaling effects in the FVC may be defined by where ω(f) is an area-averaged FVC value at the finest resolution (distributed case), with the subscript f indicating the "finest" resolution.ω(f) is defined by where ω(v) and ω(s) represent the FVC estimates for the vegetation and non-vegetation endmember spectra, expressed by with the definitions of In order to derive the maximum of ∆ω 1→f as a function of the true FVC, we consider a partial derivative of ∆ω 1→f with respect to ω, that is, The first term of Equation (39) becomes The second term of Equation (39) becomes ∆ω 1→f reaches a maximum when Solving the system of Equations ( 39), ( 40), (41), and (42) for ω yields where The coefficient C is always positive because all determinants in Equation ( 44) are positive, given that the following four relationships hold; v v > v s , vv > vs , vv > v s , and ωv > ωs .If ζ exceeds unity in Equation ( 43), the two terms are negative, and the negative sign must be selected in Equation (43) by noting the range of FVC values, (0 ≤ ω ≤ 1).Similarly, if ζ is less than unity, the first term in Equation ( 43) exceeds unity.As a result, the negative sign must be selected.In summary, the FVC value (ω max ) that provides the maximum value of ∆ω 1→f is The maximum of ∆ω 1→f , denoted ∆ω 1→f (ω max ), can be derived by substituting Equation (45) into Equation (33), which implicitly includes ω as a parameter.The result is a ratio of quadratic forms in the variable where

Numerical Validation
The results derived to express the maximum variations in FVC in terms of the scaling effects were validated by numerical experiments.The variables ω max and ∆ω 1→f (ω max ) were computed for the set of endmember spectra listed in Table 4.The non-vegetation endmember spectrum was varied along the soil line.The slope and offset of the assumed soil line are 1.2 and 0.0, respectively.The red reflectance was varied from 0.03 to 0.38 at intervals of 0.05 (total of eight cases).Table 4. True endmember spectra for the vegetation (ρ ρ ρ v ) and non-vegetation (ρ ρ ρ s ) surfaces (over a target field), and assumed endmember spectra for vegetation (ρ ρ ρ v ) and non-vegetation (ρ ρ ρ v ) surfaces in the FVC retrieval algorithm during the numerical simulations.
Red reflectance 0.05 0.10 0.05 0.03∼0.38 at 0.05 (8 cases) NIR reflectance 0.40 0.12 0.36 NIR=1.2Red Simulations were conducted according to the following steps.First, the values of ∆ω 1→f (ω) were computed by varying the true FVC from zero to unity for each of the eight cases involving soil endmember spectra assumed in the retrieval algorithms.The scaling effects were then plotted as a function of the true FVC value for the eight cases involving different soil endmembers.The results are indicated by the eight distinct lines in Figure 10.In the second step, the variables ω max and ∆ω 1→f (ω max ) were computed from Equations ( 45) and ( 46) separately based on the results from the previous step.The computed pairs of ω max and ∆ω 1→f (ω max ) (denoted by the empty circles) were then plotted along with the simulated results in the figure.The figure clearly shows the validity of the derived results (ω max and ∆ω 1→f (ω max )).
Figure 10.Numerical demonstration of ω max and ∆ω 1→f (ω max ).The values of ∆ω 1→f were computed over the domain of true FVC values (from zero to unity) using the eight endmember sets listed in Table 4.The non-vegetation endmember spectrum was varied over the red reflectance range from 0.03 to 0.38 in increments of 0.05, and the slope of the soil line was assumed to be 1.2.The red reflectance in the non-vegetation endmember are identified based on the colors and lines.The variables ω max and ∆ω 1→f (ω max ) are plotted over the results and are denoted by empty circles.

Discussion
Previous studies of the scaling effects attempted to derive appropriate error bounds and to develop algorithms for correcting scaling effects [33,39,[49][50][51].In this sense, the motivation and objectives of the present study are somewhat similar to previous studies; however, the framework of the analysis presented here differ substantially from the previous studies, which focused on the uncertainties associated with parameter retrieval at the two extreme resolutions (finest and coarsest) due to variations in the spectral measurements at each point within a pixel.The pixel-scale reflectance was held constant at all resolutions.In their work, the endmember spectral ranges (spectral variations) were parameters that eventually determined the error bounds based on the characteristics of the function (e.g., convex or concave) [39].By contrast, both the spectral and spatial contributions to the scaling effects were investigated in this study.By parameterizing both contributions for a given set of spectral and spatial variables, the monotonicity of the FVC could be determined, leading to a discussion of the error bounds.Note that the resolution transformation model introduced in our previous study enabled us to perform an analysis over the entire set of resolution cases.
This study was conducted as an extension of our previous work, in which we analyzed the scaling effects in the calculation of area-averaged NDVI (the monotonicity of the NDVI) [41].If the averaged FVC shifts monotonically within a given resolution class (ζ equals to or exceeds zero), the FVC values at either extreme resolution (the coarsest and finest resolution levels) are either the maximum or minimum values.The reasoning behind this conclusion is as follows: the extreme resolutions belong to the same resolution class, and any resolution case certainly belongs to at least one of the resolution classes; therefore, the FVC values strictly specify the error bounds on the scaling effects in an FVC estimate caused by changes in the spatial resolution.
In light of the findings associated with the error bounds discussed above, we further derived the maximum difference between FVC values at the two extreme resolution levels as a function of a true FVC value for a fixed set of true and assumed endmember spectra.We derived the expression for the maximum difference as a function of both the "true" endmember spectra and the "assumed" spectra in the algorithm.
The scope of practical applications that lend themselves to the findings of this study is limited due to the difficulties associated with accurately estimating the "true" endmember spectra in a target area.Nevertheless, the estimations of the error bounds in FVC calculations are only one type of application; any technique may be used to measure the ranges over which the true endmember spectra vary across the red-NIR reflectance space.Such approaches can lead to uncertainty in the estimates of the averaged FVC values across spectral data collected at multiple resolutions.One often encounters this type of application in the context of long-term observations of biophysical variables by multiple sensors, where the inter-sensor calibration between sensors of two different resolutions plays an important role.Further studies are needed to explore this possibility.

Conclusions
This work investigated the mechanism underlying the scaling effects in an fraction of vegetation cover (FVC) retrieval algorithm using an NDVI-isoline-based linear mixture model (LMM) in an extension of our previous analysis (which treated the scaling effects on NDVI).The analysis was performed by focusing on the monotonicity of area-averaged FVC calculations as a function of the spatial resolution.
The assumption of a two-endmember LMM facilitated the analytical treatment, which was found to be consistent with our previous investigations of NDVI scaling effects.Interestingly, the monotonic behavior of the FVC was somewhat different from that observed in NDVI calculations, even though the FVC algorithm used NDVI as a condition.The NDVI changes monotonically within a resolution class under the two-endmember LMM, whereas the FVC computed by the NDVI-isoline-based LMM does not necessarily change monotonically.This non-monotonic behavior occurs when the endmember spectra satisfies a certain condition.In other words, the NDVI and FVC may behave differently regarding their monotonic aspect, which is one of the findings of this study.
The condition of monotonicity was determined by the factor ζ, which is a function of the "true" vegetation and non-vegetation endmember spectra over a target field and the endmember spectra "assumed" in the algorithm.Remarkably, this factor was independent of the true value of the FVC over a target field.It means that distribution and fraction of the endmember components do not affect the monotonic behavior.This clearly suggests that the monotonicity was determined only by the two sets of endmember spectra.If ζ ≥ 0, then the average FVC varied monotonically with the spatial resolution within a resolution class (generated by repeated application of a simple partitioning rule).In contrast, the average FVC varied non-monotonically for ζ < 0.
If FVC varies monotonically (ζ ≥ 0), the error bounds on the FVC due to the scaling effects may be determined from the FVC values at the extreme resolution levels for a target field consisting of two endmembers (vegetation and non-vegetation).In contrast, if the FVC varies non-monotonically (ζ < 0), it becomes difficult to determine the error bounds because it is not possible to identify the positions of the maximum and minimum values during the resolution transformation.Further investigations will be needed to address this issue.The maximum scaling effects in an averaged FVC calculation (the difference between the FVC values at the extreme resolution levels), along with the actual FVC values, are then derived.These values were also found to be a function of the endmember spectra assumed in the model and over the target area.
The mechanism underlying the scaling effects observed in FVC calculations was analyzed in terms of the function monotonicity and the error bounds.The treatment developed here provides a theoretical basis for the scaling effects, which is strength of this work.The findings of this study can contribute to development of a scale-invariant algorithm for FVC retrieval under scenarios in which multiple datasets collected at different spatial resolutions are integrated in a single analysis.
Since the analyses have been performed based on two-endmember linear mixture model, the number of endmember spectra might be the major limitations of this work.For instance, increase of the number of endmember spectra would cause differences in the monotonic behavior to some extent.This point still remains unclear from this work, which should be solely investigated in the future.In addition, further validations of the findings with actual satellite data will also be needed as future efforts.

Figure 1 .
Figure 1.(a) A false-color image of a target field used to simulate the scaling effects in an area-averaged FVC calculation; (b) the reflectance spectra of the field and endmember spectra assumed in the model.One vegetation endmember spectrum assumed in the algorithm was (0.05, 0.35), denoted by a green circle.Two non-vegetation spectra were assumed: (0.12, 0.14) and (0.18, 0.22), that is, "soil endmem.-1"and "soil endmem.-2",indicated by the dark and light brown circles, respectively.

Figure 2 .
Figure 2. False-color images of the spectral data at several resolution levels used to simulate the scaling effect in the calculation of an averaged FVC.

Figure 3 .
Figure 3. Scaling effects of area-averaged FVC, calculated based on the NDVI-isoline-based LMM as a function of the spatial resolution, using two pairs of endmember spectra.The averaged FVC derived using the vegetation endmember and the soil endmem.-1spectra is indicated by the filled circles.The averaged FVC derived using the vegetation endmember and soil endmem.-2spectra is indicated by the filled squares.Increasing and decreasing trends in the averaged FVC depend on the choice of endmember spectra assumed in the model.(The trend also depends on the endmember spectra present over the target fields).

Figure 4 .
Figure 4. Illustration of the variables, the partitioning process, and the unmixing process used to implement a FVC estimate using the isoline-based LMM.

Theorem 1 Theorem 2
Within a single resolution class, the area-averaged FVC values change monotonically as a function of the spatial resolution if ζ ≥ 0. When this is the case, The area-averaged FVC values follow a trend determined by ζ such that estimation of the FVC is a non-increasing function if ζ > 1 or a non-decreasing function if 0 ≤ ζ < 1.

Figure 7 .Figure 8 .
Figure 7.An example numerical demonstrations of the relationships between the endmember spectra and the monotonic behavior of the area-averaged FVC.∆ρ ρ ρ was varied from 0 to 2π in 0.10 radian increments.The initial point and the length of the vectors were assumed to be (0, 0) and 0.25, respectively.∆ρ ρ ρ is indicated by black arrows, and the vector defined by the assumed endmember spectra, ∆ρ ρ ρ, is indicated by the colored arrows.Gray arrows indicate ρ ρ ρ v and ρ ρ ρ s .The vectors indicated in green correspond to non-decreasing trends in the averaged FVC (0 ≤ ζ < 1).The red arrows indicate that ζ falls into the class of non-increasing functions (ζ > 1).The blue arrows indicate non-monotonic case (ζ < 0).

Table 1 .
Summary of the major FVC-equivalent products.