Predicting Soil-Water and Soil-Air Transport Properties and Their Effects on Soil-Vapor Extraction Efficiency

Accurate prediction of water and air Iran sport parameters in variably saturated soil is necessary for modeling of soil-vapor extraction (SVE) at soil sites contaminated with volatile organic chemicals (VOCs). An expression for predicting saturated water permeability (k_l,s) in undisturbed soils from the soil total porosity and the field capacity soil-water content was developed by fitting a tortuous-tube fluid flow model to measured water permeability and gas diffusivity data. The new k_l,s expression gave accurate predictions when tested against independent k_l,s data. The k_l,s expression was implemented in the Campbell relative water permeability model to yield a predictive model for water permeability in variably saturated, undisturbed soil. The water permeability model, together with recently developed predictive equations for gas permeability and gas diffusivity, was used in a two-dimensional numerical SVE model that also included non-equilibrium mass transfer of VOC from a separate phase (nonaqueous phase liquid [NAPL]) to the air phase. SVE: calculations showed that gas permeability is likely the most important factor controlling VOC migration and vapor extraction efficiency. Water permeability and gas diffusivity effects became significant at water contents near and above field capacity. The NAPL-air mass transfer coefficient also had large impacts on simulated vapor extraction efficiency. The calculations suggest that realistic SVE models need to include predictive expressions for both conveciive, diffusive. and phase-partitioning processes in natural, undisturbed soils.     

Controls on the lateral channel-migration rate of braided channel systems in coarse non-cohesive sediment

Lateral movements of alluvial river channels control the extent and reworking rates of alluvial fans, floodplains, deltas, and alluvial sections of bedrock rivers. These lateral movements can occur by gradual channel migration or by sudden changes in channel position (avulsions). Whereas models exist for rates of river avulsion, we lack a detailed understanding of the rates of lateral channel migration on the scale of a channel belt. In a two-step process, we develop here an expression for the lateral migration rate of braided channel systems in coarse, non-cohesive sediment. On the basis of photographic and topographic data from laboratory experiments of braided channels performed under constant external boundary conditions, we first explore the impact of autogenic variations of the channel-system geometry (i.e. channel-bank heights, water depths, channel-system width, and channel slope) on channel-migration rates. In agreement with theoretical expectations, we find that, under such constant boundary conditions, the laterally reworked volume of sediment is constant and lateral channel-migration rates scale inversely with the channel-bank height. Furthermore, when channel-bank heights are accounted for, lateral migration rates are independent of the remaining channel geometry parameters. These constraints allow us, in a second step, to derive two alternative expressions for lateral channel-migration rates under different boundary conditions using dimensional analysis. Fits of a compilation of laboratory experiments to these expressions suggest that, for a given channel bank-height, migration rates are strongly sensitive to water discharges and more weakly sensitive to sediment discharges. In addition, external perturbations, such as changes in sediment and water discharges or base level fall, can indirectly affect lateral channel-migration rates by modulating channel-bank heights. © 2019 The Author. Earth Surface Processes and Landforms published by John Wiley & Sons, Ltd. © 2019 The Author. Earth Surface Processes and Landforms published by John Wiley & Sons, Ltd.     

Acoustic and petrophysical properties of mechanically compacted overconsolidated sands: Part 2 – Rock physics modelling and applications

