Impact of normal and shear stresses due to wheel slip on hydrological properties of an agricultural clay loam: Experimental and new computerized approach

The main purpose of this study was to evaluate the effect that mechanical stresses acting under the slipping driving wheels of agricultural equipment have on the soil’s pore system and water flow process (surface runoff generation during extreme event). The field experiment simulated low slip (1%) and high slip (27%) on a clay loam. The stress on the soil surface and changes in the amounts of water flowing from macropores were simulated using the Tires/tracks And Soil Compaction (TASC) tool and the MACRO model, respectively. Taking a 65 kW tractor on a clay loam as a reference, results showed that an increase in slip of the rear wheels from 1% to 27% caused normal stress to increase from 90.6 kPa to 104.4 kPa at the topsoil level, and the maximum shear contact stress to rise drastically from 6.0 kPa to 61.6 kPa. At 27% slip, topsoil was sheared and displaced over a distance of 0.35 m. Excessive normal and shear stress values with high slip caused severe reductions of the soil’s macroporosity, saturated hydraulic conductivity, and water quantities flowing from topsoil macropores. Assuming that, under conditions of intense rainfall on sloping land, a loss in vertical water flow would mean an increase in surface runoff, we calculated that a rainfall intensity of 100 mm h-1 and a rainfall duration of 1 h would increase the runoff coefficient to 0.79 at low slip and to 1.00 at high slip, indicating that 100% of rainwater would be transformed into surface runoff at high slip. We expect that these effects have a significant impact on soil erosion and floods in steeper terrain (slope > 15°) and across larger surface areas (> 16 m2) than those included in our study.

Relating in situ gas measurements to the surface outgassing properties of cometary nuclei

The sensitivity of the gas flow field to changes in different initial conditions has been studied for the case of a highly simplified cometary nucleus model. The nucleus model simulated a homogeneously outgassing sphere with a more active ring around an axis of symmetry. The varied initial conditions were the number density of the homogeneous region, the surface temperature, and the composition of the flow (varying amounts of H2O and CO2) from the active ring. The sensitivity analysis was performed using the Polynomial Chaos Expansion (PCE) method. Direct Simulation Monte Carlo (DSMC) was used for the flow, thereby allowing strong deviations from local thermal equilibrium. The PCE approach can be used to produce a sensitivity analysis with only four runs per modified input parameter and allows one to study and quantify non-linear responses of measurable parameters to linear changes in the input over a wide range. Hence the PCE allows one to obtain a functional relationship between the flow field properties at every point in the inner coma and the input conditions. It is for example shown that the velocity and the temperature of the background gas are not simply linear functions of the initial number density at the source. As probably expected, the main influence on the resulting flow field parameter is the corresponding initial parameter (i.e. the initial number density determines the background number density, the temperature of the surface determines the flow field temperature, etc.). However, the velocity of the flow field is also influenced by the surface temperature while the number density is not sensitive to the surface temperature at all in our model set-up. Another example is the change in the composition of the flow over the active area. Such changes can be seen in the velocity but again not in the number density. Although this study uses only a simple test case, we suggest that the approach, when applied to a real case in 3D, should assist in identifying the sensitivity of gas parameters measured in situ by, for example, the Rosetta spacecraft to the surface boundary conditions and vice versa.

Physical properties of 39 sediment cores from the Bengal Fan and northeast Indian Ocean

