Mountain wave motions determined by the Esrange MST radar

A European campaign of ground-based radar, lidar and optical measurements was carried out during the winter of 1996/1997 (28 December–2 February) to study lee waves in the northern part of Scandinavia. The participants operated ozone lidars, backscatter lidars and MST radars at ALOMAR/Andoya and Esrange/Kiruna, and an ALIS imaging system in Kiruna. The Andoya site was generally windward of the Scandinavian mountains, the Kiruna site on the leeward side. The goal of the experiment was to examine the influence of lee waves on the formation of Polar Stratospheric Clouds (PSCs). This paper studies the radar data from MST-radar ESRAD located at Esrange [68.°N, 21.°E], i.e. in the lee of the mountains. We present three cases where strong lee waves were observed: in one case they propagated upwards to the lower stratosphere and in the other two cases they were trapped or absorbed in the troposphere. We examine the local waves and the direction and strength of the local wind using the radar, the synoptic meteorological situation using weather maps (European Meteorological Bulletin) and the synoptic stratospheric temperatures using ECMWF data. We observed that waves propagate up to the stratosphere during frontal passages. When anticyclonic ridges are present, the propagation to the stratosphere is very weak. This is due to trapping of the waves at or below the tropopause. We also show that the radar data alone can be used to characterise the different weather conditions for the three cases studied (through the variation of the height of the tropopause). The synoptic stratospheric temperatures in the three cases were similar, and were above the expected threshold for PSC formation. Lidar and visual observation of PSCs and nacreous clouds, respectively, showed that these were present only in the case when the lee waves propagated up to the lower stratosphere.     

Radiative forcing by contrails

A parametric study of the instantaneous radiative impact of contrails is presented using three different radiative transfer models for a series of model atmospheres and cloud parameters. Contrails are treated as geometrically and optically thin plane parallel homogeneous cirrus layers in a static atmosphere. The ice water content is varied as a function of ambient temperature. The model atmospheres include tropical, mid-latitude, and subarctic summer and winter atmospheres. Optically thin contrails cause a positive net forcing at top of the atmosphere. At the surface the radiative forcing is negative during daytime. The forcing increases with the optical depth and the amount of contrail cover. At the top of the atmosphere, a mean contrail cover of 0.1% with average optical depth of 0.2 to 0.5 causes about 0.01 to 0.03 Wm⁻² daily mean instantaneous radiative forcing. Contrails cool the surface during the day and heat the surface during the night, and hence reduce the daily temperature amplitude. The net effect depends strongly on the daily variation of contrail cloud cover. The indirect radiative forcing due to particle changes in natural cirrus clouds may be of the same magnitude as the direct one due to additional cover.     

The role of sodium bicarbonate in the nucleation of noctilucent clouds

It is proposed that a component of meteoric smoke, sodium bicarbonate (NaHCO₃), provides particularly effective condensation nuclei for noctilucent clouds. This assertion is based on three conditions being met. The first is that NaHCO₃ is present at sufficient concentration (10⁴ cm^–3) in the upper mesosphere between 80 and 90 km. It is demonstrated that there is strong evidence for this based on recent laboratory measurements coupled with atmospheric modelling. The second condition is that the thermodynamics of NaHCO₃(H₂O)_n cluster formation allow spontaneous nucleation to occur under mesospheric conditions at temperatures below 140 K. The Gibbs free energy changes for forming clusters with n = 1 and 2 were computed from quantum calculations using hybrid density functional/Hartree-Fock (B3LYP) theory and a large basis set with added polarization and diffuse functions. The results were then extrapolated to higher n using an established dependence of the free energy on cluster size and the free energy for the sublimation of H₂O to bulk ice. A 1-dimensional model of sodium chemistry was then employed to show that spontaneous nucleation to form ice particles (n > 100) should occur between 84 and 89 km in the high-latitude summer mesosphere. The third condition is that other metallic components of meteoric smoke are less effective condensation nuclei, so that the total number of potential nuclei is small relative to the amount of available H₂O. Quantum calculations indicate that this is probably the case for major constituents such as Fe(OH)₂, FeO₃ and MgCO₃.     

