Referee 1 (RC1)

The manuscript discusses the impact of droughts on vegetation in Madagascar through the analysis of the relationships between the SPI and NDVI. Spatiotemporal characteristics of droughts were performed for the period 1981-2022, while the relationships between the indices were analysed over the time span of 2000-2022. Drought duration, frequency, severity and intensity were investigated. The topic and idea are important and interesting. In my opinion, the manuscript can be published after some corrections, which I assess as something between minor and major.

The manuscript is written in the correct language; however, in some parts, the sentences are very similar, which makes the text a bit humdrum (section discussing the spatial patterns).

The interpretation of the results in the section  “Spatial analysis of drought characteristics (Duration, Frequency, 220 Intensity, and Severity)” is poor - limited to the very general description of the figures. In the case of drought characteristics, one general sentence describes each of the characteristics and says no more than the index is higher or lower. After reading this section, no clear information remains in one's mind. The section lacks the interpretation of the charts. It can also be due to the language style, like  - Drought intensity shows the occurrence of moderate drought during the dry season. A sentence like “during dry season moderate droughts dominated” would be more direct and more transparent.

Response: Thank you for this comment. We have elaborated the interpretation of the whole spatial analysis of drought characteristics in section 3.2 and separated between the scales of SPI-3, -6, and -12 and seasonal and annual to refine the details. Line (L254-291)

The Division into the 3 regions – the justification of the Madagascar division into the three regions differing in the precipitation and vegetation conditions is convincing, however, the borders are like those made by the ruler. In nature, there are no such straight borders. Moreover, it would be informative if Fig.1 included the grids for precipitation data used in this study and some information on precipitation and air temperature in these regions to provide their differences.

Response: Thank you for this comment. We have added the annual mean precipitation and temperature over the study area (Fig.1). Concerning the division of the region and the straight borders; they are based on the characteristics of precipitation, temperature as plotted in Fig.1c-d, and vegetation of each region. It is for the sake of simplicity, of course it is simplified, but the main characteristics of the regions are coherent. Moreover, since SPI index mainly deals with precipitation data, the current study applies the division of the regions in connection to the general distribution of precipitation amount over the Island, as the amount of rainfall is well categorized from the annual mean climatology (Fig.1c). However, we have added a comment about the straight borders in line (L115-118): “Moreover, it is worth mentioning that in this study, the borders of the regions are only simple straight lines based on the distribution of annual mean precipitation and temperature over Madagascar. This is done for the sake of simplicity, but the regions are coherent mainly in terms of precipitation, which is the only input for the SPI calculation. Also, prevailing vegetation types are consistent within the three studied regions.”

Precipitation Data. The manuscript is based on model data, which usually differs from the station data. Two different data sources were used and the data from these sources probably differ. Several questions arise concerning these data. Were those data just combined – what about the differences between the sources? Why are two datasets used? It is not explained.  What is the resolution of the data? How much do they differ from station data? Randriatsara et al. (2022b) stated that the databases performed well for 1983-2015 while the manuscript uses the data for 1981-2022. What about the 1981 and 1982? Please refer to these questions.

Response: Thank you for pointing this out. Firstly, because Madagascar is facing some difficulties regarding the availability of station data (such as fewer available in situ datasets, missing values in several station datasets, and short record periods), our paper Randriatsara et al., (2022b) aimed to evaluate the performance of different gridded precipitation datasets (gauge-based, reanalysis and satellite) to find out the better performing gridded data, which could be used for any future research over the country. More information regarding the evaluation of these data is found in the paper (Randriatsara et al., 2022b, DOI: 10.1002/joc.7628)

Secondly, Randriatsara et al. (2022b) evaluated the time-period of 1983-2015 because for this period all 15 analyzed datasets were available. In present study, we focus on 1981-2022, but the choice of the dataset is based only on the evaluation in Randriatsara et al. (2022b) for the shorter time period. We consider this reasonable, as only an extra 9 years on top of 1983-2015.  

The reason for choosing these two datasets (CHIRPS and ERA5) has been added in the manuscript, please refer to the line (L127-129). The resolution of the datasets and the information of their combination are found in line (L134) and (L129-132), respectively.

Please extend the section regarding SPI. Deliver more information on how the SPIs were calculated. It will make catching the difference between the further SPI figures easier.

