Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “CLIMATE MODELS”

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 289 records · Page 16

Future changes in regional precipitation simulated by a half-degree coupled climate model: Sensitivity to horizontal resolution

The global fully coupled half-degree Community Climate System Model Version 4 (CCSM4) was integrated for a suite of climate change ensemble simulations including five historical runs, five Representative Concentration Pathway 8.5 [RCP8.5) runs, and a long Pre-Industrial control run. This study focuses on precipitation at regional scales and its sensitivity to horizontal resolution. The half-degree historical CCSM4 simulations are compared to observations, where relevant, and to the standard 1° CCSM4. Both the halfdegree and 1° resolutions are coupled to a nominal 1° ocean. North American and South Asian/Indian monsoon regimes are highlighted because these regimes demonstrate improvements due to higher resolution, primarily because of better-resolved topography. Agriculturally sensitive areas are analyzed and include Southwest, Central, and Southeast U.S., Southern Europe, and Australia. Both mean and extreme precipitation is discussed for convective and large-scale precipitation processes. Convective precipitation tends to decrease with increasing resolution and large-scale precipitation tends to increase. Improvements for the half-degree agricultural regions can be found for mean and extreme precipitation in the Southeast U.S., Southern Europe, and Australian regions. Climate change responses differ between the model resolutions for the U.S. Southwest/Central regions and are seasonally dependent in the Southeast and Australian regions. Both resolutions project a clear drying signal across Southern Europe due to increased greenhouse warming. As a result, differences between resolutions tied to the representation of convective and large-scale precipitation play an important role in the character of the climate change and depend on regional influences.

58 GEOSCIENCES↗

Global Temperature Responses to Large Tropical Volcanic Eruptions in Paleo Data Assimilation Products and Climate Model Simulations Over the Last Millennium

Abstract Large volcanic eruptions are one of the dominant perturbations to global and regional atmospheric temperatures on timescales of years to decades. Discrepancies remain, however, in the estimated magnitude and persistence of the surface temperature cooling caused by volcanic eruptions, as characterized by paleoclimatic proxies and climate models. We investigate these discrepancies in the context of large tropical eruptions over the Last Millennium using two state‐of‐the‐art data assimilation products, the Paleo Hydrodynamics Data Assimilation product (PHYDA) and the Last Millennium Reanalysis (LMR), and simulations from the National Center for Atmospheric Research Community Earth System Model‐Last Millennium Ensemble (NCAR CESM‐LME). We find that PHYDA and LMR estimate mean global and hemispheric cooling that is similar in magnitude and persistence once effects from eruptions occurring in short succession are removed. The estimates also compare well to Northern‐Hemisphere reconstructions based solely or partially on tree‐ring density, which have been proposed as the most accurate proxy estimates of surface cooling due to volcanism. All proxy‐based estimates also agree well with the magnitude of the mean cooling simulated by the CESM‐LME. Differences remain, however, in the spatial patterns of the temperature responses in the PHYDA, LMR, and the CESM‐LME. The duration of cooling anomalies also persists for several years longer in the PHYDA and LMR relative to the CESM‐LME. Our results demonstrate progress in resolving discrepancies between proxy‐ and model‐based estimates of temperature responses to volcanism, but also indicate these estimates must be further reconciled to better characterize the risks of future volcanic eruptions.

Tejedor, E.↗

ARM Data-Oriented Metrics and Diagnostics Package for Climate Model Evaluation

A Python-based metrics and diagnostics package is currently being developed by the U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) user facility Infrastructure Team at Lawrence Livermore National Laboratory to facilitate the use of long-term, high-frequency measurements from ARM in evaluating the regional climate simulation of clouds, radiation, and precipitation. This metrics and diagnostics package computes climatological means of targeted climate model simulation and generates tables and plots for comparing the model simulation with ARM observational data.

54 ENVIRONMENTAL SCIENCES↗

Assessment of Storm-Associated Precipitation and Its Extremes Using Observational Data Sets and Climate Model Short-Range Hindcasts

