Reply on RC2

<strong class="journal-contentHeaderColor">Abstract.</strong> Evapotranspiration (ET) and gross primary production (GPP) are critical fluxes contributing to the energy, water, and carbon exchanges between the atmosphere and the land surface. Land surface models such as the Community Land Model v5 (CLM5) quantify these fluxes, contribute to a better understanding of climate change's impact on ecosystems, and estimate the state of carbon budgets and water resources. Past studies have shown the ability of CLM5 to model ET and GPP magnitudes well but emphasized systematic underestimations and lower variability than in the observations. Here, we evaluate the simulated ET and GPP from CLM5 at the grid scale (CLM5_grid) and the plant functional type (PFT) scale (CLM5_PFT) with observations from eddy covariance stations from the Integrated Carbon Observation System (ICOS) over Europe. For most PFTs, CLM5_grid and CLM5_PFT compared better to ICOS than publicly available reanalysis data and estimates obtained from remote sensing. CLM5_PFT exhibited a low systematic error in simulating the ET of the ICOS measurements (average bias of -5.05 %), implying that the PFT-specific ET matches the magnitude of the observations closely. However, CLM5_PFT severely underestimates GPP, especially in deciduous forests (bias of -43.76 %). Furthermore, the simulated ET and GPP distribution moments across PFTs in CLM5_grid and CLM5_PFT, reanalyses, and remote sensing data indicate an underestimated spatiotemporal variability compared to the observations across Europe. These results are essential insights for further evaluations in CLM5 by pointing to the limitations of CLM5 in simulating the spatiotemporal variability of ET and GPP across PFTs.

Abstract. Critical soil moisture (CSM), a tipping point of soil moisture (SM) at which surface fluxes shift from the energy-limited regime to the water-limited regime, is essential for the vegetation state and the corresponding land–atmosphere coupling. However, detecting CSM and attributing water–energy limit shifts to climate and ecosystem variables are challenging as in situ observations of water, carbon fluxes, and soil moisture (SM) are sparse. In this study, CSM was assessed over China using two satellite-based methods: (i) the difference between the correlation between SM and evapotranspiration (ET) and the correlation between vapor pressure deficit (VPD) and ET and (ii) the covariance between VPD and gross primary production (GPP). ET and GPP products were based on the Penman–Monteith–Leuning (PML) ET and GPP, Global LAnd Surface Satellite (GLASS) ET and GPP, Collocation-Analyzed Multi-source Ensembled Land Evapotranspiration (CAMELE) ET, Surface Energy Balance Algorithm for Land (SEBAL) ET, two-leaf light use efficiency (TL-LUE) GPP, and solar-induced chlorophyll fluorescence (SIF)-based (GOSIF) GPP. At flux sites, ET and GPP products were evaluated by eddy-covariance-based measurements; CSM values using two satellite-based methods were assessed using the soil moisture–evaporative fraction method. Their consistency at site scales demonstrated reliable results and applicability to regional scales. Through intercomparison, the spatial pattern of CSM from multi-source ET and GPP datasets was consistent and robust in eastern and southern China. Generally, CSM decreased from south to north. The Pearl River basin and Southeastern River basin displayed a relatively high CSM for clay-rich soils (e.g., 0.39 m3 m−3 using PML ET and 10 cm depth SM) and forests (e.g., 0.35 m3 m−3 using PML ET and 10 cm depth SM). Since CSM values were higher than the average SM at four soil layers, grassland and clay were water-limited. Thus, with increased water demand, western grasslands were more susceptible to water stress. These findings highlight the variability in CSM and the primary determinants of water–energy limit shifts, offering valuable insights into the potential water limitations on ecosystems under comparable SM circumstances.

Integrating leaf functional traits improves modelled estimates of carbon and water fluxes at a subtropical evergreen conifer forest

