Evaluation of Weather Catastrophe Data for Use in Climate Change Investigations

A 1950–1994 data set of major weather losses developed by the property insurance industry was examined to assess its potential utility in climate change research and use in assessing the relevance of recent extreme losses in the United States. A process for adjusting these historical storm losses to ever-changing factors including dollar values, amount of insurance coverage per area, and the sensitivity of society to damaging storms was developed by the industry. Analysis of the temporal frequency and losses of these adjusted weather catastrophes revealed differences according to the amount of loss. Temporal changes since 1975 in the catastrophes causing $35 to $100 million in loss were strongly related to changes in U.S. population, whereas catastrophes that created insured losses greater than $100 million appear related to both shifting weather conditions and to regional population changes. This evaluation revealed that the industry's catastrophe adjustment technique did not adequately allow for changes in various demographic and social factors affecting damage; however, results suggest use of population values for normalizing the adjusted catastrophe database to allow meaningful studies of their temporal variability.     

The Cartesian method for solving partial differential equations in spherical geometry

Cartesian coordinates are used to solve the nonlinear shallow-water equations on the sphere. The two-dimensional equations, in spherical coordinates, are first embedded in a three-dimensional system in a manner that preserves solutions of the two-dimensional system. That is, solutions of the three-dimensional system, with appropriate initial conditions, also solve the two-dimensional system on the surface of the sphere. The higher dimensional system is then transformed to Cartesian coordinates. Computations are limited to the surface of the sphere by projecting the equations, gradients, and solution onto the surface. The projected gradients are approximated by a weighted sum of function values on a neighboring stencil. The weights are determined by collocation using the spherical harmonics in trivariate polynomial form. That is, the weights are computed from the requirement that the projected gradients are near exact for a small set of spherical harmonics. The method is applicable to any distribution of points and two test cases are implemented on an icosahedral geodesic grid. The method is both vectorizable and parallelizable.     

Mean Profiles of Moisture Fluxes in Snow-Filled Boundary Layers

Profiles of moisture fluxes have been examined for convective boundary layers containing clouds and snow, using data derived from aircraft measurements taken on four dates during the 1983/1984 University of Chicago lake-effect snow project. Flux profiles were derived from vertical stacks of aircraft cross-wind flight legs taken at various heights over Lake Michigan near the downwind shore. It was found that, if ice processes are taken into account, profiles of potential temperature and water content were very similar to those presented in past studies of convective boundary layers strongly heated from below. Profiles of total water content and equivalent potential temperature adjusted for ice were nearly invariant with height, except very near the top of the boundary layer, suggesting that internal boundary-layer mixing processes were rapid relative to the rates at which heat and vapour were transported into the boundary layer through entrainment and surface fluxes. Ice was found to play a significant, measurable role in boundary-layer moisture fluxes. It was estimated that 40 to 57% of the upward vapour flux was returned to the surface in the form of snow, converting about 45 to 64% of the surface latent heat flux into sensible heat in the snow-producing process. Assuming advective fluxes are relatively small (thought to be appropriate after the first few tens of km over the lake as suggested by past studies), the boundary layer was found to warm at a rate faster than could be explained by surface heat fluxes and latent heat releases alone, the remainder of the heating presumably coming from radiational processes and entrainment. Discussions of moisture phase change processes throughout the boundary layer and estimates of errors of these flux measurements are presented.     

Filtering vibroseis data in the precorrelation domain

Vibroseis data recorded at short source–receiver offsets can be swamped by direct waves from the source. The signal-to-noise ratio, where primary reflections are the signal and correlation side lobes are the noise, decreases with time and late reflection events are overwhelmed. This leads to low seismic resolution on the vibroseis correlogram. A new precorrelation filtering approach is proposed to suppress correlation noise. It is the ‘squeeze-filter-unsqueeze’ (SFU) process, a combination of ‘squeeze’ and ‘unsqueeze’ (S and U) transformations, together with the application of either an optimum least-squares filter or a linear recursive notch filter. SFU processing provides excellent direct wave removal if the onset time of the direct wave is known precisely, but when the correlation recognition method used to search for the first arrival fails, the SFU filtering will also fail. If the tapers of the source sweeps are badly distorted, a harmonic distortion will be introduced into the SFU-filtered trace. SFU appears to be more suitable for low-noise vibroseis data, and more effective when we know the sweep tapers exactly. SFU requires uncorrelated data, and is thus cpu intensive, but since it is automatic, it is not labour intensive. With non-linear sweeps, there are two approaches to the S,U transformations in SFU. The first requires the non-linear analytical sweep formula, and the second is to search and pick the zero nodes on the recorded pilot trace and then carry out the S,U transformations directly without requiring the algorithm or formula by which the sweep was generated. The latter method is also valid for vibroseis data with a linear sweep. SFU may be applied to the removal of any undesired signal, as long as the exact onset time of the unwanted signal in the precorrelation domain is known or determinable.     

