Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “CLIMATE MODELING”

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 973 records · Page 54

Hemispheric Differences in Tropical Lower Stratospheric Transport and Tracers Annual Cycle

Transport of long-lived tracers (such as O, CO, and N O) in the lower stratosphere largely determines the composition of the entire stratosphere. Stratospheric transport includes the mean residual circulation (with air rising in the tropics and sinking in the polar and middle latitudes), plus two-way isentropic (quasi-horizontal) mixing by eddies. However, the relative importance of two transport components remains uncertain. Previous studies quantified the relative role of these processes based on tropics-wide average characteristics under common assumption of well-mixed tropics. However, multiple instruments provide us with evidence that show significant differences in the seasonal cycle of ozone between the Northern (0-20N) and Southern (0-20S) tropical (NT and ST respectively) lower stratosphere. In this study we investigate these differences in tracer seasonality and quantify transport processes affecting tracers annual cycle amplitude using simulations from Goddard Earth Observing System Chemistry Climate Model (GEOSCCM) and Whole Atmosphere Community Climate Model (WACCM) and compare them to observations from the Microwave Limb Sounder (MLS) on the Aura satellite. We detect the observed contrast between the ST and NT in GEOSCCM and WACCM: annual cycle in ozone and other chemical tracers is larger in the NT than in the ST but opposite is true for the annual cycle in vertical advection. Ozone budgets in the models, analyzed based on the Transformed Eulerian Mean (TEM) framework, demonstrate a major role of quasi-horizontal mixing vertical advection in determining the NTST ozone distribution and behavior. Analysis of zonal variations in the NT and ST ozone annual cycles further suggests important role of North American and Asian Summer Monsoons (associated with strong isentropic mixing) on the lower stratospheric ozone in the NT. Furthermore, multi model comparison shows that most CCMs reproduce the observed characteristic of ozone annual cycle quite well. Thus, latitudinal variations within the tropics have to be considered in order to understand the balance between upwelling and quasi- horizontal mixing in the tropical lower stratosphere and the paradigm of well mixed tropics has to be reconsidered.

Chemistry-Climate Models↗

Multi-Frequency Analysis of Simulated versus Observed Variability in Tropospheric Temperature

Studies seeking to identify a human-caused global warming signal generally rely on climate model estimates of the “noise” of intrinsic natural variability. Assessing the reliability of these noise estimates is of critical importance. We evaluate here the statistical significance of differences between climate model and observational natural variability spectra for global-mean mid- to upper-tropospheric temperature (TMT). We use TMT information from satellites and large multimodel ensembles of forced and unforced simulations. Our main goal is to explore the sensitivity of model-versus-data spectral comparisons to a wide range of subjective decisions. These include the choice of satellite and climate model TMT datasets, the method for separating signal and noise, the frequency range considered, and the statistical model used to represent observed natural variability. Of particular interest is the amplitude of the interdecadal noise against which an anthropogenic tropospheric warming signal must be detected. We find that on time scales of 5–20 years, observed TMT variability is (on average) overestimated by the last two generations of climate models participating in the Coupled Model Intercomparison Project. This result is relatively insensitive to different plausible analyst choices, enhancing confidence in previous claims of detectable anthropogenic warming of the troposphere and indicating that these claims may be conservative. A further key finding is that two commonly used statistical models of short-term and long-term memory have deficiencies in their ability to capture the complex shape of observed TMT spectra.

54 ENVIRONMENTAL SCIENCES↗

Advancing Measurements and Understanding of the Rate and Structure of Atmospheric Warming

