Search NASA⌕ Search

SEARCH · Search NASA

Results for “Climate model”

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,099 records · Page 61

Evaluation of native Earth system model output with ESMValTool v2.6.0

Earth system models (ESMs) are state-of-the-art climate models that allow numerical simulations of the past, present-day, and future climate. To extend our understanding of the Earth system and improve climate change projections, the complexity of ESMs heavily increased over the last decades. As a consequence, the amount and volume of data provided by ESMs has increased considerably. Innovative tools for a comprehensive model evaluation and analysis are required to assess the performance of these increasingly complex ESMs against observations or reanalyses. One of these tools is the Earth System Model Evaluation Tool (ESMValTool), a community diagnostic and performance metrics tool for the evaluation of ESMs. Input data for ESMValTool needs to be formatted according to the CMOR (Climate Model Output Rewriter) standard, a process that is usually referred to as “CMORization”. While this is a quasi-standard for large model intercomparison projects like the Coupled Model Intercomparison Project (CMIP), this complicates the application of ESMValTool to non-CMOR-compliant climate model output. In this paper, we describe an extension of ESMValTool introduced in v2.6.0 that allows seamless reading and processing of “native” climate model output, i.e., operational output produced by running the climate model through the standard workflow of the corresponding modeling institute. This is achieved by an extension of ESMValTool's preprocessing pipeline that performs a CMOR-like reformatting of the native model output during runtime. Thus, the rich collection of diagnostics provided by ESMValTool is now fully available for these models. For models that use unstructured grids, a further preprocessing step required to apply many common diagnostics is regridding to a regular latitude–longitude grid. Extensions to ESMValTool's regridding functions described here allow for more flexible interpolation schemes that can be used on unstructured grids. Currently, ESMValTool supports nearest-neighbor, bilinear, and first-order conservative regridding from unstructured grids to regular grids. Example applications of this new native model support are the evaluation of new model setups against predecessor versions, assessing of the performance of different simulations against observations, CMORization of native model data for contributions to model intercomparison projects, and monitoring of running climate model simulations. For the latter, new general-purpose diagnostics have been added to ESMValTool that are able to plot a wide range of variable types. Currently, five climate models are supported: CESM2 (experimental; at the moment, only surface variables are available), EC-Earth3, EMAC, ICON, and IPSL-CM6. As the framework for the CMOR-like reformatting of native model output described here is implemented in a general way, support for other climate models can be easily added.

58 GEOSCIENCES↗

Constraints and Opportunities in GCM Model Development

Over the past 30 years climate models have evolved from relatively simple representations of a few atmospheric processes to complex multi-disciplinary system models which incorporate physics from bottom of the ocean to the mesopause and are used for seasonal to multi-million year timescales. Computer infrastructure over that period has gone from punchcard mainframes to modern parallel clusters. Constraints of working within an ever evolving research code mean that most software changes must be incremental so as not to disrupt scientific throughput. Unfortunately, programming methodologies have generally not kept pace with these challenges, and existing implementations now present a heavy and growing burden on further model development as well as limiting flexibility and reliability. Opportunely, advances in software engineering from other disciplines (e.g. the commercial software industry) as well as new generations of powerful development tools can be incorporated by the model developers to incrementally and systematically improve underlying implementations and reverse the long term trend of increasing development overhead. However, these methodologies cannot be applied blindly, but rather must be carefully tailored to the unique characteristics of scientific software development. We will discuss the need for close integration of software engineers and climate scientists to find the optimal processes for climate modeling.

Schmidt, Gavin↗

Final Report Collaborative Project. Improving the Representation of Coastal and Estuarine Processes in Earth System Models

This project aimed to improve long term global climate simulations by resolving and enhancing the representation of the processes involved in the cycling of freshwater through estuaries and coastal regions. This was a collaborative multi-institution project consisting of physical oceanographers, climate model developers, and computational scientists. It specifically targeted the DOE objectives of advancing simulation and predictive capability of climate models through improvements in resolution and physical process representation. The main computational objectives were: 1. To develop computationally efficient, but physically based, parameterizations of estuary and continental shelf mixing processes for use in an Earth System Model (CESM). 2. To develop a two-way nested regional modeling framework in order to dynamically downscale the climate response of particular coastal ocean regions and to upscale the impact of the regional coastal processes to the global climate in an Earth System Model (CESM). 3. To develop computational infrastructure to enhance the efficiency of data transfer between specific sources and destinations, i.e., a point-to-point communication capability, (used in objective 1) within POP, the ocean component of CESM.