We obtained sediment physical properties and geochemical data from 47 piston and gravity cores located in the Bay of Bengal, to study the complex history of the Late Pleistocene run-off from the Ganges and Brahmaputra rivers and its imprint on the Bengal Fan. Grain-size parameters were predicted from core logs of density and velocity to infer sediment transport energy and to distinguish different environments along the 3000-km-long transport path from the delta platform to the lower fan. On the shelf, 27 cores indicate rapidly prograding delta foresets today that contain primarily mud, whereas outer shelf sediment has 25% higher silt contents, indicative of stronger and more stable transport regime, which prevent deposition and expose a Late Pleistocene relic surface. Deposition is currently directed towards the shelf canyon 'Swatch of No Ground', where turbidites are released to the only channel-levee system that is active on the fan during the Holocene. Active growth of the channel-levee system occurred throughout sea-level rise and highstand with a distinct growth phase at the end of the Younger Dryas. Coarse-grained material bypasses the upper fan and upper parts of the middle fan, where particle flow is enhanced as a result of flow-restriction in well-defined channels. Sandier material is deposited mainly as sheet-flow deposits on turbidite-dominated plains at the lower fan. The currently most active part of the fan with 10-40 cm thick turbidites is documented for the central channel including inner levees (e.g., site 40). Site 47 from the lower fan far to the east of the active channel-levee system indicates the end of turbidite sedimentation at 300 ka for that location. That time corresponds to the sea-level lowering during late isotopic stage 9 when sediment supply to the fan increased and led to channel avulsion farther upstream, probably indicating a close relation of climate variability and fan activity. Pelagic deep-sea sites 22 and 28 contain a 630-kyear record of climate response to orbital forcing with dominant 21- and 41-kyear cycles for carbonate and magnetic susceptibility, respectively, pointing to teleconnections of low-latitude monsoonal forcing on the precession band to high-latitude obliquity forcing. Upper slope sites 115, 124, and 126 contain a record of the response to high-frequency climate change in the Dansgaard-Oeschger bands during the last glacial cycle with shared frequencies between 0.75 and 2.5 kyear. Correlation of highs in Bengal Fan physical properties to lows in the d18O record of the GISP2 ice-core suggests that times of greater sediment transport energy in the Bay of Bengal are associated with cooler air temperatures over Greenland. Teleconnections were probably established through moisture and other greenhouse-gas forcing that could have been initiated by instabilities in the methane hydrate reservoir in the oceans.

Physical properties and accumulation rates in Dronning Maud Land, Antarctica

The European Programme for Ice Coring in Antarctica includes a comprehensive pre-site survey on the inland ice plateau of Dronning Maud Land, Antarctica. The German glaciological programme during the 1997/98 field season was carried out along a 1200 km traverse on Amundsenisen and involved sampling the snow cover in pits and by shallow firn cores. This paper focuses on the accumulation studies. The cores were dated by dielectric-profiling and continuous-flow analysis. Distinct volcanogenic peaks and seasonal signals in the profiles served to establish a depth time-scale. The eruptions of Krakatoa, Tambora, an unknown volcano, Kuwae and El Chichon are well-documented in the ice. Variations of the accumulation rates over different times were inferred from the depth time-scales. A composite record of accumulation rates for the last 200 years was produced by stacking 12 annually resolved records. According to this, accumulation rates decreased in the 19th century and increased in the 20th century. The recent values are by no means extraordinary, as they do not exceed the values at the beginning of the 19th century. Variations in accumulation rates are most probably linked to temperature variations indicated in d18O records from Amundsenisen.

Summary of physical properties on cores used for petrofabric analyses triaxial experiments at DSDP Leg 67 Holes

Deformation features within the cores are studied with a view towards elucidating the structure of the Middle America Trench along the transect drilled during Leg 67. Where possible, inferences are made as to the physical environment of deformation. Extensional tectonics prevails in the area of the seaward slope and trench. Fracturing and one well-preserved normal fault are found mostly within the lower Miocene chalks, at the base of the sedimentary section. These chalks have high porosities (40%-60%) and water content (30%-190%, based on % dry wt.). Experimental triaxial compression tests conducted on both dry and water-saturated samples of chalk from Holes 495 and 499B show that only in the saturated samples is more brittle behavior observed. Brittle failure of the chalks is greatly facilitated by pore fluid pressures that lead to low effective pressures. Additional embrittlement (weakening) can take place as a result of the imposed extensional stress resulting from bending of a subducting elastic oceanic plate. The chalks exhibit, in a landward direction, an increase in density and mechanical strength and a decrease in water content. These changes are attributed to mechanical compaction that may have resulted from tectonic horizontal compression. The structure of the landward slope is not well understood because the slope sites had to be abandoned due to the presence of gas hydrate. The relationship of the chaotic, brittle deformation (observed in the cores from Hole 494A) at the base of the landward slope to tectonic processes remains unclear. The deformation observed on the slope sites (Holes 496 and 497) is mostly fracturing and near-vertical sigmoidal veinlets. These are interpreted as being the result of gas/fluid overpressurization due to the decomposition of the gas hydrate, and not due to tectonic loading of accreted sediments. Aside from four small displacement (less than 1cm) reverse faults observed in the lower Miocene chalks (which may be the product of soft-sediment deformation), there is a noticeable absence of structures reflecting a dominance of horizontal (tectonic) compression along the transect drilled. The absence of such features, the lack of continuity of sediment types across the trench-landward slope, and the normal stratigraphic sequence in Hole 494A do not support any known accretionary model.

