High-resolution seismic stratigraphy of the late Neogene of the central sector of the Colombian Pacific continental shelf: A seismic expression of an active continental margin

The sedimentary prism of the central Pacific continental shelf of Colombia was affected by regional folding and faulting, and probably later mud diapirism, from the Late Miocene to the Holocene. Interpretation of high-resolution seismic lines (2 s/dt) revealed that the prism consists of 13 high-resolution seismic units, that can be separated into 5 seismic groups.Deposition of the prism and the associated stacking pattern, are probably the response to variable uplift and subsidence in a fore-arc basin that underwent important tectonic events by the end of the Miocene. Throughout the Pliocene, the continental shelf sedimentation was affected by the growing of a dome structure probable due to mud diapirism. This fact caused peripheral faults both normal and reverse that controlled the distribution of some of the seismic units. During the Late Pleistocene (Wisconsin stage?) a eustatic sea level fall caused the shoreline to advance about 50 km westward of its present position. Because of this eustatic sea level change, a strong fluvial dissection took place and is interpreted as the probable extension of the San Juan River to the south of the present day river mouth. Within this framework it is believed that the Malaga and Buenaventura Bays were the passageways of branches of the old drainage system of the San Juan River. The inner branch circulated through the present Buenaventura Bay and runs southward leaving the mark of an apparent valley identified in the seismic information in the eastern sector of the study area. This old fluvial valley and its filling material located in the present day inner continental shelf front of Buenaventura are postulated as important targets to find placer minerals such as gold and platinum.     

Relationships between land-use and forced-pool characteristics in the Colorado Front Range

Land use practices in Colorado during the last two centuries altered the supply of sediment and water to many channels in the upper South Platte Basin. As a result of increased supply of sediment and mobility and reduced peak flows, the characteristics of pools associated with channel constrictions, referred to as forced pools, may have been altered. Increased supply of sediment and reduced transport capacity of high flows could lead to aggradation in forced pools. Channel confined by road corridors could lead to high velocities at normal flows, increased energy dissipation from riprap, or even increased pool frequency resulting from failed riprap. To assess potential alterations, four hypotheses were tested: (1) impacted streams will show significantly different mean volume of pools than a control stream; (2) mean volume of pools on streams where land-use activities increased the supply of sediment will be significantly different from streams solely affected by flow regulation; (3) the strongest change in characteristics of pools of impacted streams will be a reduced volume of pools; (4) streams affected by road corridors will show statistically lower spacing of pools than streams unaffected by roads. The downstream spacing and residual volume of twenty consecutive forced pools were surveyed on five streams in the Colorado Front Range that varied from no contemporary impact to multiple historical and contemporary impacts. ANCOVA with stepwise model selection indicated that degree of land-use (categorical), bankfull spacing of pools, upstream riffle slope and expansion ratio were all significant (α = 0.1) predictors of volume of pools (R² = 0.73). Simple linear regression of mean volume of pools and stream specific variables (gradient, drainage area and discharge) and least square means comparison of mean volume of pools indicated a need to standardize volumes of pools by slope and discharge so that the volumes of pools could be compared among different levels of land-use. Significant correlations between drainage area and volume of pools allowed volume of pools to be standardized by drainage area and thereby discharge. This dimensionless variable was also significantly correlated with channel slope, which permitted the construction of a new variable, PV_QS (volume of pools standardized by discharge and slope). Least square means comparison of mean PV_QS revealed that the control reach was significantly different from road-impacted reaches. Mean volume of pools was significantly larger in the control reach compared to all but one road-impacted stream. This was likely a function of higher wood loading in the control reach and the competence of high flows in the road-impacted reach. Streams affected by road corridors did not have significantly different bankfull spacing of pools from streams not impacted by roads. The multiple interactions among control and response variables explored in this study indicate the need to identify the most constrained and sensitive response variables when attempting to assess channel response to land use.     

Dual Kriging with Local Neighborhoods: Application to the Representation of Surfaces

Ordinary kriging, in its common formulation, is a discrete estimator in that it requires the solution of a kriging system for each point in space in which an estimate is sought. The dual formulation of ordinary kriging provides a continuous estimator since, for a given set of data, only a kriging system has to be estimated and the resulting estimate is a function continuously defined in space. The main problem with dual kriging up to now has been that its benefits can only be capitalized if a global neighborhood is used. A formulation is proposed to solve the problem of patching together dual kriging estimates obtained with data from different neighborhoods by means of a blending belt around each neighborhood. This formulation ensures continuity of the variable and, if needed, of its first derivative along neighbor borders. The final result is an analytical formulation of the interpolating surface that can be used to compute gradients, cross-sections, or volumes; or for the quick evaluation of the interpolating surface in numerous locations.     

A non-parametric automatic blending methodology to estimate rainfall fields from rain gauge and radar data

