Engineering topics

Swinbank, Richard

Publications and source records attributed to Swinbank, Richard.

13 records

NASA NTRS2006

Fibonacci Grids

Recent years have seen a resurgence of interest in a variety of non-standard computational grids for global numerical prediction. The motivation has been to reduce problems associated with the converging meridians and the polar singularities of conventional regular latitude-longitude grids. A further impetus has come from the adoption of massively parallel computers, for which it is necessary to distribute work equitably across the processors; this is more practicable for some non-standard grids. Desirable attributes of a grid for high-order spatial finite differencing are: (i) geometrical regularity; (ii) a homogeneous and approximately isotropic spatial resolution; (iii) a low proportion of the grid points where the numerical procedures require special customization (such as near coordinate singularities or grid edges). One family of grid arrangements which, to our knowledge, has never before been applied to numerical weather prediction, but which appears to offer several technical advantages, are what we shall refer to as "Fibonacci grids". They can be thought of as mathematically ideal generalizations of the patterns occurring naturally in the spiral arrangements of seeds and fruit found in sunflower heads and pineapples (to give two of the many botanical examples). These grids possess virtually uniform and highly isotropic resolution, with an equal area for each grid point. There are only two compact singular regions on a sphere that require customized numerics. We demonstrate the practicality of these grids in shallow water simulations, and discuss the prospects for efficiently using these frameworks in three-dimensional semi-implicit and semi-Lagrangian weather prediction or climate models.

Swinbank, Richard↗

NASA NTRS2005

Diagnostic Comparison of Meteorological Analyses during the 2002 Antarctic Winter

Several meteorological datasets, including U.K. Met Office (MetO), European Centre for Medium-Range Weather Forecasts (ECMWF), National Centers for Environmental Prediction (NCEP), and NASA's Goddard Earth Observation System (GEOS-4) analyses, are being used in studies of the 2002 Southern Hemisphere (SH) stratospheric winter and Antarctic major warming. Diagnostics are compared to assess how these studies may be affected by the meteorological data used. While the overall structure and evolution of temperatures, winds, and wave diagnostics in the different analyses provide a consistent picture of the large-scale dynamics of the SH 2002 winter, several significant differences may affect detailed studies. The NCEP-NCAR reanalysis (REAN) and NCEP-Department of Energy (DOE) reanalysis-2 (REAN-2) datasets are not recommended for detailed studies, especially those related to polar processing, because of lower-stratospheric temperature biases that result in underestimates of polar processing potential, and because their winds and wave diagnostics show increasing differences from other analyses between similar to 30 and 10 hPa (their top level). Southern Hemisphere polar stratospheric temperatures in the ECMWF 40-Yr Re-analysis (ERA-40) show unrealistic vertical structure, so this long-term reanalysis is also unsuited for quantitative studies. The NCEP/Climate Prediction Center (CPC) objective analyses give an inferior representation of the upper-stratospheric vortex. Polar vortex transport barriers are similar in all analyses, but there is large variation in the amount, patterns, and timing of mixing, even among the operational assimilated datasets (ECMWF, MetO, and GEOS-4). The higher-resolution GEOS-4 and ECMWF assimilations provide significantly better representation of filamentation and small-scale structure than the other analyses, even when fields gridded at reduced resolution are studied. The choice of which analysis to use is most critical for detailed transport studies (including polar process modeling) and studies involving synoptic evolution in the upper stratosphere. The operational assimilated datasets are better suited for most applications than the NCEP/CPC objective analyses and the reanalysis datasets.

stratosphere↗

NASA NTRS2005

Simulations of Dynamics and Transport during the September 2002 Antarctic Major Warming

