Shallow site archaeology: Artifact dispersal,stratigraphy, and radiocarbon dating at the Barger Gulch locality B Folsom Site,Middle Park,Colorado

Shallowly buried archaeological sites are particularly susceptible to surface and subsurface disturbance processes. Yet, because cultural deposition often operates on short time scales relative to geologic deposition, vertical artifact distributions can be used to clarify questions of site formation. In particular, patterns in artifact distributions that cannot be explained by occupation histories must be explained by natural processes that have affected sites. Buried only 10–50 cm beneath the ground surface for 10,450 ¹⁴C yr, the Folsom component at Barger Gulch Locality B (Middle Park, Colorado) exhibits many signs of post‐depositional disturbance. Through examination of variation in the vertical distribution of the artifact assemblage, we are able to establish that only a Folsom component is present. Using vertical artifact distributions, stratigraphy, and radiocarbon dating, we are able to reconstruct the series of events that have impacted the site. The Folsom occupation (˜10,450 ¹⁴C yr B.P.) was likely initially buried in a late‐Pleistocene eolian silt loam. Erosion brought the artifacts to rest on a deflation surface at some time prior to 9400 ¹⁴C yr B.P. A mollic epipedon formed in sediments that accumulated between 9400 and 7000 ¹⁴C yr B.P. Some time after 5200 ¹⁴C yr B.P., this soil was partially truncated, and artifacts that had previously dispersed upward created a secondary lag at its upper contact. This surface was buried again and artifact dispersal continued. © 2005 Wiley Periodicals, Inc.     

Use of GRACE determined secular gravity rates for glacial isostatic adjustment studies in North-America

Monthly geopotential spherical harmonic coefficients from the GRACE satellite mission are used to determine their usefulness and limitations for studying glacial isostatic adjustment (GIA) in North-America. Secular gravity rates are estimated by unweighted least-squares estimation using release 4 coefficients from August 2002 to August 2007 provided by the Center for Space Research (CSR), University of Texas. Smoothing is required to suppress short wavelength noise, in addition to filtering to diminish geographically correlated errors, as shown in previous studies. Optimal cut-off degrees and orders are determined for the destriping filter to maximize the signal to noise ratio. The halfwidth of the Gaussian filter is shown to significantly affect the sensitivity of the GRACE data (with respect to upper mantle viscosity and ice loading history). Therefore, the halfwidth should be selected based on the desired sensitivity.It is shown that increase in water storage in an area south west of Hudson Bay, from the summer of 2003 to the summer of 2006, contributes up to half of the maximum estimated gravity rate. Hydrology models differ in the predictions of the secular change in water storage, therefore even 4-year trend estimates are influenced by the uncertainty in water storage changes. Land ice melting in Greenland and Alaska has a non-negligible contribution, up to one-fourth of the maximum gravity rate.The estimated secular gravity rate shows two distinct peaks that can possibly be due to two domes in the former Pleistocene ice cover: west and south east of Hudson Bay. With a limited number of models, a better fit is obtained with models that use the ICE-3G model compared to the ICE-5G model. However, the uncertainty in interannual variations in hydrology models is too large to constrain the ice loading history with the current data span. For future work in which GRACE will be used to constrain ice loading history and the Earth's radial viscosity profile, it is important to include realistic uncertainty estimates for hydrology models and land ice melting in addition to the effects of lateral heterogeneity.     

Generalized likelihood uncertainty estimation (GLUE) using adaptive Markov Chain Monte Carlo sampling

In the last few decades hydrologists have made tremendous progress in using dynamic simulation models for the analysis and understanding of hydrologic systems. However, predictions with these models are often deterministic and as such they focus on the most probable forecast, without an explicit estimate of the associated uncertainty. This uncertainty arises from incomplete process representation, uncertainty in initial conditions, input, output and parameter error. The generalized likelihood uncertainty estimation (GLUE) framework was one of the first attempts to represent prediction uncertainty within the context of Monte Carlo (MC) analysis coupled with Bayesian estimation and propagation of uncertainty. Because of its flexibility, ease of implementation and its suitability for parallel implementation on distributed computer systems, the GLUE method has been used in a wide variety of applications. However, the MC based sampling strategy of the prior parameter space typically utilized in GLUE is not particularly efficient in finding behavioral simulations. This becomes especially problematic for high-dimensional parameter estimation problems, and in the case of complex simulation models that require significant computational time to run and produce the desired output. In this paper we improve the computational efficiency of GLUE by sampling the prior parameter space using an adaptive Markov Chain Monte Carlo scheme (the Shuffled Complex Evolution Metropolis (SCEM-UA) algorithm). Moreover, we propose an alternative strategy to determine the value of the cutoff threshold based on the appropriate coverage of the resulting uncertainty bounds. We demonstrate the superiority of this revised GLUE method with three different conceptual watershed models of increasing complexity, using both synthetic and real-world streamflow data from two catchments with different hydrologic regimes.     