Abstract. Rain gauges are unevenly spaced around the world with extremely low gauge density over developing countries. For instance, in some regions in Africa the gauge density is often less than one station per 10 000 km2. The availability of rainfall data provided by gauges is also not always guaranteed in near real time or with a timeliness suited for agricultural and water resource management applications, as gauges are also subject to malfunctions and regulations imposed by national authorities. A potential alternative is satellite-based rainfall estimates, yet comparisons with in situ data suggest they are often not optimal. In this study, we developed a short-latency (i.e. 2–3 d) rainfall product derived from the combination of the Integrated Multi-Satellite Retrievals for GPM (Global Precipitation Measurement) Early Run (IMERG-ER) with multiple-satellite soil-moisture-based rainfall products derived from ASCAT (Advanced Scatterometer), SMOS (Soil Moisture and Ocean Salinity) and SMAP (Soil Moisture Active and Passive) L3 (Level 3) satellite soil moisture (SM) retrievals. We tested the performance of this product over four regions characterized by high-quality ground-based rainfall datasets (India, the conterminous United States, Australia and Europe) and over data-scarce regions in Africa and South America by using triple-collocation (TC) analysis. We found that the integration of satellite SM observations with in situ rainfall observations is very beneficial with improvements of IMERG-ER up to 20 % and 40 % in terms of correlation and error, respectively, and a generalized enhancement in terms of categorical scores with the integrated product often outperforming reanalysis and ground-based long-latency datasets. We also found a relevant overestimation of the rainfall variability of GPM-based products (up to twice the reference value), which was significantly reduced after the integration with satellite soil-moisture-based rainfall estimates. Given the importance of a reliable and readily available rainfall product for water resource management and agricultural applications over data-scarce regions, the developed product can provide a valuable and unique source of rainfall information for these regions.

&amp;lt;p&amp;gt;Merging data from different instruments is required to construct long time data records of soil moisture (SM). This is the goal of projects such as the ESA Climate Change Initiative (CCI) for SM (Gruber et al., 2019), which uses both active and passive microwave sensors. Currently, the GLDAS v2.1 model is used as reference to re-scale active and passive time series by matching their Cumulative Density Function (CDF) to that of the model. Removing the dependency on models is important, in particular for data assimilation applications into hydrological or climate models, and it has been proposed (Van der Schalie et al., 2018) to use L-band data from one of the two instruments specifically designed to measure SM, ESA Soil Moisture and Ocean Salinity (SMOS) and NASA Soil Moisture Active Passive (SMAP) satellites, as reference to re-scale other time series.&amp;lt;br&amp;gt;To investigate this approach, AMSR-2 SM time series obtained from C1-, C2- and X-band observations using LPRM (Land Parameter Retrieval Model) were re-scaled by CDF-matching (Brocca et al., 2011) using different SMAP and SMOS official (SMAP L2 V005, SMOS L3 V300, SMOS NRT V100&amp;amp;V200) and research (SMOS IC V103) SM products as well as the SMAP and SMOS LPRM v6 SM data used by the ESA CCI. The time series re-scaled using L-band remote sensing data were compared to those re-scaled using GLDAS and were evaluated against in situ measurements at several hundred sites retrieved from the International Soil Moisture Network (Dorigo et al., 2011). The results were analyzed as a function of the land cover class and the Koppen-Geiger climate classification.&amp;lt;br&amp;gt;Overall, AMSR-2 time series re-scaled using SMAP L2, SMAP LPRM and SMOS IC data sets as reference gave the best correlations with respect to in situ measurements, similar to those obtained by the time series re-scaled using GLDAS and slightly better than those of the original AMSR-2 time series. These results imply that different SMAP and SMOS products could actually be used to replace GLDAS as reference for the re-scaling of other sensors time series within the ESA CCI. However, one must bear in mind that this study is limited to the re-scaling of AMSR-2 data at a few hundred sites.&amp;lt;br&amp;gt;For a more detailed assessment of the L-band data set to be used for a global re-scaling, it is necessary to investigate other effects such as the spatial coverage or the time series length. SMAP spatial coverage is better than that of SMOS in regions affected by radio frequency interference. In contrast, the length of SMAP time series can be too short to capture the long term SM variability for climate applications in some regions. The CDF of SMOS time series computed from the date of SMAP launch is significantly different to those of the full length SMOS time series in some regions of the Globe. Possible ways of using a coherent SMAP/SMOS L-band data set will be discussed.&amp;lt;/p&amp;gt;

