Global distribution of electron temperature and density fluctuations measured by AEROS-B

Abstract

Between July 1974 and September 1975 aboard AEROS-B the electron temperature and the electron density have been measured by the Retarding Potential Analyses and the Impedance Probe respectively in the altitude range 200–870 km. Temperature and density fluctuations along satellite paths have been analysed by digital filtering. The large amount of data allows the mean fluctuations at different heights and latitude intervals to be evaluated. The resulting two-dimensional arrays are presented for some time periods. The global distributions show different patterns in the polar and equatorial regions, and seasonal variations. Mean temperature fluctuations and fluctuations of electron density in general occur simultaneously though with different amplitudes. At high latitudes a seasonal dependence in the strength of the fluctuations of both quantities was observed.     

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.     

Greenhouse effects of aircraft emissions as calculated by a radiative transfer model

With a radiative transfer model, assessments are made of the radiative forcing in northern mid-latitudes due to aircraft emissions up to 1990. Considered are the direct climate effects from the major combustion products carbon dioxide, nitrogen dioxide, water vapor and sulphur dioxide, as well as the indirect effect of ozone production from NO_x emissions. Our study indicates a local radiative forcing at the tropopause which should be negative in summer (-0.5 to 0.0 W/m²) and either negative or positive in winter (-0.3 to 0.2 W/m²). To these values the indirect effect of contrails has to be added, which for the North Atlantic Flight Corridor covers the range -0.2 to 0.3 W/m² in summer and 0.0 to 0.3 W/m² in winter. Apart from optically dense non-aged contrails during summer, negative forcings are due to solar screening by sulphate aerosols. The major positive contributions come from contrails, stratospheric water vapor in winter and ozone in summer. The direct effect of NO₂ is negligible and the contribution of CO₂ is relatively small.     

Multiple flow direction algorithms for runoff modelling in grid based elevation models: An empirical evaluation

When pathways for groundwater flow are defined in grid based elevation models, the grid structure creates some obstacles as it only allows for eight flow directions. Earlier algorithms have used arbitrarily designed flow distributions, which have resulted in either too diverging or too converging flow patterns. This study defines a computationally simple distribution function, where the converging nature (x) can be altered. A material of reference flow distributions was produced through simulation in four elevation models. The distribution function was tested for different values on x, against the reference. The results show a stable correlation pattern between the tested function and the reference. Optimum values on x were found, which suggest that the flow distribution from a grid cell should be weighted with the slope gradient raised to the power of 4–6.     

Ground-based mobile scanning LIDAR for remote sensing of contrails

Air traffic is a source of trace gases in the upper troposphere and lower stratosphere. Contrails readily form from water vapor exhausts under favorable meteorological conditions. Since contrails are ice crystal clouds like natural cirrus clouds, they bear a greenhouse potential which has to be investigated. The IFU has built a scanning lidar system employing a pulsed Nd:YAG laser as the emitter and a 52-cm diameter telescope as the receiver. Signals are processed in several channels to investigate depolarization and wavelength dependencies of the light backscattered from ice crystals. These investigations are aimed at the formation and life cycles of contrails, their optical properties, and their climatological consequences in areas of dense air traffic. The experimental lidar setup is described and a sample measurement is shown.     

A high-resolution,three-dimensional,time dependent,nested grid model of the coupled thermosphere–ionosphere

