Floods are becoming more frequent and are hard to control due to the shortage of water conservancy projects in small and medium sized basins, especially in developing regions. Understanding the hydrologic responses and estimated flood characteristics to storm events can help to predict flood disasters and form better effective mitigation and adaptation strategies. Using observations of several representative watersheds and Artificial Neural Networks (ANN) and Principal Component Regression (PCR) models, we conducted data-driven analyses to examine the flood responding characteristics. We quantified the relative contributions of the influencing factors to the variation of each flood characteristic. Statistical analysis of the observations shows that the drainage area plays a key role in determining the distribution of lag time and peak discharge. Rainstorm variability has direct influence on floods, and typhoon-induced rainstorms with high total rainfall and rainfall intensity generate higher lag time, flood peak, unit discharge and runoff depth, but lower runoff coefficient. The ANN and PCR models accurately predicted the variations of flood features using the driving factors including physical geographical characteristics, rainstorm features, and antecedent condition. Physical geographic characteristics are key influential factors of lag time, flood peak and runoff coefficient, while the rainstorm features control the magnitude of unit discharge and runoff depth. These results indicate that floods are mainly affected by rainstorm features and physical geographic characteristics in the Yangtze River delta, and it might become more damaging with the increasing rainfall extremes and sprawling impervious surfaces in a changing environment.

Evapotranspiration (ET) plays an important role in land-atmosphere coupling of energy, water, and carbon cycles. Following deforestation, ET is typically observed to decrease substantially as a consequence of decreases in leaf area and roots and increases in runoff. Changes in ET (latent heat flux) revise the surface energy and water budgets, which further affects large-scale atmospheric dynamics and feeds back positively or negatively to long-term forest sustainability. In this study, we used observations from a recent synthesis of 29 pairs of adjacent intact and deforested FLUXNET sites to improve model parameterization of stomatal characteristics, photosynthesis, and soil water dynamics in version 1 of the Energy Exascale Earth System Model (E3SM) Land Model (ELMv1). We found that default ELMv1 predicts an increase in ET after deforestation, likely leading to incorrect estimates of the effects of deforestation on land-atmosphere coupling. The calibrated model accurately represented the FLUXNET observed deforestation effects on ET. Importantly, the search for global optimal parameters converged at values consistent with recent observational syntheses, confirming the reliability of the calibrated physical parameters. Applying this improved model parameterization to the globe scale reduced the bias of annual ET simulation by up to ~600 mm/year. Analysis on the roles of parameters suggested that future model development to improve ET simulation should focus on stomatal resistance and soil water-related parameterizations. Finally, our predicted differences in seasonal ET changes from deforestation are large enough to substantially affect land-atmosphere coupling and should be considered in such studies.

At the end of its first year of operation, we compare soil moisture retrievals from the Soil Moisture Active Passive (SMAP) mission to simulations from a land surface model with meteorological forcing downscaled from observations/reanalysis and in situ observations from sparse monitoring networks within continental United States (CONUS). The radar failure limits the duration of comparisons for the active and combined products (~3 months). Nevertheless, the passive product compares very well against in situ observations over CONUS. On average, SMAP compares to the in situ data even better than the land surface model and provides significant added value on top of the model and thus good potential for data assimilation. At large scale, SMAP is in good agreement with the model in most of CONUS with less-than-expected degradation over mountainous areas. Lower correlation between SMAP and the model is seen in the forested east CONUS and significantly lower over the Canadian boreal forests.

Accurate simulation of energy, water, and carbon fluxes exchanging between the land surface and the atmosphere is beneficial for improving terrestrial ecohydrological and climate predictions. We systematically assessed the Noah land surface model (LSM) with mutiparameterization options (Noah-MP) in simulating these fluxes and associated variations in terrestrial water storage (TWS) and snow cover fraction (SCF) against various reference products over 18 United States Geological Survey two-digital hydrological unit code regions of the continental United States (CONUS). In general, Noah-MP captures better the observed seasonal and interregional variability of net radiation, SCF, and runoff than other variables. With a dynamic vegetation model, it overestimates gross primary productivity by 40% and evapotranspiration (ET) by 22% over the whole CONUS domain; however, with a prescribed climatology of leaf area index, it greatly improves ET simulation with relative bias dropping to 4%. It accurately simulates regional TWS dynamics in most regions except those with large lakes or severely affected by irrigation and/or impoundments. Incorporating the lake water storage variations into the modeled TWS variations largely reduces the TWS simulation bias more obviously over the Great Lakes with model efficiency increasing from 0.18 to 0.76. Noah-MP simulates runoff well in most regions except an obvious overestimation (underestimation) in the Rio Grande and Lower Colorado (New England). Compared with North American Land Data Assimilation System Phase 2 (NLDAS-2) LSMs, Noah-MP shows a better ability to simulate runoff and a comparable skill in simulating Rn but a worse skill in simulating ET over most regions. This study suggests that future model developments should focus on improving the representations of vegetation dynamics, lake water storage dynamics, and human activities including irrigation and impoundments.