Recent Results from Studies of Electric Discharges in the Mesosphere

The paper reviews recent advances in studies of electric discharges in the stratosphere and mesosphere above thunderstorms, and their effects on the atmosphere. The primary focus is on the sprite discharge occurring in the mesosphere, which is the most commonly observed high altitude discharge by imaging cameras from the ground, but effects on the upper atmosphere by electromagnetic radiation from lightning are also considered. During the past few years, co-ordinated observations over Southern Europe have been made of a wide range of parameters related to sprites and their causative thunderstorms. Observations have been complemented by the modelling of processes ranging from the electric discharge to perturbations of trace gas concentrations in the upper atmosphere. Observations point to significant energy deposition by sprites in the neutral atmosphere as observed by infrasound waves detected at up to 1000 km distance, whereas elves and lightning have been shown significantly to affect ionization and heating of the lower ionosphere/mesosphere. Studies of the thunderstorm systems powering high altitude discharges show the important role of intracloud (IC) lightning in sprite generation as seen by the first simultaneous observations of IC activity, sprite activity and broadband, electromagnetic radiation in the VLF range. Simulations of sprite ignition suggest that, under certain conditions, energetic electrons in the runaway regime are generated in streamer discharges. Such electrons may be the source of X- and Gamma-rays observed in lightning, thunderstorms and the so-called Terrestrial Gamma-ray Flashes (TGFs) observed from space over thunderstorm regions. Model estimates of sprite perturbations to the global atmospheric electric circuit, trace gas concentrations and atmospheric dynamics suggest significant local perturbations, and possibly significant meso-scale effects, but negligible global effects.     

Dynamics of multiple intertidal bars over semi‐diurnal and lunar tidal cycles,North Lincolnshire,England

Multiple intertidal bars are common features of wave‐dominated sandy beaches, yet their short‐term (<1 month) and small‐scale (<1 km) morphology and dynamics remain poorly understood. This study describes the morphodynamics of multiple intertidal bars in North Lincolnshire, England, during single and lunar tidal cycles under two contrasting conditions – first when significant wave height was <0·5 m and second when significant wave height frequently exceeded 1 m. The relative importance of swash, surf and shoaling processes in determining morphological change was examined using detailed field observations and a numerical model. The beach featured four intertidal bars and both cross‐shore and longshore bar morphology evolved during the field investigation, particularly under medium to high wave‐energy conditions. Numerical modelling suggests shoaling processes are most common on the seaward two bars under calm wave conditions (H_s < 0·5 m) and that surf zone processes become more common during neap tides and under more energetic (H_s < 0·5 m) conditions. Surf processes dominate the inner two bars, though swash influence increases in a landward direction. The numerical modelling results combined with low tide survey data and high‐resolution morphological measurements strongly suggest changes in the intertidal bar morphology are accomplished by surf zone processes rather than by shoaling wave or swash processes. This is because shoaling waves do not induce significant sediment transport to have any morphological effect, whereas swash action generally does not have enough scope to act as the swash zone is much narrower than the surf zone. It was found, however, that the absolute rate of morphological change under swash action and surfzone processes are of similar magnitudes and that swash action may induce a significant amount of local morphological change when the high tide mark is located on the upper bar, making this process important for bar morphodynamics. Copyright © 2007 John Wiley & Sons, Ltd.     

Behaviour of deep immersed tunnel under combined normal fault rupture deformation and subsequent seismic shaking

