Salinograph trends as indicators of the recession characteristics of stream components

Compared to hydrograph recession analysis, which is widely applied in engineering hydrology, the quantitative assessment of stream salinity with time (i.e. the salinograph) has received significantly less attention. In particular, while in many previous hydrological studies an inverse relationship between hydrograph and salinograph responses is apparent, the concept of salinity accession (the inversely related salinity counterpart to hydrograph recession) has not been introduced nor quantitatively evaluated in previous literature. In this study, we conduct a mathematical analysis of salinograph accession, and determine new quantitative relationships between salinity accession and hydrograph recession parameters. An equation is formulated that reproduces the general trend in salinity accession. A salinity accession parameter k_c is then introduced and is shown to be the ratio of direct runoff to total stream flow recession parameters: k_r/k. The groundwater recession parameter k_g was estimated using a simple and rapid method that uses both salinograph and hydrograph data. Salinity accession type‐curves illustrate that under certain conditions, the relative steepness of individual salinographs is dependent upon the ratio of groundwater salinity to direct runoff salinity: C_g/C_r. The salinity accession algorithms are applied to two contrasting field settings: Scott Creek, South Australia and Sandy Creek, northern Queensland, Australia. It was found that k_g > k during periods of obvious stream flow recession, for the events analysed. Salinograph accession behaviour was fairly similar for both sites, despite contrasting environments. Using assumed end‐member salinities for groundwater and direct runoff based upon field observations, the behaviour of k_c from the Scott Creek site was approximately reproduced by varying the initial groundwater to runoff flow ratio: Q_g0/Q_r0, within reasonable parameter ranges. The use of salinograph information when used in addition to standard hydrograph analyses provided useful information on recession characteristics of stream components. Copyright © 2008 John Wiley & Sons, Ltd.     

Currents and sediment transport in the Mississippi Canyon and effects of Hurricane Georges

The temporal variability in currents, temperature, and particulate matter concentration were measured in the Mississippi Canyon axis where the thalweg was 300 m deep from May–July and August–November 1998 using current meters, thermographs, a light-scattering sensor, and sediment traps. Canyon sediments were sampled by coring and observed using an ROV video camera. Currents in the upper Mississippi Canyon generally oscillated up/down canyon with diurnal periodicity and were bottom-intensified. Mean current speed at 3.5 mab was approximately 8 cm s^?1 during both deployments, reaching maximum speeds of over 50 cm s^?1 under normal conditions. Based on current velocities, critical bed shear stress for resuspension of canyon-floor sediments was exceeded about 30% of the time during both deployments. In late September, Hurricane Georges passed 150 km NE of the study site, significantly intensifying current velocities, bed shear stress, resuspension, trap fluxes and temperature fluctuations. As the hurricane passed, maximum current speed reached 68 cm^?s and temperature decreased ～7 °C in less than two hours. Critical bed shear stress for sediment resuspension was exceeded approximately 50% of the time during the five days of hurricane influence. Further evidence for sediment resuspension was the five-fold (and perhaps 70–130 fold) increase in trap fluxes and compositional similarities between canyon surface sediment and material collected by traps.     

Detection Time for Plausible Changes in Annual Precipitation, Evapotranspiration, and Streamflow in Three Mississippi River Sub-Basins

We use diagnostic studies of off-line variable infiltration capacity (VIC) model simulations of terrestrial water budgets and 21st-century climate change simulations using the parallel climate model (PCM) to estimate the time required to detect predicted changes in annual precipitation (P), evapotranspiration (E), and discharge (Q) in three sub-basins of the Mississippi River Basin. Time series lengths on the order of 50–350 years are required to detect plausible P, E, and Q trends in the Missouri, Ohio, and Upper Mississippi River basins. Approximately 80–160, 50, and 140–350 years, respectively, are needed to detect the predicted P, E, and Q trends with a high degree of statistical confidence. These detection time estimates are based on conservative statistical criteria (α = 0.05 and β = 0.10) associated with low probability of both detecting a trend when it is not occurring (Type I error) and not detecting a trend when it is occurring (Type II error). The long detection times suggest that global-warming-induced changes in annual basin-wide hydro-climatic variables that may already be occurring in the three basins probably cannot yet be detected at this level of confidence. Furthermore, changes for some variables that may occur within the 21st century might not be detectable for many decades or until the following century – this may or may not be the case for individual recording station data. The long detection times for streamflow result from comparatively low signal-to-noise ratios in the annual time series. Finally, initial estimates suggest that faster detection of acceleration in the hydrological cycle may be possible using seasonal time series of appropriate hydro-climatic variables, rather than annual time series.     

Geology and eruptive history of the Rabaul Caldera area, Papua New Guinea