Letter to the editor A strange cloud in the Arctic summer stratosphere 1998 above Esrange (68°N), Sweden

When the University of Bonn lidar on the Esrange (68°N, 21°E), Sweden, was switched on in the evening of July 18, 1998, a geometrically and optically thin cloud layer was present near 14 km altitude or 400 K potential temperature, where it persisted for two hours. The tropopause altitude was 4 km below the cloud altitude. The cloud particles depolarized the lidar returns, thus must they have been aspherical and hence solid. Atmospheric temperatures near 230 K were approximately 40 K too high to support ice particles at stratospheric water vapour pressures of a few ppmv. The isentropic back trajectory on 400 K showed the air parcels to have stayed clear of active major rocket launch sites. The air parcels at 400 K had travelled from the Aleutians across Canada and the Atlantic Ocean arriving above central Europe and then turned northward to pass over above the lidar station. Parcels at levels at ±25 K from 400 K had come from the pole and joined the 400 K trajectory path above eastern Canada. Apparently the cloud existed in a filament of air with an origin different from those filaments both above and below. Possibly the 400 K level air parcels had carried soot particles from forest wild fires in northern Canada or volcanic ash from the eruption of the Korovin Volcano in the Aleutian Islands.     

Effects of geomagnetic activity on the mesospheric electric fields

The results of three series of rocket measurements of mesospheric electric fields carried out under different geomagnetic conditions at polar and high middle latitudes are analysed. The measurements show a clear dependence of the vertical electric fields on geomagnetic activity at polar and high middle latitudes. The vertical electric fields in the lower mesosphere increase with the increase of geomagnetic indexes K_p and K_p. The simultaneous increase of the vertical electric field strength and ion conductivity was observed in the mesosphere during geomagnetic disturbances. This striking phenomenon was displayed most clearly during the solar proton events of October, 1989 accompanied by very strong geomagnetic storm (K_p = 8+). A possible mechanism of generation of the vertical electric fields in the mesosphere caused by gravitational sedimentation of charged aerosol particles is discussed. Simultaneous existence in the mesosphere of both the negative and positive multiply charged aerosol particles of different sizes is assumed for explanation of the observed V/m vertical electric fields and their behaviour under geomagnetically disturbed conditions.Paper Presented at the Second IAGA/ICMA (IAMAS) Workshop on Solar Activity Forcing of the Middle Atmosphere, Prague, August 1997     

Meso-scale effects of tropical deforestation in Amazonia: preparatory LBA modelling studies

As part of the preparation for the Large-Scale Biosphere Atmosphere Experiment in Amazonia, a meso-scale modelling study was executed to highlight deficiencies in the current understanding of land surface atmosphere interaction at local to sub-continental scales in the dry season. Meso-scale models were run in 1-D and 3-D mode for the area of Rondonia State, Brazil. The important conclusions are that without calibration it is difficult to model the energy partitioning of pasture; modelling that of forest is easier due to the absence of a strong moisture deficit signal. The simulation of the boundary layer above forest is good, above deforested areas (pasture) poor. The models’ underestimate of the temperature of the boundary layer is likely to be caused by the neglect of the radiative effects of aerosols caused by biomass burning, but other factors such as lack of sufficient entrainment in the model at the mixed layer top may also contribute. The Andes generate patterns of subsidence and gravity waves, the effects of which are felt far into the Rondonian area The results show that the picture presented by GCM modelling studies may need to be balanced by an increased understanding of what happens at the meso-scale. The results are used to identify key measurements for the LBA atmospheric meso-scale campaign needed to improve the model simulations. Similar modelling studies are proposed for the wet season in Rondonia, when convection plays a major role.     

Effect of remote clouds on surface UV irradiance

