Calibration of cyclic constitutive models by oscillating functions

Most structures are subjected to more cyclic loads during their life time than static loads. These cyclic action could be a result of either natural or man-made activities and may lead to soil failure. In order to understand the response of the foundation and its interaction with these complex cyclic loadings, various researchers have over the years developed different constitutive models. Although a lot of research is being carried out on these relatively new models, little or no details exist in literature about the model-based identification of the cyclic constitutive parameters which to a large extent govern the quality of the model output. This could be attributed to the difficulties and complexities of the inverse modeling of such complex phenomena. A variety of optimisation strategies are available for the solution of the sum of least-squares problems as usually done in the field of model calibration. However, for the back analysis (calibration) of the soil response to oscillatory load functions, this article gives insight into the model calibration challenges and also puts forward a method for the inverse modeling of cyclic loaded foundation response such that high-quality solutions are obtained with minimum computational effort.     

Remediation scenarios for selenium contamination, Blackfoot watershed, southeast Idaho, USA

Extensive phosphate mining in the Blackfoot watershed of Idaho (USA) has substantially increased the selenium (Se) concentration in the river during both snowmelt and baseflow when groundwater discharge dominates. Phosphate mines create a linkage between Se-laden shale that occurs in the Phosphoria Formation and the underlying regional Wells Formation aquifer. Using a reconnaissance-level transport model, mines in the watershed were prioritized for remediation and for comparing the results of simulations of remediation scenarios with a baseline of no remediation, for which Se concentration in the river will exceed the aquatic standard along an extensive length. An accurate simulation of recharge distribution around the watershed and simulated flux to the river is essential. Remediation of mines north of the river will substantially decrease the size of the Se plume, although significant Se will continue to discharge to the river. Similarly, remediation of three mines south of the river would decrease the Se discharge to the river but allow substantial amounts to remain stored in the groundwater north and far south of the river. A lack of calibration data is not a reason to forgo remediation, but rather ongoing data collection can be used to fine-tune plans as they are implemented.     

The geology of kimberlite pipes of the Ekati property, Northwest Territories, Canada

This paper reviews key characteristics of kimberlites on the Ekati property, NWT, Canada. To date 150 kimberlites have been discovered on the property, five of which are mined for diamonds. The kimberlites intrude Archean basement of the central Slave craton. Numerous Proterozoic diabase dykes intrude the area. The Precambrian rocks are overlain by Quaternary glacial sediments. No Phanerozoic rocks are present. However, mudstone xenoliths and disaggregated sediment within the kimberlites indicate that late-Cretaceous and Tertiary cover (likely <200 m) was present at the time of emplacement. The Ekati kimberlites range in age from 45 to 75 Ma. They are mostly small pipe-like bodies (surface area mostly <3 ha but up to 20 ha) that typically extend to projected depths of 400–600 m below current surface. Pipe morphologies are strongly controlled by joints and faults. The kimberlites consist primarily of variably bedded volcaniclastic kimberlite (VK). This is dominated by juvenile constituents (olivine and lesser kimberlitic ash) and variable amounts of exotic sediment (primarily mud), with minor amounts of xenolithic wall-rock material (generally <5%). Kimberlite types include: mud-rich resedimented VK (mRVK); olivine-rich VK (oVK); sedimentary kimberlite; primary VK (PVK); tuffisitic kimberlite (TK) and magmatic kimberlite (MK). The presence and arrangement of these rock types varies widely. The majority of bodies are dominated by oVK and mRVK, but PVK is prominent in the lower portions of certain kimberlites. TK is rare. MK occurs primarily as precursor dykes but, in a few cases, forms pipe-filling intrusions. The internal geology of the kimberlites ranges from simple single-phase pipes (RVK or MK), to complex bodies with multiple, distinct units of VK. The latter include pipes infilled with steep, irregular VK blocks/wedges and at least one case in which the pipe is occupied by well-defined sub-horizontal VK phases, including a unique, 100-m-thick graded sequence. The whole-rock compositions of VK samples suggest significant loss of kimberlitic fines during eruption followed by variable dilution by surface sediment and concurrent incorporation of kimberlitic ash. Diamond distribution within the kimberlites reflects the amount and nature of mantle material sampled by individual kimberlite phases, but is modified considerably by eruption and depositional processes. The characteristics of the Ekati kimberlites are consistent with a two-stage emplacement process: (1) explosive eruption/s causing vent clearing followed by formation of a significant tephra rim/cone of highly fragmented, olivine-enriched juvenile material with varying amounts of kimberlitic ash and surface sediments (predominantly mud); and (2) infilling of the vent by direct deposition from the eruption column and/or resedimentation of crater rim materials. The presence of less fragmented, juvenile-rich PVK in the lower portions of certain pipes and the intrusion of large volumes of MK to shallow levels in some bodies suggest emplacement of relatively volatile-depleted, less explosive kimberlite in the later stages of pipe formation and/or filling. Explosive devolatilisation of CO₂-rich kimberlite magma is interpreted to have been the dominant eruption mechanism, but phreatomagmatism is thought to have played a role and, in certain cases, may have been dominant.     