Heavy precipitation, often associated with weather phenomena such as tropical cyclones, extratropical cyclones (ETCs), atmospheric rivers (ARs), and mesoscale convective systems (MCSs), can cause significant socio-economic loss. Here, in this study, we apply atmospheric feature trackers to quantify the contributions of these storm types in observational data sets and climate model short-range hindcasts. We generate a global hourly storm data set at 0.25° spatial resolution covering 2006–2020, based on the tracking results from TempestExtremes and Python FLEXible object TRacKeR. Our analyses show that these four storm types account for 67% of global annual mean precipitation and 82% of top 1% precipitation extremes, with MCSs mainly over the tropics, and ARs and ETCs over the midlatitudes. The percentage of precipitation contributions from these storms also show strong seasonality over many geographical locations. We further apply the tracking results to the Energy Exascale Earth System Model (E3SM) short-range hindcasts and evaluate how well these storms are simulated. The evaluation show that E3SM, with ∼1° resolution, significantly underestimates storm-associated precipitation totals and extremes, especially for MCSs in the tropics. Our analysis also suggests that model fails to capture the correct mean diurnal phases and amplitude of MCS precipitation. This phenomenon-based approach provides a better understanding of precipitation characteristics and can lead to enhanced model evaluation by revealing underlying problems in model physics related to precipitation processes associated with the heavy-precipitating storms.

54 ENVIRONMENTAL SCIENCES↗

Effects of Urbanization in China on the East Asian Summer Monsoon as Revealed by Two Global Climate Models

Abstract The East Asian summer monsoon (EASM), vulnerable to anthropogenic activities, exerts profound impacts on precipitation in East Asia. In recent decades, China has been experiencing rapid urbanization. It has not yet reached a consensus on its influences on EASM. Here, we introduce satellite‐derived urban impervious areas of China during 1985–2017 into two global climate models with contrasting biases of climatological EASM simulation. Both models consistently show that urbanization weakens EASM by increasing surface friction drag. The increased surface roughness reduces climatological southerly winds in the concentrated urbanization area. As a result, moisture convergence occurs south of the concentrated urbanization area with ascending motion, while moisture divergence with descending motion appears to the north, resulting in the precipitation anomalies of “southern flood and northern drought.” The dominant role of urban dynamic forcing in regulating EASM is confirmed by reducing urban roughness to the level of croplands in sensitivity experiments. As urban dynamic forcing is largely mitigated, the urbanization impact mainly driven by urban thermal forcing conversely strengthens EASM due to enhanced land‐sea thermal contrast via urbanization‐induced warming. Spatial patterns and magnitudes of the urbanization‐induced EASM weakening vary between the two models due to disparities in their simulations of urban thermodynamic forcing and atmospheric adjustments. These results underscore the robust weakening effect of urbanization on EASM by urban friction drag and document uncertainties in the related spatial pattern and magnitude resulting from model discrepancies.

Meteorology & Atmospheric Sciences↗

Developing and Testing a Novel Stochastic Ice Microphysics Parameterization for Cloud and Climate Models Using ARM Field Campaign Data (Final Progress Report)

The major goals of this project were: 1) to use field campaign measurements from DOE’s Atmospheric Radiation Measurement (ARM) program to characterize variability of important parameters describing properties of ice particles in the atmosphere; 2) based on this observational analysis, to develop a parameterization scheme for weather and climate models that stochastically varies these parameters, and implement the new scheme into a weather model called the Weather Research and Forecasting model (WRF); 3) to use WRF coupled with the new stochastic scheme to simulate ARM field campaign thunderstorm cases and analyze how accounting for this parameter variability affects the model simulations. This work was performed jointly between the National Center for Atmospheric Research, University of Oklahoma, and University of Utah. To accomplish these goals, we extended an approach previously developed to characterize the variability in the size distribution of ice particles to parameters that are explicitly represented in models (i.e., relationships between ice particle mass and size, and between particle fall velocity and size). Our project was, to our knowledge, the first to apply observationally-constrained estimates of this parameter variability describing mass-size and fall velocity-size in a modeling framework. Our results showed efficacy of the approach, evaluated using ARM observations. Similarly, to our knowledge, work in this project was the first to propose and evaluate in detail a stochastic approach for unresolved turbulent mixing in high-resolution model simulations against detailed, benchmark large eddy simulations and ARM observations. Results showed some promising behavior, particularly with increased mixing and dilution of air in thunderstorm cores with surrounding environmental air, bringing the stochastic simulations closer to the benchmark large eddy simulations; however, results were somewhat degraded using stochastic mixing compared to observations from the AMIE/DYNAMO field campaign. This project also further refined and applied a modeling methodology called “piggybacking” that can robustly separate dynamical and thermodynamic impacts of model changes, and comparison studies of different models based on cases developed from ARM observations. Finally, this project directly supported three graduate students who completed their PhDs as well as a postdoctoral research fellow.