Part one of this paper reported results from experimental compaction measurements of unconsolidated natural sand samples with different mineralogical compositions and textures. The experimental setup was designed with several cycles of stress loading and unloading applied to the samples. The setup was aimed to simulate a stress condition where sediments underwent episodes of compaction, uplift and erosion. P-wave and S-wave velocities and corresponding petrophysical (porosity and density) properties were reported. In this second part of the paper, rock physics modelling utilizing existing rock physics models to evaluate the model validity for measured data from part one were presented. The results show that a friable sand model, which was established for normally compacted sediments is also capable of describing overconsolidated sediments. The velocity–porosity data plotted along the friable sand lines not only describe sorting deterioration, as has been traditionally explained by other studies, but also variations in pre-consolidation stress or degree of stress release. The deviation of the overconsolidated sands away from the normal compaction trend on the V_P/V_S and acoustic impedance space shows that various stress paths can be predicted on this domain when utilizing rock physics templates. Fluid saturation sensitivity is found to be lower in overconsolidated sands compared to normally consolidated sands. The sensitivity decreases with increasing pre-consolidation stress. This means detectability for four-dimensional fluid saturation changes can be affected if sediments were pre-stressed and unloaded. Well log data from the Barents Sea show similar patterns to the experimental sand data. The findings allow the development of better rock physics diagnostics of unloaded sediments, and the understanding of expected 4D seismic response during time-lapse seismic monitoring of uplifted basins. The studied outcomes also reveal an insight into the friable sand model that its diagnostic value is not only for describing sorting microtextures, but also pre-consolidation stress history. The outcome extends the model application for pre-consolidation stress estimation, for any unconsolidated sands experiencing similar unloading stress conditions to this study.     

Correction of Discretization Errors Simulated at Supply Wells

Many hydrogeology problems require predictions of hydraulic heads in a supply well. In most cases, the regional hydraulic response to groundwater withdrawal is best approximated using a numerical model; however, simulated hydraulic heads at supply wells are subject to errors associated with model discretization and well loss. An approach for correcting the simulated head at a pumping node is described here. The approach corrects for errors associated with model discretization and can incorporate the user's knowledge of well loss. The approach is model independent, can be applied to finite difference or finite element models, and allows the numerical model to remain somewhat coarsely discretized and therefore numerically efficient. Because the correction is implemented external to the numerical model, one important benefit of this approach is that a response matrix, reduced model approach can be supported even when nonlinear well loss is considered.     

Can GPS-Derived Surface Loading Bridge a GRACE Mission Gap?

We investigated two ‘gap-filler’ methods based on GPS-derived low-degree surface loading variations (GPS-I and GPS-C) and a more simple method (REF-S) which extends a seasonal harmonic variation into the expected Gravity Recovery and Climate Experiment (GRACE) mission gap. We simulated two mission gaps in a reference solution (REF), which is derived from a joint inversion of GRACE (RL05) data, GPS-derived surface loading and simulated ocean bottom pressure. The GPS-I and GPS-C methods both have a new type of constraint applied to mitigate the lack of GPS station network coverage over the ocean. To obtain the GPS-C solution, the GPS-I method is adjusted such that it fits the reference solution better in a 1.5 year overlapping period outside of the gap. As can be expected, the GPS-I and GPS-C solutions contain larger errors compared to the reference solution, which is heavily constrained by GRACE. Within the simulated gaps, the GPS-C solution generally fits the reference solution better compared to the GPS-I method, both in terms of spherical harmonic loading coefficients and in terms of selected basin-averaged hydrological mass variations. Depending on the basin, the RMS-error of the water storage variations (scaled for leakage effects) ranges between 1.6 cm (Yukon) and 15.3 cm (Orinoco). In terms of noise level, the seasonal gap-filler method (REF-S) even outperforms the GPS-I and GPS-C methods, which are still affected by spatial aliasing problems. However, it must be noted that the REF-S method cannot be used beyond the study of simple harmonic seasonal variations.     

GEOCHEMICAL RECOGNITION OF SPILLED SEDIMENTS USED IN NUMERICAL MODEL VALIDATION

A fixed link (tunnel and bridge, in total 16 km) was constructed between Sweden and Denmark during 1995-2000. As part of the work, approximately 16 million tonnes of seabed materials (limestone and clay till) were dredged, and about 0.6 million tonnes of these were spilled in the water. Modelling of the spreading and sedimentation of the spilled sediments took place as part of the environmental monitoring of the construction activities. In order to verify the results of the numerical modelling of sediment spreading and sedimentation, a new method with the purpose of distinguishing between the spilled sediments and the naturally occurring sediments was developed. Because the spilled sediments tend to accumulate at the seabed in areas with natural sediments of the same size, it is difficult to separate these based purely on the physical properties. The new method is based on the geo-chemical differences between the natural sediment in the area and the spill. The basic propertiesused are the higher content of calcium carbonate material in the spill as compared to the natural sediments and the higher Ca/Sr ratio in the spill compared to shell fragments dominating the natural calcium carbonate deposition in the area. The reason for these differences is that carbonate derived from recent shell debris can be discriminated from Danien limestone, which is the material in which the majority of the dredging took place, on the basis of the Ca/Sr ratio being 488 in Danien Limestone and 237 in shell debris. The geochemical recognition of the origin of the sediments proved useful in separating the spilled from the naturally occurring sediments. Without this separation, validation of the modelling of accumulation of spilled sediments would not have been possible. The method has general validity and can be used in many situations where the origin of a given sediment is sought.     

