Search NASA⌕ Search

Engineering topics

Wan, Hui

Publications and source records attributed to Wan, Hui.

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Data for publication Wan et al. (2020, JAMES), Improving time-step convergence in an atmosphere model with simplified physics: the impacts of closure assumption and process coupling

This dataset contains simple condensation mdoel output and CAM4 model output that were used to produce figures in the following paper: Wan, H., Woodward, C. S., Zhang, S., Vogl, C. J., Stinis, P., Gardner, D. J., Rasch, P. J., Zeng, X., Larson, V. E., and Singh, B. (2020): Improving time-step convergence in an atmosphere model with simplified physics: the impacts of closure assumption and process coupling, Journal of Advances in Modeling Earth Systems, https://doi.org/10.1029/2019MS001982

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Data from Wan et al. (2020, JAMES, DOI: 10.1029/2019MS001982)

This data set includes the simple condensation model output and CAM4 model output used to generate figures presented in the following paper: Wan, H., Woodward, C. S., Zhang, S., Vogl, C. J., Stinis, P., Gardner, D. J., Rasch, P. J., Zeng, X., Larson, V. E., and Singh, B. (2020): Improving time-step convergence in an atmosphere model with simplified physics: the impacts of closure assumption and process coupling, Journal of Advances in Modeling Earth Systems, https://doi.org/10.1029/2019MS001982

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Data from Zhang et al. (2019, JAMES, DOI: 10.1029/2019MS001817)

This data set includes the EAM model output used to generate figures presented in the following paper: Zhang, S., Wan, H., Rasch, P. J., Balwinder, S., Larson, V. E., and Woodward, C. S. (2019): An objective and efficient method for assessing the impact of reduced-precision calculations on solution correctness, Journal of Advances in Modeling Earth Systems, https://doi.org/10.1029/2019MS001817

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

An Objective and Efficient Method for Assessing the Impact of Reduced-Precision Calculations On Solution Correctness

Recent studies have shown that reducing the precision of floating-point calculations in an atmospheric model can improve the models' computational performance without affecting model fidelity, but code changes are needed to accommodate lower precision, and a small number of calculations might need to retain higher precision to prevent undue round-off error. For complex and computationally expensive systems like the Energy Exascale Earth System Model (E3SM), a method is needed to objectively assess the quality of the lower-precision simulations and to quickly identify problematic code pieces This paper demonstrates that solution correctness can be evaluated through a simple and quantitative error metric based on time step convergence. The proposed test method can unambiguously detect the accumulation of precision error and objectively assess its impact on solution accuracy. Since fast physical processes are known to affect key features of the multi-year mean climate in an atmospheric model, we show that the convergence test applied to short simulations can provide useful information about the impact of reduced precision on a model’s long-term behavior. In contrast, the traditional way of climate model evaluation requires multi-year simulations and involves inspecting many physical quantities for statistical significance in the presence of natural variability. The simplicity, computational efficiency, and objectivity of the proposed method makes it very attractive for high-resolution model development. The method is also expected to be applicable to other models that numerically solve time evolution equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Impact of light-absorbing particles on snow albedo darkening and associated radiative forcing over high-mountain Asia: high-resolution WRF-Chem modeling and new satellite observations

Light-absorbing particles (LAPs), mainly dust and black carbon, can significantly impact snowmelt and regional water availability over high-mountain Asia (HMA). In this study, for the first time, online aerosol–snow interactions are enabled and a fully coupled chemistry Weather Research and Forecasting (WRF-Chem) regional model is used to simulate LAP-induced radiative forcing on snow surfaces in HMA at relatively high spatial resolution (12 km, WRF-HR) compared with previous studies. Simulated macro- and microphysical properties of the snowpack and LAP-induced snow darkening are evaluated against new spatially and temporally complete datasets of snow-covered area, grain size, and impurity-induced albedo reduction over HMA. A WRF-Chem quasi-global simulation with the same configuration as WRF-HR buta coarser spatial resolution (1°, WRF-CR) is also used to illustrate the impact of spatial resolution on simulations of snow properties and aerosol distribution over HMA. Due to a more realistic representation of terrain slopes over HMA, the higher-resolution model (WRF-HR) shows significantly better performance in simulating snow area cover, duration of snow cover, snow albedo and snow grain size over HMA, as well as an evidently better atmospheric aerosol loading and mean LAP concentration in snow. However, the differences in albedo reduction from model and satellite retrievals is large during winter due to associated overestimation in simulated snow fraction. It is noteworthy that Himalayan snow cover has high magnitudes of LAP-induced snow albedo reduction (4% - 8%) in pre-monsoon seasons (both from WRF-HR and satellite estimates), which induces a snow-mediated radiative forcing of ~30-50 W m -2 . As a result, the Himalayas (specifically the western Himalayas) hold the most vulnerable glaciers and mountain snowpack to the LAP-induced snow darkening effect within HMA. In summary, coarse spatial resolution and absence of snow–aerosol interactions over the Himalayan cryosphere will result insignificant underestimation of aerosol effects on snow melting and regional hydroclimate.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Parametric Sensitivity and Uncertainty Quantification in the Version 1 of E3SM Atmosphere Model Based on Short Perturbed Parameter Ensemble Simulations

