Deriving fault-slip histories to test for secular variation in slip,with examples from the Kunlun and Awatere faults

Although offset and age data from displaced landforms are essential for identifying earthquake clusters and thus testing whether faults slip at uniform or secularly varying rates, it is not clear how the uncertainties in such measurements should be propagated so as to yield a robust fault-slip history (i.e., record of fault displacement over time). Here we develop a Monte Carlo approach for estimating the distribution of geologically reasonable fault-slip histories that fit the offset and age data from a population of dated and displaced landforms. The model assumes that the landforms share common faulting histories, the offset and age constraints are correct, and the fault has not reversed shear sense. Analysis of the model results yields both a precise average slip rate, in the case where a linear fit is applied to the data, and a best-fit fault-slip history, in the case where the linear constraint is removed. The method can be used to test for secular variation in slip because the uncertainty on this best-fit history is quantified. By applying the method to previously published morphochronologic data from faulted late Quaternary terrace risers along the Kunlun fault in China and the Awatere fault in New Zealand, we have assessed the extent to which our modeled average slip rates match previously reported values and the data support previous interpretations of uniform slip rate. The Kunlun data set yields average slip rates of 8.7 + 3.6/?2.1 mm/yr and 5.1 + 1.6/?1.2 mm/yr (68.27% confidence), for the central and eastern reaches of the fault, respectively, both of which match previously published slip rates. Our analysis further indicates that these fault reaches have both slipped uniformly over the latest Quaternary. In contrast, analysis of data from the Saxton River site along the Awatere fault reveals a mid-Holocene deceleration in slip rate from 6.2 + 1.6/?1.4 mm/yr to 2.8 + 1.0/?0.6 mm/yr. This result contradicts previous interpretations of uniform slip along the Awatere fault. The Monte Carlo method we present here for quantifying fault-slip histories using the offset and age data from a population of faulted landforms provides an important tool for distinguishing temporally uniform from secularly varying fault slip.     

New pterosaurian remains from the Late Cretaceous (Santonian) of Hungary (Iharkút, Csehbánya Formation)

The Late Cretaceous (Santonian) continental vertebrate locality of Iharkút, western Hungary has provided numerous azhdarchid pterosaur remains including the recently described, Bakonydraco galaczi. Since the first report of these fossils, additional remains have been discovered that improve considerably our knowledge of some aspects of the anatomy of azhdarchid pterosaurs. New cranial material described here indicates an edentulous, non-crested premaxillary rostrum in Bakonydraco similar to that of Quetzalcoatlus and reveals that this rostrum was considerably thinner and more lightly built than the relatively massive, pointed mandibular symphysis. In addition, the contact surface of the upper and lower jaws of Bakonydraco at least in the symphyseal region was more similar to Tapejara wellnhoferi than to other azdarchids: the contact is irregular and the premaxilla does not fit closely the rostral part of the mandible. Among the postcranial material the atlas-axis complex possessing lateral pneumatic foramina is of great importance because it further supports the notion of interspecific variability of vertebral pneumaticity among pterosaurs.     

Coupled discrete element and finite volume solution of two classical soil mechanics problems

One dimensional solutions for the classic critical upward seepage gradient/quick condition and the time rate of consolidation problems are obtained using coupled routines for the finite volume method (FVM) and discrete element method (DEM), and the results compared with the analytical solutions. The two phase flow in a system composed of fluid and solid is simulated with the fluid phase modeled by solving the averaged Navier–Stokes equation using the FVM and the solid phase is modeled using the DEM. A framework is described for the coupling of two open source computer codes: YADE-OpenDEM for the discrete element method and OpenFOAM for the computational fluid dynamics. The particle–fluid interaction is quantified using a semi-empirical relationship proposed by Ergun [12]. The two classical verification problems are used to explore issues encountered when using coupled flow DEM codes, namely, the appropriate time step size for both the fluid and mechanical solution processes, the choice of the viscous damping coefficient, and the number of solid particles per finite fluid volume.     

Nutrient Inputs, Phytoplankton Response, and CO2 Variations in a Semi-Enclosed Subtropical Embayment, Kaneohe Bay, Hawaii

