The Role of Sediment Dynamics for Inorganic Accretion Patterns in Southern California’s Mediterranean-Climate Salt Marshes

Salt marsh resilience to sea-level rise depends on marsh plain elevation, tidal range, subsurface processes, as well as surface accretion, of which suspended-sediment concentration (SSC) is a critical component. However, spatial and temporal patterns of inorganic sedimentation are poorly quantified within and across Salicornia pacifica (pickleweed)-dominated marshes. We compared vertical accretion rates and re-examined previously published suspended-sediment patterns during dry-weather periods at Seal Beach Wetlands, which is characterized by a mix of Spartina foliosa (cordgrass) and pickleweed, and for Mugu Lagoon, where cordgrass is rare. Mugu Lagoon occurs higher in the tidal frame and receives terrigenous sediment from an adjacent creek. Feldspar marker horizons were established in winter 2013–2014 to measure accretion. Accretion rates at Seal Beach Wetlands and Mugu Lagoon were 6 ± 0.5 mm/year (mean ± SE) and 2 ± 0.3 mm/year. Also, the estimated sediment flux (g/year) across the random feldspar plots was 3.5 times higher at Seal Beach Wetlands. At Mugu Lagoon, accretion was higher near creeks, although not statistically significant. Dry-weather SSC showed similar concentrations at transect locations across sites. During wet weather, however, SSC at Mugu Lagoon increased at all locations, with concentrations decaying farther than 8 m from tidal creek edge. Based on these results from Mugu Lagoon, we conclude accretion patterns are set by infrequent large flooding events in systems where there is a watershed sediment source. Higher accretion rates at Seal Beach Wetlands may be linked to lower-marsh elevations, and thus more frequent inundation, compared with Mugu Lagoon.     

Selection of Momentum Variables for a Three-Dimensional Variational Analysis

Three choices of control variables for meteorological variational analysis (3DVAR or 4DVAR) are associated with horizontal wind: (1) streamfunction and velocity potential, (2) eastward and northward velocity, and (3) vorticity and divergence. This study shows theoretical and numerical differences of these variables in practical 3DVAR data assimilation through statistical analysis and numerical experiments. This paper demonstrates that (a) streamfunction and velocity potential could potentially introduce analysis errors; (b) A 3DVAR using velocity or vorticity and divergence provides a natural scale dependent influence radius in addition to the covariance; (c) for a regional analysis, streamfunction and velocity potential are retrieved from the background velocity field with Neumann boundary condition. Improper boundary conditions could result in further analysis errors; (d) a variational data assimilation or an inverse problem using derivatives as control variables yields smoother analyses, for example, a 3DVAR using vorticity and divergence as controls yields smoother wind analyses than those analyses obtained by a 3DVAR using either velocity or streamfunction/velocity potential as control variables; and (e) statistical errors of higher order derivatives of variables are more independent, e.g., the statistical correlation between U and V is smaller than the one between streamfunction and velocity potential, and thus the variables in higher derivatives are more appropriate for a variational system when a cross-correlation between variables is neglected for efficiency or other reasons. In summary, eastward and northward velocity, or vorticity and divergence are preferable control variables for variational systems and the former is more attractive because of its numerical efficiency. Numerical experiments are presented using analytic functions and real atmospheric observations.     

The practicalities of using CFCs and SF6 for groundwater dating and tracing

Knowledge of groundwater residence time is important in understanding key issues in the evolution of water quality, whether this occurs due to water–rock interaction or simply by mixing or contamination. The build-up in the atmosphere of the trace gases chlorofluorocarbons (CFCs) and sulphur hexafluoride (SF₆) from the middle of the last century offers a convenient way of dating waters up to ∼60 a old. The gases are well-mixed in the atmosphere so their input functions are not area-specific as is the case with ³H. While any one of these trace gases can in principle provide a groundwater age, when two or more are measured on water samples the potential exists to distinguish between different modes of flow including piston flow, exponential flow and simple end-member mixing. As with all groundwater dating methods, caveats apply. Factors such as recharge temperature and elevation must be reasonably well-constrained. Primarily for SF₆, the phenomenon of ‘excess air’ also requires consideration. Primarily for the CFCs, local sources of contamination need to be considered, as do redox conditions. For both SF₆ and the CFCs, the nature and thickness of the unsaturated zone need to be factored into residence time calculations. However, as an inexpensive dating method, the trace gases can be applied to a wide range of groundwater problems where traditional age indicators might once have been used more sparingly. Examples include tracing flowlines, detecting small modern inputs in ‘old’ waters, and pollution risk assessment. In the future, with the main CFCs already declining in the atmosphere, new anthropogenic trace gases are likely to take their place.     

Phytoplankton abundance and contributions to suspended particulate matter in the Ohio,Upper Mississippi and Missouri Rivers