The rate and geographic pattern of atmospheric warming are key indicators of historical climate change and play a role in modulating future changes in Earth’s climate. The past and future evolution of atmospheric temperature is simulated using global climate models while observational estimates of past warming are derived from satellite microwave measurements. Individual climate models simulate widely varying rates of past and future atmospheric warming and most model simulations exhibit greater tropospheric (lowest ~10 km of atmosphere) temperature change than satellite observations between 1979 and 2020. This project examined intermodel differences in the pattern of atmospheric warming and how these differences influence climate feedbacks that can amplify or damp the rate of global surface warming. A key result is that climate model representation of the current climate can influence climate feedbacks and the simulation of future changes in climate. Another focus of this project was to analyze how models respond to different input datasets, such as volcanic aerosols or sea surface temperature (in atmosphere-only simulations). We find that different prescribed inputs can affect simulated changes in atmospheric temperature, even though “input uncertainty” is often unconsidered in model-observational comparisons. A final focus of this research project was to consider the influence of natural internal climate variability on satellite era changes in climate. Although the Earth is warming substantially due to anthropogenic emissions of greenhouse gases, the observed warming rate can be modulated by natural variations in the Earth’s climate. We find that natural climate variability has slowed the rate of tropical tropospheric warming, which explains model-satellite differences in the rate of tropospheric warming.

58 GEOSCIENCES↗

Models of climate cycles recorded in Martian polar layered deposits

Mars polar deposit accumulation models are used to predict the sequences of layer thicknesses. Continuous deposition and climate-modulated deposition models are compared, and two physical models of the latter type are developed in detail. In addition to examining the sensitivity of the sequence of layer thicknesses to model type and threshold, parallels with terrestrial ocean margin deposition are drawn with a view to better understanding the effects of a prolonged interval of nondeposition on stratigraphy and the modification of the central polar stratigraphy by scarp margin effects. When the threshold for deposition lies in the midrange with respect to excursions of the climate function, deposition occurs about half of the time, and there is only a slight modulation of the layer thickness. When the threshold for deposition is high, deposition occurs only a small fraction of the total time, and groups of layers form which are separated from other groups by long periods of nondeposition.

Cutts, J. A.↗

Regionally refined test bed in E3SM atmosphere model version 1 (EAMv1) and applications for high-resolution modelling

Climate simulation with more accurate process-level representation at finer resolutions is a pressing necessity in order to provide actionable information to policy-makers regarding extreme events in a changing climate. Computational limitation is a major obstacle for building, and running high-resolution (HR, here 0.25° average grid spacing at the equator) models (HRM). A more affordable path to HRM is to use a global regionally refined model (RRM), which only simulates a portion of the globe at HR while the remaining is at low-resolution (LR, 1°). In this study, we compare the Energy Exascale Earth System Model (E3SM) atmosphere model version 1 (EAMv1) RRM with the HR mesh over the contiguous United States (CONUS) to its corresponding globally uniform LR and HR configurations, as well as to observations and reanalysis data. The RRM has a significantly reduced computational cost (roughly proportional to the HR mesh size) relative to the globally uniform HRM. Over the CONUS, we evaluate the simulation of important dynamical and physical quantities as well as various precipitation measures. Differences between the RRM and HRM over the HR region are predominantly small, demonstrating that the RRM reproduces both well- and poorly simulated behaviours of the HRM over the CONUS. Further analysis based on RRM simulations with the LR vs. HR model parameters reveals that RRM performance is greatly influenced by the different parameter choices used in the LR and HR EAMv1. This is a result of the poor scale-aware behaviour of physical parameterizations, especially for variables influencing sub-grid scale physical processes. RRM can serve as a useful framework to test physics schemes across a range of scales, leading to improved consistency in future E3SM versions. Applying nudging-to-observations techniques within the RRM framework also demonstrates significant advantages over a free-running configuration for use as a testbed, and as such represents an efficient and more robust physics testbed capability. Our conclusions provide additional confirmatory evidence that the RRM is an efficient and effective approach for HRM development and hydrologic research.

54 ENVIRONMENTAL SCIENCES↗

Achieving Realistic Runoff in the Western United States with a Land Surface Model Forced by Dynamically Downscaled Meteorology