First results are presented from a 3-D, time dependent, high resolution, nested grid model that has been developed to study mesoscale processes in the global, coupled thermosphere–ionosphere system. This new Thermosphere–Ionosphere Nested Grid (TING) model, which is an extension of the National Center for Atmospheric Researchs thermosphere–ionosphere general circulation model (NCAR–TIGCM), runs on a UNIX workstation. The TING model simultaneously calculates global (coarse resolution) and local (high resolution) distributions of neutral and plasma winds, temperature and composition. It is comprised of two coupled codes—a global TIGCM and an adjustable nested grid code which uses the same solvers as the TIGCM, but has higher spatial and temporal resolution. The size, location and level of nesting of the high resolution grid(s) are adjustable to suit the specific application. The coupling between the coarse (TIGCM) grid and the nested interior grids is via a one-way interaction scheme. In this scheme, the TIGCM output influences the nested grid model by providing initial conditions and temporally evolving boundary conditions, but the outputs from the nested grid are not permitted to influence the TIGCM. Diurnally-reproducible results of the TING model are presented for solar-maximum, winter solstice, geomagnetically-quiet conditions. The TING model successfully simulates well-known thermosphere–ionosphere features that are smeared or not modeled at the spatial resolutions used in standard TIGCMs. These include the sub-auroral electron density trough, the polar cap hole and the polar cap tongue of ionization.     

The global distribution of gravity wave energy in the lower stratosphere derived from GPS data and gravity wave modelling: Attempt and challenges

Five years of global temperatures retrieved from radio occultations measured by Champ (Challenging Minisatellite Payload) and SAC-C (Satelite de Aplicaciones Cientificas-C) are analyzed for gravity waves (GWs). In order to separate GWs from other atmospheric variations, a high-pass filter was applied on the vertical profile. Resulting temperature fluctuations correspond to vertical wavelengths between 400 m (instrumental resolution) and 10 km (limit of the high-pass filter). The temperature fluctuations can be converted into GW potential energy, but for comparison with parameterization schemes GW momentum flux is required. We therefore used representative values for the vertical and horizontal wavelength to infer GW momentum flux from the GPS measurements. The vertical wavelength value is determined by high-pass filtering, the horizontal wavelength is adopted from a latitude-dependent climatology. The obtained momentum flux distributions agree well, both in global distribution and in absolute values, with simulations using the Warner and McIntyre parameterization (WM) scheme. However, discrepancies are found in the annual cycle. Online simulations, implementing the WM scheme in the mechanistic COMMA-LIM (Cologne Model of the Middle Atmosphere—Leipzig Institute for Meteorology) general circulation model (GCM), do not converge, demonstrating that a good representation of GWs in a GCM requires both a realistic launch distribution and an adequate representation of GW breaking and momentum transfer.     

Compressive fluctuations in the solar wind and their polytropic index

Magnetohydrodynamic compressive fluctuations of the interplanetary plasma in the region from 0.3 to 1 AU have been characterized in terms of their polytropic index. Following Chandrasekhar’s approach to polytropic fluids, this index has been determined through a fit of the observed variations of density and temperature. At least three different classes of fluctuations have been identified: (1) variations at constant thermal pressure, in low-speed solar wind and without a significant dependence on distance, (2) adiabatic variations, mainly close to 1 AU and without a relevant dependence on wind speed, and (3) variations at nearly constant density, in fast wind close to 0.3 AU. Variations at constant thermal pressure are probably a subset of the ensemble of total-pressure balanced structures, corresponding to cases in which the magnetic field magnitude does not vary appreciably throughout the structure. In this case the pressure equilibrium has to be assured by its thermal component only. The variations may be related to small flow-tubes with approximately the same magnetic-field intensity, convected by the wind in conditions of pressure equilibrium. This feature is mainly observed in low-velocity solar wind, in agreement with the magnetic topology (small open flow-tubes emerging through an ensemble of closed structures) expected for the source region of slow wind. Variations of adiabatic type may be related to magnetosonic waves excited by pressure imbalances between contiguous flow-tubes. Such imbalances are probably built up by interactions between wind flows with different speeds in the spiral geometry induced by the solar rotation. This may account for the fact that they are mainly found at a large distance from the sun. Temperature variations at almost constant density are mostly found in fast flows close to the sun. These are the solar wind regions with the best examples of incompressible behaviour. They are characterized by very stable values for particle density and magnetic intensity, and by fluctuations of Alfvénic type. It is likely that temperature fluctuations in these regions are a remnant of thermal features in the low solar atmosphere. In conclusion, the polytropic index appears to be a useful tool to understand the nature of the compressive turbulence in the interplanetary plasma, as far as the frozen-in magnetic field does not play a crucial role.     