Main channel habitats of the Ohio, Missouri, and Upper Mississippi Rivers were surveyed during the summers of 2004, 2005 and 2006 using a probability-based sampling design to characterize inter-annual and inter-river variation in suspended chlorophyll a (CHLa) and related variables. Large (fivefold) differences in CHLa were observed with highest concentrations in the Upper Mississippi (32.3 ± 1.8 μg L⁻¹), intermediate values in the Missouri (19.7 ± 1.1 μg L⁻¹) and lowest concentrations in the Ohio (6.8 ± 0.5 μg L⁻¹). Inter-annual variation was small in comparison to inter-river differences suggesting that basin-specific factors exert greater control over river-wide CHLa than regional-scale processes influencing climate and discharge. The rivers were characterized by variable but generally low light conditions as indicated by depth-averaged underwater irradiance <4 E m⁻² day⁻¹ and high ratios of channel depth to euphotic depth (>3). Despite poor light conditions, regression analyses revealed that TP was the best single predictor of CHLa (R ² = 0.40), though models incorporating both light and TP performed better (R ² = 0.60). Light and nutrient conditions varied widely within rivers and were inversely related, suggesting that riverine phytoplankton may experience shifts in resource limitation during transport. Inferred grazing and sedimentation losses were large yet CHLa concentrations did not decline downriver indicating that growth and loss processes were closely coupled. The contribution by algae to suspended particulate organic matter in these rivers (mean = 41%) was similar to that of lakes (39%) but lower relative to reservoirs (61%).     

Lattice distortion in a zircon population and its effects on trace element mobility and U–Th–Pb isotope systematics: examples from the Lewisian Gneiss Complex, northwest Scotland

Zircon is a key mineral in geochemical and geochronological studies in a range of geological settings as it is mechanically and chemically robust. However, distortion of its crystal lattice can facilitate enhanced diffusion of key elements such as U and Pb. Electron backscatter diffraction (EBSD) analysis of ninety-nine zircons from the Lewisian Gneiss Complex (LGC) of northwest Scotland has revealed five zircons with lattice distortion. The distortion can take the form of gradual bending of the lattice or division of the crystal into subgrains. Zircon lattices are distorted because of either post-crystallisation plastic distortion or growth defects. Three of the five distorted zircons, along with many of the undistorted zircons in the population, were analysed by ion microprobe to measure U and Pb isotopes, Ti and REEs. Comparison of Th/U ratio, ²⁰⁷Pb/²⁰⁶Pb age, REE profile and Ti concentration between zircons with and without lattice distortion suggests that the distortion is variably affecting the concentration of these trace elements and isotopes within single crystals, within samples and between localities. REE patterns vary heterogeneously, sometimes relatively depleted in heavy REEs or lacking a Eu anomaly. Ti-in-zircon thermometry records temperatures that were either low (~700 °C) or high (>900 °C) relative to undistorted zircons. One distorted zircon records apparent ²⁰⁷Pb/²⁰⁶Pb isotopic ages (?3.0 to +0.3 % discordance) in the range of ~2,420–2,450 Ma but this does not correlate with any previously dated tectonothermal event in the LGC. Two other distorted zircons give discordant ages of 2,331 ± 22 and 2,266 ± 40 Ma, defining a discordia lower intercept within error of a late amphibolite-facies tectonothermal event. This illustrates that Pb may be mobilised in distorted zircons at lower metamorphic grade than in undistorted zircons. These differences in trace element abundances and isotope systematics in distorted zircons relative to undistorted zircons are generally interpreted to have been facilitated by subgrain walls. Trace elements and isotopes would have moved from undistorted lattice into these subgrain walls as their chemical potential is modified due to the presence of the dislocations which make up the subgrain wall. Subgrain walls provided pathways for chemical exchange between crystal and surroundings. Only five per cent of zircons in this population have lattice distortion suggesting it will not have a major impact on zircon geochronology studies, particularly as three of the five distorted zircons are from strongly deformed rocks not normally sampled in such studies. However, this does suggest there may be a case for EBSD analysis of zircons prior to geochemical analysis when zircons from highly deformed rocks are to be investigated.     

Paleolimnological evidence of the response of the central Canadian treeline zone to radiative forcing and hemispheric patterns of temperature change over the past 2000 years