54 ENVIRONMENTAL SCIENCES↗

Optical Properties and Aging of Light Absorbing Secondary Organic Aerosol

The light-absorbing organic aerosol (OA), commonly referred to as “brown carbon (BrC)”, has attracted considerable attention in recent years because of its potential to affect atmospheric radiation balance, especially in the ultraviolet region and thus impact photochemical processes. A growing amount of data has indicated that BrC is prevalent in the atmosphere, which has motivated numerous laboratory and field studies; however, our understanding of the relationship between the chemical composition and optical properties of BrC remains limited. We conducted chamber experiments to investigate the effect of various VOC precursors, NOx concentrations, photolysis time and relative humidity (RH) on the light absorption of selected secondary organic aerosols (SOA). Light absorption of chamber generated SOA samples, especially aromatic SOA, was found to increase with NOx concentration, at moderate RH, and for the shortest photolysis aging times. The highest mass absorption coefficients (MAC) value is observed from toluene SOA products formed under high NOx conditions at moderate RH, in which nitro-aromatics were previously identified as the major light absorbing compounds. BrC light absorption is observed to decrease with photolysis time, correlated with a decline of the organonitrate fraction of SOA. SOA formed from mixtures of aromatics and isoprene absorb less visible and UV light than SOA formed from aromatic precursors alone on a mass basis. However, the mixed-SOA absorption was underestimated when optical properties were predicted using a two-product SOA formation model, as done in many current climate models. Further investigation, including analysis on detailed mechanisms, are required to explain the discrepancy.

Brown Carbon↗

An experimental seasonal hydrological forecasting system over the Yellow River basin – Part 2: The added value from climate forecast models

This is the second paper of a two-part series on introducing an experimental seasonal hydrological forecasting system over the Yellow River basin in northern China. While the natural hydrological predictability in terms of initial hydrological conditions (ICs) is investigated in a companion paper, the added value from eight North American Multimodel Ensemble (NMME) climate forecast models with a grand ensemble of 99 members is assessed in this paper, with an implicit consideration of human-induced uncertainty in the hydrological models through a post-processing procedure. The forecast skill in terms of anomaly correlation (AC) for 2 m air temperature and precipitation does not necessarily decrease over leads but is dependent on the target month due to a strong seasonality for the climate over the Yellow River basin. As there is more diversity in the model performance for the temperature forecasts than the precipitation forecasts, the grand NMME ensemble mean forecast has consistently higher skill than the best single model up to 6 months for the temperature but up to 2 months for the precipitation. The NMME climate predictions are downscaled to drive the variable infiltration capacity (VIC) land surface hydrological model and a global routing model regionalized over the Yellow River basin to produce forecasts of soil moisture, runoff and streamflow. And the NMME/VIC forecasts are compared with the Ensemble Streamflow Prediction method (ESP/VIC) through 6-month hindcast experiments for each calendar month during 1982–2010. As verified by the VIC offline simulations, the NMME/VIC is comparable to the ESP/VIC for the soil moisture forecasts, and the former has higher skill than the latter only for the forecasts at long leads and for those initialized in the rainy season. The forecast skill for runoff is lower for both forecast approaches, but the added value from NMME/VIC is more obvious, with an increase of the average AC by 0.08–0.2. To compare with the observed streamflow, both the hindcasts from NMME/VIC and ESP/VIC are post-processed through a linear regression model fitted by using VIC offline-simulated streamflow. The post-processed NMME/VIC reduces the root mean squared error (RMSE) from the post-processed ESP/VIC by 5–15 %. And the reduction occurs mostly during the transition from wet to dry seasons. As a result, with the consideration of the uncertainty in the hydrological models, the added value from climate forecast models is decreased especially at short leads, suggesting the necessity of improving the large-scale hydrological models in human-intervened river basins.

54 ENVIRONMENTAL SCIENCES↗

Software Engineering Support of the Third Round of Scientific Grand Challenge Investigations: An Earth Modeling System Software Framework Strawman Design that Integrates Cactus and UCLA/UCB Distributed Data Broker