Soil moisture (SM) plays an important role in the water and energy exchanges that occur in the terrestrial surface. Soil moisture can be retrieved at a larger scale by using the visible and InfraRed (IR) bands as well as through the microwave remote sensing. Because of very high importance of SM for variety of applications, there are now two dedicated microwave satellites in the Earth’s orbit for soil moisture retrieval from space. The first, Soil Moisture and Ocean Salinity (SMOS) satellite has been launched by the European Space Agency in November 2009 and second is Soil Moisture Active and Passive (SMAP) launched by the National Aeronautics and Space Administration (NASA) in January 2015. In this review, brief background of soil moisture retrieval algorithms are presented with different applications in the area of water resources. The first section provides the introduction of the soil moisture, presents several in situ techniques for measurement of soil moisture and soil moisture retrieval algorithms from visible/IR and microwave remote sensing. Section 2 describes the satellite soil moisture applications in water resources.

The long‐term monitoring of gross primary production (GPP) is crucial to the assessment of the carbon cycle of terrestrial ecosystems. In this study, a well‐known machine learning model (random forest, RF) is established to reconstruct the global GPP data set named ECGC_GPP. The model distinguished nine functional plant types, including C3 and C4 crops, using eddy fluxes, meteorological variables, and leaf area index (LAI) as training data of RF model. Based on ERA5_Land and the corrected GEOV2 data, global monthly GPP data set at a 0.05° resolution from 1999 to 2019 was estimated. The results showed that the RF model could explain 74.81% of the monthly variation of GPP in the testing data set, of which the average contribution of LAI reached 41.73%. The average annual and standard deviation of GPP during 1999–2019 were 117.14 ± 1.51 Pg C yr−1, with an upward trend of 0.21 Pg C yr−2 (p &lt; 0.01). By using the plant functional type classification, the underestimation of cropland is improved. Therefore, ECGC_GPP provides reasonable global spatial pattern and long‐term trend of annual GPP.

This paper presents the Soil Moisture and Ocean Salinity (SMOS) mission data assimilation activities conducted at the European Centre for Medium-Range Weather Forecasts (ECMWF) to analyse soil moisture for Numerical Weather Prediction (NWP) applications. Two different approaches are presented based on SMOS brightness temperature and SMOS neural network soil moisture data assimilation, respectively. For the first approach, SMOS brightness temperature data assimilation relies on forward modelling. Long term results, spanning the SMOS period, of SMOS forward modelling, monitoring and data assimilation are presented. They emphasize the relevance of SMOS data for monitoring and to support NWP model developments. For the second approach, a SMOS soil moisture product has been produced based on a Neural Network (NN) trained on ECMWF soil moisture. So, the SMOS-ECMWF NN soil moisture product captures the SMOS signal variability in time and space, while by design its climatology is consistent with that of the ECMWF soil moisture, which makes it suitable for data assimilation purpose. This approach, initially tested for 2012 in a global scale stand alone approach, shows that SMOS NN data assimilation slightly improves the two-metre air temperature forecast in the short range at regional scale. For NWP applications this approach has been further developed with a near real time production of the SMOS-ECMWF NN soil moisture product, with the implementation of the SMOS NN data assimilation in the ECMWF Integrated Forecasting System (IFS), and with high resolution (9km) global scale testing compatible with the current ECMWF NWP system.