Quantifying groundwater discharge to a small perennial stream in southern Ontario, Canada

A study of the interaction between groundwater and surface water was undertaken within a small agricultural watershed in southern Ontario, Canada. Groundwater contributions to streamflow were measured along a section of stream during baseflow conditions and during rainfall events. Four techniques were used to estimate the contribution of groundwater to the stream along a 450 m reach (three during baseflow and one during stormflow conditions). Under baseflow conditions, streamflow measurements using the velocity–area technique indicated that the net groundwater flux to the stream during the summer months was 10 ml s⁻¹ m⁻¹. Hydrometric measurements (i.e. hydraulic gradient and hydraulic conductivity) taken using mini-piezometers installed in the sediments beneath the stream resulted in net groundwater flux estimates that were four to five times lower. Seepage meters failed to provide any measurements of water flux into or out of the stream. Therefore, based on these results, the velocity–area technique gives the best estimate of groundwater discharge. Hydrograph separations were conducted using isotopic ratios and electrical conductivity on two large rainfall events with different antecedent moisture conditions in the catchment. Both events showed that pre-event water (generally considered groundwater) dominated streamflow and tile drain flow with 64%–80% of the total discharge contributed by pre-event water. High water table conditions within the catchment resulted in greater stream discharge and a greater contribution of event water in the streamflow than that observed under low water table conditions for similar intensity storm events. The results also showed that differences in riparian zone width, vegetation and surface saturation conditions between the upper and lower catchment can influence the relative magnitude of streamflow response from the two catchment areas.     

Cleavage-fold relationships and their implications for transected folds: an example from southwest Virginia,U.S.A.

A non-coaxial deformation involving pre-folding initiation of cleavage perpendicular to bedding is proposed to explain non-axial planar cleavage associated with mesoscopic folds in part of the Appalachian foreland thrust-belt of southwest Virginia. Folds are gently plunging, asymmetric, upright to slightly inclined, sinusoidal forms with non-axial fanning cleavage. They show extreme local variations in type and degree of transection and the consistency of transection direction. These relations are further complicated by hinge migration.Cleavage-fan angles, bedding-cleavage angles and δ transection values appear influenced by fold tightness, and in part by fold flattening strain. Fold flattening increments are considered simultaneous with folding. Axial surface traces, and not cleavage traces, coincide with the principal extension direction in fold profiles. Geometric modelling of cleavage fanning and bedding-cleavage angle variations for various theoretical folding modes suggest that folding in limestone and sandstone layers was by tangential longitudinal strain. Significant shape modification and change in bedding-cleavage relations occurred after limb dips of 40 and 50° were attained in limestone and sandstone respectively. Mud-rock class 1C folds with convergent cleavage fans show features transitional between buckling and flexural flow. Initiation of ‘cleavage’ fabrics during layer-parallel shortening prior to significant folding may be important for cleavage evolution in some deformed rocks.     

The socioeconomic impacts of nuclear generating stations: An analysis of the Rancho Seco and peach bottom facilities

This paper presents the results of a comparative analysis of the socioeconomic impacts of the construction and operation of two nuclear generating stations. Although the construction of nuclear power plants is typically a multiyear process utilizing large numbers of workers and requiring large expenditures for equipment and materials, the socioeconomic changes in the areas in which the plants are located were small and temporary. The extent and size of the changes were found to be related to the size of the work force residing in the local area, the magnitude of local utility and contractor purchases, the amount of tax payments resulting from the plant, and the level of involvement of area groups over plant-related issues.     

Exposed cross-sections through the continental crust: implications for crustal structure,petrology, and evolution

The continental crust is exposed in cross-section at numerous sites on the earth's surface. These exposures, which appear to have formed by obduction along great faults during continental collision, may be recognized by exposures of deep crustal rocks exhibiting asymmetric patterns of metamorphic grade and age across the faults and by distinctive Bouguer anomaly patterns reflecting dipping basement structure and an anomalously deep mantle. From an examination of five complexes which meet these criteria, it is concluded that the most prominent layering in the crust is not compositional but metamorphic. The lower crust consists of granulite facies rocks of mafic to intermediate composition while the intermediate and shallow levels consist predominantly of amphibolite facies gneisses and greenschist facies supracrustal rocks, respectively. Post-metamorphic granitic intrusions are common at intermediate to shallow levels. Position of discontinuities in refraction velocity, where present, commonly correspond to changes in composition or metamorphic grade with depth. The continental crust is characterized by lateral and vertical heterogeneities of varying scale which are the apparent cause of the complex seismic reflections recorded by COCORP. Field observations, coupled with geochemical data, indicate a complex evolution of the lower crust which can include anatexis, multiple deformation, polymetamorphism and reworking of older crustal material. The complexity of the crust is thus the result of continuous evolution by recycling and metamorphism through time in a variety of tectonic environments.