A decision tree algorithm for investigation of model biases related to dynamical cores and physical parameterizations

An object‐based evaluation method using a pattern recognition algorithm (i.e., classification trees) is applied to the simulated orographic precipitation for idealized experimental setups using the National Center of Atmospheric Research (NCAR) Community Atmosphere Model (CAM) with the finite volume (FV) and the Eulerian spectral transform dynamical cores with varying resolutions. Daily simulations were analyzed and three different types of precipitation features were identified by the classification tree algorithm. The statistical characteristics of these features (i.e., maximum value, mean value, and variance) were calculated to quantify the difference between the dynamical cores and changing resolutions. Even with the simple and smooth topography in the idealized setups, complexity in the precipitation fields simulated by the models develops quickly. The classification tree algorithm using objective thresholding successfully detected different types of precipitation features even as the complexity of the precipitation field increased. The results show that the complexity and the bias introduced in small‐scale phenomena due to the spectral transform method of CAM Eulerian spectral dynamical core is prominent, and is an important reason for its dissimilarity from the FV dynamical core. The resolvable scales, both in horizontal and vertical dimensions, have significant effect on the simulation of precipitation. The results of this study also suggest that an efficient and informative study about the biases produced by GCMs should involve daily (or even hourly) output (rather than monthly mean) analysis over local scales.Key PointsThe complexity and the bias introduced in small‐scale phenomena due to the spectral transform method of CAM Eulerian spectral dynamical core are prominentThe classification tree algorithm with objective thresholding is successful in detecting different types of precipitation features with high spatial complexityAn efficient and informative study about the biases produced by GCMs should involve daily (or hourly) output (rather than monthly mean) analysis over local scalesPeer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/136040/1/jame20331.pdfhttp://deepblue.lib.umich.edu/bitstream/2027.42/136040/2/jame20331_am.pd

Analysis of stratospheric ozone, temperature, and minor constituent data

The objective of this research is to use available satellite measurements of temperature and constituent concentrations to test the conceptual picture of stratospheric chemistry and transport. This was originally broken down into two sub-goals: first, to use the constituent data to search for critical tests of our understanding of stratospheric chemistry and second, to examine constituent transport processes emphasizing interactions with chemistry on various time scales. A third important goal which has evolved is to use the available solar backscattered ultraviolet (SBUV) and Total Ozone Mapping Spectrometer (TOMS) data from Nimbus 7 to describe the morphology of recent changes in Antarctic and global ozone with emphasis on searching for constraints to theories. The major effort now being pursued relative to the two original goals is our effort as a theoretical team for the Arctic Airborne Stratospheric Expedition (AASE). Our effort for the AASE is based on the 3D transport and chemistry model at Goddard. Our goal is to use this model to place the results from the mission data in a regional and global context. Specifically, we set out to make model runs starting in late December and running through March of 1989, both with and without heterogeneous chemistry. The transport is to be carried out using dynamical fields from a 4D data assimilation model being developed under separate funding from this task. We have successfully carried out a series of single constituent transport experiments. One of the things demonstrated by these runs was the difficulty in obtaining observed low N2O abundances in the vortex without simultaneously obtaining very high ozone values. Because the runs start in late December, this difficulty arises in the attempt to define consistent initial conditions for the 3D model. To accomplish a consistent set of initial conditions, we are using the 2D photochemistry-transport model of Jackman and Douglass and mapping in potential temperature, potential vorticity space as developed by Schoeberl and coworkers

Software Testing and Verification in Climate Model Development

Over the past 30 years most climate models have grown from relatively simple representations of a few atmospheric processes to a complex multi-disciplinary system. Computer infrastructure over that period has gone from punch card mainframes to modem parallel clusters. Model implementations have become complex, brittle, and increasingly difficult to extend and maintain. Existing verification processes for model implementations rely almost exclusively upon some combination of detailed analysis of output from full climate simulations and system-level regression tests. In additional to being quite costly in terms of developer time and computing resources, these testing methodologies are limited in terms of the types of defects that can be detected, isolated and diagnosed. Mitigating these weaknesses of coarse-grained testing with finer-grained "unit" tests has been perceived as cumbersome and counter-productive. In the commercial software sector, recent advances in tools and methodology have led to a renaissance for systematic fine-grained testing. We discuss the availability of analogous tools for scientific software and examine benefits that similar testing methodologies could bring to climate modeling software. We describe the unique challenges faced when testing complex numerical algorithms and suggest techniques to minimize and/or eliminate the difficulties

Fundamentals of Modeling, Data Assimilation, and High-performance Computing

This lecture will introduce the concepts of modeling, data assimilation and high- performance computing as it relates to the study of atmospheric composition. The lecture will work from basic definitions and will strive to provide a framework for thinking about development and application of models and data assimilation systems. It will not provide technical or algorithmic information, leaving that to textbooks, technical reports, and ultimately scientific journals. References to a number of textbooks and papers will be provided as a gateway to the literature

Stratospheric Data Analysis System (STRATAN)

A state of the art stratospheric analyses using a coupled stratosphere/troposphere data assimilation system is produced. These analyses can be applied to stratospheric studies of all types. Of importance to this effort is the application of the Stratospheric Data Analysis System (STRATAN) to constituent transport and chemistry problems

Stratospheric General Circulation with Chemistry Model (SGCCM)

In the past two years constituent transport and chemistry experiments have been performed using both simple single constituent models and more complex reservoir species models. Winds for these experiments have been taken from the data assimilation effort, Stratospheric Data Analysis System (STRATAN)

Satellite observation and mapping of wintertime ozone variability in the lower stratosphere

Comparison is made between 30 mbar ozone fields that are generated by a transport chemistry model utilizing the winds from the Goddard Space Flight Center stratospheric data assimilation system (STRATAN), observations from the LIMS instrument on Nimbus-7, and the ozone fields that result from 'flying a mathematical simulation of LIMS observations through the transport chemistry model ozone fields. The modeled ozone fields were found to resemble the LIMS observations, but the model fields show much more temporal and spatial structure than do the LIMS observations. The 'satellite mapped' model results resemble the LIMS observations much more closely. These results are very consistent with the earlier discussions of satellite space-time sampling by Salby