Immersed tunnels are particularly sensitive to tensile and compressive deformations such as those imposed by a normal seismogenic fault rupturing underneath, and those generated by the dynamic response due to seismic waves. The paper investigates the response of a future 70 m deep immersed tunnel to the consecutive action of a major normal fault rupturing in an earthquake occurring in the basement rock underneath the tunnel, and a subsequent strong excitation from a different large-magnitude seismic event that may occur years later. Non-linear finite elements model the quasi-static fault rupture propagation through the thick soil deposit overlying the bedrock and the ensuing interaction of the rupture with the immersed tunnel. It is shown that despite imposed bedrock offset of 2 m, net tension or excessive compression between tunnel segments could be avoided with a suitable design of the joint gaskets. Then, the already deformed (“injured”) structure is subjected to strong asynchronous seismic shaking. The thick-walled tunnel is modelled as a 3-D massive flexural beam connected to the soil through properly-calibrated nonlinear interaction springs and dashpots, the supports of which are subjected to the free-field acceleration time histories. The latter, obtained with 1-D wave propagation analysis, are then modified to account for wave passage effects. The joints between tunnel segments are modeled with special non-linear hyper-elastic elements, properly accounting for their 7-bar longitudinal hydrostatic pre-stressing. Sliding is captured with special gap elements. The effect of segment length and joint properties is explored parametrically. A fascinating conclusion emerges in all analysed cases for the joints between segments that were differentially deformed after the quasi-static fault rupture: upon subsequent very strong seismic shaking, overstressed joints de-compress and understressed joints re-compress—a “healing” process that leads to a more uniform deformation profile along the tunnel. This is particularly beneficial for the precariously de-compressed joint gaskets. Hence, the safety of the immersed tunnel improves with “subsequent” strong seismic shaking!     

Cross-shelf transport into nearshore waters due to shoaling internal tides in San Pedro Bay, CA

In the summer of 2001, a coastal ocean measurement program in the southeastern portion of San Pedro Bay, CA, was designed and carried out. One aim of the program was to determine the strength and effectiveness of local cross-shelf transport processes. A particular objective was to assess the ability of semidiurnal internal tidal currents to move suspended material a net distance across the shelf. Hence, a dense array of moorings was deployed across the shelf to monitor the transport patterns associated with fluctuations in currents, temperature and salinity. An associated hydrographic program periodically monitored synoptic changes in the spatial patterns of temperature, salinity, nutrients and bacteria. This set of measurements show that a series of energetic internal tides can, but do not always, transport subthermocline water, dissolved and suspended material from the middle of the shelf into the surfzone. Effective cross-shelf transport occurs only when (1) internal tides at the shelf break are strong and (2) subtidal currents flow strongly downcoast. The subtidal downcoast flow causes isotherms to tilt upward toward the coast, which allows energetic, nonlinear internal tidal currents to carry subthermocline waters into the surfzone. During these events, which may last for several days, the transported water remains in the surfzone until the internal tidal current pulses and/or the downcoast subtidal currents disappear. This nonlinear internal tide cross-shelf transport process was capable of carrying water and the associated suspended or dissolved material from the mid-shelf into the surfzone, but there were no observation of transport from the shelf break into the surfzone. Dissolved nutrients and suspended particulates (such as phytoplankton) transported from the mid-shelf into the nearshore region by nonlinear internal tides may contribute to nearshore algal blooms, including harmful algal blooms that occur off local beaches.     

Extreme wave event along the Guyana coastline in October 2005

The objective of the research is to analyze in detail the causes and consequences of the unusual event at the coast of Guyana (South America) during October 16–19, 2005. Several sea defense structures were damaged and flooding of low-lying areas occurred. A data analysis of offshore wave and water level characteristics shows an abrupt change in wave direction from east to almost north on October 16, 2005 and a sudden increase in the offshore peak period up to extreme values. The offshore significant wave height was also relatively high, and these wave characteristics coincided with springtide conditions. The long-wave periods and the sharp transition in wave direction indicate that this event is associated with swell waves generated by a depression far away. An analysis of hurricanes and depressions reveals that a severe depression in the Northern Atlantic Ocean during October 11–15 was the origin of this swell event. Numerical computations with SWAN have been carried out to investigate the propagation of the offshore wave characteristics towards the shoreline. The SWAN model includes wave damping due to the presence of soft mud deposits. A calibration of the parameters has been carried out using joint offshore and onshore wave data from November 2006. The numerical simulations of the event in October 2005 clearly demonstrate that the mud banks damp the wave heights, but have almost no effect on the peak period. The resulting waves at the steep sea walls can be classified as surging waves causing severe runup and overtopping. The obtained insights are translated into practical recommendations for the Guyana Sea and River Defence Division in Guyana to build a sustainable management and maintenance of the sea defenses in the future.     