Abstract In this study, we calibrate a regional climate model’s (RCM) underlying land surface model (LSM). In addition to providing a realistic representation of runoff across the hydroclimatically diverse western United States, this is done to take advantage of the RCM’s ability to physically resolve meteorological forcing data in ungauged regions, and to prepare the calibrated hydrologic model for tight coupling, or the ability to represent land surface–atmosphere interactions, with the RCM. Specifically, we use a 9-km resolution meteorological forcing dataset across the western United States, from the fifth generation ECMWF Reanalysis (ERA5) downscaled by the Weather Research Forecasting (WRF) regional climate model, as an offline forcing for Noah-Multiparameterization (Noah-MP). We detail the steps involved in producing an LSM capable of accurately representing runoff, including physical parameterization selection, parameter calibration, and regionalization to ungauged basins. Based on our model evaluation from 1954 to 2021 for 586 basins with daily natural streamflow, the streamflow bias is reduced from 24.2% to 4.4%, and the median daily Nash–Sutcliffe efficiency (NSE) is improved from 0.12 to 0.36. When validating against basins with monthly natural streamflow data, we obtain a similar reduction in bias and a median monthly NSE improvement from 0.18 to 0.56. In this study, we also discover the optimal setup when using a donor-basin method to regionalize parameters to ungauged basins, which can vary by 0.06 NSE for unique designs of this regionalization method. Significance Statement This study provides useful guidance for improving a land surface model to accurately represent runoff across a spatially extensive and hydroclimatically diverse region (the western United States). The land surface model is updated to represent runoff more accurately at gauged basins, and then additionally updated for basins without observational data using a mathematical approach called the donor-basin method. We make use of a regional climate model’s reanalysis-derived meteorological data and its underlying land surface model to achieve realistic runoff. The calibrated land surface model can thus be tightly coupled in subsequent studies in a manner that should more accurately reflect runoff conditions. Findings from this study will serve as a useful reference for the atmospheric (and hydrologic) modeling communities and their ability to represent large-scale hydrology accurately.

Meteorology & Atmospheric 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↗

Plant water potential improves prediction of empirical stomatal models

Climate change is expected to lead to increases in drought frequency and severity, with deleterious effects on many ecosystems. Stomatal responses to changing environmental conditions form the backbone of all ecosystem models, but are based on empirical relationships and are not well-tested during drought conditions. Here, we use a dataset of 34 woody plant species spanning global forest biomes to examine the effect of leaf water potential on stomatal conductance and test the predictive accuracy of three major stomatal models and a recently proposed model. We find that current leaf-level empirical models have consistent biases of over-prediction of stomatal conductance during dry conditions, particularly at low soil water potentials. Furthermore, the recently proposed stomatal conductance model yields increases in predictive capability compared to current models, and with particular improvement during drought conditions. Our results reveal that including stomatal sensitivity to declining water potential and consequent impairment of plant water transport will improve predictions during drought conditions and show that many biomes contain a diversity of plant stomatal strategies that range from risky to conservative stomatal regulation during water stress. Such improvements in stomatal simulation are greatly needed to help unravel and predict the response of ecosystems to future climate extremes.

54 ENVIRONMENTAL SCIENCES↗

Relative Controls of Vapor Pressure Deficit and Soil Water Stress on Canopy Conductance in Global Simulations by an Earth System Model

Climate change impacts drought and vapor pressure deficit (VPD) that critically influence land-atmosphere water vapor and CO 2 exchange. There is a lack of model-based global evaluation of the relative roles of the two types of water stresses in limiting canopy conductance (G c ). Using the Energy Exascale Earth System Model (E3SM), we conducted four land-atmosphere coupled simulations for the historical period with and without a plant hydraulics scheme (PHS) and under present-day and quadrupling of CO 2 concentrations. By fitting an empirical model of G c to the E3SM simulations, we evaluated the sensitivity of G c to water stresses caused by VPD and soil water (SW) in the growing season. The empirical model is based on a multiplicative algorithm that adjusts G c according to VPD, soil water saturation, and their interaction. Our results show positive trends of VPD over land during the historical period, with large differences between simulations with and without PHS. Despite the differences, a majority of land pixels that experience both SW and VPD stresses are in semi-arid and tropical savanna ecosystems. G c is more sensitive to VPD than SW in less than 20% of land pixels where the empirical model provides a good fit for the simulated G c . Ignoring SW and VPD interaction may underestimate the relative role of VPD on G c by as much as 28%. Including PHS accentuates SW stress limitation on G c by simulating wetter soil and higher ET that lower VPD. Elevated CO 2 reduces G c by increasing water use efficiency and ameliorates SW stress on G c .

