Intercomparison of Stratospheric Chemistry Models under Polar Vortex Conditions

Several stratospheric chemistry modules from box, 2-D or 3-D models, have been intercompared. The intercomparison was focused on the ozone loss and associated reactive species under the conditions found in the cold, wintertime Arctic and Antarctic vortices. Comparisons of both gas phase and heterogeneous chemistry modules show excellent agreement between the models under constrained conditions for photolysis and the microphysics of polar stratospheric clouds. While the mean integral ozone loss ranges from 4–80% for different 30–50 days long air parcel trajectories, the mean scatter of model results around these values is only about ±1.5%. In a case study, where the models employed their standard photolysis and microphysical schemes, the variation around the mean percentage ozone loss increases to about ±7%. This increased scatter of model results is mainly due to the different treatment of the PSC microphysics and heterogeneous chemistry in the models, whereby the most unrealistic assumptions about PSC processes consequently lead to the least representative ozone chemistry. Furthermore, for this case study the model results for the ozone mixing ratios at different altitudes were compared with a measured ozone profile to investigate the extent to which models reproduce the stratospheric ozone losses. It was found that mainly in the height range of strong ozone depletion all models underestimate the ozone loss by about a factor of two. This finding corroborates earlier studies and implies a general deficiency in our understanding of the stratospheric ozone loss chemistry rather than a specific problem related to a particular model simulation.     

A Case Study Of An On-Ice Air Flow Over The Arctic Marginal Sea-Ice Zone

A case study of warm air advection over the Arctic marginalsea-ice zone is presented, based on aircraft observations with direct flux measurements carriedout in early spring, 1998. A shallow atmospheric boundary layer (ABL) was observed, which wasgradually cooling with distance downwind of the ice edge. This process was mainly connected with astrong stable stratification and downward turbulent heat fluxes of about 10–20 W m^-2, but wasalso due to radiative cooling. Two mesoscale models, one hydrostatic and the other non-hydrostatic,having different turbulence closures, were applied. Despite these fundamental differences betweenthe models, the results of both agreed well with the observed data. Various closure assumptions had amore crucial influence on the results than the differences between the models.Such an assumption was, for example,the parameterization of the surface roughness for momentum (z₀) and heat (z_T). This stronglyaffected the wind and temperature fields not only close to the surface but also within and abovethe temperature inversion layer. The best results were achieved using a formulation for z₀ that took intoaccount the form drag effect of sea-ice ridges together withz_T = 0.1z₀. The stability within theelevated inversion strongly depended on the minimum eddy diffusivity K_min. A simple ad hocparameterization seems applicable, where K_min is calculated as 0.005 timesthe neutral eddy diffusivity. Although the longwave radiative cooling was largest within the ABL, theapplication of a radiation scheme was less important there than above the ABL. This was related to theinteraction of the turbulent and radiative fluxes. To reproduce the strong inversion, it wasnecessary to use vertical and horizontal resolutions higher than those applied in most regional andlarge-scale atmospheric models.     

Meso-gamma scale forecasts using the nonhydrostatic model LM

Summary ?The nonhydrostatic model LM was developed for small scale operational predictions. Advances in computer development will give the possibility of operational models of a rather fine scale, which will cover the meso-gamma scale. The LM is currently applied at a scale of 7 km and an increase of the operational resolution to 2.5 km is planned for the next few years. Predictions of such high resolution require to abandon the hydrostatic assumption, which is used with most current operational weather prediction models. The LM was designed to cover all resolutions from 50 m to 50 km with an efficiency making it suitable for operational use. It is a fully elastic model, using second order centred finite differences. The time integration is done using the Klemp–Wilhelmson method, treating the slow modes by a larger time step than the fast modes. The vertical propagation of the fast waves is done implicitly. After describing the design of the LM, this paper gives examples of model predictions at the meso-γ scale. Some results of the current operational application at the resolution 7 km are presented. Deficiencies in the localisation of model generated precipitation are investigated using an idealised bell shaped mountain and applying different resolutions. In this way the convergence to the correct solution can be investigated. From these results it is concluded, that orographic filtering is necessary and the effect of such filtering on precipitation forecasts is investigated. Finally, the prediction of a squall line over northern Germany is shown in order to demonstrate the potential of the model in forecasting the meso-γ scale. Received May 15, 2001; revised September 21, 2001     