Distinguishing volcanic debris avalanche deposits from their reworked products: the Perrier sequence (French Massif Central)

Debris avalanches associated with volcanic sector collapse are usually high-volume high-mobility phenomena. Debris avalanche deposit remobilisation by cohesive debris flows and landslides is common, so they can share textural characteristics such as hummocks and jigsaw cracks. Distinguishing original deposits from reworked products is critical for geological understanding and hazard assessment because of their different origin, frequency and environmental impact. We present a methodology based on field evidence to differentiate such epiclastic breccias. Basal contact mapping constrained by accurate altitude and location data allows the reconstruction of deposit stratigraphy and geometry. Lithological analysis helps to distinguish the different units. Incorporation structures, kinematic indicators and component mingling textures are used to characterise erosion and transport mechanisms. We apply this method to the enigmatic sequence at Perrier (French Massif Central), where four units (U1–U4) have been interpreted either as debris flow or debris avalanche deposits. The sequence results from activity on the Monts Dore Volcano about 2 Ma ago. The epiclastic units are matrix supported with an almost flat top. U2 and U3 have clear debris flow deposit affinities such as rounded clasts and intact blocks (no jigsaw cracks). U1 and U4 have jigsaw cracked blocks with matrix injection and stretched sediment blocks. U1 lacks large blocks (>10 m wide) and has a homogenous matrix with an upward increase of trapped air vesicle content and size. This unit is interpreted as a cohesive debris flow deposit spawned from a debris avalanche upstream. In contrast, U4 has large mega-blocks (up to 40 m wide), sharp contacts between mixed facies zones with different colours and numerous jigsaw fit blocks (open jigsaw cracks filled by monogenic intra-clast matrix). Mega-blocks are concentrated near the deposit base and are spatially associated with major substratum erosion. This deposit has a debris avalanche distal facies with local debris flow affinities due to partial water saturation. We also identify two landslide deposits (L1 and L2) resulting from recent reworking that has produced a similar facies to U1 and U4. These are distinguishable from the original deposits, as they contain blocks of mixed U1/U4 facies, a distinctly less consolidated and more porous matrix and a fresh hummocky topography. This work shows how to differentiate epiclastic deposits with similar characteristics, but different origins. In doing so, we improve understanding of present and past instability of the Monts Dore and identify present landslide hazards at Perrier.     

A simplified model for lateral response of large diameter caisson foundations—Linear elastic formulation

The transient response of large embedded foundation elements of length-to-diameter aspect ratio D/B=2–6 is characterized by a complex stress distribution at the pier–soil interface that cannot be adequately represented by means of existing models for shallow foundations or flexible piles. On the other hand, while three-dimensional (3D) numerical solutions are feasible, they are infrequently employed in practice due to their associated cost and effort. Prompted by the scarcity of simplified models for design in current practice, we here develop an analytical model that accounts for the multitude of soil resistance mechanisms mobilized at their base and circumference, while retaining the advantages of simplified methodologies for the design of non-critical facilities. The characteristics of soil resistance mechanisms and corresponding complex spring functions are developed on the basis of finite element simulations, by equating the stiffness matrix terms and/or overall numerically computed response to the analytical expressions derived by means of the proposed Winkler model. Sensitivity analyses are performed for the optimization of the truncated numerical domain size, the optimal finite element size and the far-field dynamic boundary conditions to avoid spurious wave reflections. Numerical simulations of the transient system response to vertically propagating shear waves are next successfully compared to the analytically predicted response. Finally, the applicability of the method is assessed for soil profiles with depth-varying properties. The formulation of frequency-dependent complex spring functions including material damping is also described, while extension of the methodology to account for nonlinear soil behavior and soil–foundation interface separation is described in the conclusion and is being currently investigated.