http://people.atmos.ucla.edu/alex/ROMS/ROMSArticle2005.pdf

The regional oceanic modeling system (ROMS): a split-explicit, free-surface, topography-following-coordinate oceanic model

The fourth version of the Community Climate System Model (CCSM4) was recently completed and released to the climate community. This paper describes developments to all CCSM components, and documents fully coupled preindustrial control runs compared to the previous version, CCSM3. Using the standard atmosphere and land resolution of 1° results in the sea surface temperature biases in the major upwelling regions being comparable to the 1.4°-resolution CCSM3. Two changes to the deep convection scheme in the atmosphere component result in CCSM4 producing El Nino–Southern Oscillation variability with a much more realistic frequency distribution than in CCSM3, although the amplitude is too large compared to observations. These changes also improve the Madden–Julian oscillation and the frequency distribution of tropical precipitation. A new overflow parameterization in the ocean component leads to an improved simulation of the Gulf Stream path and the North Atlantic Ocean meridional overturning circulati...

/pdf/the-community-climate-system-model-version-4-30sj6n6rru.pdf

The Community Climate System Model Version 4

Eleven coupled climate–carbon cycle models used a common protocol to study the coupling between climate change and the carbon cycle. The models were forced by historical emissions and the Intergovernmental Panel on Climate Change (IPCC) Special Report on Emissions Scenarios (SRES) A2 anthropogenic emissions of CO2 for the 1850–2100 time period. For each model, two simulations were performed in order to isolate the impact of climate change on the land and ocean carbon cycle, and therefore the climate feedback on the atmospheric CO2 concentration growth rate. There was unanimous agreement among the models that future climate change will reduce the efficiency of the earth system to absorb the anthropogenic carbon perturbation. A larger fraction of anthropogenic CO2 will stay airborne if climate change is accounted for. By the end of the twenty-first century, this additional CO2 varied between 20 and 200 ppm for the two extreme models, the majority of the models lying between 50 and 100 ppm. The higher CO2 levels led to an additional climate warming ranging between 0.1° and 1.5°C. All models simulated a negative sensitivity for both the land and the ocean carbon cycle to future climate. However, there was still a large uncertainty on the magnitude of these sensitivities. Eight models attributed most of the changes to the land, while three attributed it to the ocean. Also, a majority of the models located the reduction of land carbon uptake in the Tropics. However, the attribution of the land sensitivity to changes in net primary productivity versus changes in respiration is still subject to debate; no consensus emerged among the models.

/pdf/climate-carbon-cycle-feedback-analysis-results-from-the-2hflmwlg71.pdf

Climate–Carbon Cycle Feedback Analysis: Results from the C4MIP Model Intercomparison

The numerical implementation of an ocean model based on the incompressible Navier Stokes equations which is designed for studies of the ocean circulation on horizontal scales less than the depth of the ocean right up to global scale is described. A "pressure correction" method is used which is solved as a Poisson equation for the pressure field with Neumann boundary conditions in a geometry as complicated as that of the ocean basins. A major objective of the study is to make this inversion, and hence nonhydrostatic ocean modeling, efficient on parallel computers. The pressure field is separated into surface, hydrostatic, and nonhydrostatic components. First, as in hydrostatic models, a two-dimensional problem is inverted for the surface pressure which is then made use of in the three-dimensional inversion for the nonhydrostatic pressure. Preconditioned conjugate-gradient iteration is used to invert symmetric elliptic operators in both two and three dimensions. Physically motivated preconditioners are designed which are efficient at reducing computation and minimizing communication between processors. Our method exploits the fact that as the horizontal scale of the motion becomes very much larger than the vertical scale, the motion becomes more and more hydrostatic and the three- dimensional Poisson operator becomes increasingly anisotropic and dominated by the vertical axis. Accordingly, a preconditioner is used which, in the hydrostatic limit, is an exact integral of the Poisson operator and so leads to a single algorithm that seamlessly moves from nonhydrostatic to hydrostatic limits. Thus in the hydrostatic limit the model is "fast," competitive with the fastest ocean climate models in use today based on the hydrostatic primitive equations. But as the resolution is increased, the model dynamics asymptote smoothly to the Navier Stokes equations and so can be used to address small- scale processes. A "finite-volume" approach is employed to discretize the model in space in which property fluxes are defined normal to faces that delineate the volumes. The method makes possible a novel treatment of the boundary in which cells abutting the bottom or coast may take on irregular shapes and be "shaved" to fit the boundary. The algorithm can conveniently exploit massively parallel computers and suggests a domain decomposition which allocates vertical columns of ocean to each processing unit. The resulting model, which can handle arbitrarily complex geometry, is efficient and scalable and has been mapped on to massively parallel multiprocessors such as the Connection Machine (CM5) using data-parallel FORTRAN and the Massachusetts Institute of Technology data-flow machine MONSOON using the implicitly parallel language Id. Details of the numerical implementation of a model which has been designed for the study of dynamical processes in the ocean from the convective, through the geostrophic eddy, up to global scale are set out. The "kernel" algorithm solves the incompressible Navier Stokes equations on the sphere, in a geometry as complicated as that of the ocean basins with ir- regular coastlines and islands. (Here we use the term "Navier Stokes" to signify that the full nonhydrostatic equations are being employed; it does not imply a particular constitutive relation. The relevant equations for modeling the full complex- ity of the ocean include, as here, active tracers such as tem- perature and salt.) It builds on ideas developed in the compu- tational fluid community. The numerical challenge is to ensure that the evolving velocity field remains nondivergent. Most

