Time-calibrated estimates of oceanographic profiles using empirical orthogonal functions and clustering

Oceanographic climatology is widely used in different applications, such as climate studies, ocean model validation and planning of naval operations. Conventional climatological estimates are based on historic measurements, typically by averaging the measurements and thereby smoothing local phenomena. Such phenomena are often local in time and space, but crucial to some applications. Here, we propose a new method to estimate time-calibrated oceanographic profiles based on combined historic and real-time measurements. The real-time measurements may, for instance, be SAR pictures or autonomous underwater vehicles providing temperature values at a limited set of depths. The method employs empirical orthogonal functions and clustering on a training data set in order to divide the ocean into climatological regions. The real-time measurements are first used to determine what climatological region is most representative. Secondly, an improved estimate is determined using an optimisation approach that minimises the difference between the real-time measurements and the final estimate.     

Wind erosion characteristics of Sahelian surface types

The assessment of wind erosion magnitudes for a given area requires knowledge of wind erosion susceptibilities of the dominant local surface types. Relative wind erosion potentials of surfaces can hardly be compared under field conditions, as each erosion event is unique in terms of duration, intensity and extent. The objective of this study was to determine and compare relative wind erosion potentials of the most representative surface types over a transect comprising most parts of southwestern Niger. For this purpose, mobile wind tunnel experiments were run on 26 dominant surface types. The effects of surface disturbance were additionally determined for 13 of these surfaces. The results, namely measurements of wind fields and mass fluxes, can be classified according to specific surface characteristics. Three basic surface groups with similar emission behaviour and aerodynamic characteristics were identified: (1) sand surfaces, (2) rough stone surfaces and (3) flat crusted surfaces. Sand surfaces feature a turbulent zone close to the surface due to the development of a saltation layer. Their surface roughness is medium to high, as a consequence of the loss of kinetic energy of the wind field to saltating particles. Sand surfaces show the highest mass fluxes due to the abundance of loose particles, but also fairly high PM10 fluxes, as potential dust particles are not contained in stable crusts or aggregates. Rough stone surfaces, due to their fragmented and irregular surface, feature the highest surface roughness and the most intense turbulence. They are among the weakest emitters but, due to their relatively high share of potential dust particles, PM10 emissions are still average. Flat crusted surfaces, in contrast, show low turbulence and the lowest surface roughness. This group of surfaces shows rather heterogeneous mass fluxes, which range from moderate to almost zero, although the share of PM10 particles is always relatively high. Topsoil disturbance always results in higher total and PM10 emissions on sand surfaces and also on flat crusted surfaces. Stone surfaces regularly exhibit a decrease in emission after disturbance, which can possibly be attributed to a reorganization which protects finer particles from entrainment. The results are comparable with field studies of natural erosion events and similar wind tunnel field campaigns. The broad range of tested surfaces and the standardized methodology are a precondition for the future regionalization of the experimental point data. Copyright © 2010 John Wiley & Sons, Ltd.     

Bank erosion of legacy sediment at the transition from vertical to lateral stream incision

Streambank erosion is a primary source of suspended sediments in many waterways of the US Atlantic Piedmont. This problem is exacerbated where banks are comprised of fine sediment produced by the intensive land use practices of early European settlers. A stream in this region, Richland Creek incises into banks comprised of three stratigraphic layers associated with historic land use: pre‐European settlement, early European agriculture and development, and water‐powered milldam operation. This study aims to identify the bank processes along a reach of Richland Creek that is eroding towards its pre‐disturbance elevation. The volume of material that has eroded along this stream since the milldam breached was calculated by differencing a reconstructed surface of the pond bed and an aerial lidar digital terrain model (DTM). Immediately downstream from the study reach, the channel is floored by bedrock and immediately upstream the rate of channel erosion approximately doubled along the longitudinal profile of Richland Creek, which indicate that the study reach spans the transition from a channel dominated by vertical incision in the upstream direction to horizontal widening in the downstream direction. The combined hydrometeorological conditions and dominant processes causing reach‐scale cut bank erosion were investigated with analyses of stream stage, precipitation, and streambank volumetric and surfaces change that was measured during nine terrestrial lidar surveys in 2010–2012. The spatial variability of erosion during a simulated precipitation event was examined in a field‐based experiment. Erosion was greatest where mill pond sediment columns detached along vertical desiccation and horizontal seepage cracks. This sediment accumulated on the bank toe throughout the study and was a source of readily‐entrained fine sediment contrary to the upper reaches where depositional accommodation space is more limited. Findings suggest that hotspots of sediment excavation progress upstream, indicating that restoration efforts should focus upon stabilizing banks at these locations. Copyright © 2015 John Wiley & Sons, Ltd.     

Sandwaves on the Southeast Vietnam Shelf recorded by high resolution seismic profiles: formation and mechanism