The marine shelf areas in subtropical and tropical regions represent only 35% of the total shelf areas globally, but receive a disproportionately large amount of water (65%) and sediment (58%) discharges that enter such environments. Small rivers and/or streams that drain the mountainous areas in these climatic zones deliver the majority of the sediment and nutrient inputs to these narrow shelf environments; such inputs often occur as discrete, episodic introductions associated with storm events. To gain insight into the linked biogeochemical behavior of subtropical/tropical mountainous watershed-coastal ocean ecosystems, this work describes the use of a buoy system to monitor autonomously water quality responses to land-derived nutrient inputs and physical forcing associated with local storm events in the coastal ocean of southern Kaneohe Bay, Oahu, Hawaii, USA. The data represent 2.5 years of near-real time observations at a fixed station, collected concurrently with spatially distributed synoptic sampling over larger sections of Kaneohe Bay. Storm events cause most of the fluvial nutrient, particulate, and dissolved organic carbon inputs to Kaneohe Bay. Nutrient loadings from direct rainfall and/or terrestrial runoff produce an immediate increase in the N:P ratio of bay waters up to values of 48 and drive phytoplankton biomass growth. Rapid uptake of such nutrient subsidies by phytoplankton causes rapid declines of N levels, return to N-limited conditions, and subsequent decline of phytoplankton biomass over timescales ranging from a few days to several weeks, depending on conditions and proximity to the sources of runoff. The enhanced productivity may promote the drawing down of pCO₂ and lowering of surface water column carbonate saturation states, and in some events, a temporary shift from N to P limitation. The productivity-driven CO₂ drawdown may temporarily lead to air-to-sea transfer of atmospheric CO₂ in a system that is on an annual basis a source of CO₂ to the atmosphere due to calcification and perhaps heterotrophy. Storms may also strongly affect proximal coastal zone pCO₂ and hence carbonate saturation state due to river runoff flushing out high pCO₂ soil and ground waters. Mixing of the CO₂-charged water with seawater causes a salting out effect that releases CO₂ to the atmosphere. Many subtropical and tropical systems throughout the Pacific region are similar to Kaneohe Bay, and our work provides an important indication of the variability and range of CO₂ dynamics that are likely to exist elsewhere. Such variability must be taken into account in any analysis of the direction and magnitude of the air?Csea CO₂ exchange for the integrated coastal ocean, proximal and distal. It cannot be overemphasized that this research illustrates several examples of how high frequency sampling by a moored autonomous system can provide details about ecosystem responses to stochastic atmospheric forcing that are commonly missed by traditional synoptic observational approaches. Finally, the work exemplifies the utility of combining synoptic sampling and real-time autonomous observations to elucidate the biogeochemical and physical responses of coastal subtropical/tropical coral reef ecosystems to climatic perturbations.     

Natural and Post-Urbanization Signatures of Hypoxia in Two Basins of Puget Sound: Historical Reconstruction of Redox Sensitive Metals and Organic Matter Inputs

Hypoxia has been observed in Hood Canal, Puget Sound, WA, USA since the 1970s. Four long sediment cores were collected in 2005 and age-dated to resolve natural and post-urbanization signatures of hypoxia and organic matter (OM) sources in two contrasting basins of Puget Sound: Main Basin and Hood Canal. Paleoecological indicators used for sediment reconstructions included pollen, stable carbon and nitrogen isotopes (??¹³C and ??¹⁵N), biomarkers of terrestrial OM (TOM), biogenic silica (BSi), and redox-sensitive metals (RSM). The sedimentary reconstructions illustrated a gradient in RSM enrichment factors as Hood Canal > Main Basin, southern > northern cores, and pre-1900s > 1900?C2005. The urbanization of Puget Sound watersheds during the 1900s was reflected as shifts in all the paleoecological signatures. Pollen distributions shifted from predominantly old growth conifer to successional alder, dominant OM signatures recorded a decrease in the proportion of marine OM (MOM) concomitant with an increase in the proportion of TOM, and the weight % of BSi decreased. However, these shifts were not coincidental with an overall increase in the enrichment of RSM or ??¹⁵N signatures indicative of cultural eutrophication. The increased percentage of TOM was independently verified by both the elemental ratios and lignin yields. In addition, isotopic signatures, BSi, and RSMs all suggest that OM shifts may be due to a reduction in primary productivity rather than an increase in OM regeneration in the water column or at the sediment/water interface. Therefore, the reconstructions suggested the Hood Canal has been under a more oxygenated ??stance?? during the twentieth century compared to prior periods. However, these 2005 cores and their resolutions do not encompass the period of high resolution water column measurements that showed short-lived hypoxia events and fish kills in Hood Canal during the early twenty-first century. The decoupling between the increased watershed-scale anthropogenic alterations recorded in the OM signatures and the relatively depleted RSM during the twentieth century suggests that physical processes, such as deep-water ventilation, may be responsible for the historical variation in oxygen levels. Specifically, climate oscillations may influence the ventilation and/or productivity of deep water in Puget Sound and particularly their least mixed regions.     