One of the most significant challenges in large-scale climate modeling, as well as in high-performance computing in other scientific fields, is that of effectively integrating many software models from multiple contributors. A software framework facilitates the integration task. both in the development and runtime stages of the simulation. Effective software frameworks reduce the programming burden for the investigators, freeing them to focus more on the science and less on the parallel communication implementation, while maintaining high performance across numerous supercomputer and workstation architectures. This document proposes a strawman framework design for the climate community based on the integration of Cactus, from the relativistic physics community, and UCLA/UCB Distributed Data Broker (DDB) from the climate community. This design is the result of an extensive survey of climate models and frameworks in the climate community as well as frameworks from many other scientific communities. The design addresses fundamental development and runtime needs using Cactus, a framework with interfaces for FORTRAN and C-based languages, and high-performance model communication needs using DDB. This document also specifically explores object-oriented design issues in the context of climate modeling as well as climate modeling issues in terms of object-oriented design.

Talbot, Bryan↗

Improving collaborations between empiricists and modelers to advance grassland community dynamics in ecosystem models

Climate change, increasing atmospheric CO 2 , and land use change have altered biogeochemical and hydrologic cycles world-wide, with grassland systems being particularly vulnerable to resulting vegetation shifts (Komatsu et al., 2019). Therefore, incorporating plant community dynamics into ecosystem models is critical for accurate forecasting of ecosystem responses to global change (Levine, 2016). Process-based ecosystem models, which simulate the biogeochemical transfers of mass and energy among biota, the subsurface, and atmosphere, require representation of dynamic composition of organisms within ecosystems. For example, these models simulate leaf and plant-level characteristics, such as electron transport rate and allometry of carbon (C) allocation, to predict how net primary productivity and other ecosystem processes respond to abiotic drivers. These models are particularly useful in scaling from organismal to ecosystem levels but are still underdeveloped in their ability to capture community change, especially in grassland ecosystems. To represent compositional changes, these models must simulate competition, mortality, establishment, and reproduction of plant populations within communities. Yet, current ecosystem modeling approaches to forecast plant community change have derived from studies of forested systems and are either too coarse to capture fine-scale community dynamics (e.g. dynamic global vegetation models (DGVMs)) or too complex to be used at large spatial scales (e.g. forest gap models).

54 ENVIRONMENTAL SCIENCES↗

Quantifying the Effect of Mixing on the Mean Age of Air in CCMVal-2 and CCMI-1 Models

The stratospheric age of air (AoA) is a useful measure of the overall capabilities of a general circulation model (GCM) to simulate stratospheric transport. Previous studies have reported a large spread in the simulation of AoA by GCMs and coupled chemistry-climate models (CCMs). Compared to observational estimates, simulated AoA is mostly too low. Here we attempt to untangle the processes that lead to the AoA differences between the models and between models and observations. AoA is influenced by both mean transport by the residual circulation and two-way mixing; we quantify the effects of these processes using data from the CCM inter-comparison projects CCMVal-2 (Chemistry-Climate Model Validation Activity 2) and CCMI-1 (Chemistry-Climate Model Initiative, phase 1). Transport along the residual circulation is measured by the residual circulation transit time (RCTT). We interpret the difference between AoA and RCTT as additional aging by mixing. Aging by mixing thus includes mixing on both the resolved and subgrid scale. We find that the spread in AoA between the models is primarily caused by differences in the effects of mixing and only to some extent by differences in residual circulation strength. These effects are quantified by the mixing efficiency, a measure of the relative increase in AoA by mixing. The mixing efficiency varies strongly between the models from 0.24 to 1.02. We show that the mixing efficiency is not only controlled by horizontal mixing, but by vertical mixing and vertical diffusion as well. Possible causes for the differences in the models' mixing efficiencies are discussed. Differences in subgrid-scale mixing (including differences in advection schemes and model resolutions) likely contribute to the differences in mixing efficiency. However, differences in the relative contribution of resolved versus parameterized wave forcing do not appear to be related to differences in mixing efficiency or AoA.

CCMI↗

Attribution of the Australian bushfire risk to anthropogenic climate change

