Electrical Imaging of Impact Structures

Electrical imaging provides important subsurface information for the construction of hypervelocity impact models. We here provide an overview and evaluation of the current electrical imaging methods used in impact cratering studies. Although apparent resistivity models are commonly used in the geoelectrical imaging of impact structures, the reliability of these models has not hitherto been determined. In order to assess these imaging approaches in impact cratering, we investigate for the first time the discrepancies between the apparent resistivity and true resistivity models of an impact structure. To this end, we present (1) a new true resistivity model of the Araguainha impact structure in central Brazil by applying L2-norm inversion to previously published data, (2) apparent resistivity model of the impact structure, and (3) models obtained from different stages of the iterative tomographic inversions. Our results show that changes in vertical resistivity gradient are significantly better defined in the true resistivity models than in the apparent resistivity model. On the basis of these results, we outline a new approach that true resistivity models can be effectively assessed by applying both L1- and L2-norm inversion schemes together with the monitoring of intermediate models from iterative inversion. The results of our study highlight the importance of tomographic inversion of resistivity data in impact cratering studies, and they provide a data modeling framework and foundation for cost-effective subsurface imaging of impact structures in the future.     

GEOCHEMICAL RECOGNITION OF SPILLED SEDIMENTS USED IN NUMERICAL MODEL VALIDATION

A fixed link (tunnel and bridge, in total 16 km) was constructed between Sweden and Denmark during 1995-2000. As part of the work, approximately 16 million tonnes of seabed materials (limestone and clay till) were dredged, and about 0.6 million tonnes of these were spilled in the water. Modelling of the spreading and sedimentation of the spilled sediments took place as part of the environmental monitoring of the construction activities. In order to verify the results of the numerical modelling of sediment spreading and sedimentation, a new method with the purpose of distinguishing between the spilled sediments and the naturally occurring sediments was developed. Because the spilled sediments tend to accumulate at the seabed in areas with natural sediments of the same size, it is difficult to separate these based purely on the physical properties. The new method is based on the geo-chemical differences between the natural sediment in the area and the spill. The basic propertiesused are the higher content of calcium carbonate material in the spill as compared to the natural sediments and the higher Ca/Sr ratio in the spill compared to shell fragments dominating the natural calcium carbonate deposition in the area. The reason for these differences is that carbonate derived from recent shell debris can be discriminated from Danien limestone, which is the material in which the majority of the dredging took place, on the basis of the Ca/Sr ratio being 488 in Danien Limestone and 237 in shell debris. The geochemical recognition of the origin of the sediments proved useful in separating the spilled from the naturally occurring sediments. Without this separation, validation of the modelling of accumulation of spilled sediments would not have been possible. The method has general validity and can be used in many situations where the origin of a given sediment is sought.     

Mean and eddy motion in the Skagerrak/northern North Sea: insight from a numerical model

Insight regarding the mean and eddy motion in the Skagerrak/northern North Sea area is gained through an analysis of model-simulated currents, hydrography, kinetic energy and relative vorticity for the 2 years 2000 and 2001. In this a -coordinate ocean model is used. Since the tidal currents are generally strong in the area, care is exercised to distinguish the mesoscale (eddy) motion from higher-frequency motion such as tides, before computing the mean and eddy kinetic energy. The model-simulated response is first compared with available knowledge of the circulation in the area, and when available, also with sea-surface temperature obtained from satellite imagery. It is concluded that the model appears to faithfully reproduce most of what is known, in particularly the upper mixed layer circulation. An analysis of the mean and eddy kinetic energy reveals that many of the mesoscale structures found in the area are recurrent. This is particularly true for the structures off the southern tip of Norway. Also in general, areas of strong mean and eddy kinetic energy are co-located. The exception is the area off the southern tip of Norway, where the eddy kinetic energy is much larger than its mean counterpart. An analysis of the relative vorticity reveals that the variability found is due to the occurrence of recurrent anticyclonic eddies. It is hypothesized that these eddies are generated due to an offshore veering of the Norwegian coastal current (NCC) as it reaches the eastern end of the Norwegian Trench plateau. Here it becomes a free jet, which is then vulnerable to either barotropic instability caused by the horizontal shear in the jet-like structure of the NCC at this point, or a baroclinic (frontal) instability. The latter may come into play when the NCC veers offshore and its relatively fresh water meets the inflowing saline water of Atlantic origin, a frontogenesis that may become strong enough for cyclogenesis to take place. Due to the depth-independent nature of the model-generated eddies, the barotropic instability is the most likely candidate. It remains to resolve the reason for the offshore veering of the NCC. The most likely candidate mechanisms are vortex squeezing or simply that the coastline curvature is large enough for the NCC to separate from the coast in a hydraulic sense.Responsible Editor: Phil Dyke     