Abstract The atmospheric component of Energy Exascale Earth System Model version 1 has included many new features in the physics parameterizations compared to its predecessors. Potential complex nonlinear interactions among the new features create a significant challenge for understanding the model behaviors and parameter tuning. Using the one‐at‐a‐time method, the benefit of tuning one parameter may offset the benefit of tuning another parameter, or improvement in one target variable may lead to degradation in another target variable. To better understand the Energy Exascale Earth System Model version 1 model behaviors and physics, we conducted a large number of short simulations (three days) in which 18 parameters carefully selected from parameterizations of deep convection, shallow convection, and cloud macrophysics and microphysics were perturbed simultaneously using the Latin hypercube sampling method. From the perturbed parameter ensemble simulations and use of different skill score functions, we identified the most sensitive parameters, quantified how the model responds to changes of the parameters for both global mean and spatial distribution, and estimated the maximum likelihood of model parameter space for a number of important fidelity metrics. Comparison of the parametric sensitivity using simulations of two different lengths suggests that perturbed parameter ensemble using short simulations has some bearing on understanding parametric sensitivity of longer simulations. Results from this analysis provide a more comprehensive picture of the Energy Exascale Earth System Model version 1 behavior. The difficulty in reducing biases in multiple variables simultaneously highlights the need of characterizing model structural uncertainty (so‐called embedded errors) to inform future development efforts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Physics-Dynamics Coupling in Weather, Climate, and Earth System Models: Challenges and Recent Progress

Geophysical models of the atmosphere and ocean invariably involve parameterizations. These represent two distinct areas: a) Subgrid processes which the model cannot (yet) resolve, due to its discrete resolution, and b) sources in the equation, due to radiation for example. Hence coupling between these physics parameterizations and the resolved fluid dynamics and also between the dynamics of the different fluids in the system (air and water) is necessary. This coupling is an important aspect of geophysical models. However, often model development is strictly segregated into either physics or dynamics. Hence, this area has many more unanswered questions than in-depth understanding. Furthermore, recent developments in the design of dynamical cores (significant increase of resolution, move to non-hydrostatic equation sets etc), extended process physics (prognostic micro physics, 3D turbulence, non-vertical radiation etc) and predicted future changes of the computational infrastructure (Exascale with its need for task parallelism, data locality and asynchronous time stepping for example) is adding even more complexity and new questions.

Gross, Markus↗
DOE OSTI2018

Role of Troposphere-Convection-Land Coupling in the Southwestern Amazon Precipitation Bias of the Community Earth System Model version 1 (CESM1)

The Amazon forest plays an important role in the tropical circulation and global water and biogeochemical cycles. However, the common Amazon precipitation biases in global climate models undermine our ability to understand future changes in the Amazon and their global consequences. The Community Earth System Model exhibits a persistent wet season dry bias over the southwestern Amazon. The bias is also characterized by inaccurate diurnal timing and intensity of rainfall events. Analysis of a diverse set of simulations suggests that the bias originates from local processes. A known land model bias in the partitioning of surface energy fluxes is found not to contribute significantly to the rainfall bias. Key model processes are identified to be the low-level divergent circulation controlling horizontal moisture flow and the sensitivity of the modeled deep convection to the lower tropospheric moisture. The lower tropospheric moisture divergence dictates the spatiotemporal distributions of convective precipitation through parameterized entrainment mixing, which generally improves the rainfall diurnal cycle. However, it also limits convective precipitation over the areas with frequent and large-scale divergence in the lower troposphere. The southwestern Amazon is such an area, particularly in the late afternoon, although uncertainty exists in the low-level flow in the simulation and reanalyses. One way to sustain rainfall in the presence of low-level divergence is a convective trigger based on both the moisture and temperature changes by large-scale processes, which is shown to help in reducing the dry bias.