The application of high resolution seismic data using boomer sound source has revealed a wide distribution of large-scale bedforms (sandwaves) on the Southeast Vietnam continental shelf. Bedforms that are a few meters high in wave height and hundreds of meters long in wavelength are primarily developed in the inner shelf (20–40 m) and considered to be formed under the present-day marine hydrodynamic conditions. Those bedforms developed in the deeper water (120 m) of the northernmost part of the continent can be interpreted as the relict morphological features formed during the latest sea-level lowstand of the late Pleistocene period. Two sediment transport paths have been identified on the basis of the bedform’s leeward orientation: northeast-southwest (along-shore) and north-south (cross-shore). A quantitative bottom current map is constructed from sandwave dimensions, surface sediments and measurement data. The strongest current velocities that gradually decrease toward the southwest are indicated by large sandwaves in the north (field B). Water depth, surficial sediment composition and bottom current are three factors that control the development of bedforms.     

Albeluvisol development with time in loamy marine sediments of southern Norway

This paper reports a study on soil development in loamy marine sediments on both sides of the Oslofjord (Vestfold, Østfold). This region has been subject to steady glacio-isostatic uplift during the whole Holocene. Hence, land surface age continuously increases from the coast inland. Several sea level curves, based on radiocarbon datings, enable estimation of land surface age for all locations. Clay illuviation starts in less than 1650 years. E horizons become lighter with age, but their lower boundary stays around 40 cm for more than 10 000 years. Albeluvic tongues develop between 4600 and 6200 years. Initially, they form along intersections of cracks. As preferential flow and leaching along the cracks continues, the tongues increase in length and width, progressively consuming the prisms between the cracks in the upper Bt horizon. The Fe_d/Fe_t ratios (weighted means of the upper meter) show a clear linear increase with soil age and may be used for “pedo-dating”, i.e. for estimating the ages of non-dated land surfaces covered with similar sediments. In contrast, the logarithmic decrease of base saturation and pH, with rapid changes in the first time but only very slight changes between 2000 and 10 000 years, makes these parameters unsuitable for “pedo-dating”.     

A late Holocene paleo-productivity record in the western Gulf of Maine,USA, inferred from growth histories of the long-lived ocean quahog (<Emphasis Type="Italic">Arctica islandica</Emphasis>)

To investigate environmental variability during the late Holocene in the western Gulf of Maine, USA, we collected a 142-year-old living bivalve (Arctica islandica) in 2004, and three fossil A. islandica shells of the Medieval Warm Period (MWP) and late MWP / Little Ice Age (LIA) period (corrected ¹⁴C_AMS = 1030 ± 78 ad; 1320 ± 45 ad; 1357 ± 40 ad) in 1996. We compared the growth record of the modern shell with continuous plankton recorder (CPR) time-series (1961–2003) from the Gulf of Maine. A significant correlation (r ² = 0.55; p < 0.0001) exists between the standardized annual growth index (SGI) of the modern shell and the relative abundance of zooplankton species Calanus finmarchicus. We therefore propose that SGI data from A. islandica is a valid proxy for paleo-productivity of at least one major zooplankton taxa. SGIs from these shells reveal significant periods of 2–6 years (NAO-like) based on wavelet analysis, multitaper method (MTM) analysis and singular spectrum analysis (SSA) during the late Holocene. Based on established physical oceanographic observation in the Gulf of Maine, we suggest that slope water variability coupled with North Atlantic Oscillation (NAO) dynamics is primarily responsible for the observed SGI variability. Special Issue: AGU OS06 special issue “Ocean’s role in climate change—a paleo perspective”.     

Weathering and the mobility of phosphorus in the catchments and forefields of the Rhône and Oberaar glaciers, central Switzerland: Implications for the global phosphorus cycle on glacial-interglacial timescales