54 ENVIRONMENTAL SCIENCES↗

Evaluation of ACCMIP Outgoing Longwave Radiation from Tropospheric Ozone Using TES Satellite Observations.

We use simultaneous observations of tropospheric ozone and outgoing longwave radiation (OLR) sensitivity to tropospheric ozone from the Tropospheric Emission Spectrometer (TES) to evaluate model tropospheric ozone and its effect on OLR simulated by a suite of chemistry-climate models that participated in the Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP). The ensemble mean of ACCMIP models show a persistent but modest tropospheric ozone low bias (5-20 ppb) in the Southern Hemisphere (SH) and modest high bias (5-10 ppb) in the Northern Hemisphere (NH) relative to TES ozone for 2005-2010. These ozone biases have a significant impact on the OLR. Using TES instantaneous radiative kernels (IRK), we show that the ACCMIP ensemble mean tropospheric ozone low bias leads up to 120mW/ sq. m OLR high bias locally but zonally compensating errors reduce the global OLR high bias to 39+/- 41mW/ sq. m relative to TES data. We show that there is a correlation (Sq. R = 0.59) between the magnitude of the ACCMIP OLR bias and the deviation of the ACCMIP preindustrial to present day (1750-2010) ozone radiative forcing (RF) from the ensemble ozone RF mean. However, this correlation is driven primarily by models whose absolute OLR bias from tropospheric ozone exceeds 100mW/ sq. m. Removing these models leads to a mean ozone radiative forcing of 394+/- 42mW/ sq. m. The mean is about the same and the standard deviation is about 30% lower than an ensemble ozone RF of 384 +/- 60mW/ sq. m derived from 14 of the 16 ACCMIP models reported in a companion ACCMIP study. These results point towards a profitable direction of combining satellite observations and chemistry-climate model simulations to reduce uncertainty in ozone radiative forcing.

troposphere↗

Recent advances in understanding secondary organic aerosol: Implications for global climate forcing: Advances in Secondary Organic Aerosol

Anthropogenic emissions and land use changes have modified atmospheric aerosol concentrations and size distributions over time. Understanding preindustrial conditions and changes in organic aerosol due to anthropogenic activities is important because these features (1) influence estimates of aerosol radiative forcing and (2) can confound estimates of the historical response of climate to increases in greenhouse gases. Secondary organic aerosol (SOA), formed in the atmosphere by oxidation of organic gases, represents a major fraction of global submicron-sized atmospheric organic aerosol. Over the past decade, significant advances in understanding SOA properties and formation mechanisms have occurred through measurements, yet current climate models typically do not comprehensively include all important processes. Our review summarizes some of the important developments during the past decade in understanding SOA formation. We also highlight the importance of some processes that influence the growth of SOA particles to sizes relevant for clouds and radiative forcing, including formation of extremely low volatility organics in the gas phase, acid-catalyzed multiphase chemistry of isoprene epoxydiols, particle-phase oligomerization, and physical properties such as volatility and viscosity. Several SOA processes highlighted in this review are complex and interdependent and have nonlinear effects on the properties, formation, and evolution of SOA. Current global models neglect this complexity and nonlinearity and thus are less likely to accurately predict the climate forcing of SOA and project future climate sensitivity to greenhouse gases. Efforts are also needed to rank the most influential processes and nonlinear process-related interactions, so that these processes can be accurately represented in atmospheric chemistry-climate models.

