Naturaleza y estructura geológica del borde costero del noroeste de Costa Rica / Geological background and structure of the northwestern coastline of Costa Rica

Introducción: El noroeste de Costa Rica representa un área con una alta riqueza geológica que evidencia rocas generadas en diferentes ambientes, como plano abisal, talud continental, plataforma continental, volcánicos continentales efusivos y explosivos, volcánicos submarinos, costeros, erosivos fluviales, erosivos gravitatorios, de intrusión y el manto terrestre. Objetivo: Describir el contexto geológico actual del noroeste de Costa Rica, mediante el análisis y recopilación de datos de campo, geoquímicos, petrológicos, estructurales, paleontológicos y radiométricos con el fin de reconstruir los eventos y etapas de deformación desde el Jurásico a la actualidad. Métodos: Se hizo una recopilación bibliográfica de estudios en diversas ramas geológicas para establecer un estado del arte del extremo noroeste de Costa Rica. Resultados: La historia geológica representa el registro en rocas desde el Jurásico al Holoceno - actualidad con once etapas definidas a partir de interpretaciones petrológicas, estructurales, geocronológicas, estratigráficas y paleontológicas agrupadas en tres etapas de depositación y afectadas por cuatro fases tectónicas compresivas. Conclusiones: Las etapas de depositación corresponden con una primera etapa magmática con afinidad oceánica entre el Jurásico y el Cretácico Superior, una segunda etapa predominantemente sedimentaria con un rango de edad entre el Cretácico Superior y el Oligoceno, y por último una etapa volcánica efusiva - explosiva ubicada desde el Plioceno hasta la actualidad. Las fases tectónicas se asocian con diferentes eventos entre placas. La primera se dio en el Cretácico Inferior y su principal resultado visible fue la emersión de la Ofiolita de Santa Elena; la segunda sucedió del Cretácico Superior temprano al Campaniano y su reconocimiento está ligado a la posición actual de la Ofiolita de Santa Elena y sus estructuras circundantes; la tercera fase aconteció desde el Eoceno Superior hasta el Mioceno y su principal producto fue la generación del tren de pliegues sedimentarios del Golfo de Santa Elena a bahía de Salinas; y por último se presenta una fase Cuaternaria que originó el eje de basculamiento de los productos piroclásticos frente a la cordillera volcánica de Guanacaste.

Introduction: The Costa Rican northwestern coastline has a wide geological uniqueness that exhibits rocks from different environments such as abyssal plane, continental slope, continental platform, volcanic eruptions (effusive, explosive, and submarine), coastal, erosional (fluvial and gravitatory), intrusions, and upper mantle rocks. Objective: To present the geologic state of the art of the northwestern coastline of Costa Rica through the bibliographic review of geochemical, petrologic, structural, paleontological and geochronological data. Methods: A bibliographic revision was done to propose a state of the art of northwestern coastline of Costa Rica. Results: The geologic record shows eleven stages from the Jurassic to the Holocene. These stages were regrouped from petrologic, structural, geochronologic, stratigraphic and paleontological interpretations in three depositional stages and four compressive tectonic phases. Conclusions: The first depositional stage is volcanic with oceanic affinity between the Jurassic and Upper Cretaceous. The second has a sedimentary predominance with an age range between the Upper Cretaceous and the Miocene. The last depositional stage is a volcanic (effusive and explosive) from the Pliocene until the present. The tectonic phases are associate with different interactions between tectonic plates. The first phase triggered the Santa Elena Ophiolite obduction during the Early Cretaceous. The second phase occurred in the Early Upper Cretaceous - Campanian and is recognized by its current position and surrounding structures. The third phase lasted from the Upper Eocene until the Miocene and its main result was the folding of the Bahía de Salinas sedimentary rocks. Finally, the Quaternary phase created a tilt axis and the aperture for the deposition of pyroclastic density currents in front of the Guanacaste volcanic ridge.

Statistical evaluation of the geochemical variability in overbank sediments in Spain.

