Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parameterization”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,009 records · Page 56

Quantifying the Contributions of Aerosol- and Snow-Produced ClNO 2 through Observations and 1D Modeling

Nitryl chloride (ClNO 2 ) is a radical reservoir that forms and accumulates in the nocturnal atmospheric boundary layer influenced by combustion emissions and chloride (e.g., sea salt and road salt). Upon sunrise, ClNO 2 rapidly photolyzes to generate highly reactive chlorine radicals (Cl • ) that affect the air quality by generating secondary air pollutants. Recent studies have shown road salt aerosols and saline snowpack to be sources of ClNO 2 in the wintertime urban environment; however, the quantitative contributions of each chloride source are not known. In this study, we examine the vertically resolved contributions of aerosol particles and saline snowpack as sources of ClNO 2 by using an observationally constrained snow–atmosphere coupled one-dimensional model applied to wintertime Kalamazoo, Michigan, U.S. Model simulations show that ClNO 2 emitted from urban snowpack can be vertically transported throughout the entire atmospheric boundary layer and can be a significant source of ClNO 2 , contributing up to ~60% of the ClNO 2 budget near the surface. Here, modeled snowpack ClNO 2 emission rates were 6 (±7) times higher than the observationally derived emission rates, suggesting that not all snow chloride is available for reaction. ClNO 2 production from both aerosol particles and snow emissions are required to best simulate the observed surface-level ClNO 2 . Using the traditional bulk parameterization for ClNO 2 produced from particles significantly overestimated ClNO 2 due to the assumption of having equivalent dinitrogen pentoxide (N 2 O 5 ) uptake and chloride availability for the entire particle population. In comparison, the chemically resolved surface area-based parameterization slightly underestimated the observations and had uncertainties deriving from ClNO 2 production from residential wood burning particles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

TWP-ICE Global Atmospheric Model Intercomparison: Convection Responsiveness and Resolution Impact

Results are presented from an intercomparison of global atmospheric model (GAM) simulations of tropical convection during the Tropical Warm Pool-International Cloud Experiment (TWP-ICE). The distinct cloud properties, precipitation, radiation, and vertical diabatic heating profiles associated with three different monsoon regimes (wet, dry, and break) from available observations are used to evaluate 9 GAM forecasts initialized daily from realistic global analyses. All models well captured the evolution of large-scale circulation and the thermodynamic fields, but cloud properties differed substantially among models. For example, liquid water path and ice water path differed by up to two orders of magnitude. Compared with the relatively well simulated top-heavy heating structures during the wet and break period, most models had difficulty in depicting the bottom-heavy heating profiles associated with cumulus congestus. The best performing models during this period were the ones whose convection scheme was most responsive to the free tropospheric humidity. Compared with the large impact of cloud and convective parameterizations on model cloud and precipitation characteristics, resolution has relatively minor impact on simulated cloud properties. However, one feature that was influence by the resolution study in several models was the diurnal cycle of precipitation. Peaking at a different time from convective precipitation, large-scale precipitation generally increases in high resolution forecasts and modulates the total precipitation diurnal cycle. Overall, the study emphasizes the importance of more environmental responsive convective parameterizations to capture various types of convection and the substantial diversity among large-scale cloud and precipitation schemes in current GAMs. This experiment has also demonstrated itself to be a very useful testbed for those developing cloud and convection schemes in these models.

58 GEOSCIENCES↗

Hydrological parameter estimations from a conservative tracer test with variable-density effects at the Boise Hydrogeophysical Research Site

Reliable predictions of groundwater flow and solute transport require an estimation of the detailed distribution of the parameters (e.g., hydraulic conductivity, effective porosity) controlling these processes. However, such parameters are difficult to estimate because of the inaccessibility and complexity of the subsurface. In this regard, developments in parameter estimation techniques and investigations of field experiments are still challenging and necessary to improve our understanding and the prediction of hydrological processes. Here we analyze a conservative tracer test conducted at the Boise Hydrogeophysical Research Site in 2001 in a heterogeneous unconfined fluvial aquifer. Some relevant characteristics of this test include: variable-density (sinking) effects because of the injection concentration of the bromide tracer, the relatively small size of the experiment, and the availability of various sources of geophysical and hydrological information. The information contained in this experiment is evaluated through several parameter estimation approaches, including a grid-search-based strategy, stochastic simulation of hydrological property distributions, and deterministic inversion using regularization and pilot-point techniques. Doing this allows us to investigate hydraulic conductivity and effective porosity distributions and to compare the effects of assumptions from several methods and parameterizations. Our results provide new insights into the understanding of variabledensity transport processes and the hydrological relevance of incorporating various sources of information in parameter estimation approaches. Among others, the variable-density effect and the effective porosity distribution, as well as their coupling with the hydraulic conductivity structure, are seen to be significant in the transport process. The results also show that assumed prior information can strongly influence the estimated distributions of hydrological properties.

