>>The analysis of effects of varying sun / sensor geometry has been done over 15 days (of which 3 have
been removed) during the peak of the growing season. This misses the highest zenith angles and times
of different vegetation conditions. I suggest to repeat the analysis for other time periods as well to gain
a full picture of sun / sensor geometry effects. Furthermore, why have only NDSIs been investigated
and not the reflectances themselves? This information would help to understand the behaviour of the
NDSIs and would support the claim in the discussion that NDSIs reduce angular effects.
>Response: The reason for not doing the analysis of the varying sun/sensor conditions at the point
in time with the highest zenith angles, is that this occurs during the dry season (two months prior
to the onset of the growing season) where there are no vegetation (herbaceous) influencing the
reflectance spectrum in the measured area. The focus of the manuscript is to investigate how
NDSI is coupled with vegetation parameters, and we hence choose to use the point in time with
most vegetation on the ground.
We agree that it would make a very interesting study to investigate how sun/sensor geometry
influences NDSI differently across the year. However, this is not a minor task and this
manuscript is long as is. We therefore feel that this is beyond the scope of this manuscript. But it
is a very good idea for a future manuscript to investigate seasonal dynamics in anisotropy of both
the reflectance spectrum on its own and on NDSI estimates. This is something that will hopefully
be possible to do in a not too distant future.
The reason for focusing on NDSI, and not on the anisotropy on the reflectance values
themselves is that it has already been done (Huber et al., 2014; Tagesson et al., 2015). The focus
of the paper by Tagesson et al. (2015) is to present all research activities at the Dahra field site.
Among them, a section of the anisotropy of the reflectance spectrum is presented. The aim of the
paper by Huber et al. (2014) is to present the ASD set-up and investigate the quality of the
measurements. A second aim is to study the effects of varying sun/sensor geometry on the
reflectance spectrum. Therefore, in order not to present the same information two times, the
effects of varying sun/sensor geometry part of this paper focus on the effects on the NDSI.
However, the comment is relevant and in the revised manuscript we have included a discussion
regarding the behaviour of the NDSI in relation to the behaviour of the reflectance spectrum and
referred to figures in Huber et al. (2014) and in Tagesson et al. (2015).
Comment: The study uses data from 15
July 2011 until 31 December 2012
. Relationships between ecosystem variables and spectral indices have been investigated for the whole time period. This means phenology from no living vegetation to max living vegetation is included. Therefore it is NOT APPROPRIATE to restrict the analysis of effects of varying sun / sensor geometry to 15 (-3) days during the peak of the growing season.
>> Why has the analysis of the relationship between reflectance / NDSI and ecosystem variables been
restricted to a linear relationship? E.g. other studies found a non-linear relationship between reflectance
and biomass due to saturation effects. Also why have only daily median reflectances / NDSIs been used
when GPP, LUE and FAPAR were daily integrals? Averages would be more appropriate in these cases.
And why have the off-nadir views not been analysed?
> Response: In case the linear relationship is strong, it indicates limited issues with saturation. For
wavelength regions where there are issues with saturations, exponential and logarithmic
regressions could fit better. However, in case the aim is to find wavelength regions which are as
sensitive as possible for investigating seasonal dynamics in an ecosystem property, wavelength
regions with saturtion issues should be avoided. Therefore linear models are better to use than
non-linear models. This was the main reason for fitting linear rather than non-linear regressions.
There is also a practical aspect to it, fitting the reduced major axis linear relationships using the
bootstraping methodology required a full month of processing for these 4 variables (GPP, LUE,
FAPAR and biomass). In case we would try several other regression models, these would require
several months of processing.
Median values were used in order to minimise the influence of errors in the analysis. Median
provides the most common model output and it is thereby more robust against outliers than
average values. This info was provided in the manuscript, but it was not mentioned the first time
that median values were used. Thank you for pointing this out to us, it has been corrected in the
revised manuscript.
We have investigated the seasonal dynamics in the off-nadir views as well, but as seen in the
figure below, there was no difference in seasonal dynamics for the different viewing angles. We
thereby choose to only use the nadir one, as it would not make any difference in the analysis.
Comment: This is not a satisfactory response. Many studies have found non-linear relationships to work much better than linear ones. Therefore, restricting the analysis to linear relationships might not yield the best wavelength combinations.
>> page 3330, lines 11-14: This is not the reason for the saturation of the NDVI. The NDVI saturates at
high biomass because the NIR reflectance is much larger than the red reflectance. NDVI therefore
reduces to R_NIR / R_NIR which equals 1.
> Response: We agree with you, and we are talking about the same thing, we are just using
different phrasing, where you consider it from an equation point of view, we consider it from a
leaf optical property point of view.
All vegetation indices using red will suffer from saturation problems. The reason for this is
related to the fact that there are only so many photons striking a plant leaf and at a certain point,
the chlorophyll absorbs nearly all the red energy to the point where no matter how much
vegetation you add, more photons cannot be absorbed because they are already all absorbed. It is
normally the red band that saturates. So any index using the red energy will suffer from the same
limitation. For example, the Enhanced Vegetation Index (EVI) is not supposed to saturate as
badly because in the equation empirical constants have been added to put more weight in the
NIR spectrum that preserves sensitivity to higher loads of biomass (more layers of leafs) because
here much more radiation is transmitted and reflected from the leaves.
Comment: No, we are not talking about the same thing!!! I say the saturation stems from the specific equation applied (i.e. normalised difference). You say the saturation stems from the red band showing no changes. R_NIR << R_RED leads to NDVI=R_NIR/R_NIR=1. If you use a different index, e.g. the simple ratio R_RED/R_NIR there are not saturation issues if R_RED is small and changes little as long as R_NIR still changes. |