The overall objective of this study is to estimate, detect and specify the main sources of variance which affect the contents of the different elements in overbank sediments across Spain. These sources of variance were assessed and compared by means of a series of analyses of variance (ANOVAs), by regarding two parameters: their significance and their contribution to the total variance. Overbank sediments, sampled in erosion banks, were studied in several locations, in basins which drain different types of geological backgrounds and land uses (urban, mining, agricultural or pristine) across the Iberian Peninsula. Forty-eight elements (mostly in the < 63 µm fraction) were analysed by ICP-OES, ICP-MS and INAA. After an isometric log ratio (ilr) transformation of the data, three ANOVA analyses were performed considering three perspectives: (1) local scale, (2) regional scale: within-profile perspective and (3) regional scale: inter-profile perspective. On a local scale, it was observed that the variability of rare earth elements (REE) depends mostly on the grain size and that heavy metals are also influenced by depth. In the analysis carried out on a regional scale, from a within-profile perspective, depth and duplicates do not influence significantly the variability of the element contents. Finally, from an inter-profile perspective, the selected sources of variance were land use and provenance, whose significance is the highest. While grain size and the selection of depth are of crucial importance in the final results, on local studies, land use and provenance are the ones that influence the most the composition of sediments in regional studies.

Geological influence on the activity concentration of terrestrial radionuclides 226Ra, 232Th and 40K in the Jos Plateau, Nigeria.

The Jos Plateau has been reported to have elevated levels of natural background radiation. A few earlier studies have measured the levels of natural radioactivity for specific locations in the area. Our interest is to investigate how geology of the study area influences the activity concentrations of natural radionuclides. Thus, the activity concentrations of terrestrial radionuclides in soil samples collected across the geological formations of the Jos Plateau were determined by gamma spectrometry technique. The mean activity concentrations of 226Ra, 232Th and 40K were found to exceed their corresponding world reference values of 35, 40 and 400 Bq kg-1, respectively. Data were compared using statistical methods, analysis of variance (ANOVA) and post hoc tests. The results revealed in some instances significant influences of geological types on the activity concentrations in the area. The spatial distribution maps of activity concentrations of 226Ra, 232Th and 40K were geostatistically interpolated by ordinary Kriging method using ArcGIS software.

Bomb-cratered coral reefs in Puerto Rico, the untold story about a novel habitat: from reef destruction to community-based ecological rehabilitation / Los arrecifes de coral con craters-bomba en Puerto Rico, la historia no contada sobre un hábitat inusual: desde la destrucción de arrecifes hasta la rehabilitación ecológica basada en la comunidad

Ecological impacts of military bombing activities in Puerto Rico have often been described as minimal, with recurrent allegations of confounding effects by hurricanes, coral diseases and local anthropogenic stressors. Reef craters, though isolated, are associated with major colony fragmentation and framework pulverization, with a net permanent loss of reef bio-construction. In contrast, adjacent non-bombarded reef sections have significantly higher benthic spatial relief and biodiversity. We compared benthic communities on 35-50 year-old bomb-cratered coral reefs at Culebra and Vieques Islands, with adjacent non-impacted sites; 2) coral recruit density and fish community structure within and outside craters; and 3) early effects of a rehabilitation effort using low-tech Staghorn coral Acropora cervicornis farming. Reef craters ranged in size from approximately 50 to 400m² and were largely dominated by heavily fragmented, flattened benthos, with coral cover usually below 2% and dominance by non-reef building taxa (i.e., filamentous algal turfs, macroalgae). Benthic spatial heterogeneity was lower within craters which also resulted in a lowered functional value as fish nursery ground. Fish species richness, abundance and biomass, and coral recruit density were lower within craters. Low-tech, community-based approaches to culture, harvest and transplant A. cervicornis into formerly bombarded grounds have proved successful in increasing percent coral cover, benthic spatial heterogeneity, and helping rehabilitate nursery ground functions.