54 ENVIRONMENTAL SCIENCES↗

The Roles of Cloud Drop Effective Radius and LWP in Determining Rain Properties in Marine Stratocumulus

Numerical simulations described in previous studies showed that adding cloud condensation nuclei to marine stratocumulus can prevent their breakup from closed into open cells. Additional analyses of the same simulations show that the suppression of rain is well described in terms of cloud drop effective radius (re). Rain is initiated when re near cloud top is around 12-14 um. Cloud water starts to get depleted when column-maximum rain intensity (Rmax) exceeds 0.1 mm h-1. This happens when cloud-top re reaches 14 um. Rmax is mostly less than 0.1 mm h-1 at re<14 um, regardless of the cloud water path, but increases rapidly when re exceeds 14 um. This is in agreement with recent aircraft observations and theoretical observations in convective clouds so that the mechanism is not limited to describing marine stratocumulus. These results support the hypothesis that the onset of significant precipitation is determined by the number of nucleated cloud drops and the height (H) above cloud base within the cloud that is required for cloud drops to reach re of 14 um. In turn, this can explain the conditions for initiation of significant drizzle and opening of closed cells providing the basis for a simple parameterization for GCMs that unifies the representation of both precipitating and non-precipitating clouds as well as the transition between them. Furthermore, satellite global observations of cloud depth (from base to top), and cloud top re can be used to derive and validate this parameterization.

54 ENVIRONMENTAL SCIENCES↗

Analysis of Cloud-resolving Simulations of a Tropical Mesoscale Convective System Observed during TWP-ICE: Vertical Fluxes and Draft Properties in Convective and Stratiform Regions

We analyze three cloud-resolving model simulations of a strong convective event observed during the TWP-ICE campaign, differing in dynamical core, microphysical scheme or both. Based on simulated and observed radar reflectivity, simulations roughly reproduce observed convective and stratiform precipitating areas. To identify the characteristics of convective and stratiform drafts that are difficult to observe but relevant to climate model parameterization, independent vertical wind speed thresholds are calculated to capture 90% of total convective and stratiform updraft and downdraft mass fluxes. Convective updrafts are fairly consistent across simulations (likely owing to fixed large-scale forcings and surface conditions), except that hydrometeor loadings differ substantially. Convective downdraft and stratiform updraft and downdraft mass fluxes vary notably below the melting level, but share similar vertically uniform draft velocities despite differing hydrometeor loadings. All identified convective and stratiform downdrafts contain precipitation below ~10 km and nearly all updrafts are cloudy above the melting level. Cold pool properties diverge substantially in a manner that is consistent with convective downdraft mass flux differences below the melting level. Despite differences in hydrometeor loadings and cold pool properties, convective updraft and downdraft mass fluxes are linearly correlated with convective area, the ratio of ice in downdrafts to that in updrafts is ~0.5 independent of species, and the ratio of downdraft to updraft mass flux is ~0.5-0.6, which may represent a minimum evaporation efficiency under moist conditions. Hydrometeor loading in stratiform regions is found to be a fraction of hydrometeor loading in convective regions that ranges from ~10% (graupel) to ~90% (cloud ice). These findings may lead to improved convection parameterizations.

analysis↗

Climate Impacts of Ice Nucleation

[1] Several different ice nucleation parameterizations in two different General Circulation Models (GCMs) are used to understand the effects of ice nucleation on the mean climate state, and the Aerosol Indirect Effects (AIE) of cirrus clouds on climate. Simulations have a range of ice microphysical states that are consistent with the spread of observations, but many simulations have higher present-day ice crystal number concentrations than in-situ observations. These different states result from different parameterizations of ice cloud nucleation processes, and feature different balances of homogeneous and heterogeneous nucleation. Black carbon aerosols have a small (-0.06 Wm -2 ) and not statistically significant AIE when included as ice nuclei, for nucleation efficiencies within the range of laboratory measurements. Indirect effects of anthropogenic aerosols on cirrus clouds occur as a consequence of increasing anthropogenic sulfur emissions with different mechanisms important in different models. In one model this is due to increases in homogeneous nucleation fraction, and in the other due to increases in heterogeneous nucleation with coated dust. The magnitude of the effect is the same however. The resulting ice AIE does not seem strongly dependent on the balance between homogeneous and heterogeneous ice nucleation. Regional effects can reach several Wm -2 . Indirect effects are slightly larger for those states with less homogeneous nucleation and lower ice number concentration in the base state. The total ice AIE is estimated at 0.27 ± 0.10 Wm -2 (1σ uncertainty). Finally, this represents a 20% offset of the simulated total shortwave AIE for ice and liquid clouds of -1.6 Wm -2 .