Abstract. Disastrous bushfires during the last months of 2019 and January 2020 affected Australia, raising the question to what extent the risk of these fires was exacerbated by anthropogenic climate change. To answer the question for southeastern Australia, where fires were particularly severe, affecting people and ecosystems, we use a physically based index of fire weather, the Fire Weather Index; long-term observations of heat and drought; and 11 large ensembles of state-of-the-art climate models. We find large trends in the Fire Weather Index in the fifth-generation European Centre for Medium-Range Weather Forecasts (ECMWF) Atmospheric Reanalysis (ERA5) since 1979 and a smaller but significant increase by at least 30 % in the models. Therefore, we find that climate change has induced a higher weather-induced risk of such an extreme fire season. This trend is mainly driven by the increase of temperature extremes. In agreement with previous analyses we find that heat extremes have become more likely by at least a factor of 2 due to the long-term warming trend. However, current climate models overestimate variability and tend to underestimate the long-term trend in these extremes, so the true change in the likelihood of extreme heat could be larger, suggesting that the attribution of the increased fire weather risk is a conservative estimate. We do not find an attributable trend in either extreme annual drought or the driest month of the fire season, September–February. The observations, however, show a weak drying trend in the annual mean. For the 2019/20 season more than half of the July–December drought was driven by record excursions of the Indian Ocean Dipole and Southern Annular Mode, factors which are included in the analysis here. The study reveals the complexity of the 2019/20 bushfire event, with some but not all drivers showing an imprint of anthropogenic climate change. Finally, the study concludes with a qualitative review of various vulnerability and exposure factors that each play a role, along with the hazard in increasing or decreasing the overall impact of the bushfires.

van Oldenborgh, Geert Jan (ORCID:0000000268989535)↗

HABITABLE ZONES AROUND MAIN-SEQUENCE STARS: NEW ESTIMATES

Identifying terrestrial planets in the habitable zones (HZs) of other stars is one of the primary goals of ongoing radial velocity (RV) and transit exoplanet surveys and proposed future space missions. Most current estimates of the boundaries of the HZ are based on one-dimensional (1D), cloud-free, climate model calculations by Kasting et al. However, this model used band models that were based on older HITRAN and HITEMP line-by-line databases. The inner edge of the HZ in the Kasting et al. model was determined by loss of water, and the outer edge was determined by the maximum greenhouse provided by a CO{sub 2} atmosphere. A conservative estimate for the width of the HZ from this model in our solar system is 0.95-1.67 AU. Here an updated 1D radiative-convective, cloud-free climate model is used to obtain new estimates for HZ widths around F, G, K, and M stars. New H{sub 2}O and CO{sub 2} absorption coefficients, derived from the HITRAN 2008 and HITEMP 2010 line-by-line databases, are important improvements to the climate model. According to the new model, the water-loss (inner HZ) and maximum greenhouse (outer HZ) limits for our solar system are at 0.99 and 1.70 AU, respectively, suggesting that the present Earth lies near the inner edge. Additional calculations are performed for stars with effective temperatures between 2600 and 7200 K, and the results are presented in parametric form, making them easy to apply to actual stars. The new model indicates that, near the inner edge of the HZ, there is no clear distinction between runaway greenhouse and water-loss limits for stars with T{sub eff} {approx}< 5000 K, which has implications for ongoing planet searches around K and M stars. To assess the potential habitability of extrasolar terrestrial planets, we propose using stellar flux incident on a planet rather than equilibrium temperature. This removes the dependence on planetary (Bond) albedo, which varies depending on the host star's spectral type. We suggest that conservative estimates of the HZ (water-loss and maximum greenhouse limits) should be used for current RV surveys and Kepler mission to obtain a lower limit on {eta}{sub Circled-Plus }, so that future flagship missions like TPF-C and Darwin are not undersized. Our model does not include the radiative effects of clouds; thus, the actual HZ boundaries may extend further in both directions than the estimates just given.

79 ASTRONOMY AND ASTROPHYSICS↗

Projected precipitation changes within the Great Lakes and Western Lake Erie Basin: a multi‐model analysis of intensity and seasonality