Heat flow values in the Sea of Okhotsk region

Tectonic structure and anomalous distributions of geophysical fields of the Sea of Okhotsk region are considered; the lack of reliable data on age of the lithosphere beneath basins of various origin in the Sea of Okhotsk is noted. Model calculations based on geological and geophysical data yielded 65 Ma (Cretaceous-Paleocene boundary) age for the Central Okhotsk rise underlain by the continental lithosphere. This estimate agrees with the age (the end of Cretaceous) derived from seismostratigraphic data. A comparative analysis of theoretical and measured heat flows in the Akademii Nauk Rise, underlain by the thinned continental crust, is performed. The analysis points to a higher (by 20%) value of the measured thermal background of the rise, which is consistent with high negative gradient of gravity anomalies in this area. Calculations yielded 36 Ma (Early Oligocene) age and lithosphere thickness of 50 km for the South Okhotsk depression, whose seafloor was formed by processes of back-arc spreading. The estimated age of the depression is supported by kinematic data on the region; the calculated thickness of the lithosphere coincides with the value estimated from data of magnetotelluric sounding here. This indicates that formation time (36 Ma) of the South Okhotsk depression was estimated correctly. Numerical modeling performed for determination of the basement age of rifting basins in the Sea of Okhotsk gave the following estimates: 18 Ma (Early Miocene) for the Deryugin Basin, 12 Ma (Middle Miocene) for the TINRO Basin, and 23 Ma (Late Oligocene) for the West Kamchatka Trough. These estimates agree with formation time (Oligocene-Quaternary) of the sedimentary cover in rifting basins of the Sea of Okhotsk derived from geological and geophysical data. Model temperature estimates are obtained for lithologic and stratigraphic boundaries of the sedimentary cover in the Deryugin and TINRO Basins and the West Kamchatka Trough; the temperature analysis indicates that the latter two structures are promising for oil and hydrocarbon gas generation; the West Kamchatka Trough possesses better reservoir properties compared to the TINRO and Deryugin Basins. The latter is promising for generation of hydrocarbon gas. Paleogeodynamic reconstructions of the Sea of Okhotsk region evolution are obtained for times of 90, 66, and 36 Ma on the base of kinematic, geomagnetic, structural, tectonic, geothermal, and other geological and geophysical data.

(Table 1) Physical properties of basalts from DSDP Legs 69 and 70

Physical properties of basalts from Deep Sea Drilling Project Sites 504 and 505, south of the Costa Rica Rift, including wet-bulk density, water content, sonic velocity, and thermal conductivity, were measured on board D/V Glomar Challenger during Legs 69 and 70. The mean wet-bulk density is 2.90±0.06 g/cm**3, porosity 5.0±2.2%, sonic velocity 5.75±0.30 km/s, and thermal conductivity 1.67±0.09 W/m°K. Basalts from this young ocean crust (5.9 m.y.) have relatively low porosity and consequently high density and sonic velocity, as compared to average DSDP basalts. Some systematic trends in depth variation of physical properties were found: down to Core 40 in Hole 504B, grain densities were lower than those deeper in the hole, whereas porosity in the upper section was higher. This can be attributed either to differences in the flow type or in the nature of alteration of basalts at the different levels in the hole.

Consolidation properties and hydraulic conductivities of sediments from ODP Leg 207 sites