Building and Testing an <Emphasis Type="Italic">a priori</Emphasis> Geophysical Model for Western Eurasia and North Africa

We construct and evaluate a new three-dimensional model of crust and upper mantle structure in Western Eurasia and North Africa (WENA) extending to 700 km depth and having 1° parameterization. The model is compiled in an a priori fashion entirely from existing geophysical literature, specifically, combining two regionalized crustal models with a high-resolution global sediment model and a global upper mantle model. The resulting WENA1.0 model consists of 24 layers: water, three sediment layers, upper, middle, and lower crust, uppermost mantle, and 16 additional upper mantle layers. Each of the layers is specified by its depth, compressional and shear velocity, density, and attenuation (quality factors, Q _P and Q _S ). The model is tested by comparing the model predictions with geophysical observations including: crustal thickness, surface wave group and phase velocities, upper mantle _n velocities, receiver functions, P-wave travel times, waveform characteristics, regional 1-D velocities, and Bouguer gravity. We find generally good agreement between WENA1.0 model predictions and empirical observations for a wide variety of independent data sets. We believe this model is representative of our current knowledge of crust and upper mantle structure in the WENA region and can successfully be used to model the propagation characteristics of regional seismic waveform data. The WENA1.0 model will continue to evolve as new data are incorporated into future validations and any new deficiencies in the model are identified. Eventually this a priori model will serve as the initial starting model for a multiple data set tomographic inversion for structure of the Eurasian continent.     

Assessing Oregon's twenty-six coastal beach areas for recreational water quality standards

Water samples from 26 Oregon beaches were analyzed for Escherichia coli (E. coli) and enterococci concentrations by the Oregon Department of Human Services (ODHS) Public Health Laboratory. Nine Oregon beaches exceeded US Environmental Protection Agency's (USEPA) single sample maximum density of 104 enterococci colony forming units (cfu) per 100 mL with levels ranging from 121 to 4325 most probable number (MPN)/100 milliliters (mL). Otter Rock at South Cove had the highest enterococci concentration at 4352 MPN/100 mL. The results from two Oregon beaches exceeded Oregon Department of Environmental Quality's (ODEQ) estuarine E. coli standard of 406 organisms/100 mL. Otter Rock at South Cove had the highest E. coli concentration at 1850 MPN/100 mL while Road's End had an E. coli density of 771 MPN/100 mL. Results of this study suggest that adopting USEPA's marine enterococci standard in lieu of the ODEQ's estuarine standard will lead to increased Oregon beach water standard failures.     

Liquefaction mitigation in silty soils using composite stone columns and dynamic compaction

The objective of this study is to develop an analytical methodology to evaluate the effectiveness ofvibro stone column (S. C.) and dynamic compaction (D.C.) techniques supplemented with wick drains to densify and mitigate liquethctionin saturated sands and non-plastic silty soils. It includes the following: (i) develop numerical models to simulate and analyze soil densification during S.C. installation and D.C. process, and (ii) identify parameters controlling post-improvement soil density in both cases, and (iii) develop design guidelines for densification of silty soils using the above techniques. An analytical procedure was developed and used to simulate soil response during S.C. and D.C. installations, and the results were compared with available case history data. Important construction design parameters and soil properties that affect the effectiveness of these techniques, and construction design choices suitable for sands and non-plastic silty soils were identified. The methodology is expected to advance the use of S.C. and DC. in silty soils reducing the reliance on expensive field trials as a design tool. The ultimate outcome of this research will be design charts and design guidelines for using composite stone columns and composite dynamic compaction techniques in liquefaction mitigation of saturated silty soils.     