Description and evaluation of the bergen climate model: ARPEGE coupled with MICOM

A new coupled atmosphere–ocean–sea ice model has been developed, named the Bergen Climate Model (BCM). It consists of the atmospheric model ARPEGE/IFS, together with a global version of the ocean model MICOM including a dynamic–thermodynamic sea ice model. The coupling between the two models uses the OASIS software package. The new model concept is described, and results from a 300-year control integration is evaluated against observational data. In BCM, both the atmosphere and the ocean components use grids which can be irregular and have non-matching coastlines. Much effort has been put into the development of optimal interpolation schemes between the models, in particular the non-trivial problem of flux conservation in the coastal areas. A flux adjustment technique has been applied to the heat and fresh-water fluxes. There is, however, a weak drift in global mean sea-surface temperature (SST) and sea-surface salinity (SSS) of respectively 0.1 °C and 0.02 psu per century. The model gives a realistic simulation of the radiation balance at the top-of-the-atmosphere, and the net surface fluxes of longwave, shortwave, and turbulent heat fluxes are within observed values. Both global and total zonal means of cloud cover and precipitation are fairly close to observations, and errors are mainly related to the strength and positioning of the Hadley cell. The mean sea-level pressure (SLP) is well simulated, and both the mean state and the interannual standard deviation show realistic features. The SST field is several degrees too cold in the equatorial upwelling area in the Pacific, and about 1 °C too warm along the eastern margins of the oceans, and in the polar regions. The deviation from Levitus salinity is typically 0.1 psu – 0.4 psu, with a tendency for positive anomalies in the Northern Hemisphere, and negative in the Southern Hemisphere. The sea-ice distribution is realistic, but with too thin ice in the Arctic Ocean and too small ice coverage in the Southern Ocean. These model deficiencies have a strong influence on the surface air temperatures in these regions. Horizontal oceanic mass transports are in the lower range of those observed. The strength of the meridional overturning in the Atlantic is 18 Sv. An analysis of the large-scale variability in the model climate reveals realistic El Niño – Southern Oscillation (ENSO) and North Atlantic–Arctic Oscillation (NAO/AO) characteristics in the SLP and surface temperatures, including spatial patterns, frequencies, and strength. While the NAO/AO spectrum is white in SLP and red in temperature, the ENSO spectrum shows an energy maximum near 3 years.     

Use of potential coefficient models for geoid undulation determinations using a spherical harmonic representation of the height anomaly/geoid undulation difference

 This paper suggests that potential coefficient models of the Earth's gravitational potential be used to calculate height anomalies which are then reduced to geoid undulations where such quantities are needed for orthometric height determination and vertical datum definition through a potential coefficient realization of the geoid. The process of the conversion of the height anomaly into a geoid undulation is represented by a height anomaly gradient term and the usual N–ζ term that is dependent on elevation and the Bouguer anomaly. Using a degree 360 expansion of 30′ elevations and the OSU91A potential coefficient model, a degree 360 representation of the correction terms was computed. The magnitude of N–ζ reached –3.4 m in the Himalaya Mountains with smaller, but still significant, magnitudes in other mountainous regions. Received: 6 May 1996; Accepted: 30 October 1996     

Results of three year observation with a superconducting gravimeter at the GeoForschungsZentrum Potsdam