54 ENVIRONMENTAL SCIENCES↗

Arctic Warming in Response to Regional Aerosol Emissions Reductions

This study examines the Arctic surface air temperature response to regional aerosol emissions reductions using three fully coupled chemistry–climate models: National Center for Atmospheric Research-Community Earth System Model version 1, Geophysical Fluid Dynamics Laboratory-Coupled Climate Model version 3 (GFDL-CM3) and Goddard Institute for Space Studies-ModelE version 2. Each of these models was used to perform a series of aerosol perturbation experiments, in which emissions of different aerosol types (sulfate, black carbon (BC), and organic carbon) in different northern mid-latitude source regions, and of biomass burning aerosol over South America and Africa, were substantially reduced or eliminated. We find that the Arctic warms in nearly every experiment, the only exceptions being the U.S. and Europe BC experiments in GFDL-CM3 in which there is a weak and insignificant cooling. The Arctic warming is generally larger than the global mean warming (i.e. Arctic amplification occurs), particularly during non-summer months. The models agree that changes in the poleward atmospheric moisture transport are the most important factor explaining the spread in Arctic warming across experiments: the largest warming tends to coincide with the largest increases in moisture transport into the Arctic. In contrast, there is an inconsistent relationship (correlation) across experiments between the local radiative forcing over the Arctic and the simulated Arctic warming, with this relationship being positive in one model (GFDL-CM3) and negative in the other two. Our results thus highlight the prominent role of poleward energy transport in driving Arctic warming and amplification, and suggest that the relative importance of poleward energy transport and local forcing/feedbacks is likely to be model dependent.

Arctic↗

Modeling Uncertainty in Climate Change: A Multi-Model Comparison

The economics of climate change involves a vast array of uncertainties, complicating both the analysis and development of climate policy. This study presents the results of the first comprehensive study of uncertainty in climate change using multiple integrated assessment models. The study looks at model and parametric uncertainties for population, total factor productivity, and climate sensitivity and estimates the pdfs of key output variables, including CO 2 concentrations, temperature, damages, and the social cost of carbon (SCC). One key finding is that parametric uncertainty is more important than uncertainty in model structure. Our resulting pdfs also provide insight on tail events.

54 ENVIRONMENTAL SCIENCES↗

The GEWEX Process Evaluation Study: GEWEX-PROES

Progress in realistically simulating the energy and water cycles in weather and climate models has been slower than desirable. Climate models submitted to the Coupled Model Intercomparison Project Phase-5 (CMIP-5), while more comprehensive than their predecessors, have shown little to no improvement in their biases in simulating key features of the climate system. As a result, uncertainties in global climate and hydrological sensitivities, and in the simulation of regional climate change, have not been reduced significantly. It is time to make use of these opportunities to significantly advance our understanding of key energy and water cycle processes at a wide range of space and time scales, and to provide a more insightful evaluation of the representation of these processes within models. This will require new ways of both analyzing the observations and diagnosing model behavior. The keys to success will be in skillfully combining different data sets and exploring relationships between them, as well as in the ability of models to reproduce those relationships correctly. In response to this challenge, a new GEWEX-wide activity, the GEWEX Process Evaluation Study, will take advantage of the opportunities that the combination of many of the existing data sets provide.

climate models↗

Final DOE-ASR Report for the Project “Advancing our Understanding and the Remote Sensing of Ice Clouds”