Foraminifers as indicators of marine pollution: a culture experiment with Rosalina leei

In order to develop a viable foraminiferal proxy for heavy metal pollutants, juvenile specimens of Rosalina leei were subjected to different mercury concentrations (0-180 ng/l). Initially considerable growth was observed in specimens kept in saline water having a mercury concentration up to 100 ng/l. But with the gradual increase in concentration of mercury the growth rate started decreasing. Total growth achieved was significantly lower in case of specimens kept at relatively higher mercury concentrations then those maintained in normal saline water. The most significant result of this experiment was the addition of abnormal chambers in the specimens kept at higher mercury concentration. Later the specimens kept at highest concentration (180 ng/l) were subjected to progressively increasing concentration of mercury to see the further effects and it was found that the specimens were still living at as high a mercury concentration as 260 ng/l although there was no growth.     

GOCE Gravity Field Processing Strategy

A processing strategy and the corresponding software architecture for the processing of GOCE (Gravity field and steady-state Ocean Circulation Explorer) observables is presented and described, with the major objective to compute a high-accuracy, high-resolution spherical harmonic model of the Earth's gravity field. The combination of two numerical solution strategies, i.e. the rigorous solution of the corresponding large normal equation systems applying parallel processing (on a PC cluster) as the core solver, and the fast semianalytic approach as a quick-look gravity field analysis (QL-GFA) tool, is proposed. Such a method fusion benefits from the advantages of the individual components: the rigorous inversion of the system providing also the full variance-covariance information, and the quickness enabling the consecutive production of intermediate gravity field solutions, for the purpose to analyse partial and incomplete data sets and to derive a diagnosis of the performance of the GOCE measurement system. The functionality and operability of the individual components are demonstrated in the framework of a closed loop simulation, which is based on a realistic mission scenario both in terms of the orbit configuration and the coloured measuring noise. Special concern is given to the accuracy of the recovered coefficients, the numerical behaviour, the required computing time, and the particular role of the individual modules within the processing chain. In the case of the core solver, it is demonstrated that the assembling and rigorous solution of large normal equation systems can be handled by using Beowulf clusters within a reasonable computing time. The application of the quick-look tool to partial data sets with short-term data gaps is demonstrated on the basis of several case studies. Additionally, the spectral analysis of the residuals of the adjustment is presented as a valuable tool for the verification of the noise characteristics of the GOCE gradiometer.     

Multicomponent ionic diffusion in porewaters: Coulombic effects revisited

The diffusion of an ion in porewaters cannot occur independently of the other ions in solution as a result of Coulombic coupling, as well as from other effects not considered here. Unfortunately, a longstanding disagreement exists about the correct form and meaning of the equations that describe Coulombic coupling in porewaters, i.e., Ben-Yaakov [Am. J. Sci. 281 (1981) 974] vs. Lasaga [Am. J. Sci. 281 (1981) 981]. This paper re-examines this controversy by reformulating the problem starting from fundamental concepts of mass and charge conservation. We show that these antagonistic formulations are both valid and, in fact, equivalent, when the different interpretations of charge balance are resolved. Most of the disagreements between Ben-Yaakov and Lasaga are then shown to result from differing methods of solution, not fundamental disparities in their models. We note, however, that the explanation for the concept of “stationary” gradients of nonreacting ions as given Ben-Yaakov is inaccurate, and such gradients do lead to diffusive fluxes that are counterbalanced by electrochemical migrational fluxes to produce no net flux (excluding advective flux). We further find that the bicarbonate diffusive flux will not balance the diffusional charge flux of sulfate during its reduction if advection is present. This latter imbalance generates compensating fluxes in the other nonreacting ions. We have applied our theory to a simplified case of sulfate reduction in a marine sediment. The results show that nonreacting ions do diffuse and that with normally expected values of porewater advection, the ratio of the bicarbonate to the sulfate flux can be far different than the ideal value of −2.