Exploring weathering and regolith transport controls on Critical Zone development with models and natural experiments

The architecture of the Critical Zone, including mobile regolith thickness and depth to the weathering front, is first order controlled by advance of a weathering front at depth and transport of sediment at the surface. Differences in conditions imposed by slope aspect in the Gordon Gulch catchment of the Boulder Creek Critical Zone Observatory present a natural experiment to explore these interactions. The weathering front is deeper and saprolite more decayed on north-facing than on south-facing slopes. Simple numerical models of weathering front advance, mobile regolith production, and regolith transport are used to test how weathering and erosion rates interact in the evolution of weathered profiles. As the processes which attempt are being made to mimic are directly tied to climate variables such as mean annual temperature, the role of Quaternary climate variation in governing the evolution of Critical Zone architecture can be explored with greater confidence.     

The Regional Ocean Modeling System (ROMS) 4-dimensional variational data assimilation systems: Part II - Performance and application to the California Current System

The Regional Ocean Modeling System (ROMS) 4-dimensional variational (4D-Var) data assimilation systems have been systematically applied to the mesoscale circulation environment of the California Current to demonstrate the performance and practical utility of the various components of ROMS 4D-Var. In particular, we present a comparison of three approaches to 4D-Var, namely: the primal formulation of the incremental strong constraint approach; the dual formulation “physical-space statistical analysis system”; and the dual formulation indirect representer approach. In agreement with theoretical considerations all three approaches converge to the same ocean circulation estimate when using the same observations and prior information. However, the rate of convergence of the dual formulation was found to be inferior to that of the primal formulation. Other aspects of the 4D-Var performance that relate to the use of multiple outer-loops, preconditioning, and the weak constraint are also explored. A systematic evaluation of the impact of the various components of the 4D-Var control vector (i.e. the initial conditions, surface forcing and open boundary conditions) is also presented. It is shown that correcting for uncertainties in the model initial conditions exerts the largest influence on the ability of the model to fit the available observations. Various important diagnostics of 4D-Var are also examined, including estimates of the posterior error, the information content of the observation array, and innovation-based consistency checks on the prior error assumptions. Using these diagnostic tools, we find that more than 90% of the observations assimilated into the model provide redundant information. This is a symptom of the large percentage of satellite data that are used and to some extent the nature of the data processing employed. This is the second in a series of three papers describing the ROMS 4D-Var systems.     

On the use of IPCC-class models to assess the impact of climate on Living Marine Resources

The study of climate impacts on Living Marine Resources (LMRs) has increased rapidly in recent years with the availability of climate model simulations contributed to the assessment reports of the Intergovernmental Panel on Climate Change (IPCC). Collaboration between climate and LMR scientists and shared understanding of critical challenges for such applications are essential for developing robust projections of climate impacts on LMRs. This paper assesses present approaches for generating projections of climate impacts on LMRs using IPCC-class climate models, recommends practices that should be followed for these applications, and identifies priority developments that could improve current projections. Understanding of the climate system and its representation within climate models has progressed to a point where many climate model outputs can now be used effectively to make LMR projections. However, uncertainty in climate model projections (particularly biases and inter-model spread at regional to local scales), coarse climate model resolution, and the uncertainty and potential complexity of the mechanisms underlying the response of LMRs to climate limit the robustness and precision of LMR projections. A variety of techniques including the analysis of multi-model ensembles, bias corrections, and statistical and dynamical downscaling can ameliorate some limitations, though the assumptions underlying these approaches and the sensitivity of results to their application must be assessed for each application. Developments in LMR science that could improve current projections of climate impacts on LMRs include improved understanding of the multi-scale mechanisms that link climate and LMRs and better representations of these mechanisms within more holistic LMR models. These developments require a strong baseline of field and laboratory observations including long time series and measurements over the broad range of spatial and temporal scales over which LMRs and climate interact. Priority developments for IPCC-class climate models include improved model accuracy (particularly at regional and local scales), inter-annual to decadal-scale predictions, and the continued development of earth system models capable of simulating the evolution of both the physical climate system and biosphere. Efforts to address these issues should occur in parallel and be informed by the continued application of existing climate and LMR models.     

Radon emanation of heterogeneous basin deposits in Kathmandu Valley,Nepal

Effective radium-226 concentration (EC_Ra) has been measured in soil samples from seven horizontal and vertical profiles of terrace scarps in the northern part of Kathmandu Valley, Nepal. The samples belong to the Thimi, Gokarna, and Tokha Formations, dated from 50 to 14 ky BP, and represent a diverse fluvio-deltaic sedimentary facies mainly consisting of gravelly to coarse sands, black, orange and brown clays. EC_Ra was measured in the laboratory by radon-222 emanation. The samples (n = 177) are placed in air-tight glass containers, from which, after an accumulation time varying from 3 to 18 days, the concentration of radon-222, radioactive decay product of radium-226 and radioactive gas with a half-life of 3.8 days, is measured using scintillation flasks. The EC_Ra values from the seven different profiles of the terrace deposits vary from 0.4 to 43 Bq kg^?1, with profile averages ranging from 12 ± 1 to 27 ± 2 Bq kg^?1. The values have a remarkable consistency along a particular horizon of sediment layers, clearly demonstrating that these values can be used for long distance correlations of the sediment horizons. Widely separated sediment profiles, representing similar stratigraphic positions, exhibit consistent EC_Ra values in corresponding stratigraphic sediment layers. EC_Ra measurements therefore appear particularly useful for lithologic and stratigraphic discriminations. For comparison, EC_Ra values of soils from different localities having various sources of origin were also obtained: 9.2 ± 0.4 Bq kg^?1 in soils of Syabru–Bensi (Central Nepal), 23 ± 1 Bq kg^?1 in red residual soils of the Bhattar-Trisuli Bazar terrace (North of Kathmandu), 17.1 ± 0.3 Bq kg^?1 in red residual soils of terrace of Kalikasthan (North of Trisuli Bazar) and 10 ± 1 Bq kg^?1 in red residual soils of a site near Nagarkot (East of Kathmandu). The knowledge of EC_Ra values for these various soils is important for modelling radon exhalation at the ground surface, in particular in the vicinity of active faults. Importantly, the study also reveals that, above numerous sediments of Kathmandu Valley, radon concentration in dwellings can potentially exceed the level of 300 Bq m^?3 for residential areas; a fact that should be seriously taken into account by the governmental and non-governmental agencies as well as building authorities.     

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.     