The stress history, permeability, and compressibility of sediments from Demerara Rise recovered during Ocean Drilling Program Leg 207 were determined using one-dimensional incremental load consolidation and low-gradient flow pump permeability tests. Relationships among void ratio, effective stress, and hydraulic conductivity are presented for sampled lithologic units and used to reconstruct effective stress, permeability, and in situ void ratio profiles for a transect of three sites across Demerara Rise. Results confirm that a significant erosional event occurred on the northeastern flank of the rise during the late Miocene, resulting in the removal of ~220 m of upper Oligocene-Miocene deposits. Although Neogene and Paleogene sediments tend to be overconsolidated, Cretaceous sediments are normally consolidated to underconsolidated, suggesting the presence of overpressure. A pronounced drop in permeability occurs at the transition from the Cretaceous black shales into the overlying Maastrichtian-upper Paleocene chalks and clays. The development of a hydraulic seal at this boundary may be responsible for overpressure in the Cretaceous deposits, leading to the lower overconsolidation ratios of these sediments. Coupled with large regional variations in sediment thickness (overburden stresses), the higher permeability overpressured Cretaceous sediments represent a regional lateral fluid conduit on Demerara Rise, possibly venting methane-rich fluids where it outcrops on the margin's northeastern flank.

Planimetric ice-flow convergence and flow-orthonormal strain rate across the Antarctic Ice Sheet, with links to model result files

Fast-flowing ice streams discharge most of the ice from the interior of the Antarctic Ice Sheet coastward. Understanding how their tributary organisation is governed and evolves is essential for developing reliable models of the ice sheet's response to climate change. Despite much research on ice-stream mechanics, this problem is unsolved, because the complexity of flow within and across the tributary networks has hardly been interrogated. Here I present the first map of planimetric flow convergence across the ice sheet, calculated from satellite measurements of ice surface velocity, and use it to explore this complexity. The convergence map of Antarctica elucidates how ice-stream tributaries draw ice from the interior. It also reveals curvilinear zones of convergence along lateral shear margins of streaming, and abundant convergence ripples associated with nonlinear ice rheology and changes in bed topography and friction. Flow convergence on ice-stream tributaries and their feeding zones is markedly uneven, and interspersed with divergence at distances of the order of kilometres. For individual drainage basins as well as the ice sheet as a whole, the range of convergence and divergence decreases systematically with flow speed, implying that fast flow cannot converge or diverge as much as slow flow. I therefore deduce that flow in ice-stream networks is subject to mechanical regulation that limits flow-orthonormal strain rates. These properties and the gridded data of convergence and flow-orthonormal strain rate in this archive provide targets for ice- sheet simulations and motivate more research into the origin and dynamics of tributarization.

Paleomagnetic properties of ODP Hole 178-1101A

Sediments recovered from a drift deposit located on the Pacific side of the Antarctic Peninsula (ODP Leg 178, Site 1101) give a physical record of a bottom current, sourced from the Weddell Sea Deep Water, for the past 3 Ma. Sediment grain size and magnetic fabric analyses indicate a contourite depositional environment and little change in the average intensity of this current. Terrigenous fluxes decreased around the time of the onset of Northern Hemisphere Glaciation, which we interpret as a freezing of the base of the Antarctic Peninsula Ice Cap. Terrigenous fluxes have increased since 1.7 Ma implying a possible return of the Antarctic Peninsula Ice Cap to a more wet-based ice sheet.

Magnetic properties of light and dark sediment layers from the Japan Sea