To prepare for the next-generation North American Land Data Assimilation System (NLDAS), three advanced land surface models [LSMs; i.e., Community Land Model, version 4.0 (CLM4.0); Noah LSM with multiphysics options (Noah-MP); and Catchment LSM-Fortuna 2.5 (CLSM-F2.5)] were run for the 1979-2014 period within the NLDAS-based framework. Unlike the LSMs currently executing in the operational NLDAS, these three advanced LSMs each include a groundwater component. In this study, the model simulations of monthly terrestrial water storage anomaly (TWSA) and its individual water storage components are evaluated against satellite-based and in situ observations, as well as against reference reanalysis products, at basinwide and statewide scales. The quality of these TWSA simulations will contribute to determining the suitability of these models for the next phase of the NLDAS. Overall, it is found that all three models are able to reasonably capture the monthly and interannual variability and magnitudes of TWSA. However, the relative contributions of the individual water storage components to TWSA are very dependent on the model and basin. A major contributor to the TWSA is the anomaly of total column soil moisture content for CLM4.0 and Noah-MP, while the groundwater storage anomaly is the major contributor for CLSM-F2.5. Other water storage components such as the anomaly of snow water equivalent also play a role in all three models. For each individual water storage component, the models are able to capture broad features such as monthly and interannual variability. However, there are large intermodel differences and quantitative uncertainties, which are motivating follow-on investigations in the NLDAS Science Testbed developed by the NASA and NCEP NLDAS teams.

Previous studies have documented that surface wind speed (u) has been increasing over the ocean but decreasing over land for the past several decades. The decreasing u at the surface over land has been referred to as terrestrial stilling. A plausible hypothesis for terrestrial stilling is an increase in surface roughness associated with changes in land surface (e.g. enhanced vegetation growth, landscape fragmentation or urbanization). One of the most widespread land surface changes is enhanced vegetation leaf area index (LAI) known as greening, particularly over the middle to high latitudes of the Northern Hemisphere where strong stilling is observed from weather station data. In this study, we examine the hypothesis that enhanced vegetation LAI is a key driver of global terrestrial stilling. We first characterized the trend in u over the ocean using long-term satellite altimeter measurements, and the trend in u over land using continuous wind records from 4305 in situ meteorological stations. We then performed initial condition ensemble Atmospheric Model Intercomparison Project-type simulations using two state-of-the-art Earth system models (IPSL-CM and CESM) to isolate the response of u to the historical increase in LAI (representing the greening) for the period 1982-2011. Both models, forced with observed sea surface temperature and sea ice and with LAI from satellite observation, captured the observed strengthening of Pacific trade winds and Southern Ocean westerly winds. However, these simulations did not reproduce the weakening of surface winds over land as significantly as it appears in the observations (-0.006 m s -1 versus -0.198 m s -1 during 1982-2011), indicating that enhanced LAI (greening) is not a dominant driver for terrestrial stilling.