54 ENVIRONMENTAL SCIENCES↗

Toward a new parameterization of hydraulic conductivity in climate models: Simulation of rapid groundwater fluctuations in N orthern C alifornia

Abstract Preferential flow through weathered bedrock leads to rapid rise of the water table after the first rainstorms and significant water storage (also known as “rock moisture”) in the fractures. We present a new parameterization of hydraulic conductivity that captures the preferential flow and is easy to implement in global climate models. To mimic the naturally varying heterogeneity with depth in the subsurface, the model represents the hydraulic conductivity as a product of the effective saturation and a background hydraulic conductivity K bkg , drawn from a lognormal distribution. The mean of the background K bkg decreases monotonically with depth, while its variance reduces with the effective saturation. Model parameters are derived by assimilating into Richards’ equation 6 years of 30 min observations of precipitation (mm) and water table depths (m), from seven wells along a steep hillslope in the Eel River watershed in Northern California. The results show that the observed rapid penetration of precipitation and the fast rise of the water table from the well locations, after the first winter rains, are well captured with the new stochastic approach in contrast to the standard van Genuchten model of hydraulic conductivity, which requires significantly higher levels of saturated soils to produce the same results. “Rock moisture,” the moisture between the soil mantle and the water table, comprises 30% of the moisture because of the great depth of the weathered bedrock layer and could be a potential source of moisture to sustain trees through extended dry periods. Furthermore, storage of moisture in the soil mantle is smaller, implying less surface runoff and less evaporation, with the proposed new model.

58 GEOSCIENCES↗

Clouds and Radiation Processes in Regional Climate Models Evaluated Using Observations Over the Ice‐free Arctic Ocean

Abstract The presence of clouds in the Arctic regulates the surface energy budget (SEB) over the sea‐ice surface and the ice‐free ocean. Following several previous field campaigns, the cloud‐radiation relationship, including cloud vertical structure and phase, has been elucidated; however, modeling of this relationship has matured slowly. In recognition of the recent decline in the Arctic sea‐ice extent, representation of the cloud system in numerical models should consider the effects of areas covered by sea ice and ice‐free areas. Using an in situ stationary meteorological observation data set obtained over the ice‐free Arctic Ocean by the Japanese Research Vessel Mirai (September 2014), coordinated evaluation of six regional climate models (RCMs) with nine model runs was performed by focusing on clouds and the SEB. The most remarkable findings were as follows: (1) reduced occurrence of unstable stratification with low‐level cloud water in all models in comparison to the observations, (2) significant differences in cloud water representations between single‐ and double‐moment cloud schemes, (3) extensive differences in partitioning of hydrometeors including solid/liquid precipitation, and (4) pronounced lower‐tropospheric air temperature biases. These issues are considered as the main sources of SEB uncertainty over ice‐free areas of the Arctic Ocean. The results from a coupled RCM imply that the SEB is constrained by both the atmosphere and the ocean (and sea ice) with considerable feedback. Coordinated improvement of both stand‐alone atmospheric and coupled RCMs would promote a more comprehensive and improved understanding of the Arctic air‐ice‐sea coupled system.

Inoue, Jun↗

Parameterizing Convective Organization Effects With a Moisture-PDF Approach in Climate Models: Concept and a Regional Case Simulation

