Dynamic river–groundwater exchange in the presence of a saline,semi‐confined aquifer

Understanding groundwater–surface water exchange in river banks is crucial for effective water management and a range of scientific disciplines. While there has been much research on bank storage, many studies assume idealized aquifer systems. This paper presents a field‐based study of the Tambo Catchment (southeast Australia) where the Tambo River interacts with both an unconfined aquifer containing relatively young and fresh groundwater (<500 μS/cm and <100 years old) and a semi‐confined artesian aquifer containing old and saline groundwater (electrical conductivity > 2500 μS/cm and >10 000 years old). Continuous groundwater elevation and electrical conductivity monitoring within the different aquifers and the river suggest that the degree of mixing between the two aquifers and the river varies significantly in response to changing hydrological conditions. Numerical modelling using MODFLOW and the solute transport package MT3DMS indicates that saline water in the river bank moves away from the river during flooding as hydraulic gradients reverse. This water then returns during flood recession as baseflow hydraulic gradients are re‐established. Modelling also indicates that the concentration of a simulated conservative groundwater solute can increase for up to ~34 days at distances of 20 and 40 m from the river in response to flood events approximately 10 m in height. For the same flood event, simulated solute concentrations within 10 m of the river increase for only ~15 days as the infiltrating low‐salinity river water drives groundwater dilution. Average groundwater fluxes to the river stretch estimated using Darcy's law were 7 m³/m/day compared with 26 and 3 m³/m/day for the same periods via mass balance using Radon (²²²Rn) and chloride (Cl), respectively. The study shows that by coupling numerical modelling with continuous groundwater–surface water monitoring, the transient nature of bank storage can be evaluated, leading to a better understanding of the hydrological system and better interpretation of hydrochemical data. Copyright © 2015 John Wiley & Sons, Ltd.     

Chromian spinel as a petrogenetic indicator in abyssal and alpine-type peridotites and spatially associated lavas

The composition of chromian spinel in alpine-type peridotites has a large reciprocal range of Cr and Al, with increasing Cr# (Cr/(Cr+Al)) reflecting increasing degrees of partial melting in the mantle. Using spinel compositions, alpine-type peridotites can be divided into three groups. Type I peridotites and associated volcanic rocks contain spinels with Cr#<0.60; Type III peridotites and associated volcanics contain spinels with Cr#>0.60, and Type II peridotites and volcanics are a transitional group and contain spinels spanning the full range of spinel compositions in Type I and Type II peridotites. Spinels in abyssal peridotites lie entirely within the Type I spinel field, making ophiolites with Type I alpine-type peridotites the most likely candidates for sections of ocean lithosphere formed at a midocean ridge. The only modern analogs for Type III peridotites and associated volcanic rocks are found in arc-related volcanic and intrusive rocks, continental intrusive assemblages, and oceanic plateau basalts. We infer a sub-volcanic arc petrogenesis for most Type III alpine-type peridotites. Type II alpine-type peridotites apparently reflect composite origins, such as the formation of an island-arc on ocean crust, resulting in large variations in the degree and provenance of melting over relatively short distances. The essential difference between Type I and Type III peridotites appears to be the presence or absence of diopside in the residue at the end of melting.Based on an examination of co-existing rock and spinel compositions in lavas, it appears that spinel is a sensitive indicator of melt composition and pressure of crystallization. The close similarity of spinel composition fields in genetically related basalts, dunites and peridotites at localities in the oceans and in ophiolite complexes indicates that its composition reflects the degree of melting in the mantle source region. Accordingly, we infer from the restricted range of spinel compositions in abyssal basalts that the degree of mantle melting beneath mid-ocean ridges is generally limited to that found in Type I alpine-type peridotites. It is apparent, therefore, that the phase boundary OL-EN-DI-SP +meltOL-EN-SP+melt has limited the degree of melting of the mantle beneath mid-ocean ridges. This was clearly not the case for many alpine-type peridotites, implying very different melting conditions in the mantle, probably involving the presence of water.     

Determination of shear strength of conglomerates using a caterpillar D9 ripper and comparison with alternative methods

Summary A method of determining shear strength in Gila Conglomerates at Inspiration Mine (now Cyprus Miami Mine), Arizona, using Caterpillar D9H and D9L rippers is described. The results are compared with data obtained using laboratory direct shear tests — both large and small scale, back analysis and boulder modelling. The significance of scale, particularly relating to boulder size is discussed.     

Geochronology and geochemistry of Grenvillian igneous suites in the northern Oaxacan Complex, southern Mexico: tectonic implications

Chemical and U–Pb isotopic analyses of metaigneous rocks in the northern Oaxacan Complex in southern Mexico indicate that they form part of two granitic–gabbroic suites intruded at 1157–1130 and 1012 Ma, which were metamorphosed under granulite facies conditions between 1004 and 980 Ma. Although the older suite has both within-plate and arc geochemical signatures, the arc characteristics (enrichment of La and Ce relative to Nb, Ta, and Th) are inferred to result from crustal contamination, a conclusion consistent with their negative _Nd signatures. The younger suite is spatially associated with anorthosites (from which we were unable to acquire a protolith age), suggesting that collectively it forms part of anorthosite–mangerite–charnockite–granite (AMCG) suites. The tholeiitic nature of the mafic rocks along with the within-plate character of the felsic rocks suggests that they were intruded during extension related to either farfield backarc rifting, rifting above a slab window, or anorogenic intercontinental rifting. Potentially correlative AMCG suites are widespread in Mexico, the Grenville Province of eastern Canada and northeastern USA, and the Andean massifs of Colombia, however, Pb isotopic data most closely resemble those in South America. These data are consistent with published hypotheses that suggest Oaxaquia represents an exotic terrane derived from Amazonia.     