Soil moisture is an important hydrologic variable due to its influence on plant growth and the surface energy budget. Most current soil moisture networks measure soil moisture at single points; these point measurements of soil moisture may not relate to the soil moisture of the larger, surrounding areas. Crop models and weather and climate models, which are sensitive to soil moisture, could be improved by having accurate soil moisture data that is relevant at scales larger than a point. Microwave C–band and X–band sensing satellites have attempted to retrieve soil moisture at relevant, large spatial scales needed for models, but these satellites are limited in their sensing depth, sensitivity to vegetation, and sensitivity to scattering by precipitation. Two recent satellite missions, the European Space Agency’s Soil Moisture and Ocean Salinity (SMOS) mission and NASA’s Soil Moisture Active Passive (SMAP) mission, retrieve soil moisture from L–band measurements of radiation, which improve upon the limitations of previous satellites. The retrieval of soil moisture from microwave satellite measurements requires knowledge of the vegetation optical thickness (τ), which is related to the vegetation water content (VWC). The multi-incidence angular design of SMOS provides its retrieval algorithm the advantage of retrieving soil moisture and τ simultaneously. SMAP will measure at a single incidence angle and will require an ancillary source of τ data. The current plan is for SMAP to use a climatology of normalized difference vegetation index (NDVI) data to generate a τ climatology. This dissertation analyzed this SMAP τ climatology by first validating SMOS τ and then by comparing SMOS τ to SMAP τ . The validation of SMOS τ was done by comparing SMOS τ during the growing season to crop yield estimates for each Iowa county. The result was that larger values of SMOS τ were correlated with larger crop yields, which suggested that SMOS τ is related to the growth of vegetation. SMOS τ was also found to be noisy and to behave strangely outside of the growing

Estimation of soil moisture at large scale has been performed using several satellite-based passive microwave sensors and a variety of retrieval methods over the past two decades. The most recent source of soil moisture is the European Space Agency Soil Moisture and Ocean Salinity (SMOS) mission. A thorough validation must be conducted to insure product quality that will, in turn, support the widespread utilization of the data. This is especially important since SMOS utilizes a new sensor technology and is the first passive L-band system in routine operation. In this paper, we contribute to the validation of SMOS using a set of four in situ soil moisture networks located in the U.S. These ground-based observations are combined with retrievals based on another satellite sensor, the Advanced Microwave Scanning Radiometer (AMSR-E). The watershed sites are highly reliable and address scaling with replicate sampling. Results of the validation analysis indicate that the SMOS soil moisture estimates are approaching the level of performance anticipated, based on comparisons with the in situ data and AMSR-E retrievals. The overall root-mean-square error of the SMOS soil moisture estimates is 0.043 m <sup xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">3</sup> /m <sup xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">3</sup> for the watershed networks (ascending). There are bias issues at some sites that need to be addressed, as well as some outlier responses. Additional statistical metrics were also considered. Analyses indicated that active or recent rainfall can contribute to interpretation problems when assessing algorithm performance, which is related to the contributing depth of the satellite sensor. Using a precipitation flag can improve the performance. An investigation of the vegetation optical depth (tau) retrievals provided by the SMOS algorithm indicated that, for the watershed sites, these are not a reliable source of information about the vegetation canopy. The SMOS algorithms will continue to be refined as feedback from validation is evaluated, and it is expected that the SMOS estimates will improve.