A Leaky-Conduit Model of Transient Flow in Karstic Aquifers

Karst Flow Model (KFM) simulates transient flow in an unconfined karstic aquifer having a well-developed conduit system. KFM treats the springshed as a two-dimensional porous matrix containing a triangulated irregular network of leaky conduits. The number and location of conduits can be specified arbitrarily, perhaps using field information as a guide, or generated automatically. Conduit networks can be tree-like or braided. Rainwater that has infiltrated down from the surface leaks into the conduits from the adjacent porous matrix at a rate dictated by Darcy’s law, then flows turbulently to the spring via the conduits. KFM is calibrated using the known steady state; geometry and recharge determine the steady fluxes in the conduits, and the head distribution determines conduit gradients and sizes. Spring flow can vary with time due to spatially and temporally variable recharge and due to prescribed variations in the elevation of the spring. KFM is illustrated by four examples run on a test aquifer consisting of 27 nodes, 42 elements, and 26 conduits. Three examples (drought, uniform rainstorm, storm-water input to one element) are simulations, while the fourth uses data from a spring-basin flooding event. The qualitative fit between the predicted and observed spring discharge in the fourth example provides support of the hypothesis that the dynamic behavior of a karst conduit system is an emergent property of a self-organized system, largely independent of the locations and properties of individual conduits.     

Petrogenesis of Sierra Nevada plutons inferred from the Sr, Nd, and O isotopic signatures of mafic igneous complexes in Yosemite Valley, California

Mafic complexes in the central Sierra Nevada batholith record valuable geochemical information regarding the role mafic magmas play in arc magmatism and the generation of continental crust. In the intrusive suite of Yosemite Valley, major and trace element compositions of the hornblende-bearing gabbroic rocks from the Rockslides mafic complex and of the mafic dikes in the North America Wall are compositionally similar to high-alumina basalt. Of these rocks, two samples have higher Ni and Cr abundances as well as higher ε_Nd values than previously recognized for the intrusive suite. Plagioclase crystals in rocks from the North America Wall and the Rockslides have prominent calcic cores and sharply defined sodic rims, a texture commonly associated with mixing of mafic and felsic magmas. In situ analyses of ⁸⁷Sr/⁸⁶Sr in plagioclase show no significant isotopic difference from the cores to the rims of these grains. We propose that the high ⁸⁷Sr/⁸⁶Sr (~0.7067) and low ε_Nd (~?3.4) of bulk rocks, the homogeneity of ⁸⁷Sr/⁸⁶Sr in plagioclase, and the high δ¹⁸O values of bulk rocks (6.6–7.3 ‰) and zircon (Lackey et al. in J Petrol 49:1397–1426, 2008) demonstrate that continental crust was assimilated into the sublithospheric mantle-derived basaltic precursors of the mafic rocks in Yosemite Valley. Contamination (20–40 %) likely occurred in the lower crust as the magma differentiated to high-alumina basalt prior to plagioclase (and zircon) crystallization. As a consequence, the isotopic signatures recorded by whole rocks, plagioclase, and zircon do not represent the composition of the underlying lithospheric mantle. We conclude that the mafic and associated felsic members of the intrusive suite of Yosemite Valley represent 60–80 % new additions to the crust and include significant quantities of recycled ancient crust.