This project has evolved during its execution, and what follows are the key project findings. This project has arguably provided the first global view of how cirrus cloud (defined as having cloud base temperature T < 235 K) nucleation physics (evaluated through satellite retrievals of ice particle number concentration Ni, effective diameter De and ice water content IWC) evolves with the seasons for a given temperature, latitude zone and surface type (e.g. ocean vs. land), based on a new satellite remote sensing method developed for this project. The retrieval method is unique in that it is very sensitive to the small ice crystals that govern the number concentration Ni, allowing Ni to be retrieved. The method currently samples single-layer cirrus clouds having visible optical depth ranging from about 0.3 to 3.0, using co-located observations from the Infrared Imaging Radiometer (IIR) and from the CALIOP (Cloud and Aerosol Lidar with Orthogonal Polarization) lidar aboard the CALIPSO (Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation) polar orbiting satellite, employing IIR channels at 10.6 μm and 12.05 μm. Retrievals of Ni are primarily used to estimate the cirrus cloud formation mechanism; that is, either homo- or heterogeneous ice nucleation (henceforth hom and het). This is possible since, in general, hom produces more than an order of magnitude more ice crystals than does het. Thus the retrievals provide insight on how these mechanisms change with the seasons for a given latitude zone or region, based on the years 2008 and 2013. Using a conservative criterion for hom cirrus, on average, the sampled cirrus clouds formed through hom occur about 43% of the time in the Arctic and 50% of the time in the Antarctic, and during winter at mid-latitudes in the Northern Hemisphere, hom cirrus occur 37% of the time. Elsewhere (and during other seasons in the Northern Hemisphere mid-latitudes), this hom cirrus fraction is lower, and it is lowest in the tropics. Thus, the microphysical properties of cirrus clouds in the Polar Regions are much different than they are in the tropics; something unknown prior to this study. Moreover, the frequency of cirrus cloud occurrence in the Polar Regions varies strongly with season, peaking during winter in the Arctic and during spring in the Antarctic. Considering these seasonal changes in microphysics and inferred cloud coverage, this leads us to speculate that the buildup of Arctic cirrus during winter may significantly contribute to tropospheric heating in that region, possibly affecting winter jet-stream dynamics and mid-latitude weather patterns through the thermal-wind balance relationship. This cirrus cloud research provides essential guidance for realistically representing cirrus clouds in climate models; guidance previously unavailable. For example, mid-latitude hom cirrus were widespread during winter over or nearby mountainous terrain, evidently due to mountain-induced waves that produce strong updrafts at cirrus cloud levels. The treatment of turbulent mountain stress and gravity waves will likely need to be improved in climate models in order to adequately represent cirrus clouds outside the tropics. Another goal of this project was to develop a ground-based 94-GHz radar retrieval for winter snowstorms, based on (1) an improved analytical framework describing the interaction of radiation from radar with snowfall and (2) the development of a steady-state snow growth model that predicts the height-evolution of the ice particle size distribution through ice particle growth by vapor diffusion, aggregation and riming (i.e. the growth of snow through collisions with supercooled cloud droplets). Although activities (1) and (2) were completed, there was insufficient time to test and finalize the radar retrieval scheme. However, activity (2) provided a new method for relating ice particle mass “m” and projected area “A” to the ice particle maximum dimension “D”. The ice cloud microphysical processes (which determine ice cloud radiative properties) in climate models are parameterized in terms of these m-D and A-D relationships. By improving these relationships, the ice cloud radiative properties in Community Atmosphere Model version 5, or CAM5 (an atmosphere global climate model, or GCM) were improved. Student funding from the University of Nevada, Reno, was combined with funds from this project to conduct some basic research on the mechanism of the North American monsoon, or NAM. Federal research on the NAM has dwindled since 2006, but atmospheric soundings taken during research vessel cruises in the Gulf of California (GC) during the North American Monsoon Experiment (NAME) were used to reveal a likely mechanism that explains the relationship between an intrusion of tropical warm water into the GC during late spring-early summer and the onset of relatively heavy NAM rainfall in northwest Mexico and the southwestern United States. These soundings, combined with reanalysis data, satellite sea surface temperatures and satellite measurements of outgoing longwave radiation were used to develop and provide evidence for a planetary-scale NAM mechanism. As far as we know, no other physical explanation has been offered for the spring-summer evolution of the NAM system.

54 ENVIRONMENTAL SCIENCES↗

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↗