Response: Thank you for pointing this out. The description of the methodology has been modified (Section 2.3.1).

Fig. 1. Colours in the legend (colour scale) should be inverted. Usually, brown is used for high altitudes, and green is used for low elevations.

Response: Thank you for this suggestion. Changed as advised

lines 102-104: due to the shift of the convergence zone (?) – please check the expression with respect to the verb. It seems that “shift” (or other similar word) is lacking in the sentence.

Response: Thank you for this comment.

The original sentence was, “During the hot-wet season, the Inter-Tropical Convergence Zone (ITCZ) covers the northern part of the country due to the convergence of the northwest monsoon wind and the trade winds (Randriamarolaza et al., 2021), leading to rainfall across the whole of Madagascar except the southern parts.”

Modified sentence: “During the hot-wet season, the Inter-Tropical Convergence Zone (ITCZ) covers the northern part of the country as the northwest monsoon wind and the trade winds converge over this area (Randriamarolaza et al., 2021). This convergence leads to rainfall across the whole of Madagascar except the southern parts.” (L105-108)

line 172: It should be “the sum of months with SPI<1” instead of “the sum of SPI values <-1” – please check it.

Response: Thank you for pointing this out. It is corrected.

Fig.2. – Please shift the shaded bels to the background of the figure so that they ....

Response: Thank you for this suggestion. We have changed it to unfilled rectangles for a better view.

Section Regional values of monthly SPI-3, SPI-6 and SPI-12 – the description of the fig. 2-3 is very scarce. In the case of Fig. 3, there is only one very general statement in the text that there were more occurrences of moderate and severe drought events over the R1 than the 195 R2 and R3, especially in the recent decade. There is much more interesting information in these figures... Please interpret these figures in more detail, not only focusing on the recent years.

Response: Thank you for emphasizing this. We have improved the description of the results shown in Fig. 2-5.

Lines 197-198: However, it is worth mentioning that based on SPI -12 all the regions experienced simultaneous drought in the years 1991, 2006, 2017, and 2021 – I can identify 7 common dry periods in this figure.

Response: Thank you for pointing this out. There are 7 simultaneous occurrences of continuous negative values. However, the simultaneous occurrence of drought (i.e., SPI≤-1) occurred only over these four years (1991, 2006, 2017, and 2021).

Lines 203-204: The temporal development of dry season SPI in Fig. 5B exhibits almost the same patterns in all three regions – I can see the similarities until 2007, but later the patterns are different  (?), aren’t they?

Response: This comment is highly appreciated; we have modified accordingly also the whole description of Fig. 5 has been adjusted. We believe the description is clear and in agreement with your comment (L232-241).

Fig 6 – I am unsure if my impression is correct because there is no clear information on how the indices were calculated. Still, It seems that the differences in the duration for SPI-3, SPI-6 and SPI-12 are due to the method of calculations. If so, should the colour scale for these figures be the same?

Response: Thank you for this comment. The explanation about the calculation of the indices (SPI-3, -6 and -12) has been improved from the methodology section 2.3.1. The computation of drought characteristics (especially for duration) for SPI-3, -6 and -12 and the seasonal and annual scales has been also better elaborated in line (L202-204): “These written equations are based on the SPI-3, -6 and -12 timescales, however for case of seasonal and annual scales, number of months with SPI less than -1 (SPI≤-1) is replaced by numbers of seasons or years with SP less than -1. Also, the number of timesteps 504 months is replaced by 42 timesteps for seasonal and annual scales.”

Moreover, we have improved the figure 6. As already mentioned, to better understand the results, we have separated the scales of SPI-3, -6, and -12 from the seasonal and annual scales.