Instrumental climate records from the central Canadian treeline zone display a pattern of variation similar to general Northern Hemisphere temperature trends. To examine whether this general correspondence extends back beyond the instrumental record, we obtained a sediment core from Lake S41, a small lake in the Northwest Territories of Canada at 63°43.11′ N, 109°19.07′ W. A radiocarbon-based chronology was developed for the core. The sediments were analyzed for organic-matter content by loss-on-ignition (LOI), biogenic-silica content (BSi), and chironomid community composition to reconstruct July air temperature and summer water temperature. The paleolimnological records were compared with records of atmospheric CO₂ concentration, solar variability, and hemispheric temperature variations over the past 2000 years. The results of the analyses suggest that widely-documented long-term variations in Northern Hemisphere temperature associated with radiative forcing, namely the cooling following the medieval period during the Little Ice Age (LIA), and twentieth century warming, are represented in the central Canadian treeline zone. There is also evidence of a brief episode of warming during the eighteenth century. As evidenced by LOI and BSi, the twentieth century warming is typified by increased lake productivity relative to the LIA. Depending upon the measure, the increased productivity of the twentieth century nearly equals or exceeds that of any other period in the past 2000 years. In contrast, the rate of chironomid head capsule accumulation decreased and remained low during the twentieth century. Although the chironomid-inferred temperature reconstructions indicate cooling during the LIA, they present no evidence of greatly increased temperatures during the twentieth century. Warming during the twentieth century might have enhanced lake stratification, and the response of the chironomid fauna to warming was attenuated by decreased oxygen and lower temperatures in the hypolimnion of the more stratification-prone lake.

Holocene lake sediment records of Arctic hydrology

Although paleoclimatic research in the Arctic has most often focused on variations in temperature, the Arctic has also experienced changes in hydrologic balance. Changes in Arctic precipitation and evaporation rates affects soils, permafrost, lakes, wetlands, rivers, ice and vegetation. Changes in Arctic soils, permafrost, runoff, and vegetation can influence global climate by changing atmospheric methane and carbon dioxide concentrations, thermohaline circulation, and high latitude albedo. Documenting past variations in Arctic hydrological conditions is important for understanding Arctic climate and the potential response and role of the Arctic in regards to future climate change. Methods for reconstructing past changes in Arctic hydrology from the stratigraphic, isotopic, geochemical and fossil records of lake sediments are being developed, refined and applied in a number of regions. These records suggest that hydrological variations in the Arctic have been regionally asynchronous, reflecting the impacts of different forcing factors including orbitally controlled insolation changes, changes in geography related to coastal emergence, ocean currents, sea ice extent, and atmospheric circulation. Despite considerable progress, much work remains to be done on the development of paleohydrological proxies and their application to the Arctic.     

Metasomatism and the crystallization of zircon megacrysts in Archaean peridotites from the Lewisian complex,NW Scotland

Zircon megacrysts are locally abundant in 1–40 cm-thick orthopyroxenite veins within peridotite host rocks in the Archaean Lewisian gneiss complex from NW Scotland. The veins formed by metasomatic interaction between the ultramafic host and Si-rich melts are derived from partial melting of the adjacent granulite-facies orthogneisses. The interaction produced abundant orthopyroxene and, within the thicker veins, phlogopite, pargasite and feldspathic bearing assemblages. Two generations of zircon are present with up to 1 cm megacrystic zircon and a later smaller equant population located around the megacryst margins. Patterns of zoning, rare earth element abundance and oxygen isotopic compositions indicate that the megacrysts crystallized from crustal melts, whereas the equant zircon represents new neocryst growth and partial replacement of the megacryst zircon within the ultramafic host. Both zircon types have U–Pb ages of ca. 2464 Ma, broadly contemporaneous with granulite-facies events in the adjacent gneisses. Zircon megacrysts locally form?>?10% of the assemblage and may be associated to zones of localized nucleation or physically concentrated during movement of the siliceous melts. Their unusual size is linked to the suppression of zircon nucleation and increased Zr solubility in the Si-undersaturated melts. The metasomatism between crustal melts and peridotite may represent an analog for processes in the mantle wedge above subducting slabs. As such, the crystallization of abundant zircon in ultramafic host rocks has implications for geochemistry of melts generated in the mantle and the widely reported depletion of high field strength elements in arc magmas.     

Dynamo-generated magnetic fields at the surface of a massive star

Spruit has shown that an astrophysical dynamo can operate in the non-convective material of a differentially rotating star as a result of a particular instability in the magnetic field (the Tayler instability). By assuming that the dynamo operates in a state of marginal instability, Spruit has obtained formulae which predict the equilibrium strengths of azimuthal and radial field components in terms of local physical quantities. Here, we apply Spruit's formulae to our previously published models of rotating massive stars in order to estimate Tayler dynamo field strengths. There are no free parameters in Spruit's formulae. In our models of 10- and  50-M_⊙  stars on the zero-age main sequence, we find internal azimuthal fields of up to 1 MG, and internal radial components of a few kG. Evolved models contain weaker fields. In order to obtain estimates of the field strength at the stellar surface, we examine the conditions under which the Tayler dynamo fields are subject to magnetic buoyancy. We find that conditions for Tayler instability overlap with those for buoyancy at intermediate to high magnetic latitudes. This suggests that fields emerge at the surface of a massive star between magnetic latitudes of about 45° and the poles. We attempt to estimate the strength of the field which emerges at the surface of a massive star. Although these estimates are very rough, we find that the surface field strengths overlap with values which have been reported recently for line-of-sight fields in several O and B stars.