/pdf/a-finite-volume-incompressible-navier-stokes-model-for-2d0vnt2c1s.pdf

A finite-volume, incompressible Navier Stokes model for studies of the ocean on parallel computers

This article describes the governing equations, computational algorithms, and other components entering into the Community Multiscale Air Quality (CMAQ) modeling system. This system has been designed to approach air quality as a whole by including state-ofthe-science capabilities for modeling multiple air quality issues, including tropospheric ozone, fine particles, acid deposition, and visibility degradation. CMAQ was also designed to have multiscale capabilities so that separate models were not needed for urban and regional scale air quality modeling. By making CMAQ a modeling system that addresses multiple pollutants and different spatial scales, it has a “one-atmosphere” perspective that combines the efforts of the scientific community. To implement multiscale capabilities in CMAQ, several issues (such as scalable atmospheric dynamics and generalized coordinates), which depend on the desired model resolution, are addressed. A set of governing equations for compressible nonhydrostatic atmospheres is available to better resolve atmospheric dynamics at smaller scales. Because CMAQ is designed to handle scale-dependent meteorological formulations and a large amount of flexibility, its governing equations are expressed in a generalized coordinate system. This approach ensures consistency between CMAQ and the meteorological modeling system. The generalized coordinate system determines the necessary grid and coordinate transformations, and it can accommodate various vertical coordinates and map projections. The CMAQ modeling system simulates various chemical and physical processes that are thought to be important for understanding atmospheric trace gas transformations and distributions. The modeling system contains three types of modeling components (Models-3): a meteorological modeling system for the description of atmospheric states and motions, emission models for man-made and natural emissions that are injected into the atmosphere, and a chemistry-transport modeling system for simulation of the chemical transformation and fate. The chemical transport model includes the following process modules: horizontal advection, vertical advection, mass conservation adjustments for advection processes, horizontal diffusion, vertical diffusion, gas-phase chemical reactions and solvers, photolytic rate computation, aqueous-phase reactions and cloud mixing, aerosol dynamics, size distributions and chemistry, plume chemistry effects, and gas and aerosol deposition velocity estimation. This paper describes the Models-3 CMAQ system, its governing equations, important science algorithms, and a few application examples. This review article cites 114 references. DOI: 10.1115/1.2128636

Review of the Governing Equations, Computational Algorithms, and Other Components of the Models-3 Community Multiscale Air Quality (CMAQ) Modeling System

Parallel ocean general circulation modeling

The numerical algorithms used to simulate the advection, diffusion, sedimentation, coagulation, and condensational growth of atmospheric aerosols are described. The model can be used in one, two, or three spatial dimensions. The continuity equation in a generalized horizontal and vertical coordinate system is developed, which allows the model to be quickly adapted to a wide variety of dynamical models of global or regional scale. Algorithms are developed to treat the various physical processes, and the results of simulations are presented, which show the strengths and weaknesses of these algorithms. Although the emphasis is on the modeling of aerosols, the work is also applicable to the simulations of the transport of gases.

A multidimensional model for aerosols - Description of computational analogs