Los impactos ecológicos de las actividades militares de bombardeos en Puerto Rico se han descrito a menudo como mínimos, con recurrentes denuncias al confundir efectos por huracanes, enfermedades de corales y estresores antropogénicos locales. Los cráteres de arrecife, aunque aislados, están relacionados con una alta fragmentación de la colonia y pulverización del contorno, con una pérdida neta permanente de arrecife en bio-construcción. En contraste, secciones adyacentes de arrecife no bombardeado tienen mayor biodiversidad y mayor relieve espacial bentónico. Comparamos las comunidades bentónicas en cráteres-bomba de arrecifes de coral con 35-50 años de antigüedad en las islas de Vieques, Puerto Rico, en comparación con los sitios adyacentes no impactados; 2) la densidad de reclutamiento de coral y estructura de la comunidad de peces dentro y fuera de los cráteres; y 3) impactos preliminares de un esfuerzo de rehabilitación basado en la comunidad arrecifal usando tecnología simple con el cultivo del coral Staghorn Acropora cervicornis. Los cráteres de arrecife se distancian en tamaño de aproximadamente 50 a 400m² y fueron dominados ampliamente por fragmentos de bentos aplanado, con una cubierta de coral generalmente por debajo de 2% y el predominio de taxones no constructores de arrecifes (es decir, tapetes de algas filamentosas, macroalgas). La heterogeneidad espacial bentónica fue significativamente menor dentro de cráteres que también resultaron en un reducido valor funcional como tierra de vivero de peces. La riqueza de especies de peces, abundancia y biomasa y densidad coral recluta fueron significativamente menores dentro de cráteres. Tecnología simple, basada en los enfoques de cultivo de comunidad, la cosecha y transplante de A. cervicornis en terrenos anteriormente bombardeados han demostrado un éxito al aumentar el porcentaje de cobertura de coral, la heterogeneidad espacial bentónica y ayudando a rehabilitar funcionalmente la tierra para vivero.

Neural network based automated interpretation algorithm for combined geophysical soundings in coastal zones.

The fresh water availability in coastal aquifers is an important problem faced by a major part of world's population dwelling there. For in situ and dynamic characterization of seawater encroachment into coastal aquifers, electrical geophysical methods are better suited. However, the resolving power of such data in distinguishing saline sands from moist clays in the subsurface is very poor. To meet this aspect and also the problem of analyzing voluminous data sets, we propose a feed forward back-propagation neural network (BPNN) based approach for the analysis of combined vertical electrical and induced polarization soundings. Our method is tested on synthetic data computed from available geo-electric sections and prevailing subsurface geological information of coastal aquifers of East Coast of India. The synthetic data are comprised of 18 combined soundings spread over five profiles. 15 out of 18 are used for training the BPNN, while 3 are used for testing. The trained BPNN (one node each in each of the input and output layers and 18 hidden nodes) showed 84.85% accuracy in testing phase for distinguishing clays from saline sands. Our method is also tested on real data concerning a shaly aquifer in Bahia, Brazil yielding an overall accuracy of 84.9%, comparable to that of synthetic case; thereby validating our approach.

[Features of cardiologic diseases prevalence among individuals engaged into geology].

The authors studied influence of work and living conditions on geologists' health state. A cohort under study demonstrated paradoxical epidemiologic situation with high mortality risk of arterial hypertension, whereas prevalence of this disease among specialists working on expeditionary shifted mode is nearly equal to that among general population. Unfavorable situation concerning cardiovascular diseases could be caused by high prevalence of cardiovascular risks and some occupational hazards. The situation is deteriorated by insufficient medical observation and psychologic traits of the specialists working on expeditionary shifted mode.

Concentration of 226Ra in rocks of the southern part of Lower Silesia (SW Poland).

The aim of the article is to present a preliminary description of rocks in the southern part of Lower Silesia and the Sudety Mountains in particular, with regard to 226Ra content. The research demonstrates that the average content of this isotope was 40.4 Bq/kg in the rocks of the southern part of Lower Silesia, and 41.7 Bq/kg in the rocks of the Sudetes. These values are slightly higher than the mean 226Ra content in the upper part of the Earth's crust, while the measured maximum content of this isotope (244 Bq/kg) is more than twice as high as the upper range of the values most frequently recorded in the upper part of the Earth's crust. The minimum values were lower than the detection limit, which was about 1 Bq/kg. These results reflect the mosaic-like geological structure of Lower Silesia, and particularly the Sudety Mountains, the occurrence of SiO2-rich igneous rocks and the products of their metamorphism, as well as numerous manifestations of uranium mineralisation or even deposital concentrations of this element. The rocks with the highest 226Ra contents include (in decreasing order): aplites, granites, gneisses and leucogranites, granite-gneisses, granodiorites and rhyolites, and, finally, mudstones. The lowest values of 226Ra content, on the other hand, were measured in sandstones, marls and conglomerates, and extremely low-in marbles and quartzites. The results show that background values of 226Ra content in the rocks of the southern part of Lower Silesia fall within a range from several to about 100 Bq/kg, which is the same as the range most frequently recorded in the upper part of the Earth's crust. Distribution of these values has log-normal character. The research demonstrates that the southern part of Lower Silesia, and the Sudetes in particular, may be marked by an increased radon potential. Particularly liable areas are: the Karkonosze granite massif, especially in its border zones, the Ladek-Snieznik and the Izera massifs, especially in their gneiss parts, as well as the Orlica massif. There is an urgent necessity to undertake more extensive and detailed research into this subject in connection with radiological hazard of radon. These statements are compliant with observed high radon concentrations in some groundwaters and tourist underground objects on these areas.