Rabaul Caldera is the most recently active (1937–1943) of four adjoining volcanic centres aligned north-south through the northern extremity of eastern New Britain. Geological mapping after the 1983–1985 Rabaul seismic and deformation crisis has partially revealed a long and complex eruption history dominated by numerous explosive eruptions, the largest accompanied by caldera collapse. The oldest exposed eruptives are the basaltic pre-caldera cone Tovanumbatir Lavas K/Ar dated at 0.5 Ma. The dacitic Rabaul Quarry Lavas exposed in the caldera wall and K/Ar dated at 0.19 Ma, are overlain by a sequence of dacitic and andesitic pyroclastic flow and fall deposits. Uplifted coral reef limestones, interbedded within the pyroclastic sequence on the northeast coast, suggest that explosive eruptions in the Rabaul area had commenced prior to the 0.125 Ma last interglacial high sea level stand. The pyroclastic sequence includes the large Boroi Ignimbrites and Malaguna Pyroclastics both ⁴⁰Ar/³⁹Ar dated at about 0.1 Ma, and the Barge Tunnel Ignimbrite ⁴⁰Ar/³⁹Ar dated at around 0.04 Ma. Few reliable ages exist for the many younger eruptives. These include Holocene ignimbrites of the latest caldera-forming eruptions—the Raluan Pyroclastics variously dated (¹⁴C) at either about 3500 or 7000 yr B.P., and the ca. 1400 yr B.P. Rabaul Pyroclastics. At least eight intracaldera eruptions have occurred since the 1400 yr B.P. collapse, building small pyroclastic and lava cones within the caldera.A major erosional episode is evident as a widespread unconformity in the upper pyroclastic stratigraphy at Rabaul. Lacking relevant radiometric ages, this episode is assumed to have occurred during last glaciation low sea levels and is here arbitarily dated at ca. ?20 ka. At least five, possibly nine, significant ignimbrite eruptions have occurred at Rabaul during the last ?20 ka. The new eruptive history differs considerably from that previously published, which considered ignimbrite eruption and caldera collapse to have first occurred at 3500 yr B.P.Rabaul volcanism has been dominated by two main types: (a) basaltic and basaltic andesite cone building eruptions; and (b) dacitic, and rarely andesitic or rhyolitic, plinian/ignimbrite eruptions of both high- and low-aspect ratio types. The 1400 yr B.P. Rabaul Ignimbrite is a type example of a low-aspect ratio, high-energy, and potentially very damaging eruption. Fine vitric ash deposits, common in the Rabaul pyroclastic sequence, demonstrate the frequent modification of eruptions by external water probably related to early caldera lakes or bays. Interbedding of these fine ashes with plinian pumice lapilli beds suggests that many early eruptions occurred from multiple vents, located in both wet and dry areas.     

Separation of PP- and PS-wave reflected seismic data using two-dimensional finite offset common-reflection-surface traveltime approximation

Recently, the interest in PS-converted waves has increased for several applications, such as sub-basalt layer imaging, impedance estimates and amplitude-versus-offset analysis. In this study, we consider the problem of separation of PP- and PS-waves from pre-stacked multicomponent seismic data in two-dimensional isotropic medium. We aim to demonstrate that the finite-offset common-reflection-surface traveltime approximation is a good alternative for separating PP- and PS-converted waves in common-offset and common shot configurations by considering a two-dimensional isotropic medium. The five parameters of the finite-offset common-reflection-surface are firstly estimated through the inversion methodology called very fast simulated annealing, which estimates all parameters simultaneously. Next, the emergence angle, one of the inverted parameters, is used to build an analytical separation function of PP and PS reflection separation based on the wave polarization equations. Once the PP- and PS-converted waves were separated, the sections are stacked to increase the signal-to-noise ratio using the special curves derived from finite-offset common-reflection-surface approximation. We applied this methodology to a synthetic dataset from simple-layered to complex-structured media. The numerical results showed that the inverted parameters of the finite offset common-reflection-surface and the separation function yield good results for separating PP- and PS-converted waves in noisy common-offset and common shot gathers.     

T -distributed Random Fields: A Parametric Model for Heavy-tailedWell-log Data1

Histograms of observations from spatial phenomena are often found to be more heavy-tailed than Gaussian distributions, which makes the Gaussian random field model unsuited. A T-distributed random field model with heavy-tailed marginal probability density functions is defined. The model is a generalization of the familiar Student-T distribution, and it may be given a Bayesian interpretation. The increased variability appears cross-realizations, contrary to in-realizations, since all realizations are Gaussian-like with varying variance between realizations. The T-distributed random field model is analytically tractable and the conditional model is developed, which provides algorithms for conditional simulation and prediction, so-called T-kriging. The model compares favourably with most previously defined random field models. The Gaussian random field model appears as a special, limiting case of the T-distributed random field model. The model is particularly useful whenever multiple, sparsely sampled realizations of the random field are available, and is clearly favourable to the Gaussian model in this case. The properties of the T-distributed random field model is demonstrated on well log observations from the Gullfaks field in the North Sea. The predictions correspond to traditional kriging predictions, while the associated prediction variances are more representative, as they are layer specific and include uncertainty caused by using variance estimates.