A mechanistic model simulation initialized on 14 September 2002, forced by 100-hPa geopotential heights from Met Office analyses, reproduced the dynamical features of the 2002 Antarctic major warming. The vortex split on approx.25 September; recovery after the warming, westward and equatorward tilting vortices, and strong baroclinic zones in temperature associated with a dipole pattern of upward and downward vertical velocities were all captured in the simulation. Model results and analyses show a pattern of strong upward wave propagation throughout the warming, with zonal wind deceleration throughout the stratosphere at high latitudes before the vortex split, continuing in the middle and upper stratosphere and spreading to lower latitudes after the split. Three-dimensional Eliassen-Palm fluxes show the largest upward and poleward wave propagation in the 0(deg)-90(deg)E sector prior to the vortex split (coincident with the location of strongest cyclogenesis at the model's lower boundary), with an additional region of strong upward propagation developing near 180(deg)-270(deg)E. These characteristics are similar to those of Arctic wave-2 major warmings, except that during this warming, the vortex did not split below approx.600 K. The effects of poleward transport and mixing dominate modeled trace gas evolution through most of the mid- to high-latitude stratosphere, with a core region in the lower-stratospheric vortex where enhanced descent dominates and the vortex remains isolated. Strongly tilted vortices led to low-latitude air overlying vortex air, resulting in highly unusual trace gas profiles. Simulations driven with several meteorological datasets reproduced the major warming, but in others, stronger latitudinal gradients at high latitudes at the model boundary resulted in simulations without a complete vortex split in the midstratosphere. Numerous tests indicate very high sensitivity to the boundary fields, especially the wave-2 amplitude. Major warmings occurred for initial fields with stronger winds and larger vortices, but not smaller vortices, consistent with the initiation of wind-deceleration by upward-propagating waves near the poleward edge of the region where wave 2 can propagate above the jet core. Thus, given the observed 100-hPa boundary forcing, stratospheric preconditioning is not needed to reproduce a major warming similar to that observed. The anomalously strong forcing in the lower stratosphere can be viewed as the primary direct cause of the major warming.

winds↗

NASA NTRS2003

The SPARC Intercomparison of Middle Atmosphere Climatologies

Our current confidence in 'observed' climatological winds and temperatures in the middle atmosphere (over altitudes approx. 10-80 km) is assessed by detailed intercomparisons of contemporary and historic data sets. These data sets include global meteorological analyses and assimilations, climatologies derived from research satellite measurements, and historical reference atmosphere circulation statistics. We also include comparisons with historical rocketsonde wind and temperature data, and with more recent lidar temperature measurements. The comparisons focus on a few basic circulation statistics, such as temperature, zonal wind, and eddy flux statistics. Special attention is focused on tropical winds and temperatures, where large differences exist among separate analyses. Assimilated data sets provide the most realistic tropical variability, but substantial differences exist among current schemes.

Randel, William↗

NASA NTRS2002

SPARC Intercomparison of Middle Atmosphere Climatologies

This atlas presents detailed incomparisons of several climatological wind and temperature data sets which cover the middle atmosphere (over altitudes approx. 10-80 km). A number of middle atmosphere climatologies have been developed in the research community based on a variety of meteorological analyses and satellite data sets. Here we present comparisons between these climatological data sets for a number of basic circulation statistics, such as zonal mean temperature, winds and eddy flux statistics. Special attention is focused on tropical winds and temperatures, where large differences exist among separate analyses. We also include comparisons between the global climatologies and historical rocketsonde wind and temperature measurements, and also with more recent lidar temperature data. These comparisons highlight differences and uncertainties in contemporary middle atmosphere data sets, and allow biases in particular analyses to be isolated. In addition, a brief atlas of zonal mean temperature and wind statistics is provided to highlight data availability and as a quick-look reference. This technical report is intended as a companion to the climatological data sets held in archive at the SPARC Data Center (http://www.sparc.sunysb.edu).

Randel, William↗

NASA NTRS2001

Lower Stratospheric Temperature Differences Between Meteorological Analyses in two cold Arctic Winters and their Impact on Polar Processing Studies

A quantitative intercomparison of six meteorological analyses is presented for the cold 1999-2000 and 1995-1996 Arctic winters. The impacts of using different analyzed temperatures in calculations of polar stratospheric cloud (PSC) formation potential, and of different winds in idealized trajectory-based temperature histories, are substantial. The area with temperatures below a PSC formation threshold commonly varies by approximately 25% among the analyses, with differences of over 50% at some times/locations. Freie University at Berlin analyses are often colder than others at T is less than or approximately 205 K. Biases between analyses vary from year to year; in January 2000. U.K. Met Office analyses were coldest and National Centers for Environmental Prediction (NCEP) analyses warmest. while NCEP analyses were usually coldest in 1995-1996 and Met Office or NCEP[National Center for Atmospheric Research Reanalysis (REAN) warmest. European Centre for Medium Range Weather Forecasting (ECMWF) temperatures agreed better with other analyses in 1999-2000, after improvements in the assimilation model. than in 1995-1996. Case-studies of temperature histories show substantial differences using Met Office, NCEP, REAN and NASA Data Assimilation Office (DAO) analyses. In January 2000 (when a large cold region was centered in the polar vortex), qualitatively similar results were obtained for all analyses. However, in February 2000 (a much warmer period) and in January and February 1996 (comparably cold to January 2000 but with large cold regions near the polar vortex edge), distributions of "potential PSC lifetimes" and total time spent below a PSC formation threshold varied significantly among the analyses. Largest peaks in "PSC lifetime" distributions in January 2000 were at 4-6 and 11-14 days. while in the 1996 periods, they were at 1-3 days. Thus different meteorological conditions in comparably cold winters had a large impact on expectations for PSC formation and on the discrepancies between different meteorological analyses. Met Office. NCEP, REAN, ECMWF and DAO analyses are commonly used for trajectory calculations and in chemical transport models; the choice of which analysis to use can strongly influence the results of such studies.

Manney, Gloria L.↗

NASA NTRS1999

Data Structure and Parallel Decomposition Considerations on a Fibonacci Grid

The Fibonacci grid, proposed by Swinbank and Purser (see companion abstract), provides attractive properties for global numerical atmospheric prediction by offering an optimally homogeneous, geometrically regular, and approximately isotropic discretization, with only the polar regions requiring special numerical treatment. It is a mathematical idealization, applied to the sphere, of the multi-spiral patterns often found in botanical structures, such as in pine cones and sunflower heads. Computationally, it is natural to organize the domain, into zones, in each of which the same pair, or triple, of "Fibonacci spirals" dominate. But the further subdivision of such zones into "tiles" of a shape and size suitable for distribution to the processors of a massively parallel computer requires very careful consideration if the subsequent spatial computations along the respective spirals, especially those computations (such as compact differencing schemes) that involve recursion, can be implemented in an efficient "load-balanced "manner without requiring excessive amounts of inter-processor communications. In this paper we show how certain "number theoretic" properties of the Fibonacci sequence (whose numbers prescribe the multiplicity of successive spirals) may be exploited in the decomposition of grid zones into tidy arrangements of triangular grid tiles, each tile possessing one side approximately parallel to the constant-latitude zone boundary. We also describe how the spatially recursive processes may be decomposed across such a tiling, and the directionality of the recursions reversed on alternate grid lines, to ensure a very high degree of load balancing throughout the execution of the computations required for one time step of a global model.

Michalakes, John↗

NASA NTRS1995

A quasi-permeable barrier in the equatorial lower stratosphere

Analyses of global isentropic maps of potential vorticity from the United Kingdom Meteorological Office (UKMO) data assimilation system reveal that there is a region of strong meridional potential vorticity (PV) gradients in the equatorial lower stratosphere. A semi-Lagrangian tracer transport model is used to show that this region of strong PV gradients acts as a quasipermeable barrier to cross-equator mass exchange on isentropic surfaces.

Chen, Ping↗

NASA NTRS1994

Evolution of the stratosphere during northern winter 1991/92 as diagnosed from U.K. Meteorological Office analyses

Meteorological analyses, produced at the U.K. Meteorological Office by data assimilation, are used to study the circulation of the stratosphere in the Northern Hemisphere during winter 1991/92. The analyses are supplemented by Lagrangian visualizations of the circulation. The main features discussed are (1) the changes in vertical structure of the circulation, (2) the merger of anticyclones that precipitated a strong stratospheric warming, (3) vortex roll up in the upper stratosphere, (4) the entrainment of air into the polar vortex in the middle and upper stratosphere, and (5) the influence of tropospheric blocking on the lower stratosphere. The study provides a meteorological basis for the interpretation of data from the Upper Atmosphere Research Satellite (UARS).

O'Neill, Alan↗

NASA NTRS1994

Spring dehydration in the Antarctic stratospheric vortex observed by HALOE

The distribution of dehydrated air in the middle and lower stratosphere during the 1992 Southern Hemisphere spring is investigated using Halogen Occultation Experiment (HALOE) observations and trajectory techniques. Comparisons between previously published Version 9 and the improved Version 16 retrievals on the 700-K isentropic surface show very slight (0.05 ppmv) increases in Version 16 CH4 relative to Version 9 within the polar vortex. Version 16 H2O mixing ratios show a reduction of 0.5 ppmv relative to Version 9 within the polar night jet and a reduction of nearly 1.0 ppmv in middle latitudes when compared to Version 9. The version 16 HALOE retrievals show low mixing ratios of total hydrogen (2CH4 + H2O) within the polar vortex on both 700 and 425 K isentropic surfaces relative to typical middle-stratospheric 2CH4 + H2O mixing ratios. The low 2CH4 + H2O mixing ratios are associated with dehydration. Slight reductions in total hydrogen, relative to typical middle-stratospheric values, are found at these levels throughout the Southern Hemisphere during this period. Trajectory calculations show that middle-latitude air masses are composed of a mixture of air from within the polar night jet and air from middle latitudes. A strong kinematic barrier to large-scale exchange is found on the poleward flank of the polar night jet at 700 K. A much weaker kinematic barrier is found at 425 K. The impact of the finite tangent pathlength of the HALOE measurements is investigated using an idealized tracer distribution. This experiment suggests that HALOE should be able to resolve the kinematic barrier, if it exists.

Pierce, R. Bradley↗

NASA NTRS1994

Mixing processes within the polar night jet

Lagrangian material line simulations are performed using U.K. Meteorological Office simulated winds and temperatures to examine mixing processes in the middle- and lower-stratospheric polar night jet during the 1992 Southern Hemisphere spring and Northern Hemisphere winter. The Lagrangian simulations are undertaken to provide insight into the effects of mixing within the polar night jet on observations of the polar vortex made by instruments onboard the Upper Atmosphere Research Satellite (UARS) during these periods. A moderate to strong kinematic barrier to large-scale isentropic exchange, similar to the barrier identified in General Circulation Model (GCM) simulations, is identified during both of these periods. Characteristic timescales for mixing by large-scale isentropic motions within the polar night jet range from 20 days in the Southern Hemisphere lower stratosphere to years in the Northern Hemisphere middle stratosphere. The long mixing timescales found in the Northern Hemisphere polar night jet do not persist. Instead, the Northern Hemisphere kinematic barriers are broken down as part of the large-scale stratospheric response to a strong tropospheric blocking event. A series of Lagrangian experiments are conducted to investigate the sensitivity of the kinematic barrier to diabatic effects and to small-scale inertial gravity wave motions. Differential diabatic descent is found to have a significant impact on mixing processes within the Southern Hemisphere middle-stratospheric jet core. The interaction between small-scale displacements by idealized, inertial gravity waves and the large-scale flow is found to have a significant impact on mixing within the polar night jet in both hemispheres. These sensitivity experiments suggest that scales of motion that are unresolved in global assimilated datasets may contribute to mass exchange across the kinematic barrier to large-scale isentropic motion.

Pierce, R. Bradley↗

NASA NTRS1994

Isentropic mass exchange between the Tropics and extratropics in the stratosphere

The isentropic mass exchange between the Tropics and extratropics in the stratosphere is investigated with a semi-Lagrangian transport model for the periods from 1 June to 31 October 1992 and from 1 December 1992 to 30 April 1993 using winds from the U.K. Meteorological Office data assimilation system. Calculations with an idealized, initially zonally symmetric tracer show that in the middle and upper stratosphere the bulk of tropical air is transported into the midlatitudes of the winter hemisphere although there exist quasi-permeable barriers in the subtropics. The transport takes place in the form of planetary-scale 'tongues' of material that are drawn poleward in association with the episodic amplification of planetary-scale waves in high latitudes of the winter hemisphere. Once air of tropical origin is transported to the midlatitudes it is irreversibly mixed with the midlatitude air in the 'surf zone.' Air of tropical origin can, however, hardly penetrate into the interior of the winter polar vortex until the breakdown of the vortex. Transport of tropical air into the midlatitudes of the summer hemisphere is strongly inhibited. In the lower stratosphere, tropical air is transported into the northern and southern midlatitudes. During the period from 1 June to 31 October 1992, the amount of tropical air transported into the Northern Hemisphere is, however, much smaller than that transported into the Southern Hemisphere, and there exist strong gradients in the tracer field in the equatorial region, indicating that there is a quasi-permeable barrier to cross-equator mass exchange. During the period from 1 December 1992 to 30 April 1993, on the other hand, roughly the same amounts of tropical air are transported into the Southern Hemisphere, and there exist strong gradients in the tracer field in the equatorial region, indicating that there is a quasi-permeable barrier to cross-equator mass exchange. During the period from 1 December 1992 to 30 April 1993, on the other hand, roughly the same amounts of tropical air are transported into the northern midlatitudes and into the southern midlatitudes, and there exist no significant transport barriers in the equatorial area.

Chen, Ping↗

NASA NTRS1994

Quasi-horizontal transport and mixing in the Antarctic stratosphre

The quasi-horizontal transport and mixing properties of the Antarctic stratosphere are investigated with a simi-Lagrangian transport model and a 'contour advection' technique for the winter and spring of 1992 using analyzed winds from the United Kingdom Meteorological Office data assimiliation system. Transport calculations show that passive tracers are well mixed inside the polar vortex as well as in the midlatitude 'surf zone.' A the vortex edge, strong radial gradients in the tracer fields are well preserved, and their evolutions follow that of the potential vorticity until some time after the breakdown of the polar vortex. In the middle stratosphere there is little tracer exchange across the vortex edge in August and September. Some vortex air is eroded into the surf zone in filamentary form in October, and very strong exchange of air occurs between high and middle latitudes in November. In the lower stratosphere the vortex is not so isolated from the midlatitudes as in the middle stratosphere, and there is more mass exchange across the vortex edge. Calculations of the lengthening of material contours using the contour advection technique show that in the middle stratosphere, strong stirring (i.e., stretching and folding of material elements) occurs in the inner vortex, while the strongest stirring occurs in the midlatitude surf zone and the weakest occurs at the vortex edge. In the lower strtosphere, strong stirring occurs in the inner vortex. Stirring is moderate at the vortex edge and in the midlatitudes.

Chen, Ping↗