A study of cross-shore maximum upwelling intensity along the Northwest Africa coast

Satellite images of sea surface temperature (SST) show that the location of cross-shore SST minimum (LCSM) stretches along the isobaths in the Northwest Africa Upwelling System. To understand and interpret these observations better, we set up a two-dimensional analytical model that takes into account the surface and bottom Ekman transport and the alongshore geostrophic current, as well as bottom friction and variations in bottom topography. The structure of vertical velocity with a realistic topography clearly illustrates the variations of SST drop in a sample cross-shore section. Some idealized theoretical model experiments are carried out to examine the effects of eddy viscosity, Coriolis force, and cross-shore wind on the location of the cross-shore maximum upwelling intensity. The results show that the cross-shore wind largely impacts on the location where the coldest water outcrops to the surface through an adjustment of the cross-shore pressure gradient. This is also verified by the remotely sensed data, which indicate that the maximum correlation coefficient between cross-shore wind stress and the depth of LCSM is ?0.65 with a lag of approximately 1 day.     

Development of a Vector-Based Method for Coastal Bluffline Mapping Using LiDAR Data and a Comparison Study in the Area of Lake Erie

This paper introduces a new method for bluffline extraction based on analysis of surface normal vectors derived from a 2.5-dimensional Delaunay triangulation network generated using a set of LiDAR points. The developed method has been implemented and tested for the extraction of blufflines in a section of the Lake Erie coastline near Painesville, Ohio. Comparison of the research results with those generated by a method for bluffline extraction based on slope information along elevation profiles on transects shows an improvement in both horizontal and vertical accuracies in situations where the bluff has a relatively sharp edge formation.     

Climatic changes in Eurasia and Africa at the last glacial maximum and mid-Holocene: reconstruction from pollen data using inverse vegetation modelling

In order to improve the reliability of climate reconstruction, especially the climatologies outside the modern observed climate space, an improved inverse vegetation model using a recent version of BIOME4 has been designed to quantitatively reconstruct past climates, based on pollen biome scores from the BIOME6000 project. The method has been validated with surface pollen spectra from Eurasia and Africa, and applied to palaeoclimate reconstruction. At 6 cal ka BP (calendar years), the climate was generally wetter than today in southern Europe and northern Africa, especially in the summer. Winter temperatures were higher (1–5°C) than present in southern Scandinavia, northeastern Europe, and southern Africa, but cooler in southern Eurasia and in tropical Africa, especially in Mediterranean regions. Summer temperatures were generally higher than today in most of Eurasia and Africa, with a significant warming from ∼3 to 5°C over northwestern and southern Europe, southern Africa, and eastern Africa. In contrast, summers were 1–3°C cooler than present in the Mediterranean lowlands and in a band from the eastern Black Sea to Siberia. At 21 cal ka BP, a marked hydrological change can be seen in the tropical zone, where annual precipitation was ∼200–1,000 mm/year lower than today in equatorial East Africa compared to the present. A robust inverse relationship is shown between precipitation change and elevation in Africa. This relationship indicates that precipitation likely had an important role in controlling equilibrium-line altitudes (ELA) changes in the tropics during the LGM period. In Eurasia, hydrological decreases follow a longitudinal gradient from Europe to Siberia. Winter temperatures were ∼10–17°C lower than today in Eurasia with a more significant decrease in northern regions. In Africa, winter temperature was ∼10–15°C lower than present in the south, while it was only reduced by ∼0–3°C in the tropical zone. Comparison of palaeoclimate reconstructions using LGM and modern CO₂ concentrations reveals that the effect of CO₂ on pollen-based LGM reconstructions differs by vegetation type. Reconstructions for pollen sites in steppic vegetation in Europe show warmer winter temperatures under LGM CO₂ concentrations than under modern concentrations, and reconstructions for sites in xerophytic woods/scrub in tropical high altitude regions of Africa are wetter for LGM CO₂ concentrations than for modern concentrations, because our reconstructions account for decreased plant water use efficiency.     

Choosing Between Two Kind of Sampling Patterns Using Geostatistical Simulation: Regularly Spaced or at High Uncertainty Locations?

Data from a mineral deposit are commonly obtained by core drilling. This kind of sampling involves high costs, limiting the number of drill holes. Additional holes should be located to bring the maximum benefit. The benefit can be evaluated by various ways and must take into account the goals of sampling. This article presents a case study where the deposit has been sampled and the new drillings must be added to reduce the uncertainty about a transfer function, Net Present Value (NPV) of the mining project. There are basically two ways to choose locations where new drillings should be placed for cases where the aim of sampling is to reduce uncertainty about a global function: the addition of new drillings outlining a quasi regular grid with previously collected drillings or the addition of new drillings on the locations of high uncertainty about the attribute of interest (or the attribute that is considered most influential in the transfer function). The performances of these patterns on reducing the uncertainty measured by the function selected are compared. The results point out that the most efficient pattern relates to the distribution (histogram) of the uncertainty about the attribute of interest. Thus, the choice of which sampling pattern should be adopted varies depending on data distribution and its influence on the transfer function.     