Sakaguchi, Koichi↗
DOE OSTI2018

Methodology for Building Granular Testing in Multicomponent Scientific Software

Computational science and engineering communities develop complex application software with multiple mathematical models that need to interact with one another. Partly due to complexity of verifying scientific software, and partly because of the way incentives work in science, there has been insufficient testing of these codes. With a spotlight on the results produced with scientific software, and increasing awareness of software testing and verification as a critical contributor to the reliability of these results, testing is gaining more attention by the developing teams. However, many science teams struggle to find a good solution for themselves due either to lack of training or lack of resources within the team. In this experience paper, we describe test development methodologies utilized in two different scenarios: one explains a methodology for building granular tests where none existed before, while the second demonstrates a methodology for selecting test cases that build confidence in the software through a process similar to scaffolding. The common insight from both the experiences is that testing should be a part of software design from the beginning for better software and scientific productivity.

test building, granular testing, Earth system mode↗
DOE OSTI2018

Impact of numerical choices on water conservation in the E3SM Atmosphere Model version 1 (EAMv1)

Abstract. The conservation of total water is an important numerical feature for global Earth system models. Even small conservation problems in the water budget can lead to systematic errors in century-long simulations. This study quantifies and reduces various sources of water conservation error in the atmosphere component of the Energy Exascale Earth System Model. Several sources of water conservation error have been identified during the development of the version 1 (V1) model. The largest errors result from the numerical coupling between the resolved dynamics and the parameterized sub-grid physics. A hybrid coupling using different methods for fluid dynamics and tracer transport provides a reduction of water conservation error by a factor of 50 at 1∘ horizontal resolution as well as consistent improvements at other resolutions. The second largest error source is the use of an overly simplified relationship between the surface moisture flux and latent heat flux at the interface between the host model and the turbulence parameterization. This error can be prevented by applying the same (correct) relationship throughout the entire model. Two additional types of conservation error that result from correcting the surface moisture flux and clipping negative water concentrations can be avoided by using mass-conserving fixers. With all four error sources addressed, the water conservation error in the V1 model becomes negligible and insensitive to the horizontal resolution. The associated changes in the long-term statistics of the main atmospheric features are small. A sensitivity analysis is carried out to show that the magnitudes of the conservation errors in early V1 versions decrease strongly with temporal resolution but increase with horizontal resolution. The increased vertical resolution in V1 results in a very thin model layer at the Earth's surface, which amplifies the conservation error associated with the surface moisture flux correction. We note that for some of the identified error sources, the proposed fixers are remedies rather than solutions to the problems at their roots. Future improvements in time integration would be beneficial for V1.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Investigation of short-term effective radiative forcing of fire aerosols over North America using nudged hindcast ensembles

Aerosols from fire emissions can potentially have large impact on clouds and radiation. However, fire aerosol sources are often intermittent, and their effect on weather and climate is difficult to quantify. Here we investigated the short-term effective radiative forcing of fire aerosols using the global aerosol–climate model Community Atmosphere Model version 5 (CAM5). Different from previous studies, we used nudged hindcast ensembles to quantify the forcing uncertainty due to the chaotic response to small perturbations in the atmosphere state. Daily mean emissions from three fire inventories were used to consider the uncertainty in emission strength and injection heights. The simulated aerosol optical depth (AOD) and mass concentrations were evaluated against in situ measurements and reanalysis data. Overall, the results show the model has reasonably good predicting skills. Short (10-day) nudged ensemble simulations were then performed with and without fire emissions to estimate the effective radiative forcing. Results show fire aerosols have large effects on both liquid and ice clouds over the two selected regions in April 2009. Ensemble mean results show strong negative shortwave cloud radiative effect (SCRE) over almost the entirety of southern Mexico, with a 10-day regional mean value of −3.0 W m −2 . Over the central US, the SCRE is positive in the north but negative in the south, and the regional mean SCRE is small (−0.56 W m −2 ). For the 10-day average, we found a large ensemble spread of regional mean shortwave cloud radiative effect over southern Mexico (15.6 % of the corresponding ensemble mean) and the central US (64.3 %), despite the regional mean AOD time series being almost indistinguishable during the 10-day period. Moreover, the ensemble spread is much larger when using daily averages instead of 10-day averages. This demonstrates the importance of using a large ensemble of simulations to estimate the short-term aerosol effective radiative forcing.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Impact of numerical choices on water conservation in the E3SM Atmosphere Model Version 1 (EAM V1)