Inclusions in diamonds from the K14 and K10 kimberlites, Buffalo Hills, Alberta, Canada: diamond growth in a plume?

Analyses of mineral inclusions, carbon isotopes, nitrogen contents and nitrogen aggregation states in 29 diamonds from two Buffalo Hills kimberlites in northern Alberta, Canada were conducted. From 25 inclusion bearing diamonds, the following paragenetic abundances were found: peridotitic (48%), eclogitic (32%), eclogitic/websteritic (8%), websteritic (4%), ultradeep? (4%) and unknown (4%). Diamonds containing mineral inclusions of ferropericlase, and mixed eclogitic-asthenospheric-websteritic and eclogitic-websteritic mineral associations suggests the possibility of diamond growth over a range of depths and in a variety of mantle environments (lithosphere, asthenosphere and possibly lower mantle).

Eclogitic diamonds have a broad range of C-isotopic composition (δ¹³C=−21‰ to −5‰). Peridotitic, websteritic and ultradeep diamonds have typical mantle C-isotope values (δ¹³C=−4.9‰ av.), except for two ¹³C-depleted peridotitic (δ¹³C=−11.8‰, −14.6‰) and one ¹³C-depleted websteritic diamond (δ¹³C=−11.9‰). Infrared spectra from 29 diamonds identified two diamond groups: 75% are nitrogen-free (Type II) or have fully aggregated nitrogen defects (Type IaB) with platelet degradation and low to moderate nitrogen contents (av. 330 ppm-N); 25% have lower nitrogen aggregation states and higher nitrogen contents (30% IaB; <1600 ppm-N).

The combined evidence suggests two generations of diamond growth. Type II and Type IaB diamonds with ultradeep, peridotitic, eclogitic and websteritic inclusions crystallised from eclogitic and peridotitic rocks while moving in a dynamic environment from the asthenosphere and possibly the lower mantle to the base of the lithosphere. Mechanisms for diamond movement through the mantle could be by mantle convection, or an ascending plume. The interaction of partial melts with eclogitic and peridotitic lithologies may have produced the intermediate websteritic inclusion compositions, and can explain diamonds of mixed parageneses, and the overlap in C-isotope values between parageneses. Strong deformation and extremely high nitrogen aggregation states in some diamonds may indicate high mantle storage temperatures and strain in the diamond growth environment. A second diamond group, with Type IaA–IaB nitrogen aggregation and peridotitic inclusions, crystallised at the base of the cratonic lithosphere. All diamonds were subsequently sampled by kimberlites and transported to the Earth's surface.     

Pattern,style and timing of British–Irish Ice Sheet advance and retreat over the last 45 000 years: evidence from NW Scotland and the adjacent continental shelf