Precise temperature measurement of groundwater for earthquake-prediction study

Precise measurements of groundwater temperature, with a resolution of 0.0001°C, have been made for the purpose of earthquake prediction at four water wells in the Tokai district in Japan. The characteristics of temporal variations in groundwater temperature differ from well to well. From records of a well located in an area where use of groundwater is heavy, movement of groundwater was sensed by temperature variations. By contrast, water temperatures in a 500 m well far from cultural areas were extremely steady, with fluctuations of less than 0.0005°C. This well appears to be suitable for monitoring possible temperature changes related to the occurrence of an impending large earthquake.     

中长周期数字化面波记录与中国东南地区地壳结构

本文使用适配滤波频时分析技术首次对中国数字地震台网的中长周期面波记录进行处理,获得穿过东南地区的82条勒夫波频散数据.使用随机反演理论,获得了东南沿海地区4°×4°网格的纯路径频散数据.这些频散的周期为1.95-68.27s,弥补了长周期面波所不能分辨的浅层结构.在网格反演的基础上,使用Harkrider的面波及演程序得出了中国东南地区的地壳和上地幔结构,浅部可分辨到1km,深部可达80km.在分辨率保证的前提下得出东南地区深至80km的三维剪切波速度结构.     

On the regional climatic impact of contrails: microphysical and radiative properties of contrails and natural cirrus clouds

The impact of contrail-induced cirrus clouds on regional climate is estimated for mean atmospheric conditions of southern Germany in the months of July and October. This is done by use of a regionalized one-dimensional radiative convective model (RCM). The influence of an increased ice cloud cover is studied by comparing RCM results representing climatological values with a modified case. In order to study the sensitivity of this effect on the radiative characteristics of the ice cloud, two types of additional ice clouds were modelled: cirrus and contrails, the latter cloud type containing a higher number of smaller and less of the larger cloud particles. Ice cloud parameters are calculated on the basis of a particle size distribution which covers the range from 2 to 2000 m, taking into consideration recent measurements which show a remarkable amount of particles smaller than 20 m. It turns out that a 10% increase in ice cloud cover leads to a surface temperature increase in the order of 1K, ranging from 1.1 to 1.2K in July and from 0.8 to 0.9K in October depending on the radiative characteristics of the air-traffic-induced ice clouds. Modelling the current contrail cloud cover which is near 0.5% over Europe yields a surface temperature increase in the order of 0.05 K.     

Deciphering seismic and normal-force fluctuation signatures of debris flows: An experimental assessment of effects of flow composition and dynamics

Debris flows are gravity-driven mass movements that are common natural hazards in mountain regions worldwide. Previous work has shown that measurements of ground vibrations are capable of detecting the timing, speed, and location of debris flows. A remaining question is to what extent additional flow properties, such as grain-size distribution and flow depth can be inferred reliably from seismic data. Here, we experimentally explore the relation of seismic vibrations and normal-force fluctuations with debris-flow composition and dynamics. We use a 5.4 m long and 0.3 m wide channel inclined at 20°, equipped with a geophone plate and force plate. We show that seismic vibrations and normal-force fluctuations induced by debris flows are strongly correlated, and that both are affected by debris-flow composition. We find that the effects of the large-particle distribution on seismic vibrations and normal-force fluctuations are substantially more pronounced than the effects of water fraction, clay fraction, and flow volume, especially when normalized by flow depth. We further show that for flows with similar coarse-particle distributions seismic vibrations and normal-force fluctuations can be reasonably well related to flow depth, even if total flow volume, water fraction, and the size distribution of fines varies. Our experimental results shed light on how changes in large-particle, clay, and water fractions affect the seismic and force-fluctuation signatures of debris flows, and provide important guidelines for their interpretation.     