Anthropogenic influences on shoreline and nearshore evolution in the San Francisco Bay coastal system

Analysis of four historical bathymetric surveys over a 132-year period has revealed significant changes to the morphology of the San Francisco Bar, an ebb-tidal delta at the mouth of San Francisco Bay estuary. From 1873 to 2005 the San Francisco Bar vertically-eroded an average of 80 cm over a 125 km² area, which equates to a total volume loss of 100 ± 52 million m³ of fine- to coarse-grained sand. Comparison of the surveys indicates the entire ebb-tidal delta contracted radially, with the crest moving landward an average of 1 km. Long-term erosion of the ebb-tidal delta is hypothesized to be due to a reduction in the tidal prism of San Francisco Bay and a decrease in coastal sediment supply, both as a result of anthropogenic activities. Prior research indicates that the tidal prism of the estuary was reduced by 9% from filling, diking, and sedimentation. Compilation of historical records dating back to 1900 reveals that a minimum of 200 million m³ of sediment has been permanently removed from the San Francisco Bay coastal system through dredging, aggregate mining, and borrow pit mining. Of this total, ∼54 million m³ of sand-sized or coarser sediment was removed from central San Francisco Bay. With grain sizes comparable to the ebb-tidal delta, and its direct connection to the bay mouth, removal of sediments from central San Francisco Bay may limit the sand supply to the delta and open coast beaches.     

An amalgamated meter-thick sedimentary package enabled by the 2011 Tohoku tsunami in El Garrapatero,Galapagos Islands

Tsunamis and storms instigate sedimentological and geomorphological changes to the coastal system, both long-term and ephemeral. To accurately predict future coastal hazards, one must identify the records that are generated by the processes associated with these hazards and recognize what will be preserved. Using eyewitness accounts, photographs, and sedimentology, this study documents pre- and post-tsunami conditions and constrains the timing and process of depositional events during and following the 11 March 2011 Tohoku tsunami in the coastal system at El Garrapatero, Galapagos Islands. While the tsunami acted as both an erosional and depositional agent, the thick, fan-like sand sheet in El Garrapatero was primarily emplaced by overwash deposition during high tide from swell waves occurring between 19–25 March and 17–22 April 2011. The swell waves were only able to access the terrestrial coastal system via a channel carved by the 2011 Tohoku tsunami through the barrier sand dune. This combined deposit could result in an overestimation of the hazard if interpreted to be the result of only one event (either tsunami or wind-generated waves). An analogous sand layer, younger than 1390–1530 cal yr BP, may record a similar, prior event.     

Environmental Determinants of Emergent Macrophyte Vegetation in Pacific Northwest Estuarine Tidal Wetlands

We investigated whether within wetland environmental conditions or surrounding land cover measured at multiple scales were more influential in structuring regional vegetation patterns in estuarine tidal wetlands in the Pacific Northwest, USA. Surrounding land cover was characterized at the 100, 250, and 1,000 m, and watershed buffer scales. Vegetation communities were characterized by high species richness, lack of monotypic zonation, and paucity of invasive species. The number of species per site ranged between 4 and 20 (mean?±?standard deviation?=?10.2?±?3.1). Sites supported a high richness (mean richness of native species 8.7?±?2.8) and abundance of native macrophytes (mean relative abundance 85 %?±?19 %). Vegetation assemblages were dominated by a mix of grasses, sedges, and herbs with Sarcocornia pacifica and Distichlis spicata being common at sites in the oceanic zone of the estuary and Carex lyngbyei and Agrostis stolonifera being common at the fresher sites throughout the study area. The vegetation community was most strongly correlated with salinity and land cover within close proximity to the study site and less so with land cover variables at the watershed scale. Total species richness and richness of native species were negatively correlated with the amount of wetland in the buffer at all scales, while abundance of invasive species was significantly correlated to within wetland factors, including salinity and dissolved phosphorus concentrations. Landscape factors related to anthropogenic disturbances were only important at the 100-m buffer scale, with anthropogenic disturbances further from the wetland not being influential in shaping the vegetation assemblage. Our research suggests that the traditional paradigms of tidal wetland vegetation structure and environmental determinants developed in east coast US tidal wetlands might not hold true for Pacific Northwest wetlands due to their unique chemical and physical factors, necessitating further detailed study of these systems.