We propose a parameterization scheme of convective organization effects based on a moisture-distribution approach, which can reflect aggregation of convective cells within a model grid as well as the interaction between convection and spatial heterogeneity in free-troposphere moisture. With this concept, convective cells in an aggregated state are surrounded by air that is moister than the grid-mean condition, which provides a “shielding” effect favorable for deeper convective updrafts. Such effects are represented quantitatively via utilizing a schematic diagram about the geometric interpretation of convective cluster, dry area, and their contact area, in which both the convective cluster size and dry area fraction increase during the development of convection. Our sensitivity analyses indicate that the new scheme performs well in capturing the mean precipitation features. Excluding the convective organization effects in the model leads to a considerable reduction in the simulated precipitation magnitude. Notably, there exists an inverse relationship between background mean moisture condition and sub-grid moisture variability for a given rain rate, suggesting that with a high spatial heterogeneity in free-troposphere moisture, convection can still develop and maintain its strength under relatively dry background conditions, consistent with the convection-permitting model simulation. Overall, our new parameterization of convective organization effects can successfully reproduce the relationship between precipitation and sub-grid moisture variability, which is a missing element in traditional convection parameterization schemes but important for the simulations of precipitation variability at various scales in climate models.

54 ENVIRONMENTAL SCIENCES↗

Accurate assessment of land–atmosphere coupling in climate models requires high-frequency data output

Land–atmosphere (L–A) interactions are important for understanding convective processes, climate feedbacks, the development and perpetuation of droughts, heatwaves, pluvials, and other land-centered climate anomalies. Local L–A coupling (LoCo) metrics capture relevant L–A processes, highlighting the impact of soil and vegetation states on surface flux partitioning and the impact of surface fluxes on boundary layer (BL) growth and development and the entrainment of air above the BL. A primary goal of the Climate Process Team in the Coupling Land and Atmospheric Subgrid Parameterizations (CLASP) project is parameterizing and characterizing the impact of subgrid heterogeneity in global and regional Earth system models (ESMs) to improve the connection between land and atmospheric states and processes. A critical step in achieving that aim is the incorporation of L–A metrics, especially LoCo metrics, into climate model diagnostic process streams. However, because land–atmosphere interactions span timescales of minutes (e.g., turbulent fluxes), hours (e.g., BL growth and decay), days (e.g., soil moisture memory), and seasons (e.g., variability in behavioral regimes between soil moisture and latent heat flux), with multiple processes of interest happening in different geographic regions at different times of year, there is not a single metric that captures all the modes, means, and methods of interaction between the land and the atmosphere. And while monthly means of most of the LoCo-relevant variables are routinely saved from ESM simulations, data storage constraints typically preclude routine archival of the hourly data that would enable the calculation of all LoCo metrics. Here, we outline a reasonable data request that would allow for adequate characterization of sub-daily coupling processes between the land and the atmosphere, preserving enough sub-daily output to describe, analyze, and better understand L–A coupling in modern climate models. A secondary request involves embedding calculations within the models to determine mean properties in and above the BL to further improve characterization of model behavior. Higher-frequency model output will (i) allow for more direct comparison with observational field campaigns on process-relevant timescales, (ii) enable demonstration of inter-model spread in L–A coupling processes, and (iii) aid in targeted identification of sources of deficiencies and opportunities for improvement of the models.

54 ENVIRONMENTAL SCIENCES↗

Using the GoAmazon-CHUVA measurements to understand what causes the biases in the onset of the rainy season in Amazonia in climate models

The onset of the Amazon rainy season has a large temporal and spatial variability, and strong impact on aerosols, ecosystems, fire, carbon fluxes, dry season length, agriculture, and hydropower. Two major droughts occurred in the region in the last 20 years, and the dry season has increased in length by about one month. These events highlight the urgency for improving our understanding and capability to project the rainy season onset and drought variability. However, global climate models (CMIP3 and CMIP5) appear to underestimate the past variability, and also project virtually no future change of the onset of rainy season over the Amazon even when they are forced by strong greenhouse forcing under the RCP8.5 scenario. Why these models underestimate the variability of the rainy season onset, and whether this bias implies an underestimate of sensitivity of their dry season length to anthropogenic radiative forcing remain unclear.