Aviation-Produced Aerosols and Contrails

Liquid and solid particles in the plumes of jet aircraft cruising in the upper troposphere and lower stratosphere lead to the formation of ice clouds (contrails), modify the microphysical properties of existing cirrus clouds, and provide sites for heterogeneous chemical reactions. Characterization of aviation-produced particles in terms of physico-chemical properties is an important step in assessing the global impact of aircraft emissions upon atmospheric chemistry and climate parameters. Chemistry and microphysics of the gas-aerosol system in aircraft plumes and its evolution in the atmosphere is a field of intense research. This paper reviews the current knowledge (mid-1998) and outlines possible atmospheric implications.     

Propagation of seasonal temperature signals into an aquifer upon bank infiltration

Infiltrating river water carries the temperature signal of the river into the adjacent aquifer. While the diurnal temperature fluctuations are strongly dampened, the seasonal fluctuations are much less attenuated and can be followed into the aquifer over longer distances. In one-dimensional model with uniform properties, this signal is propagated with a retarded velocity, and its amplitude decreases exponentially with distance. Therefore, time shifts in seasonal temperature signals between rivers and groundwater observation points may be used to estimate infiltration rates and near-river groundwater velocities. As demonstrated in this study, however, the interpretation is nonunique under realistic conditions. We analyze a synthetic test case of a two-dimensional cross section perpendicular to a losing stream, accounting for multi-dimensional flow due to a partially penetrating channel, convective-conductive heat transport within the aquifer, and heat exchange with the underlying aquitard and the land surface. We compare different conceptual simplifications of the domain in order to elaborate on the importance of different system elements. We find that temperature propagation within the shallow aquifer can be highly influenced by conduction through the unsaturated zone and into the underlying aquitard. In contrast, regional groundwater recharge has no major effect on the simulated results. In our setup, multi-dimensionality of the flow field is important only close to the river. We conclude that over-simplistic analytical models can introduce substantial errors if vertical heat exchange at the aquifer boundaries is not accounted for. This has to be considered when using seasonal temperature fluctuations as a natural tracer for bank infiltration.     

Modelling of meteorological conditions at an urban scale for the PUMA winter campaign

For the winter 2000 campaign of the Pollution of Urban Midlands Atmosphere project, observation and numerical modelling of meteorological conditions over the West Midlands conurbation, UK, was undertaken. Modelling was performed using the regional atmospheric meteorological system (RAMS). This paper presents a comparison of modelled and observed wind and temperature for 25 and 26 January 2000. The RAMS model uses two nested grids with a mesh size of 2 km for the inner grid which is embedded in the outer grid with a mesh size of 8 km. Statistical evaluation of the model results against the observational data of wind speed, direction and temperature at 10 m was conducted. In general, the modelling results are in a reasonable agreement with observation. The statistical evaluation suggests that model performance is poorer for the inner grid than the outer grid as the model uncertainties (mainly mean bias) transfer from the outer to inner one. The low indices of agreement of temperature and wind are mainly associated with the systematic root-mean-square-difference values. For temperature, the systematic bias may also be affected by representation of cloud amount by the model. For wind, the model tends to have a poor performance for calm conditions, as under a stable anti-cyclonic situation local wind patterns associated with topography may develop, although the topography of the region is relatively flat. The results for the inner grid reveal some subtle spatial patterns at a scale smaller than 10 km near hills and valleys with differences in elevation of a few hundred metres.     

Analyzing bank filtration by deconvoluting time series of electric conductivity