[Spatial heterogeneity of Ulmus pumila open forest ecosystem in Otindag sandy land].

Ulmus pumila open forest ecosystem plays an important role in Otindag Sandy Land. Using biostatistics (U test and Pearson correlation analysis) and geostatistics (Semivariogram analysis and Kriging interpolation), the spatial heterogeneity of soil organic matter (SOM), soil total nitrogen (STN), soil pH (SpH), soil water content (SWC), and the correlations between them, as well as their relationships with herbaceous coverage (HC) and Ulmus pumila were studied. The results showed that spatial heterogeneity and spatial dependence were apparent in SOM, STN, SpH, SWC and HC. With Kriging interpolation in the sampling area, the contour maps for these properties were drawn. Using the contour maps, the relationship between Ulmus pumila and spatial heterogeneity of the studied properties was investigated. The possible reason for the formation of fertility island under the Ulmus pumila trees in the ecosystem was discussed.

Multi-scale variation in spatial heterogeneity for microbial community structure in an eastern Virginia agricultural field.

To better understand the distribution of soil microbial communities at multiple spatial scales, a survey was conducted to examine the spatial organization of community structure in a wheat field in eastern Virginia (USA). Nearly 200 soil samples were collected at a variety of separation distances ranging from 2.5 cm to 11 m. Whole-community DNA was extracted from each sample, and community structure was compared using amplified fragment length polymorphism (AFLP) DNA fingerprinting. Relative similarity was calculated between each pair of samples and compared using geostatistical variogram analysis to study autocorrelation as a function of separation distance. Spatial autocorrelation was found at scales ranging from 30 cm to more than 6 m, depending on the sampling extent considered. In some locations, up to four different correlation length scales were detected. The presence of nested scales of variability suggests that the environmental factors regulating the development of the communities in this soil may operate at different scales. Kriging was used to generate maps of the spatial organization of communities across the plot, and the results demonstrated that bacterial distributions can be highly structured, even within a habitat that appears relatively homogeneous at the plot and field scale. Different subsets of the microbial community were distributed differently across the plot, and this is thought to be due to the variable response of individual populations to spatial heterogeneity associated with soil properties.

Including geophysical data in ground water model inverse calibration.

A nonlinear regression method is developed that can be used to estimate parameters of a ground waterflow model from a combination of observations of hydrological variables and observations of geophysical properties that are functionally related with the hydraulic conductivity. The procedure estimates: parameters characterizing the hydraulic conductivity field (e.g., zonal or pilot point values); geophysical properties that have been observed and that are functionally related with the hydraulic conductivity parameters; and a few parameters of the function that relates the hydraulic conductivity parameters with the geophysical properties (the type of function is assumed known). A fidelity factor, sigma(r)2, of a term of the minimized objective function reflects the faith one has in the validity of this functional relationship. The estimation methodology has been tested by means of synthetic models. The experimental results demonstrate that the number of estimated hydraulic conductivity parameters can be increased by adding geophysical observations to the set of hydrological observations that are traditionally used for model calibration. The improvement of the estimated hydraulic conductivity field and the simulated hydraulic head field can be significant but is dependent on the number, the locations, and the uncertainty of geophysical observations. The sensitivity of the estimation results to the value of sigma(r) is small for the studied problems except when the uncertainty of geophysical observations is high. In the latter case, a large sigma(r) value was found to be optimal to avoid that hydraulic conductivity estimates are closely tied to corresponding but highly uncertain geophysical observations.

Cokriging optimization of monitoring network configuration based on fuzzy and non-fuzzy variogram evaluation.