We present here results and analyses of a series of numerical experiments performed with a spectral general circulation model (GCM). The purpose of the GCM experiments is to examine the role of radiation/cloud processes in the general circulation of the troposphere and stratosphere. The experiments were primarily motivated by the significant improvements in the GCM zonal mean simulation as refinements were made in the model treatment of clear-sky radiation and cloud-radiative interactions. The GCM with the improved cloud/radiation model is able to reproduce many observed features, such as: a clear separation between the wintertime tropospheric jet and the polar night jet; winter polar stratospheric temperatures of about 200 K; interhemispheric and seasonal asymmetries in the zonal winds. In a set of sensitivity experiments, we have stripped the cloud/radiation model of its improvements, the result being a significant degradation of the zonal mean simulations by the GCM. Through these experiments ...

The Response of a Spectral General Circulation Model to Refinements in Radiative Processes

Results are presented from a high-resolution global ocean model that is driven through three decadal cycles of increasingly realistic prescribed atmospheric forcing from the period 1985–1995. The model used (the Parallel Ocean Program) is a z level primitive equation model with active thermohaline dynamics based on the formulation of Bryan [1969] rewritten for massively parallel computers. Improvements to the model include an implicit free-surface formulation of the barotropic mode [Dukowicz and Smith, 1994] and the use of pressure averaging for increasing the numerical time step. This study extends earlier 0.5° simulations of Semtner and Chervin [1992] to higher horizontal resolution with improved treatments of ocean geometry and surface forcing. The computational grid is a Mercator projection covering the global ocean from 77°N to 77°S and has 20 vertical levels. Three successive simulations have been performed on the CM-5 Connection Machine system at Los Alamos using forcing fields from the European Centre for Medium-Range Weather Forecasts (ECMWF). The first run uses monthly wind stresses for 1985–1995 and restoring of surface temperature and salinity to the Levitus [1982] seasonal climatology. The second run is the same but with 3 day-averaged rather than monthly averaged wind stress fields, and the third is the same as the second but uses the monthly climatological ECMWF heat fluxes of Barnier et al. [1995] instead of restoring to climatological sea surface temperatures. Many features of the wind-driven circulation are well represented in the model solutions, such as the overall current patterns, the numerous regions of hydrodynamic instability which correspond to those observed by satellite altimetry, and the filamented structure of the Antarctic Circumpolar Current. However, some features such as the separation points of the Gulf Stream and Kuroshio and the transport through narrow passages such as the Florida Straits are clearly inaccurate and indicate that still higher resolution may be required to correct these deficiencies. Water mass properties and some aspects of the thermohaline circulation are also not always well reproduced, which is partly due to the relatively short length of the integrations. The use of the ECMWF heat fluxes, rather than restoring to climatological surface temperatures, leads to stronger and more realistic surface and deep western boundary currents (primarily in the Atlantic) as well as more realistic meridional heat transport; this is primarily because the equilibrium meridional heat transport implied by the ECMWF surface fluxes is quite large. The ECMWF heat fluxes also produce improved seasonal cycles of sea surface temperature and height in both the northern and southern hemispheres. The 3-day wind forcing gives rise to modes of model variability that are clearly seen in synoptic observations, such as the large-scale 20–100-day oscillations seen in the TOPEX/POSEIDON data, which are barotropic oscillations induced by the high-frequency wind forcing. Additional studies on other aspects of the simulations described here are being conducted by a variety of investigators, and some of these are briefly described.

Global eddy-resolving ocean simulations driven by 1985-1995 atmospheric winds

We describe the results of January and July simulations carded out with a nine-level spectral model, employing a rhomboidal truncation at wavenumber 15. Sea-surface temperature, sea-ice distribution and solar zenith angle are held constant in each simulation. The model includes interactive clouds and radiative processes after Ramanathan et al. (1983). Selected fields are shown which highlight the model's strengths and weaknesses. The latitude-height distribution of the zonal wind is successfully simulated. The model captures the separation between the wintertime westerly jets in the troposphere and stratosphere and thus simulates the sign reversal in the vertical wind shear across the jet axis in the upper troposphere. In addition to the zonal wind, we show also the zonally averaged temperature, meridional wind and vertical velocity. Regional distributions of sea-level pressure, surface air temperature, precipitation and a number of other fields defined at various pressure levels are compared in ...

Robert C. Malone

Papers

Parallel ocean general circulation modeling

A multidimensional model for aerosols - Description of computational analogs

The Response of a Spectral General Circulation Model to Refinements in Radiative Processes

Global eddy-resolving ocean simulations driven by 1985-1995 atmospheric winds

January and July Simulations with a Spectral General Circulation Model