54 ENVIRONMENTAL SCIENCES↗

Responses in the Subpolar North Atlantic in Two Climate Model Sensitivity Experiments with Increased Stratospheric Aerosols

Abstract The subpolar North Atlantic (SPNA) shows contrasting responses in two sensitivity experiments with increased stratospheric aerosols, offering insight into the physical processes that may impact the Atlantic meridional overturning circulation (AMOC) in a warmer climate. In one, the upper ocean becomes warm and salty, but in the other it becomes cold and fresh. The changes are accompanied by diverging AMOC responses. The first experiment strengthens the AMOC, opposing the weakening trend in the reference simulation. The second experiment shows a much smaller impact. Both simulations use the Community Earth System Model with the Whole Atmosphere Community Climate Model component (CESM-WACCM) but differ in model versions and stratospheric aerosol specifications. Despite both experiments using similar approaches to increase stratospheric aerosols to counteract the rising global temperature, the contrasting SPNA and AMOC responses indicate a considerable dependency on model physics, climate states, and model responses to forcings. This study focuses on examining the physical processes involved with the impact of stratospheric aerosols on the SPNA salinity changes and their potential connections with the AMOC and the Arctic. We find that in both cases, increased stratospheric aerosols act to enhance the SPNA upper-ocean salinity by reducing freshwater export from the Arctic, which is closely tied to the Arctic sea ice changes. The impact on AMOC is primarily through the thermal component of the surface buoyancy fluxes, with negligible contributions from the freshwater component. These experiments shed light on the physical processes that dictate the important connections between the SPNA, the Arctic, the AMOC, and their subsequent feedbacks on the climate system.

54 ENVIRONMENTAL SCIENCES↗

A climate model projection weighting scheme accounting for performance and interdependence: Model Projection Weighting Scheme

Uncertainties of climate projections are routinely assessed by considering simulations from different models. Observations are used to evaluate models, yet there is a debate about whether and how to explicitly weight model projections by agreement with observations. Here we present a straightforward weighting scheme that accounts both for the large differences in model performance and for model interdependencies, and we test reliability in a perfect model setup. We provide weighted multimodel projections of Arctic sea ice and temperature as a case study to demonstrate that, for some questions at least, it is meaningless to treat all models equally. The constrained ensemble shows reduced spread and a more rapid sea ice decline than the unweighted ensemble. We argue that the growing number of models with different characteristics and considerable interdependence finally justifies abandoning strict model democracy, and we provide guidance on when and how this can be achieved robustly.

54 ENVIRONMENTAL SCIENCES↗

Projected increase of the East Asian summer monsoon ( Meiyu ) in Taiwan by climate models with variable performance

Abstract The active phase of the East Asian summer monsoon (EASM) in Taiwan during May and June, known as Meiyu , produces substantial precipitation for water uses in all sectors of society. Following a companion study that analysed the historical increase in the Meiyu precipitation, the present study conducted model evaluation and diagnosis based on the EASM lifecycle over Taiwan. Higher and lower skill groups were identified from 17 Couple Model Intercomparison Project Phase 5 (CMIP5) models, with five models in each group. Despite the difference in model performance, both groups projected a substantial increase in the Meiyu precipitation over Taiwan. In the higher skill group, weak circulation changes and reduced low‐level convergence point to a synoptically unfavourable condition for precipitation. In the lower skill group, intensified low‐level southwesterly winds associated with a deepened upper level trough enhance moisture pooling. Thus, the projected increase in Meiyu precipitation will likely occur through the combined effects of (1) the extension of a strengthened North Pacific anticyclone enhancing southwesterlies; and (2) more systematically, the Clausius–Clapeyron relationship that increases precipitation intensity in a warmer climate. The overall increase in the Meiyu precipitation projected by climate models of variable performance supports the observed tendency toward more intense rainfall in Taiwan and puts its early June 2017 extreme precipitation events into perspective.

Tung, Yu‐Shiang↗

Climate models show colorado drying sooner and with greater certainty east of the continental divide