Clouds affect local surface UV irradiance, even if the horizontal distance from the radiation observation site amounts to several kilometers. In order to investigate this effect, which we call remote clouds effect, a 3-dimensional radiative transfer model is applied. Assuming the atmosphere is subdivided into a quadratic based sector and its surrounding, we quantify the influence of changing cloud coverage within this surrounding from 0% to 100% on surface UV irradiance at the sector center. To work out this remote clouds influence as a function of sector base size, we made some calculations for different sizes between 10 km × 10 km and 100 km × 100 km. It appears that in the case of small sectors (base size 20 km × 20 km) the remote clouds effect is highly variable: Depending on cloud structure, solar zenith angle and wavelength, the surface UV irradiance may be enhanced up to 15% as well as reduced by more than 50%. In contrast, for larger sectors it is always the case that enhancements become smaller by 5% if sector base size exceeds 60 km × 60 km. However, these values are upper estimates of the remote cloud effects and they are found only for special cloud structures. Since these structures might occur but cannot be regarded as typical, different satellite observed cloud formations (horizontal resolution about 1 km × 1 km) have also been investigated. For these more common cloud distributions we find remote cloud effects to be distinctly smaller than the corresponding upper estimates, e.g., for a sector with base size of 25 km × 25 km the surface UV irradiance error due to ignoring the actual remote clouds and replacing their influence with periodic horizontal boundary conditions is less than 3%, whereas the upper estimate of remote clouds effect would suggest an error close to 10%.     

The present-day and future impact of NOx emissions from subsonic aircraft on the atmosphere in relation to the impact of NOx surface sources

The effect of present-day and future NO_x emissions from aircraft on the NO_x and ozone concentrations in the atmosphere and the corresponding radiative forcing were studied using a three-dimensional chemistry transport model (CTM) and a radiative model. The effects of the aircraft emissions were compared with the effects of the three most important anthropogenic NO_x surface sources: road traffic, electricity generation and industrial combustion. From the model results, NO_x emissions from aircraft are seen to cause an increase in the NO_x and ozone concentrations in the upper troposphere and lower stratosphere, and a positive radiative forcing. For the reference year 1990, the aircraft emissions result in an increase in the NO_x concentration at 250 hPa of about 20 ppt in January and 50 ppt in July over the eastern USA, the North Atlantic Flight Corridor and Western Europe, corresponding to a relative increase of about 50%. The maximum increase in the ozone concentrations due to the aircraft emissions is about 3-4 ppb in July over the northern mid-latitudes, corresponding to a relative increase of about 3-4%. The aircraft-induced ozone changes cause a global average radiative forcing of 0.025 W/m² in July. According to the ANCAT projection for the year 2015, the aircraft NO_x emissions in that year will be 90% higher than in the year 1990. As a consequence of this, the calculated NO_x perturbation by aircraft emissions increases by about 90% between 1990 and 2015, and the ozone perturbation by about 50-70%. The global average radiative forcing due to the aircraft-induced ozone changes increases by about 50% between 1990 and 2015. In the year 2015, the effects of the aircraft emissions on the ozone burden and radiative forcing are clearly larger than the individual effects of the NO_x surface sources. Taking chemical conversion in the aircraft plume into account in the CTM explicitly, by means of modified aircraft NO_x emissions, a significant reduction of the aircraft-induced NO_x and ozone perturbations is realised. The NO_x perturbation decreases by about 40% and the ozone perturbation by about 30% in July over Western Europe, the eastern USA and the North Atlantic Flight Corridor.     

Use of coupled ozone fields in a 3-D circulation model of the middle atmosphere