Predicting the future response of ice sheets to climate warming and rising global sea level is important but difficult. This is especially so when fast-flowing glaciers or ice streams, buffered by ice shelves, are grounded on beds below sea level. What happens when these ice shelves are removed? And how do the ice stream and the surrounding ice sheet respond to the abruptly altered boundary conditions? To address these questions and others we present new geological, geomorphological, geophysical and geochronological data from the ice-stream-dominated NW sector of the last British–Irish Ice Sheet (BIIS). The study area covers around 45 000 km² of NW Scotland and the surrounding continental shelf. Alongside seabed geomorphological mapping and Quaternary sediment analysis, we use a suite of over 100 new absolute ages (including cosmogenic-nuclide exposure ages, optically stimulated luminescence ages and radiocarbon dates) collected from onshore and offshore, to build a sector-wide ice-sheet reconstruction combining all available evidence with Bayesian chronosequence modelling. Using this information we present a detailed assessment of ice-sheet advance/retreat history, and the glaciological connections between different areas of the NW BIIS sector, at different times during the last glacial cycle. The results show a highly dynamic, partly marine, partly terrestrial, ice-sheet sector undergoing large size variations in response to sub-millennial-scale climatic (Dansgaard–Oeschger) cycles over the last 45 000 years. Superimposed on these trends we identify internally driven instabilities, operating at higher frequency, conditioned by local topographic factors, tidewater dynamics and glaciological feedbacks during deglaciation. Specifically, our new evidence indicates extensive marine-terminating ice-sheet glaciation of the NW BIIS sector during Greenland Stadials 12 to 9 – prior to the main ‘Late Weichselian’ ice-sheet glaciation. After a period of restricted glaciation, in Greenland Interstadials 8 to 6, we find good evidence for rapid renewed ice-sheet build-up in NW Scotland, with the Minch ice-stream terminus reaching the continental shelf edge in Greenland Stadial 5, perhaps only briefly. Deglaciation of the NW sector took place in numerous stages. Several grounding-zone wedges and moraines on the mid- and inner continental shelf attest to significant stabilizations of the ice-sheet grounding line, or ice margin, during overall retreat in Greenland Stadials 3 and 2, and to the development of ice shelves. NW Lewis was the first substantial present-day land area to deglaciate, in the first half of Greenland Stadial 3 at a time of globally reduced sea-level c. 26 kabp , followed by Cape Wrath at c. 24 kabp. The topographic confinement of the Minch straits probably promoted ice-shelf development in early Greenland Stadial 2, providing the ice stream with additional support and buffering it somewhat from external drivers. However, c. 20–19 kabp , as the grounding-line migrated into shoreward deepening water, coinciding with a marked change in marine geology and bed strength, the ice stream became unstable. We find that, once underway, grounding-line retreat proceeded in an uninterrupted fashion with the rapid loss of fronting ice shelves – first in the west, then the east troughs – before eventual glacier stabilization at fjord mouths in NW Scotland by ~17 kabp. Around the same time, ~19–17 kabp , ice-sheet lobes readvanced into the East Minch – possibly a glaciological response to the marine-instability-triggered loss of adjacent ice stream (and/or ice shelf) support in the Minch trough. An independent ice cap on Lewis also experienced margin oscillations during mid-Greenland Stadial 2, with an ice-accumulation centre in West Lewis existing into the latter part of Heinrich Stadial 1. Final ice-sheet deglaciation of NW mainland Scotland was punctuated by at least one other coherent readvance at c. 15.5 kabp , before significant ice-mass losses thereafter. At the glacial termination, c. 14.5 kabp , glaciers fed outwash sediment to now-abandoned coastal deltas in NW mainland Scotland around the time of global Meltwater Pulse 1A. Overall, this work on the BIIS NW sector reconstructs a highly dynamic ice-sheet oscillating in extent and volume for much of the last 45 000 years. Periods of expansive ice-sheet glaciation dominated by ice-streaming were interspersed with periods of much more restricted ice-cap or tidewater/fjordic glaciation. Finally, this work indicates that the role of ice streams in ice-sheet evolution is complex but mechanistically important throughout the lifetime of an ice sheet – with ice streams contributing to the regulation of ice-sheet health but also to the acceleration of ice-sheet demise via marine ice-sheet instabilities.     

Vertical Stacking Statistics of Multi-facies Object-Based Models

Mathematical Geosciences - Equations describing facies proportions and amalgamation ratios are derived for randomly placed objects belonging to two or three foreground facies embedded in a...     

Technical, financial and location criteria for the design of land application system treatment

Conventional water treatments when downsized to cover the demands of small and scattered populations are costly and ineffective. In fact, many small towns cannot afford the high average treatment costs nor provide the skill and complex management requirements. Even when these facilities are available, waste water might not be properly treated. Land application system treatments (LAST) have been shown to be preferable on the basis of their technical effectiveness and financial viability. The recently approved Spanish regulations (RDL 1620/2007, pp 50369–50661) represent an opportunity for the implementation of well-designed LASTs. This paper presents an assessment methodology for the design of LASTs by combining technical, financial and location criteria. The method is applied to the design of LASTs to cover the water disposal demand of twelve municipalities located within the protected natural landscape, and surrounding areas of ‘El Rebollar’, Salamanca, Spain.     

Modeling nonlinear seagrass clonal growth: Assessing the efficiency of space occupation across the seagrass flora

The clonal growth of 9 seagrass species was modeled using a simulation model based on observed clonal growth rules (i.e., spacer length, rhizome elongation rates, branching rates, branching angle) and shoot mortality rates for seagrass species. The results of the model confirmed the occurrence of complex, nonlinear growth of seagrass clones derived from internal dynamics of space occupation. The modeled clones progressed from a diffuse-limited aggregation (DLA), dendritic growth, identified with a guerrilla strategy of space occupation, to a compact (Eden) growth, comparable to the phalanx strategy of space occupation, once internal recolonization of gaps, left by dead shoots within the clone, begins. The time at which seagrass clones shifted from diffuse limited to compact growth was predictable from the branching angle and frequency of the species and varied from 1 yr to several decades among species. As a consequence the growth behavior and the apparent growth strategy of the species changes with the development of the clones. The results of the model demonstrate that the emergent complexity of seagrass clonal growth is contained within the simple set of growth rules that can be used to represent clonal growth.     