Accurate global mapping of soil moisture is the goal of the Soil Moisture Active Passive (SMAP) mission, which is expected to improve the estimation of water, energy, and carbon exchanges between the land and the atmosphere. Like other satellite products, the SMAP soil moisture retrievals need to be validated, with the validation relying heavily on in situ measurements. However, a one-to-one comparison is ill advised due to the spatial mismatch of the large SMAP footprint (∼40 km) and the point scale in situ measurements. This study uses a recently developed hyper-resolution land surface model—HydroBlocks—as a tool to upscale in situ soil moisture measurements for the SMAPVEX15 (SMAP Validation Experiment 2015) field campaign during 2–18 August 2015. Calibrated against in situ observation, HydroBlocks shows a satisfactory Kling-Gupta efficiency (KGE) of 0.817 and RMSE of 0.019 m3/m3 for the calibration period. These results indicate that HydroBlocks can be used to upscale in situ measurements for this site. Different from previous studies, here in situ measurements are upscaled using a land surface model without bias correction. The upscaled soil moisture is then used to evaluate SMAP (passive) soil moisture products. The comparison of the upscaled network to SMAP shows that the retrievals are generally able to capture the areal-averaged soil moisture temporal variations. However, SMAP appears to be oversensitive to summer precipitation. We expect these findings can be used to improve the SMAP soil moisture product and thus facilitate its usage in studying the water, energy, and carbon cycles.

The diversity in hydrologic models has historically led to great controversy on the "correct" approach to process-based hydrologic modeling, with debates centered on the adequacy of process parameterizations, data limitations and uncertainty, and computational constraints on model analysis. In this paper, we revisit key modeling challenges on requirements to (1) define suitable model equations, (2) define adequate model parameters, and (3) cope with limitations in computing power. We outline the historical modeling challenges, provide examples of modeling advances that address these challenges, and define outstanding research needs. We illustrate how modeling advances have been made by groups using models of different type and complexity, and we argue for the need to more effectively use our diversity of modeling approaches in order to advance our collective quest for physically realistic hydrologic models.

Land-use and land-cover change significantly modify local land-surface characteristics and water/energy exchanges, which can lead to atmospheric circulation and regional climate changes. In particular, deforestation accounts for a large portion of global land-use changes, which transforms forests into other land cover types, such as croplands and grazing lands. Many previous efforts have focused on observing and modeling land–atmosphere–water/energy fluxes to investigate land–atmosphere coupling induced by deforestation. However, interpreting land–atmosphere–water/energy-flux responses to deforestation is often complicated by the concurrent impacts from shifts in land-surface properties versus background atmospheric forcings. In this study, we used 29 paired FLUXNET sites, to improve understanding of how deforested land surfaces drive changes in surface-energy-flux partitioning. Each paired sites included an intact forested and non-forested site that had similar background climate. We employed transfer entropy, a method based on information theory, to diagnose directional controls between coupling variables, and identify nonlinear cause–effect relationships. Transfer entropy is a powerful tool to detective causal relationships in nonlinear and asynchronous systems. The paired eddy covariance flux measurements showed consistent and strong information flows from vegetation activity (gross primary productivity (GPP)) and physical climate (e.g. shortwave radiation, air temperature) to evaporative fraction (EF) over both non-forested and forested land surfaces. More importantly, the information transfers from radiation, precipitation, and GPP to EF were significantly reduced at non-forested sites, compared to forested sites. We then applied these observationally constrained metrics as benchmarks to evaluate the Energy Exascale Earth System Model (E3SM) land model (ELM). ELM predicted vegetation controls on EF relatively well, but underpredicted climate factors on EF, indicating model deficiencies in describing the relationships between atmospheric state and surface fluxes. Moreover, changes in controls on surface energy flux partitioning due to deforestation were not detected in the model. We highlight the need for benchmarking model simulated surface-energy fluxes and the corresponding causal relationships against those of observations, to improve our understanding of model predictability on how deforestation reshapes land surface energy fluxes.

Land surface models range in complexity of terrestrial evapotranspiration, yet it is unknown how model complexity translates to accuracy of modeled evapotranspiration estimates. Here, we use the International Land Model Benchmarking system to assess ET estimates from three models of varying complexity driven by the same forcing datasets: an earth system model, a terrestrial biosphere model, and a stand-alone ET model. The performance assessment includes both temporal and spatial evaluation, and different plant functional types across China. Our results indicate that the most complex model, an earth system model, performed best against the benchmarking datasets and metrics. Terrestrial biosphere model performed best in simulating inter-annual variability of ET, while earth system model performed best in simulating the seasonal cycle. The more complex models (earth system model and terrestrial biosphere model) perform better in forest, shrub and crop ecosystems, while the simpler model (stand-alone ET model) perform better in grass ecosystems. Our study demonstrates the impact of model complexity on ET estimates and highlights directions for future ET model improvements.