COMPUTATIONAL GEOMETRY AND GEOGRAPHY *

New developments in computational geometry are applicable to geographic data processing. Computationally efficient procedures incorporating logarithmic searches over appropriate data structures, and other shortcuts, can yield important savings in computation time and data storage space. There are “fast”methods for two-dimensional problems such as point inclusion, nearest-neighbor and two-nearest-neighbors determination, and grid overlay. Worst-case analysis of algorithms provides useful performance bounds for the developers of geographic production software.     

Scattering in the Atmosphere of Venus. I. Line profiles and curves of growth for isotropic scattering

We have computed line profiles and curves of growth for both reflected and transmitted radiation for typical lines in CO₂ bands (in the photographic infrared) which occur in the spectrum of Venus. In our model the pressure variation with altitude was considered and the base of the cloud deck was set at the 2 bar level. The temperature was held constant at 250K and a Voigt profile was used for the lineshape. We also assumed that the scale height of the cloud particles was equal to the scale height of the gas. The calculations were made for four values of the scattering optical thickness (τ_c = 0.1, 1.0, 10, and 100) using a continuum single scattering albedo $\text{ω}{}_{\mathrm{<mtext>c</mtext>}}\text{= 0.9975}$ (which gives a Bond albedo of 0.896 for τ_c = 100, the value observed for Venus at these wavelengths). Curves of growth are also presented for reflected radiation which has been averaged over the visible disk for three values of the Venus phase angle (0, 86, and 166°).     

Physics of an active region loop system

The structure of the active region loops is investigated by the study of a loop complex which undergoes a dramatic evacuation of most of the mass it contains. The need for continual energy deposition in loops is emphasized by the apparent cessation of energy input to the loops studied and their subsequent behavior. Estimates are made of the energy necessary to form and to maintain the loops, and of the relative importance of radiation and thermal conduction as energy loss mechanisms. Models based on the observed EUV emission are used to place limits on the size of loops seen in various lines and on the density and temperature structure. We find that the cool cores of active region loops are likely to be no more than a few hundred kilometers in radius and that several such cool threads may be imbedded in a common hot outer sheath. The primary energy loss on a large scale is radiation with thermal conduction contributing to local disturbances. There is a tendency for the development of apparently unstable condensations or knots along the length of a loop. Higher resolution observations will be necessary to confirm some of our predictions.     

Hill stability and distance curves for the general three-body problem

The notion of Hill stability is extended from the circular restricted 3-body problem to the general three-body problem; it is even extended to systems of positive energy and the Hill's curves with their corresponding forbidden zones are generalized.Hill stable systems of negative energy present a hierarchy: they have a close binary that can be neither approached nor disrupted by the third body. This phenomenon becomes particularly clear with the distance curves presentation.The three limiting cases, restricted, planetary and lunar are analysed as well as some real stellar cases.     

Zero velocity surfaces for the general planar three-body problem

The angular momentum and the energy integral of the planar three-body problem are used to establish regions of the physical space where motion is allowed to take place. Although forbidden regions exist for both negative and positive values of the energy of the system, the known integrals of the motion always allow for at least one of the three bodies to escape.     

Multi‐variable relations for soil effects on seismic ground motion

Soil effects on peak ground acceleration, velocity and elastic response spectra (5% damping) are expressed by simple approximate relations in terms of five key parameters: (a) the fundamental vibration period of the non‐linear soil, T_S, (b) the period of a bedrock site of equal thickness, T_b, (c) the predominant excitation period, Te, (d) the peak seismic acceleration at outcropping bedrock, a, and (e) the number of significant excitation cycles, n. Furthermore, another relation is proposed for the estimation of T_S in terms of the soil thickness H, the average shear wave velocity of the soil V?_S,o and a. The aforementioned parameters were first identified through a simplified analytical simulation of the site excitation. The multivariable approximate relations were then formulated via a statistical analysis of relevant data from more than 700 one‐dimensional equivalent‐linear seismic ground response analyses, for actual seismic excitations and natural soil conditions. Use of these relations to back‐calculate the numerical results in the database gives an estimate of their error margin, which is found to be relatively small and unbiased. The proposed relations are also independently verified through a detailed comparison with strong motion recordings from seven well‐documented case studies: (a) two sites in the San Fernando valley during the Northridge earthquake, and (b) five different seismic events recorded at the SMART‐1 accelerometer array in Taiwan. It is deduced that the accuracy of the relations is comparable to that of the equivalent‐linear method. Hence, they can be readily used as a quick alternative for routine applications, as well as for spreadsheet computations (e.g. GIS‐aided seismic microzonation studies) where numerical methods are cumbersome to implement. Copyright © 2003 John Wiley & Sons, Ltd.     