With a detailed chemistry scheme for the middle atmosphere up to 70 km which has been added to the 3-D Karlsruhe simulation model of the middle atmosphere (KASIMA), the effects of coupling chemistry and dynamics through ozone are studied for the middle atmosphere. An uncoupled version using an ozone climatology for determining heating rates and a coupled version using on-line ozone are compared in a 10-month integration with meteorological analyses for the winter 1992/93 as the lower boundary condition. Both versions simulate the meteorological situation satisfactorily, but exhibit a too cold lower stratosphere. The on-line ozone differs from the climatological data between 20 and 40 km by exhibiting too high ozone values, whereas in the lower mesosphere the ozone values are too low. The coupled model version is stable and differs only above 40 km significantly from the uncoupled version. Direct heating effects are identified to cause most of the differences. The well-known negative correlation between temperature and ozone is reproduced in the model. As a result, the coupled version slightly approaches the climatological ozone field. Further feedback effects are studied by using the on-line ozone field as a basis for an artificial climatology. For non-disturbed ozone conditions realistic monthly and zonally averaged ozone data are sufficient to determine the heating rates for modelling the middle atmosphere.     

The role of dynamics in total ozone deviations from their long-term mean over the Northern Hemisphere

Total ozone anomalies (deviation from the long-term mean) are created by anomalous circulation patterns. The dynamically produced ozone anomalies can be estimated from known circulation parameters in the layer between the tropopause and the middle stratosphere by means of statistics. Satellite observations of ozone anomalies can be compared with those expected from dynamics. Residual negative anomalies may be due to chemical ozone destruction. The statistics are derived from a 14 year data set of TOMS (Total Ozone Mapping Spectrometer January 1979-Dec. 1992) and corresponding 300 hPa geopotential (for the tropopause height) together with 30 hPa temperature (for stratospheric waves) at 60°N. The correlation coefficient for the linear multiple regression between total ozone (dependent variable) and the dynamical parameters (independent variables) is 0.88 for the zonal deviations in the winter of the Northern Hemisphere. Zonal means are also significantly dependent on circulation parameters, besides showing the known negative trend function of total ozone observed by TOMS. The significant linear trend for 60°N is 3 DU/year in the winter months taking into account the dependence on the dynamics between the tropopause region and the mid-stratosphere. The highest correlation coefficient for the monthly mean total ozone anomalies is reached in November with 0.94.     

Mesoscale energetics and flows induced by sea-land and mountain-valley contrasts

We study the relative importance of sea-land and mountain-valley thermal contrasts in determining the development of thermally forced mesoscale circulations (TFMCs) over a mountainous peninsula. We first analyse the energetics of the problem, and using this theory, we interprete the numerical simulations over Calabria, a mountainous peninsula in southern Italy. The CSU 3-D nonlinear numerical model is utilised to simulate the dynamics and the thermodynamics of the atmospheric fields over Calabria. Results show the importance of orography in determining the pattern of the flow and the local climate in a region as complex as Calabria. Analysis of the results shows that the energetics due to the sea-land interactions are more efficient when the peninsula is flat. The importance of the energy due to the sea-land decreases as the mountain height of the peninsula increases. The energy stored over the mountain gains in importance, untill it is released by the readjustment of the warm mountain air as it prevails over the energy released by the inland penetration of the sea breeze front. For instance, our results show that over a peninsula 100 km wide the energy over the mountain and the energy in the sea-land contrast are of the same order when the height of the mountain is about 700 m, for a 1500 m convective boundary layer (CBL) depth. Over the Calabrian peninsula, the energy released by the hot air in the CBL of the mountain prevails over the energy released by the inland penetration of the sea air. Calabria is about 1500 m high and about 50 km wide, and the CBL is of the order of 1500 m. The energy over the mountain is about four time larger than the energy contained in the sea-land contrast. Furthermore, the energetics increase with the patch width of the peninsula, and when its half width is much less than the Rossby radius, the MAPE of the sea breeze is negligible. When its half width is much larger than the Rossby radius, the breezes from the two opposing coastlines do not interact. Over Calabria peninsula, numerical simulations show that the flow is highly ageostrophic, and that the flow intensity increases from sunrise to reach its maximum in the afternoon but before sunset, which suggests that, in the late part of the day, the conversion of potential energy into kinetic energy is balanced by the dissipation.     

The effect of coastal upwelling on the sea-breeze circulation at Cabo Frio, Brazil: a numerical experiment