Impact angle influence in high velocity dust collisions during planetesimal formation

We have examined the influence of impact angle in collisions between small dust aggregates and larger dust targets through laboratory experiments. Targets consisted of μm-sized quartz dust and had a porosity of about 67%; the projectiles, between 1 and 5 mm in diameter, were slightly more compact (64% porosity). The collision velocity was centered at 20 m/s and impact angles range from 0° to 45°. At a given impact angle, the target gained mass for projectiles smaller than a threshold size, which decreases with increasing angle from about 3 mm to 1 mm. The fact that growth is possible up to the largest angles studied supports the idea of planetesimal formation by sweep-up of small dust aggregates.     

Assembly and concept of a web-based GIS within the paleolimnological project CONTINENT (Lake Baikal,Russia)

Web-based Geographical Information Systems (GIS) are excellent tools within interdisciplinary and multi-national geoscience projects to exchange and visualize project data. The web-based GIS presented in this paper was designed for the paleolimnological project ‘High-resolution CONTINENTal paleoclimate record in Lake Baikal’ (CONTINENT) (Lake Baikal, Siberia, Russia) to allow the interactive handling of spatial data. The GIS database combines project data (core positions, sample positions, thematic maps) with auxiliary spatial data sets that were downloaded from freely available data sources on the world wide web. The reliability of the external data was evaluated and suitable new spatial datasets were processed according to the scientific questions of the project. GIS analysis of the data was used to assist studies on sediment provenance in Lake Baikal, or to help answer questions such as whether the visualization of present-day vegetation distribution and pollen distribution supports the conclusions derived from palynological analyses. The refined geodata are returned back to the scientific community by using online data publication portals. Data were made citeable by assigning persistent identifiers (DOI) and were published through the German National Library for Science and Technology (TIB Hannover, Hannover, Germany).     

Mudstone compaction curves in basin modelling: a study of Mesozoic and Cenozoic Sediments in the northern North Sea

Basin modelling studies are carried out in order to understand the basin evolution and palaeotemperature history of sedimentary basins. The results of basin modelling are sensitive to changes in the physical properties of the rocks in the sedimentary sequences. The rate of basin subsidence depends, to a large extent, on the density of the sedimentary column, which is largely dependent on the porosity and therefore on the rate of compaction. This study has tested the sensitivity of varying porosity/depth curves and related thermal conductivities for the Cenozoic succession along a cross‐section in the northern North Sea basin, offshore Norway. End‐member porosity/depth curves, assuming clay with smectite and kaolinite properties, are compared with a standard compaction curve for shale normally applied to the North Sea. Using these alternate relationships, basin geometries of the Cenozoic succession may vary up to 15% from those predicted using the standard compaction curve. Isostatic subsidence along the cross‐section varies 2.3–4.6% between the two end‐member cases. This leads to a 3–8% difference in tectonic subsidence, with maximum values in the basin centre. Owing to this, the estimated stretching factors vary up to 7.8%, which further gives rise to a maximum difference in heat flow of more than 8.5% in the basin centre. The modelled temperatures for an Upper Jurassic source rock show a deviation of more than 20 °C at present dependent on the thermal conductivity properties in the post‐rift succession. This will influence the modelled hydrocarbon generation history of the basin, which is an essential output from basin modelling analysis. Results from the northern North Sea have shown that varying compaction trends in sediments with varying thermal properties are important parameters to constrain when analysing sedimentary basins.