Knowing the travel-time distributions from infiltrating rivers to pumping wells is important in the management of alluvial aquifers. Commonly, travel-time distributions are determined by releasing a tracer pulse into the river and measuring the breakthrough curve in the wells. As an alternative, one may measure signals of a time-varying natural tracer in the river and in adjacent wells and infer the travel-time distributions by deconvolution. Traditionally this is done by fitting a parametric function such as the solution of the one-dimensional advection-dispersion equation to the data. By choosing a certain parameterization, it is impossible to determine features of the travel-time distribution that do not follow the general shape of the parameterization, i.e., multiple peaks. We present a method to determine travel-time distributions by nonparametric deconvolution of electric-conductivity time series. Smoothness of the inferred transfer function is achieved by a geostatistical approach, in which the transfer function is assumed as a second-order intrinsic random time variable. Nonnegativity is enforced by the method of Lagrange multipliers. We present an approach to directly compute the best nonnegative estimate and to generate sets of plausible solutions. We show how the smoothness of the transfer function can be estimated from the data. The approach is applied to electric-conductivity measurements taken at River Thur, Switzerland, and five wells in the adjacent aquifer, but the method can also be applied to other time-varying natural tracers such as temperature. At our field site, electric-conductivity fluctuations appear to be an excellent natural tracer.     

Using the pseudospectral technique on curved grids for 2D acoustic forward modelling1

The Fourier pseudospectral method has been widely accepted for seismic forward modelling because of its high accuracy compared to other numerical techniques. Conventionally, the modelling is performed on Cartesian grids. This means that curved interfaces are represented in a ‘staircase fashion‘causing spurious diffractions. It is the aim of this work to eliminate these non-physical diffractions by using curved grids that generally follow the interfaces. A further advantage of using curved grids is that the local grid density can be adjusted according to the velocity of the individual layers, i.e. the overall grid density is not restricted by the lowest velocity in the subsurface. This means that considerable savings in computer storage can be obtained and thus larger computational models can be handled. One of the major problems in using the curved grid approach has been the generation of a suitable grid that fits all the interfaces. However, as a new approach, we adopt techniques originally developed for computational fluid dynamics (CFD) applications. This allows us to put the curved grid technique into a general framework, enabling the grid to follow all interfaces. In principle, a separate grid is generated for each geological layer, patching the grid lines across the interfaces to obtain a globally continuous grid (the so-called multiblock strategy). The curved grid is taken to constitute a generalised curvilinear coordinate system, where each grid line corresponds to a constant value of one of the curvilinear coordinates. That means that the forward modelling equations have to be written in curvilinear coordinates, resulting in additional terms in the equations. However, the subsurface geometry is much simpler in the curvilinear space. The advantages of the curved grid technique are demonstrated for the 2D acoustic wave equation. This includes a verification of the method against an analytic reference solution for wedge diffraction and a comparison with the pseudospectral method on Cartesian grids. The results demonstrate that high accuracies are obtained with few grid points and without extra computational costs as compared with Cartesian methods.     

Grid resolution study of ground water flow and transport

Three-dimensional grids representing a heterogeneous, ground water system are generated at 10 different resolutions in support of a site-scale flow and transport modeling effort. These grids represent hydrostratigraphy near Yucca Mountain, Nevada, consisting of 18 stratigraphic units with contrasting fluid flow and transport properties. The grid generation method allows the stratigraphy to be modeled by numerical grids of different resolution so that comparison studies can be performed to test for grid quality and determine the resolution required to resolve geologic structure and physical processes such as fluid flow and solute transport. The process of generating numerical grids with appropriate property distributions from geologic conceptual models is automated, thus making the entire process easy to implement with fewer user-induced errors. The series of grids of various resolutions are used to assess the level at which increasing resolution no longer influences the flow and solute transport results. Grid resolution is found to be a critical issue for ground water flow and solute transport. The resolution required in a particular instance is a function of the feature size of the model, the intrinsic properties of materials, the specific physics of the problem, and boundary conditions. The asymptotic nature of results related to flow and transport indicate that for a hydrologic model of the heterogeneous hydrostratigraphy under Yucca Mountain, a horizontal grid spacing of 600 m and vertical grid spacing of 40 m resolve the hydrostratigraphic model with sufficient precision to accurately model the hypothetical flow and solute transport to within 5% of the value that would be obtained with much higher resolution.