The stability of a remediated bed in Hamilton Harbour,Lake Ontario,Canada

ABSTRACT In situ measurements of lakebed sediment erodibility were made on three sites in Hamilton Harbour, Lake Ontario, using the benthic flume Sea Carousel. Three methods of estimating the surface erosion threshold (τ_c(0)) from a Carousel time series were evaluated: the first method fits measures of bed strength to eroded depth (the failure envelope) and evaluates threshold as the surface intercept; the second method regresses mean erosion rate (E_m) with bed shear stress and solves for the floc erosion rate (E_f) to derive the threshold for E_m = E_f = 1 × 10^?5 kg m^?2 s^?1; the third method extrapolates a regression of suspended sediment concentration (S) and fluid transmitted bed shear stress (τ₀) to ambient concentrations. The first field site was undisturbed (C) and acted as a control; the second (W) was disturbed through ploughing and water injection as part of lakebed treatment, whereas the third site (OIP) was disturbed and injected with an oxidant used for remediation of contaminated sediment. The main objectives of this study were: (1) to evaluate the three different methods of deriving erosion threshold; (2) to compare the physical behaviour of lacustrine sediments with their marine estuarine counterparts; and (3) to examine the effects of ploughing and chemical treatment of contaminated sediment on bed stability. Five deployments of Sea Carousel were carried out at the control site. Mean erosion thresholds for the three methods were: τ_c(0) = 0·5 (±0·06), 0·27 (±0·01) and 0·34 (±0·03) Pa respectively. Method 1 overpredicted bed strength as it was insensitive to effects in the surface 1–2 mm, and the fit of the failure envelope was also highly subjective. Method 2 exhibited a wide scatter in the data (low correlation coefficients), and definition of the baseline erosion rate (E_f) is largely arbitrary in the literature. Method 3 yielded stable (high correlation coefficients), reproducible and objective results and is thus recommended for evaluation of the erosion threshold. The results of this method correlated well with sediment bulk density and followed the same trend as marine counterparts from widely varying sites. Mass settling rates, expressed as a decay constant, k, of S(t), were strongly related to the maximum turbidity at the onset of settling (S_max) and were also in continuity with marine counterparts. Thus, it appears that differences in salinity had little effect on mass settling rates in the examples presented, and that biological activity dominated any effects normally attributable to changes in salinity. Bedload transport of eroded aggregates (2–4 mm in diameter) took place by rolling below a mean tangential flow velocity (U_y) of 0·32 ms^?1 and by saltation at higher velocities. Mass transport as bedload was a maximum at U_y = 0·4 ms^?1, although bedload never exceeded 1% of the suspended load. The proportion of material moving as bedload was greatest at the onset of erosion but decreased as flow competence increased. Given the low bulk density and strength of the lakebed sediment, the presence of a bedload component is notable. Bedload transport over eroding cohesive substrates should be greater in estuaries, where both sediment density and strength are usually higher. Significant differences between the ploughed and control sites were apparent in both the erosion rate and the friction coefficient (φ), and suggest that bed recovery after disruption is rapid (< 24 h). τ_c(0) increased linearly with time after ploughing and recovered to the control mean value within 3 days. The friction coefficient was reduced to zero by ploughing (diagnostic of fluidization), but increased linearly with time, regaining control values within 6 days. No long‐term reduction in bed strength due to remediation was apparent.     

Recurrent mesoproterozoic continental magmatism in South-Central Norway

We report U–Pb dates and Lu–Hf isotope data, obtained by LAM-ICPMS, for zircons from metamorphic rocks of the Setesdalen valley, situated in the Telemark block south of the classic Telemark region of southern Norway. The samples include infracrustal rocks from the metamorphic basement, metaigneous rocks and metasediments from the Byglandsfjorden supracrustal cover sequence, and metaigneous rocks which intruded the whole succession. The main crustal evolution took place from 1,550–1,020 Ma, beginning with the emplacement of juvenile tonalitic melts; the contribution of older crustal material increased with time. Around 1,320 Ma, further addition of juvenile material occurred, involving both mafic and felsic melts, metamorphism and deformation. Acid magmas with high FeO*/MgO were intruded at 1,215 Ma, coinciding with underplating elsewhere in South Norway. The period starting at 1,215 Ma is represented by supracrustal rocks, principally metarhyolites with minor mafic material and immature sediments of the Byglandsfjorden Group. The crust generation processes ended with the intrusion of diorites and granodiorites at 1,030 Ma, late in the Sveconorwegian orogeny. Regional processes of metamorphism and deformation (around 1,290 and 1,000 Ma) can be related to the assembly of Rodinia. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.