Many studies have examined the aridity of the Colorado River Basin and the possible impacts of climate change which could further strain already over-allocated water resources in the region. Fewer studies have examined the multiple Colorado Rocky Mountain headwater regions specifically. This is especially true of areas East of the Continental Divide, despite water originating there being critical to cities and agriculture in Eastern Colorado and further downstream. This paper explores and compares drying trends in the Eastern and Western Colorado Rocky Mountains using single-model initial-condition large ensembles from ten global climate models. The use of multiple models allows us to identify signals that are consistent across different physics parameterizations, model grids, and other model intricacies. The large ensembles also allow us to quantify the time of emergence of these climate change signals--that is, when did (or when will) the long term change due to anthropogenic greenhouse gasses exceed the internal variability of the climate system. Consistent with previous studies, we find evidence of drying on both sides of the Continental Divide. That drying is more pronounced, occurs sooner, and is more consistent across global climate models in the East, however, highlighting the region’s importance despite generally receiving less attention than the West.

Rugg, Allyson↗

Understanding the Varied Influence of Midlatitude Jet Position on Clouds and Cloud Radiative Effects in Observations and Global Climate Models

This study examines the dynamical mechanisms responsible for changes in midlatitude clouds and cloud radiative effects (CRE) that occur in conjunction with meridional shifts in the jet streams over the North Atlantic, North Pacific, and Southern Oceans. When the midlatitude jet shifts poleward, extratropical cyclones and their associated upward vertical velocity anomalies closely follow. As a result, a poleward jet shift contributes to a poleward shift in high-topped storm-track clouds and their associated longwave CRE. However, when the jet shifts poleward, downward vertical velocity anomalies increase equatorward of the jet, contributing to an enhancement of the boundary layer estimated inversion strength (EIS) and an increase in low cloud amount there. Because shortwave CRE depends on the reflection of solar radiation by clouds in all layers, the shortwave cooling effects of midlatitude clouds increase with both upward vertical velocity anomalies and positive EIS anomalies. Over midlatitude oceans where a poleward jet shift contributes to positive EIS anomalies but downward vertical velocity anomalies, the two effects cancel, and net observed changes in shortwave CRE are small. Global climate models generally capture the observed anomalies associated with midlatitude jet shifts. However, there is large intermodel spread in the shortwave CRE anomalies, with a subset of models showing a large shortwave cloud radiative warming over midlatitude oceans with a poleward jet shift. In these models, midlatitude shortwave CRE is sensitive to vertical velocity perturbations, but the observed sensitivity to EIS perturbations is underestimated. Consequently, these models might incorrectly estimate future midlatitude cloud feedbacks in regions where appreciable changes in both vertical velocity and EIS are projected

54 ENVIRONMENTAL SCIENCES↗

Simulation of African Easterly Waves in a Global Climate Model

Abstract African easterly waves (AEWs) exert significant influence on local and downstream high-impact weather including tropical cyclone (TC) genesis over the Atlantic. Accurate representation of AEWs in climate and weather prediction models therefore is necessary for skillful predictions. In this study, we examine simulated AEWs, including their evolution, vertical structure, and linkage to tropical cyclone genesis, in the NASA Goddard Earth Observing System Model, version 5 (GEOS-5) atmospheric global climate model. Identified by the leading empirical orthogonal function mode of time-filtered precipitation, the observed westward propagating AEWs along the southern track over the Atlantic are largely captured in GEOS-5, but with a slower phase speed and significantly weaker amplitude downstream off the West Africa coast. The weak downstream development of AEWs in GEOS-5 is accompanied by much reduced TC genesis over the main development region. Further analyses suggest that the slow westward propagation and weaker AEW amplitude downstream can be ascribed to a weak African easterly jet, while overestimated negative (positive) meridional potential vorticity (PV) gradients appear to the north (south) of 10°N in GEOS-5. The greatly overestimated positive meridional PV gradient to the south of 10°N is expected to generate strong horizontal stretching in the AEW wave pattern in the model, which hinders organization of convection and its feedback to sustain the AEW development. Persistent and vigorous AEW precipitation in the Guinea Highlands of the West Africa coast could also be responsible for reduced westward propagation of AEWs in the model.