A number of optimization approaches regarding monitoring network design and sampling optimization procedures have been reported in the literature. Cokriging Estimation Variance (CEV) is a useful optimization tool to determine the influence of the spatial configuration of monitoring networks on parameter estimations. It was used in order to derive a reduced configuration of a nitrate concentration monitoring well network. The reliability of the reduced monitoring configuration suffers from the uncertainties caused by the variographer's choices and several inherent assumptions. These uncertainties can be described considering the variogram parameters as fuzzy numbers and the uncertainties by means of membership functions. Fuzzy and non-fuzzy approaches were used to evaluate differences among well network configurations. Both approaches permitted estimates of acceptable levels of information loss for nitrate concentrations in the monitoring network of the aquifer of the Plain of Modena, Northern Italy. The fuzzy approach was found to require considerably more computational time and numbers of wells at comparable level of information loss.

[Application of geostatistical analysis in study of macrophytes community's pattern].

In this paper, the Log-log semivariogram and theoretical models of Geostatistics and fractal theory were applied to study the mosaic spatial pattern of macrophytes community and its different characters on different scales in littoral zone of the Baoan Lake, Hubei Province. The macrophytes in the littoral zone with a square area of 25,600 m2, which was divided in to 1,024 grids, were investigated, and the macrophytes species in every grid were recorded. If the species was discovered in the grid, the variable of the grid had a value of 1, otherwise, it was 0. The community had a value of the total species discovered in the grid. By the method, the sampling results of Vallisneria spiralis, Ceratophyllum oryzetorum, Myriophyllum spicatum and the community were made into four worksheets. With help of the software GS+ for Windows to treat these worksheets, we got their log-log semivariogram, isotropic model, anisotropic model, isotropic fractal dimension and anisotropic fractal dimension. According to the results of isotropic models, the exponential models of Vallisneria spiralis and Myriophyllum spicatum indicated their contagious distribution with a high proportion of spatial structure, which indicated that their distribution changed lesser at microscale. In their range of 48.6 m and 34.2 m respectively, they had a spatial autocorrelation. Isotropic models of Ceratophyllum oryzetorum and the community were linear, and their proportions of spatial structure were low, indicating their great change at microscale and their moderate contagious distribution. The isotropic fractal dimension and anisotropic fractal dimension results indicated the dominant pattern of Vallisneria spiralis in the direction of NE-SW, the dominant patterns of Myriophyllum spicatum and Ceratophyllum oryzetorum in the direction of N-S, and the community's dominant pattern in the direction of NW-SE. Compared with other analysis tools, the method had advantages in deciding whether the community components had a spatial autocorrelation and in quantitatively classifying the aggregation density. It provides us a new way to quantitatively describe the pattern of the macrophytes community.

Geostatistics and Bayesian updating for transmissivity estimation in a multiaquifer system in Manitoba, Canada.

In ground water flow and transport modeling, the heterogeneous nature of porous media has a considerable effect on the resulting flow and solute transport. Some method of generating the heterogeneous field from a limited dataset of uncertain measurements is required. Bayesian updating is one method that interpolates from an uncertain dataset using the statistics of the underlying probability distribution function. In this paper, Bayesian updating was used to determine the heterogeneous natural log transmissivity field for a carbonate and a sandstone aquifer in southern Manitoba. It was determined that the transmissivity in m2/sec followed a natural log normal distribution for both aquifers with a mean of -7.2 and - 8.0 for the carbonate and sandstone aquifers, respectively. The variograms were calculated using an estimator developed by Li and Lake (1994). Fractal nature was not evident in the variogram from either aquifer. The Bayesian updating heterogeneous field provided good results even in cases where little data was available. A large transmissivity zone in the sandstone aquifer was created by the Bayesian procedure, which is not a reflection of any deterministic consideration, but is a natural outcome of updating a prior probability distribution function with observations. The statistical model returns a result that is very reasonable; that is homogeneous in regions where little or no information is available to alter an initial state. No long range correlation trends or fractal behavior of the log-transmissivity field was observed in either aquifer over a distance of about 300 km.

Utility of semivariogram for spatial variation of soil nutrients and the robust analysis of semivariogram.

The spatial variation of soil nutrients in topsoil (0-20 cm) was analyzed using semivariogram in the Zunhua County of Hebei Province, China. The effect on semivariogram with randomly deleted data and kriged estimates using various reduced sample sizes was also analyzed. The semivariograms of available N, total N, available P, organic matter were best described by a spherical model, except for available K, which best fitted a complex structure of exponential model and linear with sill model. The ratio of nugget to total sample variance ranged from 34.4% to 68.4%, indicating the spatial correlation of tested soil nutrients on a large scale was moderately dependent. Among five soil nutrients, available nitrogen and available phosphorus had the shortest spatial correlation range (5 km and 5.5 km), available K had the longest range (25.5 km), whereas total nitrogen and organic matter had intermediate spatial correlation range (14.5 km and 8.5 km). The semivariograms of available N, total N, available P, and organic matter were insensitive to a 50%-60% reduction in original sampling density, while for available K, it is up to 70%. The estimated spatial distributions of total N by kriging, under various reduced sample sizes, all correlated significantly (P = 0.001) with those obtained from original data. The results showed that the semivariogram was a relatively robust tool when used in a large region and sufficient spatial variation information could be retained regardless of a higher deletion proportion of the original data. The original sample data could be reduced by kriging and the estimates showed no loss of spatial information, however, the results may be unreliable unless a clearly identified semivariogram model could be obtained. The results may provide useful information for determining the appropriate sampling densities for these scales of soil survey.

Correlations in fossil extinction and origination rates through geological time.

Recent analyses have suggested that extinction and origination rates exhibit long-range correlations, implying that the fossil record may be controlled by self-organized criticality or other scale-free internal dynamics of the biosphere. Here we directly test for correlations in the fossil record by calculating the autocorrelation of extinction [corrected] and origination rates through time. Our results show that extinction rates are uncorrelated beyond the average duration of a stratigraphic interval. Thus, they lack the long-range correlations predicted by the self-organized criticality hypothesis. In contrast, origination rates show strong autocorrelations due to long-term trends. After detrending, origination rates generally show weak positive correlations at lags of 5-10 million years (Myr) and weak negative correlations at lags of 10-30 Myr, consistent with aperiodic oscillations around their long-term trends. We hypothesize that origination rates are more correlated than extinction rates because originations of new taxa create new ecological niches and new evolutionary pathways for reaching them, thus creating conditions that favour further diversification.

Cycle anatomy and variability in the storm-dominated type cincinnatian (Upper Ordovician): coming to grips with cycle delineation and genesis.

Although parasequence and sequence are scale-independent terms, they are frequently applied only to specific scales of cycles. For example, meter-scale cycles are commonly assumed to be parasequences or PACs. In the Upper Ordovician Kope and Fairview Formations of northern Kentucky, we examined a succession of 50 meter-scale cycles that have been variously interpreted as deepening-upward, shallowing-upward, or showing no relationship with water depth. Our analysis shows that these cycles, characterized by shifts in storm-bed proximality, are highly variable in their thickness and internal construction. Most cycles are best considered high-frequency sequences, because deepening-upward intervals are common, and many cycles contain evidence of abrupt basinward shifts in facies as expected at sequence boundaries. A minority fit the parasequence model of shallowing-upward cycles bounded by flooding surfaces. Larger, 20 m scale cycles are defined by systematic thickening and thinning trends of meter-scale cycles. However, meter-scale cycles do not display any systematic trends in cycle anatomy as a function of position within the 20 m cycles or position within the Kope and Fairview Formations. The high cycle variability and the lack of systematic stratigraphic organization with respect to longer-term cyclicity reflect either the irregularity of relative sea-level changes, the poor recording of sea-level changes in this deep-water setting, or the generation of these cycles by climate-induced cyclicity in storm intensity. These three mechanisms would generate similar patterns at the outcrop scale, so it is not possible at the present to distinguish between them.

Organic geochemistry of sediments from the continental margin off southern New England, U.S.A.--Part I. Amino acids, carbohydrates and lignin.

Total organic carbon (TOC), lignin, amino acids, sugars and amino sugars were measured in recent sediments for the continental margin off southern New England. The various organic carbon fractions decreased in concentration with increasing distance from shore. The fraction of the TOC that was accounted for by these major components also decreased with increasing distance from shore. The concentration of lignin indicated that only about 3-5% of the organic carbon in the nearshore sediment was of terrestrial origin. The various fractions were highly correlated, which was consistent with a simple linear mixing model of shelf organic matter with material form the slope and rise and indicated a significant transport of sediment from the continental shelf to the continental slope and rise.

Numerical simulations and the problem of periodicity in the cratering record.

Ages of craters in the record of impacts on earth may be uniformly period, totally random, or a mixture of the two. These alternatives are studied through numerical simulation wherein time-series analysis is performed on real and simulated sequences to which random noise has been added to represent age-dating uncertainty. We conclude that the real record is most likely to have been generated by a mixture of random and periodic impacts, with the random events constituting the majority.