The conservation of total water is an important numerical feature for global Earth system models. Even small conservation problems in the water budget can lead to systematic errors in century-long simulations for sea level rise projection. This study quantifies and reduces various sources of water conservation error in the atmosphere component of the Energy Exascale Earth System Model. Several sources of water conservation error have been identified during the development of the version 1 (V1) model. The largest errors result from the numerical coupling between the resolved dynamics and the parameterized sub-grid physics. A hybrid coupling using different methods for fluid dynamics and tracer transport provides a reduction of water conservation error by a factor of 50 at 1° horizontal resolution as well as consistent improvements at other resolutions. The second largest error source is the use of an overly simplified relationship between the surface moisture flux and latent heat flux at the interface between the host model and the turbulence parameterization. This error can be prevented by applying the same (correct) relationship throughout the entire model. Two additional types of conservation error that result from correcting the surface moisture flux and clipping negative water concentrations can be avoided by using mass-conserving fixers. With all four error sources addressed, the water conservation error in the V1 model is negligible and insensitive to the horizontal resolution. The associated changes in the long-term statistics of the main atmospheric features are small. A sensitivity analysis is carried out to show that the magnitudes of the conservation errors decrease strongly with temporal resolution but increase with horizontal resolution. The increased vertical resolution in the new model results in a very thin model layer at the Earth’s surface, which amplifies the conservation error associated with the surface moisture flux correction. We note that for some of the identified error sources, the proposed fixers are remedies rather than solutions to the problems at their roots. Future improvements in time integration would be beneficial for this model.

Zhang, Kai↗
DOE OSTI2017

A new and inexpensive non-bit-for-bit solution reproducibility test based on time step convergence (TSC1.0)

A test procedure is proposed for identifying numerically significant solution changes in evolution equations used in atmospheric models. The test issues a fail signal when any code modifications or computing environment changes lead to solution differences that exceed the known time step sensitivity of the reference model. Initial evidence is provided using the Community Atmosphere Model (CAM) version 5.3 that the proposed procedure can be used to distinguish rounding-level solution changes from impacts of compiler optimization or parameter perturbation, which are known to cause substantial differences in the simulated climate. The test is not exhaustive since it does not detect issues associated with diagnostic calculations that do not feedback to the model state variables. Nevertheless, it provides a practical and objective way to assess the significance of solution changes. The short simulation length implies low computational cost. The independence between ensemble members allows for parallel execution of all simulations, thus facilitating fast turnaround. The new method is simple to implement since it does not require any code modifications. We expect that the same methodology can be used for any geophysical model to which the concept of time step convergence is applicable.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Can nudging be used to quantify model sensitivities in precipitation and cloud forcing?

Efficient simulation strategies are crucial for the development and evaluation of high resolution climate models. This paper evaluates simulations with constrained meteorology for the quantification of parametric sensitivities in the Community Atmosphere Model version 5 (CAM5). Two parameters are perturbed as illustrating examples: the convection relaxation time scale (TAU), and the threshold relative humidity for the formation of low-level stratiform clouds (rhminl). Results suggest that the fidelity and computational efficiency of the constrained simulations depend strongly on 3 factors: the detailed implementation of nudging, the mechanism through which the perturbed parameter affects precipitation and cloud, and the magnitude of the parameter perturbation. In the case of a strong perturbation in convection, temperature and/or wind nudging with a 6-hour relaxation time scale leads to non-negligible side effects due to the distorted interactions between resolved dynamics and parameterized convection, while a 1-year free running simulation can satisfactorily capture the annual mean precipitation sensitivity in terms of both global average and geographical distribution. In the case of a relatively weak perturbation the large-scale condensation scheme, results from 1-year free-running simulations are strongly affected by noise associated with internal variability, while nudging winds effectively reduces the noise, and reasonably reproduces the response of precipitation and cloud forcing to parameter perturbation. These results indicate that caution is needed when using nudged simulations to assess precipitation and cloud forcing sensitivities to parameter changes in general circulation models. We also demonstrate that ensembles of short simulations are useful for understanding the evolution of model sensitivities.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Quantifying the impact of sub-grid surface wind variability on sea salt and dust emissions in CAM5