The effect of coastal upwelling on sea-breeze circulation in Cabo Frio (Brazil) and the feedback of sea-breeze on the upwelling signal in this region are investigated. In order to study the effect of coastal upwelling on sea-breeze a non-linear, three-dimensional, primitive equation atmospheric model is employed. The model considers only dry air and employs boundary layer formulation. The surface temperature is determined by a forcing function applied to the Earths surface. In order to investigate the seasonal variations of the circulation, numerical experiments considering three-month means are conducted: January-February-March (JFM), April-May-June (AMJ), July-August-September (JAS) and October-November-December (OND). The model results show that the sea-breeze is most intense near the coast at all the seasons. The sea-breeze is stronger in OND and JFM, when the upwelling occurs, and weaker in AMJ and JAS, when there is no upwelling. Numerical simulations also show that when the upwelling occurs the sea-breeze develops and attains maximum intensity earlier than when it does not occur. Observations show a similar behavior. In order to verify the effect of the sea-breeze surface wind on the upwelling, a two-layer finite element ocean model is also implemented. The results of simulations using this model, forced by the wind generated in the sea-breeze model, show that the sea-breeze effectively enhances the upwelling signal.     

Two-dimensional simulations of katabatic layers observed during the GIMEX experiment

The hydrostatic model SALSA is used to simulate a particular event observed during the Greenland Ice Margin EXperiment “GIMEX” (on July 12th, 1991). The time evolution of the large-scale flow was incorporated in the model through time dependent boundary conditions which were updated using the closest upwind sounding. A turbulent scheme for the stable boundary layer and an appropriate parametrization of the surface fluxes implemented in the same model, are used for this study. The simulation results are discussed and compared to the available observations. The computed turbulent fluxes are correctly estimated. The model predicts a mixing zone of about 1500 m high which is in good agreement with tundra site observations. Over the ice cap, the katabatic layer is correctly simulated by the model. Its height of 80–300 m is well estimated. The comparison between the simulation and observations taken at ice cap sites is reasonably valid. The ablation computed along the ice cap corresponds well to the values reconstructed of observations at sites 4 and 9. Finally, a sensibility study to a specified westward geostrophic wind (2 ms⁻¹) shows that the consideration of this latter improves the simulated tundra wind evolution.     

Measurements and modelling of atmospheric pollution over the Paris area: an overview of the ESQUIF Project

The Étude et Simulation de la QUalité de lair en Ile de France (ESQUIF) project is the first integrated project dedicated to the study of the processes leading to air pollution events over the Paris area. The project was carried out over two years (summer 1998 to winter 2000) to document all types of meteorological conditions favourable to air quality degradation, and in particular to photo oxydant formation. The goals of ESQUIF are (1) to improve our understanding of the relevant chemical and dynamical processes and, in turn, improve their parametrizations in numerical models, and (2) to improve and validate existing models dedicated to pollution analysis, scenarios and/or forecasting, by establishing a comprehensive and thorough database. We present the rationale of the ESQUIF project and we describe the experimental set-up. We also report on the first experiments which took place during the summer of 1998 involving surface networks, and remote sensing instruments as well as several aircraft. Focusing on three days of August 1998, the relative contributions of long-range transported and locally-produced ozone to the elevated ozone concentrations observed during this period are discussed and chemistry-transport model preliminary results on this period are compared to measurements.     

Preliminary laboratory studies of the optical scattering properties of the crystal clouds

Ice crystal clouds have an influence on the radiative budget of the earth; however, the exact size and nature of this influence has yet to be determined. A laboratory cloud chamber experiment has been set up to provide data on the optical scattering behaviour of ice crystals at a visible wavelength in order to gain information which can be used in climate models concerning the radiative characteristics of cirrus clouds. A PMS grey-scale probe is used to monitor simultaneously the cloud microphysical properties in order to correlate these closely with the observed radiative properties. Preliminary results show that ice crystals scatter considerably more at 90° than do water droplets, and that the halo effects are visible in a laboratorygenerated cloud when the ice crystal concentration is sufficiently small to prevent masking from multiple scattering.     

