Characterizing uncertainty in Community Land Model version 5 hydrological applications in the United States

Land surface models such as the Community Land Model Version 5 (CLM5) are essential tools for simulating the behavior of the terrestrial system. Despite the extensive application of CLM5, limited attention has been paid to the underlying uncertainties associated with its hydrological parameters and how these uncertainties affect water resource applications. To address this long-standing issue, we use five meteorological datasets to conduct a comprehensive hydrological parameter uncertainty characterization of CLM5 over the hydroclimatic gradients of the conterminous United States. Key datasets produced from the uncertainty characterization experiment include: a benchmark dataset of CLM5 default hydrological performance, parameter sensitivities for 28 hydrological metrics, and large-ensemble outputs for CLM5 hydrological predictions. The presented datasets will assist CLM5 calibration and support broad applications, such as evaluating drought and flood vulnerabilities. The datasets can be used to identify the hydroclimatological conditions under which parametric uncertainties demonstrate substantial effects on hydrological predictions and clarify where further investigations are needed to understand how hydrological prediction uncertainties interact with other Earth system processes.

Inverse Modeling of Hydrologic Parameters in CLM4 via Generalized Polynomial Chaos in the Bayesian Framework

In this work, generalized polynomial chaos (gPC) expansion for land surface model parameter estimation is evaluated. We perform inverse modeling and compute the posterior distribution of the critical hydrological parameters that are subject to great uncertainty in the Community Land Model (CLM) for a given value of the output LH. The unknown parameters include those that have been identified as the most influential factors on the simulations of surface and subsurface runoff, latent and sensible heat fluxes, and soil moisture in CLM4.0. We set up the inversion problem in the Bayesian framework in two steps: (i) building a surrogate model expressing the input–output mapping, and (ii) performing inverse modeling and computing the posterior distributions of the input parameters using observation data for a given value of the output LH. The development of the surrogate model is carried out with a Bayesian procedure based on the variable selection methods that use gPC expansions. Our approach accounts for bases selection uncertainty and quantifies the importance of the gPC terms, and, hence, all of the input parameters, via the associated posterior probabilities.

Improved Land Evapotranspiration Simulation of the Community Land Model Using a Surrogate-Based Automatic Parameter Optimization Method

Land surface evapotranspiration (ET) is important in land-atmosphere interactions of water and energy cycles. However, regional ET simulation has a great uncertainty. In this study, a highly-efficient parameter optimization framework was applied to improve ET simulations of the Community Land Model version 4.0 (CLM4) in China. The CLM4 is a model at land scale, and therefore, the monthly ET observation was used to evaluate the simulation results. The optimization framework consisted of a parameter sensitivity analysis (also called parameter screening) by the multivariate adaptive regression spline (MARS) method and sensitivity parameter optimization by the adaptive surrogate modeling-based optimization (ASMO) method. The results show that seven sensitive parameters were screened from 38 adjustable parameters in CLM4 using the MARS sensitivity analysis method. Then, using only 133 model runs, the optimal values of the seven parameters were found by the ASMO method, demonstrating the high efficiency of the method. For the optimal parameters, the ET simulations of CLM4 were improved by 7.27%. The most significant improvement occurred in the Tibetan Plateau region. Additional ET simulations from the validation years were also improved by 5.34%, demonstrating the robustness of the optimal parameters. Overall, the ASMO method was found to be efficient for conducting parameter optimization for CLM4, and the optimal parameters effectively improved ET simulation of CLM4 in China.

Combinatorial Optimization for WRF Physical Parameterization Schemes: A Case Study of Three-Day Typhoon Simulations over the Northwest Pacific Ocean

Quantifying a set of suitable physics parameterization schemes for the Weather Research and Forecasting (WRF) model is essential for obtaining highly accurate typhoon forecasts. In this study, a systematic Tukey-based combinatorial optimization method was proposed to determine the optimal physics schemes of the WRF model for 15 typhoon simulations over the Northwest Pacific Ocean, covering all available schemes of microphysics (MP), cumulus (CU), and planetary boundary layer (PBL) physical processes. Results showed that 284 scheme combination searches were sufficient to find the optimal scheme combinations for simulations of track (km), central sea level pressure (CSLP, hPa), and 10 m maximum surface wind (10-m wind, m s−1), compared with the 700 sets of full combinations (i.e., 10 MP × 7 CU × 10 PBL). The decrease in the typhoon simulation error (i.e., root mean square error between simulation and observations) with this optimal scheme combination was 34%, 33.92%, and 25.67% for the track, CSLP, and 10-m wind, respectively. Overall, the results demonstrated that the optimal scheme combination yields reasonable results, and the Tukey-based optimization method is very effective and efficient in terms of computational resources.