This paper evaluates the impact of sub-grid variability of surface wind on sea salt and dust emissions in the Community Atmosphere Model version 5 (CAM5). The basic strategy is to calculate emission fluxes multiple times, using different wind speed samples of a Weibull probability distribution derived from model-predicted grid-box mean quantities. In order to derive the Weibull distribution, the sub-grid standard deviation of surface wind speed is estimated by taking into account four mechanisms: turbulence under neutral and stable conditions, dry convective eddies, moist convective eddies over the ocean, and air motions induced by mesoscale systems and fine-scale topography over land. The contributions of turbulence and dry convective eddy are parameterized using schemes from the literature. Wind variabilities caused by moist convective eddies and fine-scale topography are estimated using empirical relationships derived from an operational weather analysis data set at 15 km resolution. The estimated sub-grid standard deviations of surface wind speed agree well with reference results derived from 1 year of global weather analysis at 15 km resolution and from two regional model simulations with 3 km grid spacing. The wind-distribution-based emission calculations are implemented in CAM5. In terms of computational cost, the increase in total simulation time turns out to be less than 3 %. Simulations at 2° resolution indicate that sub-grid wind variability has relatively small impacts (about 7 % increase) on the global annual mean emission of sea salt aerosols, but considerable influence on the emission of dust. Among the considered mechanisms, dry convective eddies and mesoscale flows associated with topography are major causes of dust emission enhancement. With all the four mechanisms included and without additional adjustment of uncertain parameters in the model, the simulated global and annual mean dust emission increase by about 50 % compared to the default model. By tuning the globally constant dust emission scale factor, the global annual mean dust emission, aerosol optical depth, and top-of-atmosphere radiative fluxes can be adjusted to the level of the default model, but the frequency distribution of dust emission changes, with more contribution from weaker wind events and less contribution from stronger wind events. Lastly, in Africa and Asia, the overall frequencies of occurrence of dust emissions increase, and the seasonal variations are enhanced, while the geographical patterns of the emission frequency show little change.

58 GEOSCIENCES↗
DOE OSTI2015

Resolution‐dependent behavior of subgrid‐scale vertical transport in the Z hang‐ M c F arlane convection parameterization

Abstract To better understand the behavior of quasi‐equilibrium‐based convection parameterizations at higher resolution, we use a diagnostic framework to examine the resolution‐dependence of subgrid‐scale vertical transport of moist static energy as parameterized by the Zhang‐McFarlane convection parameterization (ZM). Grid‐scale input to ZM is supplied by coarsening output from cloud‐resolving model (CRM) simulations onto subdomains ranging in size from 8 × 8 to 256 × 256 km 2 . Then the ZM‐based parameterization of vertical transport of moist static energy for scales smaller than the subdomain size ( ) are compared to those directly calculated from the CRM simulations ( ) for different subdomain sizes. The ensemble mean decreases by more than half as the subdomain size decreases from 128 to 8 km across while decreases with subdomain size only for strong convection cases and increases for weaker cases. The resolution dependence of is determined by the positive‐definite grid‐scale tendency of convective available potential energy (CAPE) in the convective quasi‐equilibrium (QE) closure. Further analysis shows the actual grid‐scale tendency of CAPE (before taking the positive definite value) and behave very similarly as the subdomain size changes because they are both tied to grid‐scale advective tendencies. We can improve the resolution dependence of significantly by averaging the grid‐scale tendency of CAPE over an appropriately large area surrounding each subdomain before taking its positive definite value. Even though the ensemble mean decreases with increasing resolution, its variability increases dramatically. cannot capture such increase in the variability, suggesting the need for stochastic treatment of convection at relatively high spatial resolution (8 or 16 km).

54 ENVIRONMENTAL SCIENCES↗