line 221-222: “... (up to 26 months in the southern R1 region (Fig. 6c) ...” – bracket is missing.

Response: Thank you. Added.

line 236-237: “Drought intensity shows the occurrence of moderate drought during the dry season ...” Does it mean moderate drought dominated? Please make the sentence more informative, particularly “show the occurrence”.

 Response: Thank you for this comment. It is changed as “Meanwhile, moderate drought (i.e. -1.0 ≥ SPI ≥ -1.49) dominates during the dry season (Fig. 6Bk) with values between 1.1 and 1.5.” (L286-287) 

line 233: “Drought intensity shows the occurrence of moderate drought during the dry season...” – Please deliver some numbers ...

Response: Added as advised. “Meanwhile, moderate drought (i.e. -1.0 ≥ SPI ≥ -1.49) dominates during the dry season (Fig. 6Bk) with values between 1.1 and 1.5.” (L286-287)

Does Fig. 6 show the average values of duration, frequency, severity and intensity? Please include this information in the capture.

Response: Thank you for pointing this out. Yes, Fig.6 shows the average values of duration, frequency and intensity, but for severity, it is the accumulation of all values ≤-1. Added as “The figure shows the averaged values of duration, frequency and intensity, but for severity, it represents the accumulation of all values ≤-1.” L(335-336)

The citations in lines 247-252 (copied below) indicate that no changes in droughts occurred in the study area. It is unclear that the authors' intention was to focus on the drought-driving factors.  Please rewrite the sentences to drive the reader’s attention to the factors. 

“According to Harrington et al. (2022), based on a combination of observations and climate modeling, the likelihood of poor rains experienced over the southern part of Madagascar was not found significantly increased due to human-caused climate change since it is overwhelmed by natural variability. However, a perceptible change in drought will emerge over the region if anthropogenic activities increase the global mean temperatures by more than 2 °C above preindustrial levels (IPCC 2021; Harrington et al., 2022).”

Response: Thank you for this comment. We apologize for the confusion. The sentence was not intended to indicate that there are no changes in the occurrence of drought but instead to show that the poor rains over south Madagascar was not found perceptibly linked to the human-caused climate change because of the domination of the natural variability.

So, we have paraphrased the sentence as: “According to Harrington et al. (2022), based on a combination of observations and climate modeling, the increase of poor rains experienced over the southern part of Madagascar was not significantly linked to the anthropogenic climate change because of the overwhelming of the natural variability. However, if the anthropogenic activities increase the global mean temperatures by more than 2 °C above preindustrial levels, the changes in drought will amplify (IPCC 2021; Harrington et al., 2022). This 305 confirms that even though the anthropogenic activities have not been significantly identified as the main cause of changes in drought due to the domination of natural variability, the increase of such activities will expose clearly the amplification of drought events.” (L302-308)

Please explain the abbreviations of IOD, SIOD also, SST, and ENSO, although the last two are commonly used.

Response: Thank you for pointing this out. Added accordingly in line (L309-310).

line 285-286: “The anomaly of NDVI at the end of the Event-I (Fig. 7d) shows that the area with negative values has extended with an increased amplitude of up to -0.2 compared to the beginning of the event (Fig. 7a).” Difficult to understand, particularly when it goes to amplitude. What is meant by amplitude here? The same in line 352 – negative amplitude? How can amplitude be negative? Amplitude has its particular meaning.

Response: Thank you for pointing this out. The sentence is changed as “The results agree with intuitive expectations in the case of Event-I: NDVI anomalies at the end of the Event-I (Fig. 7d) show that the areas with negative values (between -0.03 and -0.06) expand to the south and some parts of the north| compared to the beginning of the event (Fig. 7a).” (L348-350)

line 371: Is the decrease in NDVI for R1 statistically significant?

Response: yes, the linear trend for R1 is statistically significant at 95% confidence level with p-value = 0.00006877. Added in Line (L440-441)

lines 391-392: The correlation coefficients are very low. Regarding the large number of pixels are these relationships really so valid? In the further part of the manuscript, it is extended, but the statement of the existence of the relationships should be weaker in this sentence.

Response: Thank you for this comment.

The mentioned sentence in lines (391-392) is “However, it is worth mentioning that the southern region (R1) exhibits higher correlation coefficient values of about 0.26 to 0.29 compared to the regions R2 and R3 (between 0.05 and 0.19).”

The sentence is based on Fig.S6 which displays scatterplots correlations between the time-series of NDVI and SPIs averaged over the three regions. However, Fig.11 is the one that displays the correlations between NDVI and SPIs in each grid point throughout the whole time-period.

Regarding the low correlation coefficients found, we have paraphrased and adjusted some sentences as “It is worth mentioning that, even though the correlation coefficient values are generally low (i.e. less than 0.5), they are statistically significant at the 95% confidence level. This indicates that drought occurrences are indeed among the factors contributing to the deterioration of Madagascar’s vegetation. However, it is not only drought that caused the changes in vegetation, but there are also other factors, such as human-induced deforestation as the population relies heavily on fuelwood.” (L466-470)

Conclusions also include a kind of discussion. This discussion is interesting and valuable; however, it does not agree with the title of the section. I have some trouble with the structure because the discussion is also included in the particular sections... My advice is just to skip the discussion in the conclusions and include the results arising from your analysis. In such a case, shift the discussion parts to the results. The second possibility is to modify the section title to Discussion and Conclusions, which is also not the best solution (discussion is also included in the results). If none of the other reviewers raised that problem, skip this comment.

 Response: Thank you for this comment. However, as the second reviewer has not commented on this point, we decided to skip it.

Referee 2:

Review of ‘Historical changes in drought characteristics and its impact on vegetation cover over Madagascar’

I read the manuscript of Randriatsara et al with interest. Here, the author’s analyse spatio-temporal characteristics of droughts over Madagascar, together with MODIS based NDVI. They show that the  occurrence and intensity of droughts increased over the southern region, and that, especially over this region, there is a relationship between the SPI and NDVI.

The paper is clearly written and the impact of droughts on vegetation is a very relevant topic. I do have a few concerns regarding the analyses, and I believe that some of the conclusions are not supported by the current methodology or the (non-quantitative) description of the results. I will list my concerns below with suggestions to improve the methodology.

I have two concerns regarding the analyses of NDVI differences and NDVI time series.

 

1) The NDVI difference is calculated against a reference year (2001). Although the text mentions that 2001 was a relatively normal period, it appears to follow a severely dry period in R2 and R3. And any non-normal month during 2001 would lead to false conclusions about the NDVI difference later in the manuscript. Therefore, I recommend to calculate the NDVI anomaly relative to the monthly mean NDVI of the whole time series.