Rock magnetic/paleoclimatic/diagenetic relationships of sediments spanning the last 0.78 Ma have been investigated using samples collected from light and dark layers recovered at ODP Sites 794 (Yamato Basin) and 795 (Japan Basin). Rock-magnetic parameters (K, Kfd, ARM, SIRM, S-ratio) are shown to reflect diagenetic processes and climate-related variations in the concentration, mineralogy and grain-size of the magnetic minerals contained within the sediments. The magnetic mineralogy is dominated by ferrimagnetic (magnetite-type) minerals with a small contribution made by hematite and iron sulphides such as pyrrhotite and/or greigite. Magnetic mineral concentration and grain size vary between light and dark layers with the former characterized by a higher magnetic content and a finer magnetic grain size. Magnetite dissolution, related to sulfate reduction due to bacterial degradation of organic matter, is the process responsible for the magnetic characteristics observed in the dark layers, testifying to the reducing conditions in the basin. Variations in the rock magnetic properties of the sediments are strongly correlated with global oxygen isotope fluctuations, with glacial stages characterized by a lower magnetic mineral content and a coarser magnetic grain size relative to interglacial stages. Major downcore changes in the magnetic properties observed at Site 794 can be related to changes in the oceanographic conditions of the basin associated with the flow of the warm Tsushima Current into the Japan Sea at about 0.35-0.40 Ma ago.

Bottom water properties observed during Nathanial B. Palmer cruises NBP00-08 and NBP04-08 to the Southern Ocean

We report on observations of dense shelf water overflows and Antarctic Bottom Water (AABW) formation along the continental margin of the Adelie and George V Land coast between 140°E and 149°E. Vertical sections and bottom layer water mass properties sampled during two RVIB Nathaniel B Palmer hydrographic surveys (NBP00-08, December 2000/January 2001 and NBP04-08, October 2004) describe the spreading of cold, dense shelf water on the continental slope and rise from two independent source regions. The primary source region is the Adelie Depression, exporting high-salinity dense shelf water through the Adelie Sill at 143°E. An additional eastern source region of lower-salinity dense shelf water from the Mertz Depression is identified for the first time from bottom layer properties northwest of the Mertz Sill and Mertz Bank (146°E-148°E) that extend as far as the Buffon Channel (144.75°E) in summer. Regional analysis of satellite-derived ice production estimates over the entire region from 1992 to 2005 suggests that up to 40% of the total ice production for the region occurs over the Mertz Depression and therefore this area is likely to make a significant contribution to the total dense shelf water export. Concurrent time series from bottom-mounted Microcats and ADCP instruments from the Mertz Polynya Experiment (April 1998 to May 1999) near the Adelie Sill and on the upper continental slope (1150 m) and lower continental rise (3250 m) to the north describe the seasonal variability in downslope events and their interaction with the ambient water masses. The critical density for shelf water to produce AABW is examined and found to be 27.85 kg/m**3 from the Adelie Depression and as low as 27.80 kg/m**3 from the Mertz Depression. This study suggests previous dense shelf water export estimates based on the flow through the Adelie Sill alone are conservative and that other regions around East Antarctica with similar ice production to the Mertz Depression could be contributing to the total AABW in the Australian-Antarctic Basin.

Debris flow magnitudes: a global catalog

The episodic occurrence of debris flow events in response to stochastic precipitation and wildfire events makes hazard prediction challenging. Previous work has shown that frequency-magnitude distributions of non-fire-related debris flows follow a power law, but less is known about the distribution of post-fire debris flows. As a first step in parameterizing hazard models, we use frequency-magnitude distributions and cumulative distribution functions to compare volumes of post-fire debris flows to non-fire-related debris flows. Due to the large number of events required to parameterize frequency-magnitude distributions, and the relatively small number of post-fire event magnitudes recorded in the literature, we collected data on 73 recent post-fire events in the field. The resulting catalog of 988 debris flow events is presented as an appendix to this article. We found that the empirical cumulative distribution function of post-fire debris flow volumes is composed of smaller events than that of non-fire-related debris flows. In addition, the slope of the frequency-magnitude distribution of post-fire debris flows is steeper than that of non-fire-related debris flows, evidence that differences in the post-fire environment tend to produce a higher proportion of small events. We propose two possible explanations: 1) post-fire events occur on shorter return intervals than debris flows in similar basins that do not experience fire, causing their distribution to shift toward smaller events due to limitations in sediment supply, or 2) fire causes changes in resisting and driving forces on a package of sediment, such that a smaller perturbation of the system is required in order for a debris flow to occur, resulting in smaller event volumes.