Numerical simulation of a South China Sea typhoon Leo (1999)

Summary ?A South China Sea typhoon, Leo (1999), was simulated using the Penn State/NCAR mesoscale model MM5 with the Betts-Miller convective parameterization scheme (BMEX). The simulation had two nested domains with resolutions at 54 and 18 km, and the forecast duration was 36 hours. The model was quite successful in predicting the track, the rapid deepening, the central pressure, and the maximum wind speed of typhoon Leo as verified with reports from the Hong Kong Observatory (HKO). The structure of the eye, the eye wall, and the spiral convective cloud band simulated in the model are found to be comparable to corresponding features identified in satellite images for the storm, and also with those reported by other authors. A trajectory analysis was performed. Three kinds of trajectory were found: (1) spirally rising trajectories near the eye wall; (2) spirally rising/descending trajectories in the convective/cloud free belt; (3) straight and fast rising trajectories in a heavy convection zone along one of the cloud bands on the periphery of the tropical cyclone. Both the HKO and the U.S. Joint Typhoon Warning Center (JTWC) reported the rapid deepening of Leo started around 00 UTC 29 April. In the model, the eye was first formed in the lower troposphere, and it extended to the upper troposphere within a few hours. We speculate that the spin-up of cyclonic rotation in the low-level eye enhanced the positive vorticity along the low-level eye wall. The positive vorticity was then transported to the upper troposphere by convection, leading to an extension and growth of the eye into the upper troposphere. To examine the impact of convective parameterization scheme (CPS) on the simulation, the Grell scheme (GLEX) was also tested. The GLEX predicted a weaker typhoon with a wilder eye that extended not as high up in the upper troposphere as BMEX. The different structures of the eye between the BMEX and GLEX suggest that the mesoscale features of the eye are dependent on the convection. In other words, the vertical and horizontal distribution of convective heating is essential to the development and structure of the eye. Received December 18, 2001; accepted May 7, 2002 Published online: March 20, 2003     

A Laboratory and Theoretical Study on the Uptake of SO2 Gas by Large and Small Water Drops Containing Heavy Metal Ions

Laboratory experiments have been carried out to investigate the uptake of sulfur dioxide by water drops containing heavy metal ions where the metal ions serve as catalysts to oxidise the taken up S(IV) into S(VI). During the gas uptake the drops were freely suspended at their terminal velocity in the airstream of the Mainz vertical wind tunnel. Two series of experiments were carried out, one with large millimeter size water drops containing manganese or iron ions, and the other with small water drops containing manganese ions and having radii in hundreds of micron size range. The experimental results were compared against model computations using the Kronig–Brink model and the fully mixed model, modified for the case that heavy metal ions present in the liquid phase act as catalysts for the oxidising process. The results of the model calculations show that there are only small differences between the predicted gas uptake according to the two models. In addition it was found that the experimental obtained results from the uptake of SO₂ by water drops containing heavy metal ions for both, large and small water drops did agree with the model results.     

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.     

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.