Response: We appreciate this comment. Regarding the choice of reference year 2001, it was based on the averaged SPI for the whole Madagascar from Fig.S1 (SPI-3, -6, and -12 scales) and Fig.S2 (seasonal and annual scales), instead of SPI for each region.

However, by checking SPI value for each region in 2001 from Fig. 2-5, there is an occurrence of wetter conditions over the Southern region (R1). This could be the reason for the intensified negative values over the south Madagascar during all the analysis while using 2001 as a reference.

Therefore, we have computed the difference against the NDVI mean throughout the study period and adjusted all the interpretations accordingly. Line (L330-340)

 

2) There is a trend in NDVI values for R1 for certain months. The manuscript mentions that this is related to climate change and natural variability (could rising CO2 concentrations also be one?). This means that the results presented in for example Fig. 9, L354 – 359, are (at least partly) the result of this decreasing trend, rather than the drought events itself. Detrending of the NDVI time series could possibly be used to separate the drought related effects from the negative trend effects. How does this negative trend impact the correlation results presented in for example Fig. 11 and Fig. S6?

 Response: Thank you for this comment. Yes, there is an NDVI trend for some months (from November to March), as displayed in Fig.S5. And the sentence which states about the natural variability and anthropogenic climate change from (L442-445) mentions that this decline of vegetation in southern Madagascar could be attributed to the delaying of the wet season rainfall, which is caused by natural variability and anthropogenic climate change. The sentence rather points out that the delayed wet season rainfall contributes to the insufficient rainfall from November to January, which results in vegetation decline from November to January (Fig. S5).

Obviously, it is related to the line (old reference: L354 – 359, new reference: (L417-423) because when the wet season rainfall is delayed, a negative drought index occurs for that season since SPI is calculated based on only rainfall data. But the purpose of including the natural variability and the anthropogenic climate change in the sentence is to mention what has caused this delayed wet season rainfall.

Moreover, we have performed and changed the correlation analyses to the correlations between SPIs and detrended NDVI anomaly (Fig.11 and S6). The results are very similar to the previous analysis when we used raw NDVI data without the calculation of detrended anomaly. However, the correlation coefficients have slightly decreased after the computation of detrended NDVI anomaly.

The NDVI anomalies in Fig. 7 – 9 are difficult to interpret without a reference value. What is the normal variation (or change) that one finds? Would it be an option to compare the values and changes to the 5 – 95% interval from all years? Furthermore, for Fig 7, the beginning and ending of the dry period might not reflect the worse drought period or the strongest vegetation response. I would recommend to compare the mean NDVI anomaly for the selected drought period with the normal anomaly that is found during this period.

Response: Thank you for this comment.

We believe that the confusion occurs due to the fact that we used the year 2001 as the reference. However, according to question 1), we have recomputed the differences against the NDVI mean as suggested. So, in this case, all values are regarded as anomalies relative to the NDVI mean.