Gamma-Ray Data Processing and Integration for Lode-Au Deposits Exploration

A large number of mineral deposits are associated with hydrothermal processes, especially auriferous deposits. In such processes, studies on percolating fluids may indicate the presence of potash (K), among other elements. In this study, aerogammaspectrometric data-processing methodologies are evaluated, especially those methods based on the suppression of the primary contribution of potassium, the result of lithological and soil variations, and to environmental conditions. Resulting maps point out the contribution of hydrothermal K. This processing procedure was used because of the association of hydrothermal K and auriferous mineralizations according to the deposit model defined for the studied region. Intensity maps locate the areas with great influence of hydrothermal K. Data integration required to improve a change in the gammaspectrometric data processing in order to positively correlate hydrothermalised areas. Data integration could distinguish high and medium favorable targets for mineral exploration of lode-Au deposits in the studied region.     

Long-term TNT and DNT contamination: 1-D modeling of natural attenuation in the vadose zone: case study,Portugal

The vadose zone of a trinitrotoluene (TNT) and dinitrotoluene (DNT) contaminated site was investigated to assess the mobility of those explosives under natural conditions. Located in the left margin of the River Tejo Basin, Portugal, the site is located on unconsolidated sediments. Wastewaters associated with the 50-year explosives production were disposed in excavated ponds, from where water would infiltrate and pollute the unsaturated and saturated parts of the local aquifers. Two boreholes were drilled to 9 m depth in such a former waste pond to investigate the contaminant’s fate in the vadose zone. Sediment samples were taken every 1–2 m for analysis of the poly-nitroaromatics (p-NACs) and organic volatile compounds, pH, organic carbon content, cation exchange capacity and grain size analysis. The main contaminant was TNT representing >70 % of the total p-NACs concentration that peaked approximately 7 mg/kg in one borehole, even if the median in both boreholes was of ~1 mg/kg. DNT was 4–30 % of the total p-NACs and nitrotoluene (NT), up to 5 %. No other (volatile) organic compound was detected. The predominance of TNT as the main contaminant implies that any natural mass reduction has been inefficient to clean the site. Several 1-D model simulations of p-NACs cleaning of the vadose zone under natural conditions indicated that the most probable scenario of combined advection and partitioning will only remove TNT after 10’s of years, whereas DNT and NT will hardly be removed. Such low concentrations and long times for the p-NACs removal, suggest that by now those compounds have been washed-out to a level below standard limits.     

Observing the critical zone on a critical budget: The Peri Lake experimental catchment

The Peri Lake Experimental Catchment, in Southern Brazil, is a small coastal lake (5 km², 7 m depth) with a total catchment area of 20 km². The catchment is mainly covered by subtropical Atlantic Forest. Peri Lake is a recreational destination, an important ecosystem for biodiversity preservation, and a major water supply source. Even though there is a clear social and economic relevance, the information on this ecosystem is scarce, especially regarding the dynamics of water balance and water quality. We built this observatory using a low-cost, low-technology monitoring system to estimate the major components of the water budget on what we called a critical budget. We monitor meteorology, rainfall interception in two plots, overland flow and groundwater connectivity in a representative hillslope; and streamflow and velocity are measured in several small headwater rivers mostly starting during 2015. Geographic information on topography, land cover, geology, soil characteristics and the location of all the equipment installed in the field are also provided. The data set can be used to understand rainfall interception in Atlantic forests, catchment connectivity and streamflow permanence in coastal areas, spatial patterns of baseflow, and the modelling of complex processes in the critical zone involving the interaction between surface and groundwater, that are important in coastal lake ecosystems. Besides being used for research and education, we hope this observatory built on a shoestring budget will encourage fieldwork particularly in underrepresented and underfunded regions of the world.     

Rates of silicate dissolution in deep-sea sediment: In situ measurement using U/U of pore fluids