Hydrogen and helium isotope inner radiation belts in the Earth’s magnetosphere

Radial transport theory for inner radiation zone MeV ions has been extended by combining radial diffusive transport and losses due to Coulomb friction with local generation of D, T and ³He ions from nuclear reactions taking place on the inner edge of the inner radiation zone. Based on interactions between high energy trapped protons and upper atmospheric constituents we have included a nuclear reaction yield D, T and ³He flux source that was numerically derived from a nuclear reaction model code originally developed at the Institute of Nuclear Researches in Moscow, Russia. Magnetospheric transport computations have been made covering the L-shell range L=1.0–1.6. The resulting MeV energy D, T and ³He ion flux distributions show a strong influence of the local nuclear source mechanism on the inner zone energetic D, T and ³He ion content.     

On the application of the classic Kessler and Berry schemes in Large Eddy Simulation models with a particular emphasis on cloud autoconversion, the onset time of precipitation and droplet evaporation

Many Large Eddy Simulation (LES) models use the classic Kessler parameterisation either as it is or in a modified form to model the process of cloud water autoconversion into precipitation. The Kessler scheme, being linear, is particularly useful and is computationally straightforward to implement. However, a major limitation with this scheme lies in its inability to predict different autoconversion rates for maritime and continental clouds. In contrast, the Berry formulation overcomes this difficulty, although it is cubic. Due to their different forms, it is difficult to match the two solutions to each other. In this paper we single out the processes of cloud conversion and accretion operating in a deep model cloud and neglect the advection terms for simplicity. This facilitates exact analytical integration and we are able to derive new expressions for the time of onset of precipitation using both the Kessler and Berry formulations. We then discuss the conditions when the two schemes are equivalent. Finally, we also critically examine the process of droplet evaporation within the framework of the classic Kessler scheme. We improve the existing parameterisation with an accurate estimation of the diffusional mass transport of water vapour. We then demonstrate the overall robustness of our calculations by comparing our results with the experimental observations of Beard and Pruppacher, and find excellent agreement.     

Statistical interpolation of ozone measurements from satellite data (TOMS, SBUV and SAGE II) using the kriging method

This study demonstrates that ordinary kriging in spherical coordinates using experimental semi-variograms provides highly usable results, especially near the pole in winter and/or where there could be data missing over large areas. In addition, kriging allows display of the spatial variability of daily ozone measurements at different pressure levels. Three satellite data sets were used: Total Ozone Mapping Spectrometer (TOMS) data, Solar Backscattered UltraViolet (SBUV), and the Stratospheric Aerosol and Gas Experiment (SAGE II) ozone profiles. Since SBUV is a nadir-viewing instrument, measurements are only taken along the sun-synchronous polar orbits of the satellite. SAGE II is a limb-viewing solar occultation instrument, and measurements have high vertical resolution but poor daily coverage. TOMS has wider coverage with equidistant distribution of data (resolution 1° × 1.25°) but provides no vertical information. Comparisons of the resulting SBUV-interpolated (column-integrated) ozone field with TOMS data are strongly in agreement, with a global correlation of close to 98%. Comparisons of SBUV-interpolated ozone profiles with daily SAGE II profiles are relatively good, and comparable to those found in the literature. The interpolated ozone layers at different pressure levels are shown.     

Statistical distributions of water content and sizes for clouds above Europe

The statistics of quantities involved in the synthesis of cloud scenes have been investigated from an original data base. Frequency distributions of ice and water content (IWC), horizontal and vertical sizes (L and H), and top temperatures (T) of clouds above Europe have been derived for nine types of clouds (As, Cb, Ci, Cg, LwCg, OrCg, Cs, Ns, Sc). It appears that the cumulated frequency plots can be well fitted with log-normal or Weibull profiles, and that for IWC and T cloud types can be split into two or three classes according to slopes in logarithmic coordinates. Cross-correlation coefficients between IWC, L, H and T have been also derived. Implications for the physics of the cloud build-up processes are briefly outlined. Critical analysis and comparison of other published results are proposed.