Meteorology & Atmospheric Sciences↗

Characterization of 3D Cirrus Cloud and Radiation Fields Using ARS/AIRS/MODIS data and its Application to Climate Model

During the report period, we have made the following research accomplishments. First, we performed analysis for a number of MODIS scenes comprising of heavy dust events and ice clouds, covering regions of frequent dust outbreaks in East Asia, Middle East, and West Africa, as well as areas associated with long-range dust transports over the Equatorial Tropical Atlantic Ocean. These scenes contain both dust/aerosols and clouds. We collected suitable aerosol/ice-cloud data, correlated ice cloud and aerosol parameters by means of statistical analysis, and interpreted resulting correlation trends based on the physical principles governing cloud microphysics. Aerosol and cloud optical depths and cloud effective particle size inferred from MODIS for selected domains were analyzed from which the parameters including dust aerosol number concentration, ice cloud water path, and ice particle number concentration were subsequently derived. We illustrated that the Twomey (solar albedo) effect can be statistically quantified based on the slope of best-fit straight lines in the correlation study. Analysis of aerosol and cloud retrieval products revealed that for all cases, the region with a larger dust aerosol optical depth is always characterized by a smaller cloud particle size, consistent with the Twomey hypothesis for aerosol-cloud interactions. Second, we developed mean correlation curves with uncertainties associated with small ice-crystal concentration observations for the mean effective ice crystal size (De) and ice water content (IWC) by dividing the atmosphere into three characteristic regions: Tropics cirrus, Midlatitude cirrus, including a temperature classification to improve correlation, and Arctic ice clouds. We illustrated that De has a high correlation with IWC based on theoretical consideration and analysis of thousands of observed ice crystal data obtained from a number of ARM-DOE field campaigns and other experiments. The correlation has the form: ln(De) = a + b ln(IWC) + c ((ln(IWC))2, where a, b, and c are fitting coefficients and are functions of three regions. We demonstrated that this correlation can be effectively incorporated in GCMs and climate models that predict IWC - a significant advance in ice microphysics parameterization for interactive cloud-radiation analysis and feedback. Substantial July mean differences are shown in the OLR (W/m2) and precipitation (mm/day) patterns between UCLA GCM simulations based on Des determined from the De-IWC correlations and the control run using a fixed ice crystal size. Third, in order to improve the computation of spectral radiative transfer processes in the WRF model, we developed a consistent and efficient radiation scheme that can better resolve the spectral bands, determine the cloud optical properties, and provide more reliable and accurate radiative heating fields. In the newly developed radiation module, we have implemented in WRF a modified and improved version referred to as the Fu-Liou-Gu scheme, which includes a combination of delta-four-stream and delta-two-stream approximations for solar and IR flux calculations, respectively. This combination has been proven to be computationally efficient and at the same time to produce a high degree of accuracy. The incorporation of nongray gaseous absorption in multiple scattering atmospheres was based on the correlated k-distribution method. The solar and IR spectra are divided into 6 and 12 bands, respectively, according to the location of absorption bands of H2O, CO2, O3, CH4, N2O, and CFCs. We further included absorption by the water vapor continuum and a number of minor absorbers in the solar spectrum leading to an additional absorption of solar flux in a clear atmosphere on the order of 1-3 W/m2. Additionally, we incorporated the ice microphysics parameterization that includes an interactive mean effective ice crystal size in association with radiation parameterizations. The Fu-Liou-Gu scheme is an ideal tool for the simulation of radiative transfer and ice microphysics within the domain of WRF. It is particularly useful for studying direct and indirect aerosol radiative effects associated with ice cloud formation. The newly implemented radiation module has been demonstrated to work well in WRF and can be effectively used for studies related to cirrus cloud formation and evolution as well as aerosol-cloud-radiation interactions. With the newly implemented radiation scheme, the simulations of cloud cover and ice water path have been improved for cirrus clouds, with a more consistent comparison with the corresponding MODIS observations, especially for optically thin cirrus with an improvement of about 20% in the simulated mean ice water path.

97 MATHEMATICS AND COMPUTING↗