Sunlit and shaded leaf separation proposed by Norman (1982) is an effective way to upscale from leaf to canopy in modeling vegetation photosynthesis. The Boreal Ecosystem Productivity Simulator (BEPS) makes use of this methodology, and has been shown to be reliable in modeling the gross primary productivity (GPP) derived from CO2flux and tree ring measurements. In this study, we use BEPS to investigate the effect of canopy architecture on the global distribution of GPP. For this purpose, we use not only leaf area index (LAI) but also the first ever global map of the foliage clumping index derived from the multiangle satellite sensor POLDER at 6 km resolution. The clumping index, which characterizes the degree of the deviation of 3‐dimensional leaf spatial distributions from the random case, is used to separate sunlit and shaded LAI values for a given LAI. Our model results show that global GPP in 2003 was 132 ± 22 Pg C. Relative to this baseline case, our results also show: (1) global GPP is overestimated by 12% when accurate LAI is available but clumping is ignored, and (2) global GPP is underestimated by 9% when the effective LAI is available and clumping is ignored. The clumping effects in both cases are statistically significant (p &lt; 0.001). The effective LAI is often derived from remote sensing by inverting the measured canopy gap fraction to LAI without considering the clumping. Global GPP would therefore be generally underestimated when remotely sensed LAI (actually effective LAI by our definition) is used. This is due to the underestimation of the shaded LAI and therefore the contribution of shaded leaves to GPP. We found that shaded leaves contribute 50%, 38%, 37%, 39%, 26%, 29% and 21% to the total GPP for broadleaf evergreen forest, broadleaf deciduous forest, evergreen conifer forest, deciduous conifer forest, shrub, C4 vegetation, and other vegetation, respectively. The global average of this ratio is 35%.