Bulk dissolution rates for sediment from ODP Site 984A in the North Atlantic are determined using the ²³⁴U/²³⁸U activity ratios of pore water, bulk sediment, and leachates. Site 984A is one of only several sites where closely spaced pore water samples were obtained from the upper 60 meters of the core; the sedimentation rate is high (11-15 cm/ka), hence the sediments in the upper 60 meters are less than 500 ka old. The sediment is clayey silt and composed mostly of detritus derived from Iceland with a significant component of biogenic carbonate (up to 30%).The pore water ²³⁴U/²³⁸U activity ratios are higher than seawater values, in the range of 1.2 to 1.6, while the bulk sediment ²³⁴U/²³⁸U activity ratios are close to 1.0. The ²³⁴U/²³⁸U of the pore water reflects a balance between the mineral dissolution rate and the supply rate of excess ²³⁴U to the pore fluid by α-recoil injection of ²³⁴Th. The fraction of ²³⁸U decays that result in α-recoil injection of ²³⁴U to pore fluid is estimated to be 0.10 to 0.20 based on the ²³⁴U/²³⁸U of insoluble residue fractions. The calculated bulk dissolution rates, in units of g/g/yr are in the range of 4 × 10⁻⁷ to 2 × 10⁻⁶ yr⁻¹. There is significant down-hole variability in pore water ²³⁴U/²³⁸U activity ratios (and hence dissolution rates) on a scale of ca. 10 m. The inferred bulk dissolution rate constants are 100 to 10⁴ times slower than laboratory-determined rates, 100 times faster than rates inferred for older sediments based on Sr isotopes, and similar to weathering rates determined for terrestrial soils of similar age. The results of this study suggest that U isotopes can be used to measure in situ dissolution rates in fine-grained clastic materials.The rate estimates for sediments from ODP Site 984 confirm the strong dependence of reactivity on the age of the solid material: the bulk dissolution rate (R_d) of soils and deep-sea sediments can be approximately described by the expression R_d ≈ 0.1 Age⁻¹ for ages spanning 1000 to 5 × 10⁸ yr. The age of the material, which encompasses the grain size, surface area, and other chemical factors that contribute to the rate of dissolution, appears to be a much stronger determinant of dissolution rate than any single physical or chemical property of the system.     

Transition from debris flow to hyperconcentrated flow in a submarine channel (the Cretaceous Cerro Toro Formation, southern Chile)

ABSTRACT It is important to understand the exact process whereby very large amounts of sediment are transported. This paper reports peculiar conglomerate beds reflecting the transition of submarine debris flows into hyperconcentrated flows, something that has been well documented only in subaerial debris-flow events until now. Voluminous debris flows generated along a Cretaceous submarine channel, southern Chile, transformed immediately into multiphase flows. Their deposits overlie fluted or grooved surfaces and comprise a lower division of clast-supported and imbricated pebble–cobble conglomerate with basal inverse grading and an upper division of clast- to matrix-supported, disorganized conglomerate with abundant intraformational clasts. The conglomerate beds suggest temporal succession of turbidity current, gravelly hyperconcentrated flow, and mud-rich debris flow phases. The multiphase flows resulted from progressive dilution of gravelly but cohesive debris flows that could hydroplane, in contrast to the flow transitions in subaerial environments, which involve mostly non-cohesive debris flows. This finding has significant implications for the definition, classification, and hazard assessment of submarine mass-movement processes and characterization of submarine reservoir rocks.     

Robust Rule-Based Aggradational Lobe Reservoir Models

Stratigraphic rule-based reservoir models approximate sedimentary dynamics to generate numerical models of reservoir architecture with realistic spatial distributions of petrophysical properties for reservoir forecasting and to support development decision making. A few intuitive rules for the sequential placement of surfaces bounding reservoir units render realistic reservoir heterogeneity, continuity, and spatial organization of petrophysical property distributions that are difficult to obtain using conventional geostatistical pixel- and object-based subsurface methods. While these methods are emerging in applications specifically for deepwater and fluvial clastic reservoirs, there are some remaining obstacles to broad application, such as selection of rule parameters and addressing emergent non-stationarities over the sequence of the placed surfaces. Firstly, there is a need to tune rule parameters to ensure the models honor natural heterogeneities. We demonstrate this for the compensational rule. Secondly, invariants over model sequence (from the base to the top of the model) may occur with respect to shape, volume, undulation, and gradients of surfaces. For example, for a stack of compensational lobes, the volume of individual lobes may decrease due to the onlapping of previous bathymetry and also increasing undulation over model sequence may occur. In addition, for stacking of compensational lobes, the interfacial width and average gradient of the composite surface may initially increase, but then saturate and stabilize. Such non-stationarities represent numerical artifacts that may bias the results from these rule-based models. It is essential that these features are quantified and mitigated as a prerequisite for robust application of rule-based aggradational lobe methods for reservoir modeling.