In this study we evaluate the dynamics of the biophile element phosphorus (P) in the catchment and proglacial areas of the Rhône and Oberaar glaciers (central Switzerland). We analysed erosion and dissolution rates of P-containing minerals in the subglacial environment by sampling water and suspended sediment in glacier outlets during three ablation and two accumulation seasons. We also quantified biogeochemical weathering rates of detrital P in proglacial sedimentary deposits using two chronosequences of samples of fresh, suspended, material obtained from the Oberaar and Rhône water outlets, Little-Ice-Age (LIA) moraines and Younger Dryas (YD) tills in each catchment. Subglacial P weathering is mainly a physical process and detrital P represents more than 99% of the precipitation-corrected total P denudation flux (234 and 540 kg km⁻² yr⁻¹ for the Rhône and Oberaar catchments, respectively). The calculated detrital P flux rates are three to almost five times higher than the world average flux. The precipitation-corrected soluble reactive P (SRP) flux corresponds to 1.88-1.99 kg km⁻² yr⁻¹ (Rhône) and 2.12-2.44 kg km⁻² yr⁻¹ (Oberaar), respectively. These fluxes are comparable to those of tropical rivers draining transport-limited, tectonically inactive weathering areas.In order to evaluate the efficiency of detrital P weathering in the Rhône and Oberaar proglacial areas, we systematically graded apatite grains extracted from the chronosequence in each catchment relative to weathering-induced changes in their surface morphologies (grades 1-4). Fresh apatite grains are heavily indented and dissolution rounded (grade 1). LIA grains from two 0-10 cm deep moraine samples show extensive dissolution etching, similar to surface grains from the YD profile (mean grades 2.7, 3.5 and 3.5, respectively). In these proglacial deposits, the weathering front deepens progressively as a function of time due to biocorrosion in the evolving acidic pedosphere, with mechanical indentations on grains acting as sites of preferential dissolution. We also measured iron-bound, organic and detrital P concentrations in the chronosequence and show that organic and iron-bound P has almost completely replaced detrital P in the top layers of the YD profiles. Detrital P weathering rates are calculated as 310 and 280 kg km⁻² yr⁻¹ for LIA moraines and 10 kg km⁻² yr⁻¹ for YD tills. During the first 300 years of glacial sediment exposure P dissolution rates are shown to be approximately 70 times higher than the mean global dissolved P flux from ice-free continents. After 11.6 kyr the flux is 2.5 times the global mean. These data strengthen the argument for substantial changes in the global dissolved P flux on glacial-interglacial timescales. A crude extrapolation from the data described here suggests that the global dissolved P flux may increase by 40-45% during the first few hundred years of a deglaciation phase.     

Fluvial channels on Titan: Initial Cassini RADAR observations

Cassini radar images show a variety of fluvial channels on Titan's surface, often several hundreds of kilometers in length. Some (predominantly at low- and mid-latitude) are radar-bright and braided, resembling desert washes where fines have been removed by energetic surface liquid flow, presumably from methane rainstorms. Others (predominantly at high latitudes) are radar-dark and meandering and drain into or connect polar lakes, suggesting slower-moving flow depositing fine-grained sediments. A third type, seen predominantly at mid- and high latitudes, have radar brightness patterns indicating topographic incision, with valley widths of up to 3 km across and depth of several hundred meters. These observations show that fluvial activity occurs at least occasionally at all latitudes, not only at the Huygens landing site, and can produce channels much larger in scale than those observed there. The areas in which channels are prominent so far amount to about 1% of Titan's surface, of which only a fraction is actually occupied by channels. The corresponding global sediment volume inferred is not enough to account for the extensive sand seas. Channels observed so far have a consistent large-scale flow pattern, tending to flow polewards and eastwards.     

Burning molten metallic spheres: One class of ball lightning?

Abrahamson and Dinniss [2000. Ball lightning caused by oxidation of nanoparticle networks from normal lightning strikes on soil. Nature 403, 519–521] proposed a theory of ball lighting in which silicon nanoparticles undergo slow oxidation and emit light. Paiva et al. [2007. Production of ball-lightning-like luminous balls by electrical discharges in silicon. Physical Review Letters 98, 048501] reported that an electric arc to silicon produced long-lasting luminous white spheres showing many characteristics of ball lightning. We show experimentally that these consist of burning molten silicon spheres with diameters in the 0.1–1 mm range. The evidence of our experiments leads us to propose that a subset of ball lightning events may consist of macro-scale molten spheres of burning metallic materials likely to be ejected from a conventional lightning strike to earth.     

Seismic scale saturation relations in turbidite reservoirs undergoing waterflood

Estimates of the effective fluid modulus from seismic cannot be directly converted to the true pore-volume weighted mean saturation     determined from fluid flow principles by using the saturation laws currently in use. One of the reasons is that seismic waves sample the reservoir geology and production induced saturation heterogeneity in a different way from the fluids. This mismatch prevents accurate quantitative evaluation of saturation changes from 4D seismic analysis. To tackle this problem, a reservoir-related saturation law is developed for a turbidite reservoir – this geology being chosen because the architecture for a single sand package can be modelled as a stack of horizontal beds. An effective medium and perturbation theory are applied to the determination of the seismic properties of this model. This calculation provides a relationship that connects the true saturation     to the effective fluid modulus from seismic via statistical measures of the vertical spread of the porosity and saturation variations in the reservoir. These statistics can be extracted from the simulation model and if known, enable the new saturation law to deliver a significant improvement in accuracy when estimating     compared to other well-known laws. The relationship that has been developed also captures the effect of inter-bedded shales and can therefore be used to estimate true saturation in regions of the reservoir with moderate to low net-to-gross, provided the fraction of the shale component is known. In practice, the final choice of saturation law depends upon the reservoir information available, the assumptions that can be tolerated and the accuracy required in any particular reservoir characterization study.