The superconducting gravimeter (SG) TT70 has been continuously recording gravity data at the GeoForschungsZentrum (GFZ) Potsdam absolute gravity site since July 1992. The recorded data are edited and preprocessed by spike and step detection and elimination and gap filling. An atmospheric pressure correction is carried out on gravity data in the time domain with a complex admittance before tidal fitting. The atmospheric pressure admittance is calculated from tide free output of SG data and local atmospheric pressure using the cross spectral method. The ground water level admittance is determined by a single coefficient. Improvements with these corrections are shown in analysis results. New tidal parameters for Potsdam site are determined and compared with recordings of an Askania spring gravimeter at a nearby site. Deviations against the Wahr-Dehant-model are shown. Polar motion data of the IERS (International Earth Rotation Service, Paris) are used to calculate variations of centrifugal acceleration caused by polar motion (pole tide). These are compared with the corrected tide free output of SG series. For drift determination the polar motion correction is applied on SG data. The nutation period equivalent to the Earth's Nearly Diurnal Free Wobble is calculated from the SG data with a value of T_FCN = (437.4 ± 1.5) sidereal days. This result is compared with those obtained from other SG stations. Received 19 December 1995; Accepted 13 September 1996     

The hierarchical tessellation model and its use in spatial analysis

The hierarchid tessellation model belongs to a class of spatial data models based on the recursive decomposition of space. The quadtree is one such tessellation and is characterized by square cells and a 1:4 decomposition ratio. To relax these constraints in the tessellation, a generalized hierarchical tessellation data model, called Adaptive Recursive Tessellations (ART), has been proposed. ART increases flexibility in the tessellation by the use of rectangular cells and variable decomposition ratios. In ART, users can specify cell sizes which are intuitively meaningful to their applications, or which can reflect the scales of data. ART is implemented in a data structure called Adaptive Recursive Run-Encoding (ARRE), which is a variant of two-dimensional run-encoding whose running path can vary with the different tessellation structures incorporated in an ART model. Given the recognition of the benefits of implementing statistical spatial analysis in GIS, the use of hierarchical tessellation models such as ART in spatial analysis is discussed. Three examples are introduced to show how ART can: (1) be applied to solve the quadrat size problem in quadrat analysis of point patterns; (2) act as the data model in the variable resolution block kriging technique for geostatistical data to reduce variation in kriging error; and (3) facilitate the evaluation of spatial autocorrelation for area data at multiple map resolutions via the construction of a connectivity matrix for calculating spatial autocorrelation indices based on ARRE.     

Assessing the Resolving Power of Photographic Systems Used forTerrestrial Remote Sensing

It is often necessary to measure the resolving power of an imaging system comprising "off the shelf" components by a method that can also describe resolution in the object space. Furthermore, the results should allow useful comparisons to be made with alternative systems, where these have undergone a similar evaluation procedure. An example is given of measuring the resolving power of two imaging systems employing catadioptric objectives by a method that does not necessitate extensive laboratory tests, but utilizes the photographic product provided through field trials, or in the course of operational use. While the method reported does not claim to be a rigorous evaluation of the resolving power, it has been found to provide a convenient and practical guide to the image quality associated with a complete camera system, and its relation to the object space.     

Remote sensing of horticultural plantations in Kumarsain tehsil in Shimla district, Himachal Pradesh

The favourable agroclimatic conditions for orchards especially apples have increased the acreages in Himachal Pradesh (HP) which has significantly contributed in the growth of state economy. Realizing the importance of horticulture in HP and its changing scenario of the land use/land cover, a study was conducted to identify and map apple and almond plantations in the Kumarsain tehsil of Shimla district using Remote Sensing (RS) techniques. IRS-IB LISS-II False Colour Composite (FCC) diapositives of October 27, March 30 and April 20, 1992 were visually analysed for mapping apple and almond plantations. The results indicate that IRS LISS-II data of April 20 on 1∶50,000 scale was found very useful for identification and mapping of apple and almond plantations in this region. Accuracy of interpretation was also tested on sample basis assuming a binomial distribution for the probability of success/failure of sample points. The overall interpretation accuracy assessed based on 40 sample points was found to be 87 per cent at 90 per cent confidence limits.