ABSTRACT The Great Lakes region encompasses the largest freshwater lake network in the world and supports a diverse network of agriculture, transportation, and tourism. Recently, Lake Erie has experienced increased hypoxia events, which have been attributed to agricultural practices and changes in run‐off. Here we examine the projected changes in extreme precipitation events to address concerns regarding regional agriculture, surface run‐off, and subsequent water quality. Precipitation projections within the overall Great Lakes Basin and the Western Lake Erie Basin subregion are examined using climate model simulations of varying spatial resolutions to understand historical precipitation and projected future precipitation. We develop three model ensembles for the historical period (1980–1999) and the mid‐century (2041–2060) that cover a range of spatial resolutions and future emissions scenarios, including: (1) 12 global model members from the fifth Climate Model Intercomparison Project ( CMIP5 ) using Representative Concentration Pathway ( RCP ) 8.5, (2) ten regional climate model ( RCM ) members from the North American Regional Climate Change Assessment Program driven by CMIP3 global models using the A2 emissions scenario, and (3) two high resolution RCM simulations ( RCM4 ) driven by CMIP5 global models using the RCP 8.5 scenario. For the historical period, all model ensembles overestimate winter and spring precipitation, and many of the models simulate a summer drying that is not observed. At mid‐century, most of the models predict a 10–20% increase in precipitation depending on the time of year. Daily probability distribution functions from three model ensembles reveal spring seasonal increases in high precipitation event probabilities when compared to the historical period, suggesting an increase in the frequency of high intensity precipitation at mid‐century. Overall, the presence of lakes or higher spatial resolution does not ensure improved representation of historical processes, and more complex interactions between large‐scale dynamics, local feedbacks, and physical parameterizations drive the model spread.

Basile, Samantha J.↗

Final Report for “Simulating the Arctic Winter Longwave Indirect Effects. A New Parameterization for Frost Flower Aerosol Salt Emissions” (DESC0006679) for 9/15/2011 through 9/14/2015

Description of the Project: This project has improved the aerosol formulation in a global climate model by using innovative new field and laboratory observations to develop and implement a novel wind-driven sea ice aerosol flux parameterization. This work fills a critical gap in the understanding of clouds, aerosol, and radiation in polar regions by addressing one of the largest missing particle sources in aerosol-climate modeling. Recent measurements of Arctic organic and inorganic aerosol indicate that the largest source of natural aerosol during the Arctic winter is emitted from crystal structures, known as frost flowers, formed on a newly frozen sea ice surface [Shaw et al., 2010]. We have implemented the new parameterization in an updated climate model making it the first capable of investigating how polar natural aerosol-cloud indirect effects relate to this important and previously unrecognized sea ice source. The parameterization is constrained by Arctic ARM in situ cloud and radiation data. The modified climate model has been used to quantify the potential pan-Arctic radiative forcing and aerosol indirect effects due to this missing source. This research supported the work of one postdoc (Li Xu) for two years and contributed to the training and research of an undergraduate student. This research allowed us to establish a collaboration between SIO and PNNL in order to contribute the frost flower parameterization to the new ACME model. One peer-reviewed publications has already resulted from this work, and a manuscript for a second publication has been completed. Additional publications from the PNNL collaboration are expected to follow.

54 ENVIRONMENTAL SCIENCES↗

Improving tide-estuary representation in MPAS-Ocean (Final Report)

The overarching goal of this research is to produce a skill-assessed tide-estuary representation, including BGC (biogeochemistry) components of the ecosystem, into the MPAS-O global climate modeling system. Specifically, in the original proposal we proposed one-way nesting, into global MPAS-O, of a high-resolution model for TES (tidal and estuary systems) including the coupling of the BGC components. We have completed the proposed tasks during the past 4 years. Specifically, after extensive consultation with the MPAS-O group during the first year, we have decided to go beyond the original proposal and to two-way nest SCHISM into MPAS-O. As a result, SCHISM-CoSiNE (with CoSiNE being the BGC model) is now fully integrated into MPAS-O as a subroutine (more precisely, an analysis member) and can be invoked as MPAS-O is running with two-way exchange of information between the TES and global models (with proper spatial interpolation done on the fly). As originally proposed, we have also developed the new finite-volume (FV) model MPAS-OI (aka ELCIRC-sub) for TES. We have conducted several tests using the coupled MPAS-SCHISM- CoSiNE model and also MPAS-OI for Columbia River, San Francisco Bay and Chesapeake Bay. Four journal publications have been produced so far as the result of the project and 1 is under review (cf. Publication section). With these achievements, we will be ready to conduct fully coupled simulations of TES and global ocean anywhere, as soon as some missing capabilities in MPAS-O (e.g. tides, improved ability to handle the boundary conditions, and streamlining of the global BGC model) are worked out. The significance of the outcome of this project cannot be overstated. It marks the first time, as far as we know, that a global climate model is tightly coupled to a bona fide nearshore and terrestrial model (SCHISM system); in fact SCHISM is the first bona fide compound flooding model that seamlessly covers the coastal ocean, estuaries, rivers and watersheds (Ye et al. 2020; Zhang et al. 2020). The coupled model will reduce the errors and uncertainties of the MPAS-O global climate model system. It also complements and extends on-going E3SM project, as the latter seeks in the long-term to include components required to simulate impacts of sea-level change and storm surge on coastal regions, including wave models and focusing resolution in coastal and storm-track regions, as an important component of the Cryosphere-ocean system. Note that a spectral wave model (the unstructured-grid version of Wave Watch III, aka Wind Wave Model) is already incorporated inside the SCHISM system and is therefore directly accessible in the coupled system.