Estimating land surface variables and sensitivity analysis for CLM and VIC simulations using remote sensing products

Assessment of Land Surface Models (LSMs) at heterogeneous terrain and climate regimes is essential for understanding complex hydrological and biophysical parameterization. This study utilized the two LSMs, Community Land Model (CLM 4.0) and three layer Variable Infiltration Capacity (VIC-3L), to estimate the interaction between land surface and atmosphere by means of energy fluxes including net radiation (R_N), sensible heat flux (H), latent heat flux (LE), and ground heat flux (G). The modeled energy fluxes were analyzed at two sites: Freeman Ranch-2 (FR2) located in the lowland region of Texas (272m), and Providence 301 (P301) located on the mountains of Sierra Nevada in California (2015m) from 2003 to 2013. R_N was underestimated by CLM with bias -25.06Wm^-2 due to its snow hydrology scheme at P301. LE was overestimated by the VIC during summer precipitation and had a positive bias of 5.51Wm^-2, whereas CLM showed a negative bias of -6.58Wm^-2 at the FR2 site. G was considered as a residual term in CLM, which caused weak performance at P301, while VIC calculated G as a function of soil temperature, depth, and hydraulic conductivity. In addition, The MOD16 showed similar results with models at FR2; however, at P301, they yielded a correlation value of 0.85 and 0.21 for LSMs and MOD16, respectively. The later has lower correlation with in situ specifically in summer season caused by erroneous biophysical or meteorological inputs to the algorithms. The sensitivity analysis between soil moisture and turbulent fluxes, exhibited negative trend (especially for LE at P301) due to topography and snow cover. The results from this study are conducive to improvements in models and satellite based characterization of water and energy fluxes, especially at rugged terrain with high elevation, where observational experiments are difficult to conduct.

Multi-Resolution Climate Ensemble Parameter Analysis with Nested Parallel Coordinates Plots

Due to the uncertain nature of weather prediction, climate simulations are usually performed multiple times with different spatial resolutions. The outputs of simulations are multi-resolution spatial temporal ensembles. Each simulation run uses a unique set of values for multiple convective parameters. Distinct parameter settings from different simulation runs in different resolutions constitute a multi-resolution high-dimensional parameter space. Understanding the correlation between the different convective parameters, and establishing a connection between the parameter settings and the ensemble outputs are crucial to domain scientists. The multi-resolution high-dimensional parameter space, however, presents a unique challenge to the existing correlation visualization techniques. We present Nested Parallel Coordinates Plot (NPCP), a new type of parallel coordinates plots that enables visualization of intra-resolution and inter-resolution parameter correlations. With flexible user control, NPCP integrates superimposition, juxtaposition and explicit encodings in a single view for comparative data visualization and analysis. We develop an integrated visual analytics system to help domain scientists understand the connection between multi-resolution convective parameters and the large spatial temporal ensembles. Our system presents intricate climate ensembles with a comprehensive overview and on-demand geographic details. We demonstrate NPCP, along with the climate ensemble visualization system, based on real-world use-cases from our collaborators in computational and predictive science.

Visualization of Time-Varying Weather Ensembles across Multiple Resolutions

Uncertainty quantification in climate ensembles is an important topic for the domain scientists, especially for decision making in the real-world scenarios. With powerful computers, simulations now produce time-varying and multi-resolution ensemble data sets. It is of extreme importance to understand the model sensitivity given the input parameters such that more computation power can be allocated to the parameters with higher influence on the output. Also, when ensemble data is produced at different resolutions, understanding the accuracy of different resolutions helps the total time required to produce a desired quality solution with improved storage and computation cost. In this work, we propose to tackle these non-trivial problems on the Weather Research and Forecasting (WRF) model output. We employ a moment independent sensitivity measure to quantify and analyze parameter sensitivity across spatial regions and time domain. A comparison of clustering structures across three resolutions enables the users to investigate the sensitivity variation over the spatial regions of the five input parameters. The temporal trend in the sensitivity values is explored via an MDS view linked with a line chart for interactive brushing. The spatial and temporal views are connected to provide a full exploration system for complete spatio-temporal sensitivity analysis. To analyze the accuracy across varying resolutions, we formulate a Bayesian approach to identify which regions are better predicted at which resolutions compared to the observed precipitation. This information is aggregated over the time domain and finally encoded in an output image through a custom color map that guides the domain experts towards an adaptive grid implementation given a cost model. Users can select and further analyze the spatial and temporal error patterns for multi-resolution accuracy analysis via brushing and linking on the produced image. In this work, we collaborate with a domain expert whose feedback shows the effectiveness of our proposed exploration work-flow.