Characterizing natural variability in complex hydrological systems using Passive Microwave based Climate Data Records: a case study for the Okavango Delta
- 1VanderSat B.V., Water and Climate Unit, Haarlem, Netherlands
- 2European Space Agency Climate Office, ECSAT, Harwell Campus, Didcot, Oxfordshire, UK
- 3CLIMERS, TU Wien, Department of Geodesy and Geoinformation, Vienna, Austria
- 4Climate System Analysis Group, University of Cape Town, Cape Town, South Africa
- 1VanderSat B.V., Water and Climate Unit, Haarlem, Netherlands
- 2European Space Agency Climate Office, ECSAT, Harwell Campus, Didcot, Oxfordshire, UK
- 3CLIMERS, TU Wien, Department of Geodesy and Geoinformation, Vienna, Austria
- 4Climate System Analysis Group, University of Cape Town, Cape Town, South Africa
Abstract. The Okavango river system in southern Africa is known for its strong interannual variability of hydrological conditions. Here we present how this is exposed in surface soil moisture, land surface temperature, and vegetation optical depth as derived from the Land Parameter Retrieval Model using an inter-calibrated, long term, multi-sensor passive microwave satellite data record (1998–2020). We also investigate how these interannual variations relate to state-of-the-art climate reanalysis data from ERA5-Land. We analyzed both the upstream river catchment and the Okavango Delta, supported by independent data records of discharge measurements, precipitation and vegetation dynamics observed by optical satellites. The seasonal vegetation optical depth anomalies have a strong correspondence with MODIS Leaf Area Index (correlation catchment: 0.74, Delta: 0.88). Land surface temperature anomalies derived from passive microwave observations match best with those of ERA5-Land (catchment: 0.88, Delta: 0.81), as compared to MODIS nighttime LST (catchment: 0.70, Delta: 0.65). Although surface soil moisture anomalies from passive microwave observations and ERA5-Land correlate reasonably well (catchment: 0.72, Delta: 0.69), an in-depth evaluation over the Delta uncovered situations where passive microwave satellites record strong fluctuations, while ERA5-Land does not. This is further analyzed using information on inundated area, river discharge and precipitation. The passive microwave soil moisture signal demonstrates a response to both the inundated area and precipitation. ERA5-Land however, which by default does not account for any lateral influx from rivers, only shows a response to the precipitation information that is used as forcing. This also causes the reanalysis model to miss record low land surface temperature values as it underestimates the latent heat flux in certain years. These findings demonstrate the complexity of this hydrological system and suggest that future land surface model generations should also include lateral land surface exchange. Also, our study highlights the importance of maintaining and improving climate data records of soil moisture, vegetation and land surface temperature from passive microwave observations and other observation systems.
Robin van der Schalie et al.
Status: final response (author comments only)
-
RC1: 'Comment on hess-2021-637', Anonymous Referee #1, 15 Feb 2022
The paper investigates the behaviour of different hydrological and eco-physiological variables in the Okavango River region of Southern Africa. This river has the advantage of being an endoreic basin with a large floodplain at the end of the river. The paper shows a number of differences, particularly between soil moisture products (satellite vs. ERA-5) which are then analysed to show the value of satellite measurements compared to modelling data. The paper is well written and well organised. It does not provide any fundamentally new results but it shows that the lack of lateral influx in the ECMWF (ERA-5) model results in an underestimation of soil moisture and consequently of the ERA-5 evapotranspiration fluxes. The other interest of this paper is to show the evolution over 20 years of the river flow (just before the delta) and the flooded area in relation to satellite measurements of soil moisture, vegetation characteristics and soil temperature. I propose to accept the article with minor revisions.
Minor comments:
Line 143 : indicate millions of m3 instead of Mm3 (or E6 m3). Best would be to indicate also the mean annual river discharge. I guess this is around 300 m3/s.
Line 177: Explain why only descending TB are used in this study.
Line 184: Indicate the AMSR-2 orbital hours (Asc/Desc), same than AMSR-E ?
Line 192: Isn’t any local overpass of TRMM closer than AMSR (1:30 pm/am)? Explain better why 10:30 pm and 4:30 am are the best orbits.
Line 217-220: What does MD means ?
Line 224: Can authors give more details about the “E-type gauge plates” ?
Line 451: It is interesting to introduce a second set of precipitation data (IMERG). However, it would be interesting to show its co-evolution with ERA-PR for example in figure 5. How does the better correlation obtained with PR-E5 indicate a better rainfall product than IMERG?
Line 464: indicate in the figure caption 6 that this is ROI2
Line 533: indicate “(not shown)” after “This for example could also cause the difference in LST-MW and LST-E5 in 2010 and 2011”
Line 554: the end of the sentence is missing.
-
AC1: 'Reply on RC1', Robin van der Schalie, 14 Apr 2022
Dear anonymous referee #1,
Many thanks for your efforts in reviewing the paper, we are pleased to read the positive and constructive response. Below we will reply to the minor comments, either by clarifying our choices or by providing a plan to resolve issues in the revision.
Kind Regards,
The authors
Minor comments:
Line 143 : indicate millions of m3 instead of Mm3 (or E6 m3). Best would be to indicate also the mean annual river discharge. I guess this is around 300 m3/s.
Thank you for the advice, in the revision we will adjust the text according to state “millions of m3” and add that this is per year to put it in the right temporal perspective. That would indeed mean ~300 m3/s.
Line 177: Explain why only descending TB are used in this study.
For passive microwave retrievals (e.g. the soil moisture, land surface temperature and vegetation optical depth) the night time retrievals are of higher quality (e.g. Owe et al., 2008; Van der Schalie et al., 2021). This is caused by the thermal equilibrium during nighttime and the model assumption that the vegetation temperature is equal to that of the land surface temperature. This is not achieved during daytime, which in reality has a much higher variability over time (e.g minutes to hours). This is reflected back in the data quality, being more noisy. Especially as we include TRMM in the analysis, which does not have a stable local overpass time due to its non-polar orbit. Also, because of the timescale of the evaluation being monthly/seasonal, the daytime is not necessary to reach the goal of the study. We will clarify this better in the revised manuscript.
Line 184: Indicate the AMSR-2 orbital hours (Asc/Desc), same than AMSR-E ?
That is correct, AMSR2 and AMSR-E have similar overpass times. This will be clarified in the new revision of the manuscript.
Line 192: Isn’t any local overpass of TRMM closer than AMSR (1:30 pm/am)? Explain better why 10:30 pm and 4:30 am are the best orbits.
The orbital characteristics of TRMM cause the local overpass time to vary over time. So if you would only choose overpasses close to 1:30, there would be large temporal gaps in the dataset. Therefore, we have chosen to loosen this time constraint to include observations between 10:30 pm and 4:30 am. Of course, when available, we always choose the one that is closest to the AMSR2 time of overpass. One of the assumptions, supported by the paper of Van der Schalie et al. (2021) on using a similarly merged dataset for L-band retrieval input, is that especially for the soil moisture and vegetation optical depth these time differences have little impact. For the land surface temperature, in the short term (e.g. a few days) this can have an impact, however this is assumed to be smoothed out when looking at the timescales we look at within this paper. We will make sure to highlight these choices in the revision.
Line 217-220: What does MD means ?
We noticed that we indeed did not properly introduce that abbreviation. It is for MODIS and will be properly introduced in the revised version.
Line 224: Can authors give more details about the “E-type gauge plates” ?
Unfortunately at the moment this is all the information that can be found on the instrument used to do the measurements. We have reached out to the Okavango Research Institute of the University of Botswana for extra information and will add it as soon as we know more. (http://okavangodata.ub.bw/ori/monitoring/water/#)
Line 451: It is interesting to introduce a second set of precipitation data (IMERG). However, it would be interesting to show its co-evolution with ERA-PR for example in figure 5. How does the better correlation obtained with PR-E5 indicate a better rainfall product than IMERG?
Perhaps “better rainfall product” is a bit too general of a statement made. This relates to the results in the table that show that with a combination of PRE5 and the ODIAMD or ORD, we find a stronger fit to SSMMW than when we use PRIM. So here we assume that the best rainfall product would logically support a better fit with SSM. We understand that it might not be as straightforward as this, therefore we will clarify this assumption in the revised paper.
Line 464: indicate in the figure caption 6 that this is ROI2
Thank you for noticing this, we will adjust this in the revision.
Line 533: indicate “(not shown)” after “This for example could also cause the difference in LST-MW and LST-E5 in 2010 and 2011”
We will adjust this accordingly, or add these images in an appendix (under discussion).
Line 554: the end of the sentence is missing.
This mistake will be corrected in the new revision of the paper.
-
AC1: 'Reply on RC1', Robin van der Schalie, 14 Apr 2022
-
RC2: 'Comment on hess-2021-637', Anonymous Referee #2, 24 Mar 2022
-
AC2: 'Reply on RC2', Robin van der Schalie, 14 Apr 2022
Dear anonymous referee #2,
Many thanks for your efforts in reviewing the paper, we are pleased to read the positive and constructive response. Below we will reply to the individual comments, either by clarifying our choices or by providing a plan to resolve issues in the revision. With this, we hope we can take away the last concerns, especially concerning the figures and analysis.
Kind regards,
The authors
Majors comments:
1) Concerning the PMW CDR:
- Too little information is given on the inter-calibration and retrieval model. I acknowledge that these two compounds are supported by published materials but the author could add information in an annex on the inter-calibration (such as the cost function used in the optimization during inter-calibration). Also a figure with raw (original) CDR time series over the two regions of interest (ROIs) before retrieving the climatological mean.
Thank you for your advice. We thought the published material in Van der Schalie et al. (2021) would be sufficient to cover the intercalibration. However, as it is perceived as too little, we will take this into account for the revision and extend the section to include more information on the intercalibration.
It is possible to add a figure with the original (non-anomaly) data in the appendix. However, we explicitly made the choice during the paper writing to not include this because we want the main focus to be on their ability to properly detect anomalies. As we feel adding this might be distracting to the main message. Secondly, for the ROI2 we show the underlying climatology of different datasets in Figure 6.
- In Table 1, the resolution of the pixel measured by PMW at the surface is needed as it is said that Brightness Temperatures (BT) have been aggregated at 0.25° to obtain gridded products.
Thank you for this advice, we will make sure that in the revision we are more specific in distinguishing between the resolution of the gridded product and the resolution of the passive microwave footprints.
- If descending orbit only have been processed such that it can be compared with MODIS, such information is needed p7-L192
The descending orbits were not only chosen because of the MODIS comparison. For passive microwave retrievals (e.g. the soil moisture, land surface temperature and vegetation optical depth) the night time retrievals are also of higher quality. This is caused by the thermal equilibrium during nighttime and the model assumption that the vegetation temperature is equal to that of the land surface temperature. This is not achieved during daytime, which in reality has a much higher variability over time (e.g minutes to hours). This is reflected back in the data quality, being more noisy. Especially as we include TRMM in the analysis, which does not have a stable local overpass time due to its non-polar orbit. Also, because of the timescale of the evaluation being monthly/seasonal, the daytime is not necessary to reach the goal of the study. We will clarify this better in the revised manuscript to avoid confusion.
2) Concerning the LST analysis:
- Land surface temperature from ERA5 (LSTe5) has been extracted only for the first layer (0-7cm), is there any information on the penetration depth from the PMW observation? The infrared MODIS-based LSTmd is used for comparison, as infrared LST has no penetration depth, how this could impact the analysis. Please comment on p12-L315.
Penetration depth for Ka-band observations are about 1 mm, slightly varying with soil wetness (Holmes et al., 2013; Ulaby et al., 1986), so between the IR and ERA5 depths. This information will be added to the revision.
The mismatch in depth is also a reason why we choose for night time comparisons, as there is much more thermal stability expected. Especially when looking over longer periods (e.g. weeks, months) we assume that the slightly different definitions of soil temperature should still show a similarity in underlying anomalies. Not in the absolute sense, but relatively, as for example can be seen in the comparison of the anomalies in Figure 4. We will make sure that this is clear for the reader in the revision.
- No Time series is plotted for LSTe5 and LSTmd. For a systematic analysis, these two must be added in Figure 4. It should better support the author's statement on LSTe5 through the manuscript (p22-L533; p23-L564) and in the abstract. This is not shown in the analysis yet.
We agree that it would be better to show those images. We previously chose to leave those out to have manageable number of figures. For the revision we plan on either including it in Figure 4 or in an appendix.
- I would suggest adding LSTmd climatology in Figure 6
For both improved interpretability and after Section 4.3 showed a reduced skill of LSTMD as compared to LSTE5 and LSTMW, we made the decision to exclude LSTMD for this image. If requested, we can add it to the appendix.
3) Concerning the VOD:
- Climatology for ROi1 could be added in an annex to see if the LAI and VOD seasons are less correlated over catchment as it is stated p23-L595.
The statement relates to Figure 3, which is not on different seasonal dynamics within a single season, but more about the intraseasonal (season to season) comparison. Here you can see a more sustained increase in the anomalies for VOD in 2008-2012, as compared to the LAI. This is highlighted in the discussion and linked to the ability of VOD to detect the build-up of biomass in this longer wet period, which is not perceived in the LAI.
-Xband is less sensitive to leaves over dense forest, any experiment has been conducted in using/not using Xband for VOD in the omega-tau model?
The retrievals of soil moisture and vegetation optical depth are entangled, therefore you cannot replace the VOD within the algorithm used by something else in the Land Parameter Retrieval Model. X-band is not less sensitive to the leaves over dense forest, the issue is that the signal gets saturated with the vegetation signal and shows less variability therefore. Choosing another frequency, e.g. L-band, reduces the time coverage of the total record and therefore also the value of the anomalies. For example L-band only goes back to halfway 2010 with SMOS.
4) Concerning the Figures:
- All scatterplots must have Xlabel and Ylabel for clearer reading. I suggest introducing correlation numbers inside the figures.
In the revised version we will add the X and Y labels where missing.
- All correlation numbers must have at most 3 digits as the 4th is not meaningful.
We agree with this statement, and therefore will adjust where necessary to have a maximum of 3 digits (like 0.77) for correlations.
- "Absolute anomalies" in the title is misleading as "absolute" has another mathematical meaning. Replace by “raw”?
We agree that it is currently not completely clear what “absolute” anomalies mean in the context of this study. However, we think “raw” is also confusing. Therefore we will highlight the meaning of the absolute anomalies in the context of this paper. Concerning the z-score, we will explain its meaning as “standardized anomalies”.
- Add PRim in Figure 5 as well as PRe5 and PRim for RO1 (can be in annexes) as it is stated that PRe5 has high positive anomalies over the catchment (p21-L478) with no supporting information.
To make sure all the information is included in the figures, we will include the PRE5 (ROI1) and PRIM (ROI1/2)in Figure 5 so it includes information on all support data.
5) Concerning the Linear regression experiment:
- RMSE for Z-score is difficult to analyze, pleased replace by bias and std metrics in table2
Due to the completely different underlying values for river discharge, inundated area, precipitation, and soil moisture, we decided to use the z-score anomaly as a tool to normalize this difference and improve the interpretability of the balance between their contributions in the regression activity. The advantage of using z-scores or normalized values in regression is that regression coefficients are directly comparable or interpretable in terms of strength of relationship between dependent and different explanatory variables. In the revision we will better clarify this choice. On top of that, in a regression exercise the bias does not have much meaning.
- Please consider doing the Linear regression experiment for the catchment ROI1 to see if the SSM is more related to the precipitation upstream (as stated p22-L517).
A linear regression activity does not add the same value as in the Delta. Because ERA5-Land is mostly driven by its precipitation forcing, it will always match best there compared to anything else, without providing any insight into the true quality. The unique opportunity in the Delta is the contribution of the Okavango River inflow and the inundation, so the water is not coming from an individual source.
- In the Table specify the considered ROI.
Thank you for your advice, we will specify the ROI.
- It said that OIAD show some lagging from ORD, could you find optimal lag with cross-correlation between ORD-SSM and SSM-OIAD. This might lead to finding some buffering effect in SSM between ORD and OIAD.
That is indeed an interesting suggestion. We will have a look to see if this is feasible and potentially include in the revision.
Minors comments :
-p2-L36, miss-record
-p2-L51, “the The”
-p3-L63, BAMS, acronym is not defined
-p3-L93, PMW is not defined,
-p4-L111, use Section instead of “Chapter”
-p4-L114, LPRM is not defined yet E5L should be E5
-p7-L195, The sentence is misleading since not only the Xband is used
-p9-L224, what is an E-type gauge?
-p11-L289, VODCA is not defined
-add Z-score equation
-Caption of figure 3 seems misplaced (not attached to the figure on the same page)
-p21-L474, what “memory” replaces to buffer effect?
-p21-L489, SSM not SM-In Fig6: +-15d is used for visualization only or to compute anomalies also? If it has been used for
computing anomalies, this could lead to over-smooth the anomalies with the 90 days moving average
window.
-p23-L534: could you be more specific. The increase of SMM with available solar energy, increases ET and
avoids a false increase of LST but how is LSTe5 between 2011-2014?
-p23-L553 verb is missing
We thank you for the thorough review and will correct or clarify all minor points in the revision.
-
AC2: 'Reply on RC2', Robin van der Schalie, 14 Apr 2022
Robin van der Schalie et al.
Robin van der Schalie et al.
Viewed
HTML | XML | Total | BibTeX | EndNote | |
---|---|---|---|---|---|
562 | 99 | 18 | 679 | 5 | 5 |
- HTML: 562
- PDF: 99
- XML: 18
- Total: 679
- BibTeX: 5
- EndNote: 5
Viewed (geographical distribution)
Country | # | Views | % |
---|
Total: | 0 |
HTML: | 0 |
PDF: | 0 |
XML: | 0 |
- 1