climate↗

Objective Quantification of Convective Clustering Observed During the AMIE/DYNAMO Two-Day Rain Episodes

One critical bottleneck in developing and evaluating ways to represent the mesoscale organization of convection in cumulus parameterization schemes is that there is no single accepted method of objectively quantifying the degree of convective organization or clustering from observations. This study addresses this need using high-quality S-PolKa radar data from the Atmospheric Radiation Measurement Madden-Julian Oscillation Investigation Experiment/Dynamics of the Madden-Julian Oscillation (AMIE/DYNAMO) field campaign. We first identify convective elements (contiguous convective echoes [CCEs]) from radar reflectivity observations using the rain type classification algorithm of Powell et al. (2016, https://doi.org/10.1175/JTECH-D-15-0135.1). Then we apply scalar clustering metrics, including the organization index (I org ) of Tompkins and Semie, to the radar CCEs to test their ability of quantifying convective clustering during the observed two-day rain episodes. Our results show two distinct phases of convective clustering during the two-day rain episodes, with each phase covering about 10 hr before (Phase 1) and after (Phase 2) the time of peak rain rate. In Phase 1 clustering, the number of CCEs increases and convective cells cluster as new cells form preferentially near existing convective entities. The number of CCEs decreases as the environment stabilizes in Phase 2 clustering, during which already clustered cells with associated stratiform clouds are preferred over the isolated ones. I org is capable of capturing convective clustering in both phases. The possible mechanisms for convective clustering are discussed, including cold pool-updraft feedback, moisture-convection interaction, and mesoscale circulations. Our results suggest that parameterizations of convective organization should represent the feedback processes that are responsible for the convective clustering during both phases.

54 ENVIRONMENTAL SCIENCES↗

Improving Land Surface Temperature Simulation in CoLM Over the Tibetan Plateau through Fractional Vegetation Cover Derived from a Remotely Sensed Clumping Index and Model-Simulated Leaf Area Index

Parameterizations of fractional vegetation cover (FVC) in land surface models have important effects on simulations of surface energy budget, especially in arid and semiarid regions. This study uses a FVC scheme in which FVC is derived from leaf area index and a remotely sensed clumping index. The performance of the new scheme (SMFVC) is evaluated against Moderate Resolution Imaging Spectroradiometer (MODIS) land surface temperature (LST) and in situ soil temperature observations, together with two other FVC schemes, a general FVC scheme (CTL) based on land cover map and a climatology-based FVC scheme (RSFVC) that uses long-term remotely sensed Normalized Difference Vegetation Index of MODIS. The three FVC schemes were implemented in the Common Land Model (CoLM) and applied in the Tibetan Plateau using the same forcing data and default parameters. Our results demonstrate that FVC schemes have significant influence on the CoLM performance: (1) the RSFVC and SMFVC schemes significantly reduce the LST biases found in CTL, particularly in grassland and during summer; (2) soil temperature evaluation by in situ observations from three networks on the Tibetan Plateau corroborates the LST results; and (3) the improvements are mainly related to representing temporal (seasonal) variability and subgrid heterogeneity of FVC, which improves surface albedo and surface energy balance. In other words, by including more vegetation characteristics, such as using a clumping index, land surface models may better simulate surface vegetation condition and further better represent the land surface energy budget over the Tibetan Plateau.

Li, Chengwei↗

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↗

Improved Spatiotemporal Representativeness and Bias Reduction of Satellite-Based Evapotranspiration Retrievals via Use of In Situ Meteorology and Constrained Canopy Surface Resistance

Evapotranspiration (ET) is a key component of the atmospheric and terrestrial water and energy budgets. Satellite-based vegetation index approaches have used remotely-sensed vegetation and reanalysis meteorological properties with surface energy balance models to estimate global ET (MOD16 ET). We reconstructed satellite retrievals using in situ meteorology (Argonne-ET) and evaluated them using a dense network of surface turbulent heat flux measurements. Argonne-ET resolves spatial heterogeneity of ET across the U.S. Southern Great Plains that is not well characterized by MOD16; MOD16 ET exhibits spatial autocorrelation (Moran's I (MI) significantly > 0 in May – Aug.), whereas in situ and Argonne-ET exhibit characteristics of a random spatial process (MI not significantly different from 0). The skill in resolving ET temporal variability is not significantly different between MOD16 and Argonne-ET (correlation coefficient = 0.75 and 0.72, respectively). However, the root mean square errors were significantly lower for Argonne-ET (36 W m -2 ) than MOD16 (43 W m -2 ), and MOD16 exhibits substantial bias in annual ET relative to in situ measurements (-38%). This is attributed to overestimation in the dry canopy surface resistance (rs) parameterization. Using rs constrained to the range of typical measured values, Argonne-ET substantially reduces the bias in the annual ET (+1%). The improved ET estimates are critical for regional water budget analyses. The methodology presented herein also demonstrates the ability to retrieve high temporally resolved (30-min; cf. 8-day MOD16) ET that can be used for development of processed-based diagnostics of model biases, and to elucidate avenues to improve ET model parameterizations.

54 ENVIRONMENTAL SCIENCES↗

The Impact of Surface Heterogeneities and Land‐Atmosphere Interactions on Shallow Clouds Over ARM SGP Site

Abstract Continental shallow clouds are an important component of the land‐atmosphere coupled climate system because of their role in modulating energy and water budgets. The parameterization of these clouds in climate models presents a significant challenge. We evaluate the potential impact of subgrid‐scale (for climate models) land‐atmosphere interactions on shallow clouds using nested large‐eddy simulations (LESs). We compare LESs that allow for land‐atmosphere interactions to ones that artificially suppress all or part of them by smoothing out surface heterogeneities in surface heat, moisture and radiation fluxes within the domain. The LES domains are coupled to an interactive land surface model and nested inside a mesoscale domain. Three summertime shallow convection cases over the Department of Energy Atmospheric Radiation Measurement Facility's Southern Great Plains site are examined. A consistent increase in cloud water content is observed in all cases when surface heterogeneities and the land‐atmosphere interactions they induce in the LES domain are removed. By comparing experiments where only surface flux heterogeneities induced by cloud shadows are removed to ones where all surface flux heterogeneities are removed, we find that cloud shading‐induced surface heterogeneities can have a larger impact on shallow convection than static land surface heterogeneities in our cases. Furthermore, the impact of cloud shading on cloud water content and cloud size is found to vary significantly with the solar incidence angle. Our results suggest that the impact of cloud shadows needs to be considered when parameterizing land‐atmosphere interactions in the presence of shallow clouds for regional or global climate models.

54 ENVIRONMENTAL SCIENCES↗

Development and Evaluation of an Explicit Treatment of Aerosol Processes at Cloud Scale Within a Multi‐Scale Modeling Framework (MMF)

Abstract Modeling the aerosol lifecycle in traditional global climate models (GCM) is challenging for a variety of reasons, not the least of which is the coarse grid. The multiscale modeling framework (MMF), in which a cloud resolving model replaces conventional parameterizations of cloud processes within each GCM grid column, provides a promising framework to address this challenge. Here we develop a new version of MMF that for the first time treats aerosol processes at cloud scale to improve the aerosol‐cloud interaction representation in the model. We demonstrate that the model with the explicit aerosol treatments shows significant improvements of many aspects of the simulated aerosols compared to the previous version of MMF with aerosols parameterized at the GCM grid scale. The explicit aerosol treatments produce a significant increase of the column burdens of black carbon (BC), primary organic aerosol, and sulfate by up to 40% in many remote regions, a decrease of the sea‐salt aerosol burdens by 40% in remote regions. These differences are caused by the differences in aerosol convective transport and wet removal between these two models. The new model also shows reduced bias of BC surface concentration in North America and BC vertical profiles in the high latitudes. However, the biased‐high BC concentrations in the upper troposphere over the remote Pacific regions remain, requiring further improvements on other process representations (e.g., secondary activation neglected in the model).

54 ENVIRONMENTAL SCIENCES↗

Structure and Evolution of Mesoscale Convective Systems: Sensitivity to Cloud Microphysics in Convection‐Permitting Simulations Over the United States

Abstract Regional climate simulations over the continental United States were conducted for the 2011 warm season using the Weather Research and Forecasting model at convection‐permitting resolution (4 km) with two commonly used microphysics parameterizations (Thompson and Morrison). Sensitivities of the simulated mesoscale convective system (MCS) properties and feedbacks to large‐scale environments are systematically examined against high‐resolution geostationary satellite and 3‐D mosaic radar observations. MCS precipitation including precipitation amount, diurnal cycle, and distribution of hourly precipitation intensity are reasonably captured by the two simulations despite significant differences in their simulated MCS properties. In general, the Thompson simulation produces better agreement with observations for MCS upper level cloud shield and precipitation area, convective feature horizontal and vertical extents, and partitioning between convective and stratiform precipitation. More importantly, Thompson simulates more stratiform rainfall, which agrees better with observations and results in top‐heavier heating profiles from robust MCSs compared to Morrison. A stronger dynamical feedback to the large‐scale environment is therefore seen in Thompson, wherein an enhanced mesoscale vortex behind the MCS strengthens the synoptic‐scale trough and promotes advection of cool and dry air into the rear of the MCS region. The latter prolongs the MCS lifetimes in the Thompson relative to the Morrison simulations. Hence, different treatment of cloud microphysics not only alters MCS convective‐scale dynamics but also has significant impacts on their macrophysical properties such as lifetime and precipitation. As long‐lived MCSs produced 2–3 times the amount of rainfall compared to short‐lived ones, cloud microphysics parameterizations have profound impact in simulating extreme precipitation and the hydrologic cycle.

54 ENVIRONMENTAL SCIENCES↗

Insensitivity of the Cloud Response to Surface Warming Under Radical Changes to Boundary Layer Turbulence and Cloud Microphysics: Results From the Ultraparameterized CAM

We study the cloud response to a +4K surface warming in a new multi- scale climate model that uses enough interior resolution to begin explicitly resolving boundary layer turbulence (i.e. Ultra-Parameterization or UP), and the associated low clouds. UP's predictions are compared against those from standard Super-Parameterization (SP). The mean global cloud feedback turns out to be remarkably neutral across all of our simulations, despite some radical changes in both cloud microphysical parameter settings and cloud resolving model grid resolution. The overall cloud response to warming is the positive low-cloud feedbacks over land, negative feedbacks at high latitudes, and weak feedbacks over the low-latitude oceans. The most distinct effects of UP appear to be the result of tuning decisions impacting high latitude cloud feedback. UP's microphysics are tuned to optimize the model present-day, top-of-atmosphere radiation fluxes against CERES observations, by lowering the cloud ice-liquid phase-shift temperature ramp, adjusting the ice/liquid autoconversion rate, and increasing the ice fall speed. This reduces baseline high-latitude low cloud amounts, damping the phase change cloud feedback at high latitudes, leading to a slightly more positive global cloud feedback compared to SP. A sensitivity test that isolates these microphysical impacts from UP's grid resolution conrms that the microphysical settings are mostly responsible for the dierences between SP and UP cloud feedback.

54 ENVIRONMENTAL SCIENCES↗

The Vertical Momentum Budget of Shallow Cumulus Convection: Insights From a Lagrangian Perspective

A Lagrangian Particle Dispersion Model is embedded into large eddy simulations to diagnose the responses of shallow cumulus convection to a small-amplitude large-scale temperature perturbation. The Lagrangian framework allows for a decomposition of the vertical momentum budget and diagnosis of the forces that regulate cloudy updrafts. The results are used to shed light on the parameterization of vertical velocity in convective schemes, where the treatment of the effects of entrainment as well as buoyancy-induced and mechanically induced pressure gradients remains highly uncertain. We show that both buoyancy-induced and mechanically induced pressure gradients are important for the vertical momentum budget of cloudy updrafts, whereas the entrainment dilution term is relatively less important. Based on the analysis of the dominant force balance, we propose a simple model to derive the perturbation pressure gradient forces. We further illustrate that the effective buoyancy and dynamic perturbation pressure can be approximated to a good extent using a simple cylindrical updraft model given the cloud radius. This finding has the potential for improving the parameterization of vertical velocity in convective schemes and the development of a unified scheme for cumulus convection.

54 ENVIRONMENTAL SCIENCES↗

Prognostic Power of Extreme Rainfall Scaling Formulas Across Space and Time Scales

Some studies documenting changes in extreme precipitation use scaling formulas to approximate the large percentiles of the rainfall distribution from average dynamical and thermodynamical variables, called predictors. Here we instead assess the performance of these formulas as approximations to the rain rates in individual events. We evaluate the accuracies of the scaling relationships as functions of spatial and temporal scales by analyzing tropical rainfall in a superparameterized model. Relationships using full vertical profiles of the predictors are more accurate than those using their values at specific vertical levels because they better characterize the specific dynamics of each event. Both types of scaling relationships perform well over a range of length scales from 200 to 2,000 km and time scales from an hour to a week, and their precision is higher in the case of simulations with superparameterization than with parameterized convection. Uncertainties emerging from the local use of the scaling relationships suggest that they may only characterize the intensification of individual extremes for a warming of 4–5 K or larger. Finally, we argue that these formulas can be used to reconstruct the tail of the rainfall distribution directly from its predictors without prior information on P. While scalings have been used as diagnostic equations conditioned on the occurrence of extreme rainfall, they are actually able to mimic the prognostic behavior of climate model parameterizations on a variety of scales when estimating the intensity, frequency, and spatial patterns of extremes.

54 ENVIRONMENTAL SCIENCES↗

The Community Land Model Version 5: Description of New Features, Benchmarking, and Impact of Forcing Uncertainty

The Community Land Model (CLM) is the land component of the Community Earth System Model (CESM) and is used in several global and regional modeling systems. In this paper, we introduce model developments included in CLM version 5 (CLM5), which is the default land component for CESM2. We assess an ensemble of simulations, including prescribed and prognostic vegetation state, multiple forcing data sets, and CLM4, CLM4.5, and CLM5, against a range of metrics including from the International Land Model Benchmarking (ILAMBv2) package. CLM5 includes new and updated processes and parameterizations: (1) dynamic land units, (2) updated parameterizations and structure for hydrology and snow (spatially explicit soil depth, dry surface layer, revised groundwater scheme, revised canopy interception and canopy snow processes, updated fresh snow density, simple firn model, and Model for Scale Adaptive River Transport), (3) plant hydraulics and hydraulic redistribution, (4) revised nitrogen cycling (flexible leaf stoichiometry, leaf N optimization for photosynthesis, and carbon costs for plant nitrogen uptake), (5) global crop model with six crop types and time-evolving irrigated areas and fertilization rates, (6) updated urban building energy, (7) carbon isotopes, and (8) updated stomatal physiology. New optional features include demographically structured dynamic vegetation model (Functionally Assembled Terrestrial Ecosystem Simulator), ozone damage to plants, and fire trace gas emissions coupling to the atmosphere. Conclusive establishment of improvement or degradation of individual variables or metrics is challenged by forcing uncertainty, parametric uncertainty, and model structural complexity, but the multivariate metrics presented here suggest a general broad improvement from CLM4 to CLM5.

58 GEOSCIENCES↗

Improving Precipitation Estimation Using Convolutional Neural Network

Precipitation process is generally considered to be poorly represented in numerical weather/climate models. Statistical downscaling (SD) methods, which relate precipitation with model resolved dynamics, often provide more accurate precipitation estimates compared to model's raw precipitation products. Herein we introduce the convolutional neural network model to foster this aspect of SD for daily precipitation prediction. Specifically, we restrict the predictors to the variables that are directly resolved by discretizing the atmospheric dynamics equations. In this sense, our model works as an alternative to the existing precipitation-related parameterization schemes for numerical precipitation estimation. We train the model to learn precipitation-related dynamical features from the surrounding dynamical fields by optimizing a hierarchical set of spatial convolution kernels. We test the model at 14 geogrid points across the contiguous United States. Results show that provided with enough data, precipitation estimates from the convolutional neural network model outperform the reanalysis precipitation products, as well as SD products using linear regression, nearest neighbor, random forest, or fully connected deep neural network. Evaluation for the test set suggests that the improvements can be seamlessly transferred to numerical weather modeling for improving precipitation prediction. Based on the default network, we examine the impact of the network architectures on model performance. Also, we offer simple visualization and analyzing approaches to interpret the models and their results. Our study contributes to the following two aspects: First, we offer a novel approach to enhance numerical precipitation estimation; second, the proposed model provides important implications for improving precipitation-related parameterization schemes using a data-driven approach.

54 ENVIRONMENTAL SCIENCES↗