Development of ST:REAM: a reach‐based stream power balance approach for predicting alluvial river channel adjustment

River channel sediment dynamics are important in integrated catchment management because changes in channel morphology resulting from sediment transfer have important implications for many river functions. However, application of existing approaches that account for catchment‐scale sediment dynamics has been limited, largely due to the difficulty in obtaining data necessary to support them. It is within this context that this study develops a new, reach‐based, stream power balance approach for predicting river channel adjustment. The new approach, named ST:REAM (sediment transport: reach equilibrium assessment method), is based upon calculations of unit bed area stream power (ω) derived from remotely sensed slope, width and discharge datasets. ST:REAM applies a zonation algorithm to values of ω that are spaced every 50 m along the catchment network in order to divide the branches of the network up into relatively homogenous reaches. ST:REAM then compares each reach's ω value with the ω of its upstream neighbour in order to predict whether or not the reach is likely to be either erosion dominated or deposition dominated. The paper describes the application of ST:REAM to the River Taff in South Wales, UK. This test study demonstrated that ST:REAM can be rapidly applied using remotely sensed data that are available across many river catchments and that ST:REAM correctly predicted the status of 87.5% of sites within the Taff catchment that field observations had defined as being either erosion or deposition dominated. However, there are currently a number of factors that limit the usefulness of ST:REAM, including inconsistent performance and the need for additional, resource intensive, data to be collected to both calibrate the model and aid interpretation of its results. Copyright © 2014 John Wiley & Sons, Ltd.     

A determination of source properties of large intraplate earthquakes in Alaska

Historically, large and potentially hazardous earthquakes have occurred within the interior of Alaska. However, most have not been adequately studied using modern methods of waveform modeling. The 22 July 1937, 16 October 1947, and 7 April 1958 earthquakes are three of the largest events known to have occurred within central Alaska (M _s =7.3,M _s =7.2 andM _s =7.3, respectively). We analyzed teleseismic body waves to gain information about the focal parameters of these events. In order to deconvolve the source time functions from teleseismic records, we first attempted to improve upon the published focal mechanisms for each event. Synthetic seismograms were computed for different source parameters, using the reflectivity method. A search was completed which compared the hand-digitized data with a suite of synthetic traces covering the complete parameter space of strike, dip, and slip direction. In this way, the focal mechanism showing the maximum correlation between the observed and calculated traces was found. Source time functions, i.e., the moment release as a function of time, were then deconvolved from teleseismic records for the three historical earthquakes, using the focal mechanisms which best fit the data. From these deconvolutions, we also recovered the depth of the events and their seismic moments. The earthquakes were all found to have a shallow foci, with depths of less than 10 km.The 1937 earthquake occurred within a northeast-southwest band of seismicity termed the Salcha seismic zone (SSZ). We confirm the previously published focal mechanism, indicating strike-slip faulting, with one focal plane parallel to the SSZ which was interpreted as the fault plane. Assuming a unilateral fault model and a reasonable rupture velocity of between 2 and 3 km/s, the 21 second rupture duration for this event indicates that all of the 65 km long SSZ may have ruptured during this event. The 1947 event, located to the south of the northwest-southeast trending Fairbanks seismic zone, was found to have a duration of about 11 seconds, thus indicating a rupture length of up to 30 km. The rupture duration of the 1958 earthquake, which occurred near the town of Huslia, approximately 400 km ENE of Fairbanks, was found to be about 9 seconds. This gives a rupture length consistent with the observed damage, an area of 16 km by 64 km.     

Detection of the old stellar component of the major Galactic bar

We present near-infrared colour–magnitude diagrams and star counts for a number of regions along the Galactic plane. It is shown that along the l =27°, b =0° line of sight there is a feature at 5.7±0.7 kpc with a density of stars at least a factor of 2 and probably more than a factor of 5 times that of the disc at the same position. This feature forms a distinct clump on an H versus J − H diagram and is seen at all longitudes from the bulge to about l =28°, but at no longitude greater than this. The distance to the feature at l =20° is about 0.5 kpc further than at l =27°, and by l =10° it has merged with, or has become, the bulge. Given that at l =27° and l =21° there is also a clustering of very young stars, the only component that can reasonably explain what is seen is a bar with half-length of around 4 kpc and a position angle of about 43°±7°.     

Improved seismic hazard model with application to probabilistic seismic demand analysis

An improved seismic hazard model for use in performance‐based earthquake engineering is presented. The model is an improved approximation from the so‐called ‘power law’ model, which is linear in log–log space. The mathematics of the model and uncertainty incorporation is briefly discussed. Various means of fitting the approximation to hazard data derived from probabilistic seismic hazard analysis are discussed, including the limitations of the model. Based on these ‘exact’ hazard data for major centres in New Zealand, the parameters for the proposed model are calibrated. To illustrate the significance of the proposed model, a performance‐based assessment is conducted on a typical bridge, via probabilistic seismic demand analysis. The new hazard model is compared to the current power law relationship to illustrate its effects on the risk assessment. The propagation of epistemic uncertainty in the seismic hazard is also considered. To allow further use of the model in conceptual calculations, a semi‐analytical method is proposed to calculate the demand hazard in closed form. For the case study shown, the resulting semi‐analytical closed form solution is shown to be significantly more accurate than the analytical closed‐form solution using the power law hazard model, capturing the ‘exact’ numerical integration solution to within 7% accuracy over the entire range of exceedance rate. Copyright © 2007 John Wiley & Sons, Ltd.