Quantitative estimation of rainfall fields has been a crucial objective from early studies of the hydrological applications of weather radar. Previous studies have suggested that flow estimations are improved when radar and rain gauge data are combined to estimate input rainfall fields. This paper reports new research carried out in this field. Classical approaches for the selection and fitting of a theoretical correlogram (or semivariogram) model (needed to apply geostatistical estimators) are avoided in this study. Instead, a non-parametric technique based on FFT is used to obtain two-dimensional positive-definite correlograms directly from radar observations, dealing with both the natural anisotropy and the temporal variation of the spatial structure of the rainfall in the estimated fields. Because these correlation maps can be automatically obtained at each time step of a given rainfall event, this technique might easily be used in operational (real-time) applications. This paper describes the development of the non-parametric estimator exploiting the advantages of FFT for the automatic computation of correlograms and provides examples of its application on a case study using six rainfall events. This methodology is applied to three different alternatives to incorporate the radar information (as a secondary variable), and a comparison of performances is provided. In particular, their ability to reproduce in estimated rainfall fields (i) the rain gauge observations (in a cross-validation analysis) and (ii) the spatial patterns of radar fields are analyzed. Results seem to indicate that the methodology of kriging with external drift [KED], in combination with the technique of automatically computing 2-D spatial correlograms, provides merged rainfall fields with good agreement with rain gauges and with the most accurate approach to the spatial tendencies observed in the radar rainfall fields, when compared with other alternatives analyzed.     

Validation of MODIS Vegetation Continuous Fields for monitoring deforestation and forest degradation: two cases in Mexico

This study assesses whether MODIS Vegetation Continuous Fields percent tree cover (PTC) data can detect deforestation and forest degradation. To assess the usefulness of PTC for detecting deforestation, we used a data set consisting of eight forest and seven non-forest categories. To evaluate forest degradation, we used data from two temperate forest types in three conservation states: primary (dense), secondary (moderately degraded) and open (heavily degraded) forest. Our results show that PTC can differentiate temperate forest from non-forest categories (p = 0.05) and thus suggests PTC can adequately detect deforestation in temperate forests. In contrast, single-date PTC data does not appear to be adequate to detect forest degradation in temperate forests. As for tropical forest, PTC can partially discriminate between forest and non-forest categories.     

Coupled climate impacts of the Drake Passage and the Panama Seaway

Tectonically-active gateways between ocean basins have modified ocean circulation over Earth history. Today, the Atlantic and Pacific are directly connected via the Drake Passage, which forms a barrier to the time-mean geostrophic transport between the subtropics and Antarctica. In contrast, during the warm early Cenozoic era, when Antarctica was ice-free, the Drake Passage was closed. Instead, at that time, the separation of North and South America provided a tropical seaway between the Atlantic and Pacific that remained open until the Isthmus of Panama formed in the relatively recent geological past. Ocean circulation models have previously been used to explore the individual impacts of the Drake Passage and the Panama Seaway, but rarely have the two gateways been considered together, and most explorations have used very simple atmospheric models. Here we use a coupled ocean–ice–atmosphere model (GFDL’s CM2Mc), to simulate the impacts of a closed Drake Passage both with and without a Panama Seaway. We find that the climate response to a closed Drake Passage is relatively small when the Panama Seaway is absent, similar to prior studies, although the coupling to a dynamical atmosphere does increase the temperature change. However, with a Panama Seaway, closing Drake Passage has a much larger effect, due to the cessation of deep water formation in the northern hemisphere. Both gateways alter the transport of salt by ocean circulation, with the Panama Seaway allowing fresh Pacific water to be imported to the North Atlantic, and the Drake Passage preventing the flow of saline subtropical water to the circum-Antarctic, a flow that is particularly strong when the Panama Seaway is open. Thus, with a Panama Seaway and a closed Drake Passage, the Southern Ocean tends to be relatively salty, while the North Atlantic tends to be relatively fresh, such that the deep ocean is ventilated from the circum-Antarctic. Ensuing changes in the ocean heat transport drive a bi-polar shift of surface ocean temperatures, and the Intertropical Convergence Zone migrates toward the warmer southern hemisphere. The response of clouds to changes in surface ocean temperatures amplifies the climate response, resulting in temperature changes of up to 9 °C over Antarctica, even in the absence of land-ice feedbacks. These results emphasize the importance of tectonic gateways to the climate history of the Cenozoic, and support a role for ocean circulation changes in the glaciation of Antarctica.     

Yucatán karst features and the size of Chicxulub crater