Remote sensing-based global carbon flux and water products have gained worldwide prominence due to their feasibility to capture spatial patterns, interannual variability and long-term trend of ecosystem carbon and water dynamics from landscape to global scales. The Breathing Earth System Simulator (BESS), a sophisticated coupled process-based diagnostic terrestrial model, integrates multiple physical and bio-geological processes, including atmospheric transfer, canopy radiative transfer, and soil system. BESS has demonstrated robust performance across temporal and spatial scales from 1982 to 2019 when equipped with MODIS Atmosphere and Land products. This study aims to produce global gross primary productivity (GPP) and evapotranspiration (ET) products from 2018 to 2024 using the BESS model, utilizing datasets such as land surface temperature (LST), albedo, leaf area index (LAI), and shortwave radiation (SWR) from the GCOM-C SGLI satellite. Operated by JAXA, Japan, GCOM-C SGLI provides medium spatial resolution (250&amp;#8239;m to 1&amp;#8239;km) and an 8-day temporal resolution, offering improved sensor stability compared to BESS&amp;#8217;s default inputs from MODIS datasets. The main results of this study are: (1) The initial GPP and ET outputs showed fine showed fine linear relations with measurements of GPP (R2 = 0.56), and ET (R2 = 0.45) from 300 collected global flux sites from Fluxnet Ameriflux, Ozflux, Asiaflux; (2) Interannual variations from in response to climate extremes were captured, such as the 2020 Russian heatwave and the 2022 extreme dry summer in southwestern China. (3) From a sensitivity analysis of 21 BESS input parameters, 8 were selected and optimized using the Distributed Alternating Direction Method of Multipliers (D-ADMM), resulting in improved performance with GPP (R&amp;#178;&amp;#8239;=&amp;#8239;0.81) and ET (R&amp;#178;&amp;#8239;=&amp;#8239;0.72). Our GPP and ET products (250m and 1km) will be soon released on the JAXA website, aiming to improve better understanding of global climate and carbon cycle changes.

The Soil Moisture and Ocean Salinity (SMOS) mission delivers global surface soil moisture fields at high temporal resolution which is of major relevance for water management and climate predictions. Between April 26 and May 9, 2010, an airborne campaign with the L-band radiometer EMIRAD-2 was carried out within one SMOS pixel (44 × 44 km) in the Skjern River Catchment, Denmark. Concurrently, ground sampling was conducted within three 2 × 2 km patches (EMIRAD footprint size) of differing land cover. By means of this data set, the objective of this study is to present the validation of SMOS L1C brightness temperatures <i xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">TB</i> of the selected node. Data is stepwise compared from point via EMIRAD to SMOS scale. From ground soil moisture samples, <i xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">TB</i> 's are pointwise estimated through the L-band microwave emission of the biosphere model using land cover specific model settings. These <i xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">TB</i> 's are patchwise averaged and compared with EMIRAD <i xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">TB</i> 's. A simple uncertainty assessment by means of a set of model runs with the most influencing parameters varied within a most likely interval results in a considerable spread of <i xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">TB</i> 's (5-20 K). However, for each land cover class, a combination of parameters could be selected to bring modeled and EMIRAD data in good agreement. Thereby, replacing the Dobson dielectric mixing model with the Mironov model decreases the overall RMSE from 11.5 K to 3.8 K. Similarly, EMIRAD data averaged at SMOS scale and corresponding SMOS <i xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">TB</i> 's show good accordance on the single day where comparison is not prevented by strong radio-frequency interference (RFI) (May 2, avg. RMSE = 9.7 K). While the advantages of solid data sets of high spatial coverage and density throughout spatial scales for SMOS validation could be clearly demonstrated, small temporal variability in soil moisture conditions and RFI contamination throughout the campaign limited the extent of the validation work. Further attempts over longer time frames are planned by means of soil moisture network data as well as studies on the impacts of organic layers under natural vegetation and higher open water fractions at surrounding grid nodes.

Real surface soil moisture retrieved from the Soil Moisture and Ocean Salinity (SMOS) satellite is downscaled and assimilated in a fully integrated hydrological and soil‐vegetation‐atmosphere transfer (MIKE SHE SW‐ET) model using a bias aware ensemble transform Kalman filter (Bias‐ETKF). Satellite‐derived soil moisture assimilation in a catchment scale model is typically restricted by two challenges: (1) passive microwave is too coarse for direct assimilation and (2) the data tend to be biased. The solution proposed in this study is to disaggregate the SMOS bias using a higher resolution land cover classification map that was derived from Landsat thermal images. Using known correlations between SMOS bias and vegetation type, the assimilation filter is adapted to calculate biases online, using an initial bias estimate. Real SMOS‐derived soil moisture is assimilated in a precalibrated catchment model in Denmark. The objective is to determine if any additional gains can be achieved by SMOS surface soil moisture assimilation beyond the optimized model. A series of assimilation experiments were designed to (1) determine how effectively soil moisture corrections propagate downward in the soil column, (2) compare the efficacy of in situ versus SMOS assimilation, and (3) determine how soil moisture assimilation affects fluxes and discharge in the catchment. We find that assimilation of SMOS improved R 2 soil moisture correlations in the upper 5 cm compared to a network of 30 in situ sensors for most land cover classes. Assimilation also brought modest gains in R 2 at 25 cm depth but slightly degraded the correlation at 50 cm depth. Assimilation overcorrected discharge peaks.

Hudson Bay (HB) is the largest semi-inland sea in the Northern Hemisphere, connecting with the Arctic Ocean through the Foxe Basin and the northern Atlantic Ocean through the Hudson Strait. HB is covered by ice and snow in winter, which completely melts in summer. For about six months each year, satellite remote sensing of sea surface salinity (SSS) is possible over open water. SSS links freshwater contributions from river discharge, sea ice melt/freeze, and surface precipitation/evaporation. Given the strategic importance of HB, SSS has great potential in monitoring the HB freshwater cycle and studying its relationship with climate change. However, SSS retrieved in polar regions (poleward of 50°) from currently operational space-based L-band microwave instruments has large uncertainty (~ 1 psu) mainly due to sensitivity degradation in cold water (&lt;5°C) and sea ice contamination. This study analyzes SSS from NASA Soil Moisture Active and Passive (SMAP) and European Space Agency (ESA) Soil Moisture and Ocean Salinity(SMOS) missions in the context of HB freshwater contents. We found that the main source of the year-to-year SSS variability is sea ice melting, in particular, the onset time and places of ice melt in the first couple of months of open water season. The freshwater contribution from surface forcing P-E is smaller in magnitude comparing with sea ice contribution but lasts on longer time scale through the whole open water season. River discharge is comparable with P-E in magnitude but peaks before ice melt. The spatial and temporal variations of freshwater contents largely exceed the remote sensed SSS uncertainty. This fact justifies the use of remote sensed SSS for monitoring the HB freshwater cycle.