58 GEOSCIENCES↗

Subsurface Cooling and Sea Surface Temperature Pattern Formation Over the Equatorial Pacific

Abstract The equatorial cold tongue region has not warmed up in response to historical radiative forcing in the real world, contrary to the strong warming often simulated by climate models. Here we demonstrate that climate models fail to represent one or both of the key processes driving observed sea surface temperature (SST) pattern formation: a realistic surface wind stress pattern shaping subsurface cooling through wind‐driven circulation changes, and effective connectivity between subsurface and surface temperatures via upwelling and mixing. Consequently, none of the models approximate the observed lack of cold tongue SST warming and strengthening of zonal SST gradient across the equatorial Pacific. Furthermore, those that come closest achieve this due to interhemispheric warming differences rather than equatorial dynamics as observed. Addressing different origins of subsurface cooling in observations and simulations, and how they connect to SST, will lead to improved understanding of tropical Pacific SST changes to date and how they will evolve in the future.

Jiang, Feng [Lamont‐Doherty Earth Observatory Colu↗

A note on the maintenance of the atmospheric kinetic energy

The winter simulations of the GLAS climate model and the NCAR community climate model are used to examine the maintenance of the atmospheric kinetic energy. It is found that the kinetic energy is generated in the lower latitudes south of the maximum westerlies, transported northward and then, destroyed in the midlatitudes north of the maximum westerlies. Therefore, the atmospheric kinetic energy is maintained by the counterbalance between the divergence (convergence) of kinetic energy flux and generation (destruction) of kinetic energy in lower (middle) latitudes.

Chen, T.-C.↗

Hierarchical sensitivity analysis for simulating barrier island geomorphologic responses to future storms and sea-level rise

This paper presents a new application of an advanced hierarchical sensitivity analysis of a new climate model of barrier island geomorphological evolution. The implemented sensitivity analysis in this study integrates a hierarchical uncertainty framework with a variance-based global sensitivity analysis to decompose the different model input uncertainties. The analysis can provide quantitative and accurate measurements for the relative importance of uncertain model input factors while considering their dependence relationships. The climate model used in this research was the barrier island profile (BIP) model, which is a new computer code developed to simulate barrier island morphological evolution over periods ranging from years to decades under the impacts of accelerated future sea-level rise and long-term changes in the storm climate. In the application of the model, the BIP model was used to evaluate the responses of a series of barrier island cross-sections derived for Santa Rosa Island, Florida, to random storm events and five potential accelerated rates of sea-level rise projected over the next century. The uncertain model input factors thus include the scenario uncertainty caused by alternative future sea-level rise scenarios and the parametric uncertainties of random storm parameters and dune characteristics. The study results reveal that the occurrence of storms is the most important factor for the evolution of sand dunes on the barrier island and the impact of sea-level rise is essential to the morphological change of the island backshore environment. The analysis can provide helpful insights for coastal management and planning. This hierarchical sensitivity analysis is mathematically general and rigorous and can be applied to a wide range of climate models.

Dai, Heng↗

Evaluating Large-Eddy Simulations of Boundary-Layer Cloud Response to Climate Change using MAGIC Observations

Over the cool waters off California, a persistent and extensive deck of low-lying clouds forms atop a shallow, turbulent layer of moist air extending 500-2000 meters above the ocean, above which the air is much warmer and drier. Its west edge varies with weather and season and is typically half way to Hawaii. This cloud deck and similar cloud decks elsewhere reflect sunlight back to space, cooling the underlying ocean and indeed the entire planet. They must be accurately represented by a climate model to correctly simulate our current climate or future global warming. Most climate models suggest these cloud decks will thin and break up in a warmer CO 2 -enriched climate, so extra sunlight would reach the ocean surface and amplify warming worldwide. But these models have difficulty accurately simulating the cloud thickness and extent because the turbulent cloud-forming air currents are much too small to be simulated by their coarse computational grid. In 2012-2013, the DOE Atmospheric Radiation Measurement MAGIC program instrumented a container ship transiting between California and Hawaii to observe clouds and rain, small particles (aerosols), and weather conditions. Data were obtained from 20 4-6 day cruises in each direction. This study tested whether a sophisticated ‘large-eddy simulation’ (LES) of the turbulent air motions that create the clouds in this region accurately predicts the variations in cloud properties seen across the different cruises. If so, that model is more believable for predicting cloud changes in this region in a future warmer climate. Our LES used a numerical mesh over a 6.4x6.4 km region with 50 m horizontal and 5 m vertical grid spacing to simulate the cloud-forming eddies and the sharp overlying transition between moist and dry air. To better compare with MAGIC observations, we configured our model to simulate a region that moved along with the ship. We used a weather forecast model to specify effects of the large-scale air flow. We compared our results with balloon ascents, a radar and other instruments pointed upward at the clouds, measurements of downward sunlight and infrared radiation, aerosols, and weather variables for all suitable cruise periods. Our LES was remarkably skillful in predicting the cloudiness variations measured by the ship. MAGIC also sampled mesoscale cells, regions of cloud thickening and thinning over distances of 10-50 km. An LES run over a larger region 100 km on a side shows these cells form spontaneously, enhanced by strong radiative heat loss at the cloud tops.

54 ENVIRONMENTAL SCIENCES↗

Final Technical Report for "Ice nuclei relation to aerosol properties: Data analysis and model parameterization for IN in mixed-phase clouds" (DOE/SC00002354)

Clouds play an important role in weather and climate. In addition to their key role in the hydrologic cycle, clouds scatter incoming solar radiation and trap infrared radiation from the surface and lower atmosphere. Despite their importance, feedbacks involving clouds remain as one of the largest sources of uncertainty in climate models. To better simulate cloud processes requires better characterization of cloud microphysical processes, which can affect the spatial extent, optical depth and lifetime of clouds. To this end, we developed a new parameterization to be used in numerical models that describes the variation of ice nuclei (IN) number concentrations active to form ice crystals in mixed-phase (water droplets and ice crystals co-existing) cloud conditions as these depend on existing aerosol properties and temperature. The parameterization is based on data collected using the Colorado State University continuous flow diffusion chamber in aircraft and ground-based campaigns over a 14-year period, including data from the DOE-supported Mixed-Phase Arctic Cloud Experiment. The resulting relationship is shown to more accurately represent the variability of ice nuclei distributions in the atmosphere compared to currently used parameterizations based on temperature alone. When implemented in one global climate model, the new parameterization predicted more realistic annually averaged cloud water and ice distributions, and cloud radiative properties, especially for sensitive higher latitude mixed-phase cloud regions. As a test of the new global IN scheme, it was compared to independent data collected during the 2008 DOE-sponsored Indirect and Semi-Direct Aerosol Campaign (ISDAC). Good agreement with this new data set suggests the broad applicability of the new scheme for describing general (non-chemically specific) aerosol influences on IN number concentrations feeding mixed-phase Arctic stratus clouds. Finally, the parameterization was implemented into a regional cloud-resolving model to compare predictions of ice crystal concentrations and other cloud properties to those observed in two intensive case studies of Arctic stratus during ISDAC. Our implementation included development of a prognostic scheme of ice activation using the IN parameterization so that the most realistic treatment of ice nuclei, including their budget (gains and losses), was achieved. Many cloud microphysical properties and cloud persistence were faithfully reproduced, despite a tendency to under-predict (by a few to several times) ice crystal number concentrations and cloud ice mass, in agreement with some other studies. This work serves generally as the basis for improving predictive schemes for cloud ice crystal activation in cloud and climate models, and more specifically as the basis for such a scheme to be used in a Multi-scale Modeling Format (MMF) that utilizes a connected system of cloud-resolving models on a global grid in an effort to better resolve cloud processes and their influence on climate.

54 ENVIRONMENTAL SCIENCES↗