Regarding Fig. 7, we have mentioned in the manuscript that whenever two different months are compared, the results will be irrelevant due to the seasonal cycle of vegetation regeneration. As we explained in this line (L364-369) “This suggests that it is not appropriate to assess vegetation loss based on a single month by referring to the starting and ending months of the continuous occurrence of negative SPI values. This is because whenever the starting month falls within the negative peak months (from August to November, Fig.S3 and S4) and the ending month falls in other months, the intensity of vegetation loss is always greater during these negative peak months than in other months. Therefore, the comparison of anomalies should be performed between the same month from different years. For instance, more decline in vegetation cover is found in January 2022 (Fig. S4a) compared to January 2016 (Fig. S3a). This is relevant without the influence of seasonal rainfall and can be used to examine the impact of negative SPI values on vegetation changes for a specific month from different years.”

 

L280 – 283: How are the start and end of the drought events determined? And how sensitive are the results in Fig 7 to the selection of the start and end of the drought event?

Response: Thank you for this comment. The start and the end were selected based on the simultaneous occurrence of negative values of SPI for each of the regions (Fig. 2-4), as well as the values for the whole island (Fig. S1). It is mentioned in line (L342-345): “The selection was based both on the SPI values for each of the regions from Fig. 2-4 (as marked with green rectangles in Fig. 2) and the SPI values averaged over the whole island (Fig. S1). The episodes will be further denoted as "Event-I" (spanning October 2005 - October 2006), "Event-II" (January 2016 - April 2017), and "Event-III" (September 2020 - December 2022).”

Regarding Fig. 7, we have mentioned in the text that for the case of Fig. 7 dealing with individual month is not appropriate due to the influence of seasonal rainfall. The sentence is from line (L364-369): “This suggests that it is not appropriate to assess vegetation loss based on a single month by referring to the starting and ending months of the continuous occurrence of negative SPI values. This is because whenever the starting month falls within the negative peak months (from August to November, Fig.S3 and S4) and the ending month falls in other months, the intensity of vegetation loss is always greater during these negative peak months than in other months. Therefore, the comparison of anomalies should be performed between the same month from different years. For instance, more decline in vegetation cover is found in January 2022 (Fig. S4a) compared to January 2016 (Fig. S3a). This is relevant without the influence of seasonal rainfall and can be used to examine the impact of negative SPI values on vegetation changes for a specific month from different years.”

Generally, the manuscript would benefit from more quantitative statements. I give three examples. 1. L190 – L191 “The SPI-3 (Fig. 2) shows that the occurrence of moderate to severe drought events (i.e. SPI values between -1 and -1.99) in the recent decade are more frequent over the southern region (R1) than over the western (R2) and eastern (R3) regions”. What is more frequent? Is this a significant change? Because I count a similar number of drought events in each region. 2. L228 – L229: “It is worth mentioning that for all the drought characteristics for all the three timescales, the magnitudes are higher over the southern region, especially for SPI-12, compared to the rest of the area”. Based on the figure, I do not agree with this statement, but if this is relevant, it would be good to have the corresponding values in a table. 3. L233 – 234: “Drought appears less frequent during the dry season, compared to the wet season”. What is less frequent? And ‘appears’; does this mean that there is a statistically significant difference?

Response: Thank you for pointing out this. We have improved the explanation of Fig.2-5 as well as all sentences from the section 3.2 (Spatial analysis of drought characteristics).

From the old reference in L190 – L191, the occurrence of moderate to severe drought events is the one that was mentioned to become more frequent. We changed the sentence to “The occurrence of moderate (i.e. SPI values between -1 and -1.49) to severe drought events (i.e. SPI values between -1 and -1.99) in the recent decade are more frequent over the southern region (R1) than over the western (R2) and eastern (R3) regions as shown in SPI-3 (Fig.2).  (L219-221)”.

Regarding the old reference to L228 – L229, as mentioned above, we have adjusted the figure and modified the description of the results. The interpretation is adjusted as “It is worth mentioning that overall, for all three timescales (Fig.6A), the drought characteristics’ magnitudes are higher over the eastern and southern regions, especially with prominent values for SPI-12 over the southern part, while some of the western and the central parts display lower values. (L266-269)”