Hydraulic Effects of Shales in Fluvial-Deltaic Deposits: Ground-Penetrating Radar, Outcrop Observations, Geostatistics, and Three-Dimensional Flow Modeling for the Ferron Sandstone, Utah

Ground-penetrating radar (GPR) surveys, outcrop measurements, and cores provide a high-resolution 3D geologic model to investigate the hydraulic effects of shales in marine-influenced lower delta-plain distributary channel deposits within the Cretaceous-age Ferron Sandstone at Corbula Gulch in central Utah, USA. Shale statistics are computed from outcrop observations. Although slight anisotropy was observed in mean length and variogram ranges parallel and perpendicular to pale of low , the anisotropy is not statistically significant and the estimated mean length is 5.4 m. Truncated Gaussian simulation was used to create maps of shales that are placed on variably dipping stratigraphic surfaces interpreted from high-resolution 3D GPR surveys, outcrop interpretations, and boreholes. Sandstone permeability is estimated from radar responses calibrated to permeability measurements from core samples. Experimentally designed flow simulations examine the effects of variogram range, shale coverage fraction, and trends in shale coverage on predicted upscaled permeability, breakthrough time, and sweep efficiency. Approximately 1500 flow simulations examine three different geologic models, flow in the 3 coordinate directions, 16 geostatistical parameter combinations, and 10 realizations for each model. ANOVA and response models computed from the flow simulations demonstrate that shales decrease sweep, recovery, and permeability, especially in the vertical direction. The effect on horizontal flow is smaller. Flow predictions for ideal tracer displacements at Corbula Gulch are sensitive to shale-coverage fraction, but are relatively insensitive to twofold variations in variogram range or to vertical trends in shale coverage. Although the hydraulic effects of shale are statistically significant, the changes in flow responses rarely exceed 20%. As a result, it may be reasonable to use simple models when incorporating analogous shales into models of reservoirs or aquifers.     

On the Dynamics of the Jovian Ionosphere and Thermosphere: III. The Modelling of Auroral Conductivity

Recent work has been concerned with calculating the three-dimensional ion concentrations and Pedersen and Hall conductivities within the auroral region of Jupiter for varying conditions of incident electron precipitation. Using the jovian ionospheric model, we present results that show the auroral ionospheric response to changing the incoming flux of precipitating electrons (for constant initial energy) and also the response to changing the initial energy (for both constant flux and constant energy flux). The results show that, for expected energy fluxes of precipitating particles, the average auroral integrated Pedersen conductivity attains values in excess of 1 mho. In addition, it is shown that electrons with an initial energy of around 60 keV are particularly effective at generating auroral conductivity: Particles of this energy penetrate most effectively to the layer of the jovian ionosphere at which the auroral conductivity is at a maximum.     

Contents to Volume 8

Volume Contents

Contents to Volume 8     

Yielding oscillator subjected to simple pulse waveforms: numerical analysis & closed‐form solutions

Numerical and analytical solutions are presented for the elastic and inelastic response of single‐degree‐of‐freedom yielding oscillators to idealized ground acceleration pulses. These motions are typical of near‐fault earthquake recordings generated by forward rupture directivity and may inflict damage in the absence of substantial structural strength and ductility capacity. Four basic pulse waveforms are examined: (1) triangular; (2) sinusoidal; (3) exponential; and (4) rectangular. In the first part of the article, a numerical study is presented of the effect of oscillator period, strength, damping, post‐yielding stiffness and number of excitation cycles, on inelastic response. Results are presented in the form of dimensionless graphs and regression formulas that elucidate the salient features of the problem. It is shown that conventional R–µ relations may significantly underestimate ductility demand imposed by near‐fault motions. The second part of the article concentrates on elastic‐perfectly plastic oscillators. Closed‐form solutions are derived for post‐yielding response and associated ductility demand. It is shown that all three ground motion histories (i.e. acceleration, velocity, and displacement) control oscillator response—contrary to the widespread view that ground velocity alone is of leading importance. The derived solutions provide insight on the physics of inelastic response, which is often obscured by the complexity of numerical algorithms and actual earthquake motions. The model is evaluated against numerical results from near‐field recordings. A case study is presented. Copyright © 2006 John Wiley & Sons, Ltd.