The buried Chicxulub impact structure is marked by a dramatic ring of sinkholes (called cenotes if containing water), and adjacent less prominent partial rings, which have been shown to coincide with maxima in horizontal gravity gradients and a topographic depression. These observations, along with the discreteness and spacing of the features, suggest a formation mechanism involving faulting in the outer slump zone of the crater, which would thus have a diameter of approximately 180 km.
An opposing view, based primarily on the interpretation of gravity data, is that (he crater is much larger than the cenote ring implies. Given the association of the known cenote ring with faults, we here examine northern Yucatan for similar rings in gravity, surface features and elevation, which we might expect to be associated with outer concentric faults in the case of a larger, possibly multiring, structure.
No such outer rings have been found, although definite patterns are seen in the distribution of karst features outside the crater rim. We explain these patterns as resulting mainly from deformation related to the block fault zone that parallels tbe shelf edge of eastern Yucatan.     

Patterns of biomass-size spectra from oligotrophic waters of the Northwest Atlantic

The study of the distribution of biomass by size provides an ataxonomic approach for analyzing the structure of the pelagic ecosystem. However, empirical data regarding planktonic size-structure in offshore areas are scarce. Here, we report the results of a study of the planktonic biomass size-distribution at several stations located within two regions (Sargasso Sea and New England Seamounts Area) of the Northwest Atlantic. The biomass size-spectra covered a body-size range from bacteria to mesozooplankton and a depth range from the surface to 400 m. It is shown that the slope of the normalized biomass size-spectrum (NBS-spectrum) varies depending on whether volume or carbon units are used. The transformation from volume to carbon units makes the slope of the NBS-spectrum approximately 0.15 units more negative. The distribution of normalized-biomass by size was linear (plotted on a log–log scale) at all stations. The slopes of the NBS-spectra (volume scale) ranged from −0.96 to −1.01. There were no significant differences among the slopes of the NBS-spectra within either of the two areas studied. In addition, no significant differences were detected between the stations in the Sargasso Sea and those located in the New England Seamounts area. Apart from a tendency towards a decrease in the intercept of the normalized-biomass axis of the size-spectra in deeper waters, the NBS-spectra were also very similar through depth. The slopes of the NBS-spectra in biovolume units are not significantly different from −1.0 (p<0.01), and therefore are in agreement with Sheldon’s Linear Biomass Hypothesis. By contrast, the slopes of the NBS-spectra in carbon units are significantly different (p<0.01) from −1.0 (range=–1.09 to –1.17; all stations TOGETHER=–1.14) and their numerical values are in close agreement with Platt and Denman’s model and with the findings of [Rodríguez and Mullin, 1986a]) for the North Pacific Central Gyre. The results of this study support the hypothesis that the planktonic size-structure of offshore oligotrophic systems is a conservative property.     

Geophysical evidence for a large impact structure on the Falkland (Malvinas) Plateau

A large, roughly circular structural basin is recognised on the Falkland (Malvinas) Plateau to the NW of West Falkland (Gran Malvina) Island (S 51°00′, W 62°00′). The basin, seen in seismic‐reflection profiles and evident as a large negative gravity anomaly, has a diameter of ~250 km. The age of the basin is estimated to be Late Palaeozoic. It is completely buried by younger sediments and has no topographic expression on the sea floor. We propose that the basin and geophysical anomalies, especially the combination of a large circular negative gravity anomaly with a rim of positive anomalies, and a marked circular series of positive magnetic anomalies in the same area, may be best explained by the presence of a large buried impact structure.     

Foreshocks and aftershocks of the M W = 7.1, 1992, earthquake in the Atrato region, Colombia

We study the October 18, M _W = 7.1, 1992 Atrato earthquake, and its foreshocks and aftershocks, which occurred in the Atrato valley, northwestern Colombia. The main shock was preceded by several foreshocksof which the M _W = 6.6, October 17 earthquacke was the largest. Inparticular, we examine foreshocks and aftershocks performing joint-hypocenter relocations using high quality Pn and Sn wave readingsfrom permanent regional networks. We observed a few hours prior to the main shock a sudden increase of foreshocks. Maybe this could be used as a predictor since foreshocks have been known for other major events in the region. Our locations align for 90 km with a trend of 5^° ±4^° in agreement with the Harvard CMT solution showing the faultplane trending 9^° to be the plane of rupture. In relation to theepicenter of the main shock, maximum intensities were located to thesouth, consistent with a rupture that traveled from north to south witha larger energy release in the south as suggested by an empirical Green'sfunction study (Li and Toksöz, 1993; Ammon et al., 1994). The boundarybetween the Panama and North Andes blocks has been placed close to thePanama-Colombia border as either a sharp boundary or a diffuse zone. TheAtrato earthquake, however, shows that the plate boundary between thePanama and North Andes microblocks is a diffuse deformation zone. Thiszone has a width of at least 2^° stretching from 78^°W to 76^°W. Quantification of earthquake moment release (during the past30 years) in this zone shows a similar amount of moment release in thewestern and eastern parts of this zone.