Regarding the old reference L233 – 234, the sentence basically describes drought frequency; either drought occurs more frequently or less frequently. We, therefore, have paraphrased the sentence to better clarify it as: “Drought occurrence is more frequent during the wet season (Fig. 6Bd, frequency up to 23%), specifically over southern Madagascar, than during the dry season (Fig. 6Be).” (L275-277)

Fig 11 and Fig S6 represents the correlation between the NDVI and SPI. Are these monthly mean values? I have the impression that this figure is comparing seasonality in NDVI with anomalies in precipitation; most of the variation along the x-axis is the effect of seasonality, while the magnitude of the drought effect is much smaller. I would recommend to do this analyses with NDVI anomaly values, rather than absolute NDVI values.

 Response: Thank you for this comment. This is already done. As mentioned above, we have changed the correlations analyses in Fig.11 and S6 to correlations between SPIs and detrended NDVI anomaly as suggested.

Minor points

The study links the SPI (index for precipitation drought) to NDVI (agricultural droughts). However, vegetation is not affected by a lack of rainfall, it ‘cares about’ the vapour pressure deficit and the amount of water in the root zone. Agricultural droughts have a stronger link with soil moisture availability. However, precipitation data has a better availability, compared to soil moisture data, and therefore, the SPI (or SPEI) is a frequently used index for studying agricultural droughts. I recommend to add a bit of introduction or discussion on the processes linking drought, SPI, vegetation, and NDVI.

Response: We added the following sentence to the Introduction, when discussing the choice of indices: “Even though incorporating an index based on soil moisture would be beneficial for analyzing drought impacts on vegetation, SPI is frequently for studying agricultural droughts since it requires only precipitation data, which has better availability.” (L76-78)

I find the overview of drought characteristics in 2.3.3 very clear. It would be helpful for the reader to add that these drought characteristics are calculated for the three considered timescales (SPI3, SPI6, SPI12, and for the seasonal and annual). I do have two questions regarding these calculations. Is the Number_Timesteps lower for the seasonal and annual timescale? For Fig 6d-f: I do not understand what the drought duration of 4 or 9 months means for these seasonal and annual time scale. Could this be explained in the text?

Response: Thank you for pointing out this comment.

We have added explanations concerning this matter to the methodology section, in line (L201-204) as “These drought characteristics are calculated for the three considered timescales (SPI3, SPI6, SPI12, and for the seasonal and annual. These written equations are based on the SPI-3, -6 and -12 timescales, however for case of seasonal and annual scales, number of months with SPI less than -1 (SPI≤-1) is replaced by numbers of seasons or years with SP less than -1. Also, the number of timesteps 504 months is replaced by 42 timesteps for seasonal and annual scales.”

As well as, for the old Fig.6d-f, we have inserted more information in the caption of Fig.6, that for the seasonal and annual scales, the units are [seasons] and [years], respectively, in line (L330-336): “Figure 6: Spatial patterns of drought characteristics described in Section 2.3.3 (A) for SPI-3, -6, and -12 timescales and (B) for seasonal and annual scales. For (A), the units of Duration are [months], of Frequency [%] (the percentage of months with SPI ≤-1 relative to the number of all timesteps), of Severity the sum of SPI ≤-1, and of Intensity it is the average of the severity during the months with SPI ≤-1. For seasonal and annual scales (panel B), the units of the Duration are [seasons] for the case of season and [years] for the annual scale. Frequency, severity and intensity in (B) have the same units as in (A), only the timesteps are seasons and years instead of months. The figure shows the averaged values of duration, frequency and intensity, but for severity, it represents the accumulation of all values ≤-1”

For Eq. 1: the NDVI is usually calculated with specifically the Red part of the visible spectrum. Did you also use the wavelength of Red light to calculate the NDVI or all visible light?

Response: Thank you for pointing this out. We apologize for the confusion. The data we used is based on the Red band. So, we have adjusted for the right formula and information in line (L157)

L148: The website refers to Version 6, while you mention that version 6.1 was used.

Response: Thank you for this comment. Corrected to the right website in line (L164): “https://lpdaac.usgs.gov/products/mod13c2v061/”