Variability in bedform morphology and kinematics with transport stage

Bedform geometry is widely recognized to be a function of transport stage. Bedform aspect ratio (height/length) increases with transport stage, reaches a maximum, then decreases as bedforms washout to a plane bed. Bedform migration rates are also linked to bedform geometry, in so far as smaller bedforms in coarser sediment tend to migrate faster than larger bedforms in finer sediment. However, how bedform morphology (height, length and shape) and kinematics (translation and deformation) change with transport stage and suspension have not been examined. A series of experiments is presented where initial flow depth and grain size were held constant and the transport stage was varied to produce bedload dominated, mixed‐load dominated and suspended‐load dominated conditions. The results show that the commonly observed pattern in bedform aspect ratio occurs because bedform height increases then decreases with transport stage, against a continuously increasing bedform length. Bedform size variability increased with transport stage, leading to less uniform bedform fields at higher transport stage. Total translation‐related and deformation‐related sediment fluxes all increased with transport stage. However, the relative contribution to the total flux changed. At the bedload dominated stage, translation‐related and deformation‐related flux contributed equally to the total flux. As the transport stage increased, the fraction of the total load contributed by translation increased and the fraction contributed by deformation declined because the bedforms got bigger and moved faster. At the suspended‐load dominated transport stage, the deformation flux increased and the translation flux decreased as a fraction of the total load, approaching one and zero, respectively, as bedforms washed out to a plane bed.     

Response of low‐angle dunes to variable flow

Current understanding of bedform dynamics is largely based on field and laboratory observations of bedforms in steady flow environments. There are relatively few investigations of bedforms in flows dominated by unsteadiness associated with rapidly changing flows or tides. As a consequence, the ability to predict bedform response to variable flow is rudimentary. Using high‐resolution multibeam bathymetric data, this study explores the dynamics of a dune field developed by tidally modulated, fluvially dominated flow in the Fraser River Estuary, British Columbia, Canada. The dunes were dominantly low lee angle features characteristic of large, deep river channels. Data were collected over a field ca 1·0 km long and 0·5 km wide through a complete diurnal tidal cycle during the rising limb of the hydrograph immediately prior to peak freshet, yielding the most comprehensive characterization of low‐angle dunes ever reported. The data show that bedform height and lee angle slope respond to variable flow by declining as the tide ebbs, then increasing as the tide rises and the flow velocities decrease. Bedform lengths do not appear to respond to the changes in velocity caused by the tides. Changes in the bedform height and lee angle have a counterclockwise hysteresis with mean flow velocity, indicating that changes in the bedform geometry lag changes in the flow. The data reveal that lee angle slope responds directly to suspended sediment concentration, supporting previous speculation that low‐angle dune morphology is maintained by erosion of the dune stoss and crest at high flow, and deposition of that material in the dune trough.     

Multi-species chemical data assimilation with the Danish Eulerian hemispheric model: system description and verification

Satellite retrievals of atmospheric composition provide a wealth of data on a global scale. These complement results from atmospheric chemistry-transport models (CTMs), and can be combined using data assimilation. We present two assimilation schemes coupled to the Danish Eulerian Hemispheric Model (DEHM), a three-dimensional, off-line CTM with full photochemistry: a variant on the ensemble Kalman filter and the three-dimensional variational scheme. The aim of this paper is to describe the two schemes and present an initial assessment of their impacts on model skill. Retrievals of multiple atmospheric trace gases are assimilated, namely: NO₂ tropospheric column densities, CH₄ total column densities, and partial column concentrations of O₃, CO and CH₄; these data are retrieved from four satellite sensors. Data for each species are assimilated independently of one another, and other species are only adjusted indirectly via the model’s chemistry and dynamics. Assimilation results are compared with measurements from surface monitoring stations and other satellite retrievals, and preliminary validation results are presented.Reference simulations (without assimilation) grossly underestimate surface CO concentrations, and both assimilation schemes eliminate this large and systematic model bias. The assimilation improves the spatial correlation of modelled CO with surface observations, and improves the spatial correlation between forecasts and retrievals for CO, NO₂ and O₃. Results for CH₄ show a loss of skill due to a mismatch in model bias between two assimilated CH₄ data-sets. Finally, we discuss differences in methodology and results between this paper and a recent study on multi-species chemical data assimilation. Joint optimisation of initial conditions and emission rates offers a promising direction for improving modelled boundary-layer concentrations.     

Initiation of large-volume silicic centers in the Yellowstone hotspot track: insights from H<Subscript>2</Subscript>O- and F-rich quartz-hosted rhyolitic melt inclusions in the Arbon Valley Tuff of the Snake River Plain

During the onset of caldera cluster volcanism at a new location in the Snake River Plain (SRP), there is an increase in basalt fluxing into the crust and diverse silicic volcanic products are generated. The SRP contains abundant and compositionally diverse hot, dry, and often low-δ¹⁸O silicic volcanic rocks produced through time during the formation of individual caldera clusters, but more H₂O-rich eruptive products are rare. We report analyses of quartz-hosted melt inclusions from pumice clasts from the upper and lower Arbon Valley Tuff (AVT) to gain insight into the initiation of caldera cluster volcanism. The AVT, a voluminous, caldera-forming rhyolite, represents the commencement of volcanism (10.44 Ma) at the Picabo volcanic field of the Yellowstone hotspot track. This is a normal δ¹⁸O rhyolite consisting of early and late erupted members (lower and upper AVT, respectively) with extremely radiogenic Sr isotopes and unradiogenic Nd isotopes, requiring that ~50 % of the mass of these elements is derived from melts of Archean upper crust. Our data reveal distinctive features of the early erupted lower AVT melt including: variable F concentrations up to 1.4 wt%, homogenous and low Cl concentrations (~0.08 wt%), H₂O contents ranging from 2.3 to 6.4 wt%, CO₂ contents ranging from 79 to 410 ppm, and enrichment of incompatible elements compared to the late erupted AVT, subsequent Picabo rhyolites, SRP rhyolites, and melt inclusions from other metaluminous rhyolites (e.g., Bishop Tuff, Mesa Falls Tuff). We couple melt inclusion data with Ti measurements and cathodoluminescence (CL) imaging of the host quartz phenocrysts to elucidate the petrogenetic evolution of the AVT rhyolitic magma. We observe complex and multistage CL zoning patterns, the most critical being multiple truncations indicative of several dissolution–reprecipitation episodes with bright CL cores (higher Ti) and occasional bright CL rims (higher Ti). We interpret the high H₂O, F, F/Cl, and incompatible trace element concentrations in the context of a model involving melting of Archean crust and mixing of the crustal melt with basaltic differentiates, followed by multiple stages of fractional crystallization, remelting, and melt extraction. This multistage process, which we refer to as distillation, is further supported by the complex CL zoning patterns in quartz. We interpret new Δ¹⁸O(Qz-Mt) isotope measurements, demonstrating a 0.4 ‰ or ~180 °C temperature difference, and strong Sr isotopic and chemical differences between the upper and lower AVT to represent two separate eruptions. Similarities between the AVT and the first caldera-forming eruptions of other caldera clusters in the SRP (Yellowstone, Heise and Bruneau Jarbidge) suggest that the more evolved, lower-temperature, more H₂O-rich rhyolites of the SRP are important in the initiation of a caldera cluster during the onset of plume impingement.     

Reefs at the edge: coral community structure around Rapa,southernmost French Polynesia

Rapa (27°36′ S, 144°20′ W) is a small (~40 km²) volcanic island isolated in the Southern Austral Archipelago, where direct anthropogenic stressors are extremely limited. Here, we present the results of the first quantitative survey of coral community structure across habitats and depths around the island. Despite its geographical isolation in the depauperate South Central Pacific, its small size and unfavourable environmental conditions (competition with macroalgae, low sea surface temperatures, reduced reef accretion), the diversity of scleractinian corals at Rapa is particularly high (112 species from 32 genera, including 37 species of Acropora) in comparison to other French Polynesian islands and subtropical Pacific locations. Our results indicate that the abundance (>100 colonies per 10 m² recorded at nine of the 17 sampling stations) and cover (>40% at four stations) of corals are relatively high for a marginal reef location. Strong spatial heterogeneity was found, with high variation in diversity, abundance, cover and community composition among stations. Variation in community composition was related to habitat types, with distinct assemblages among fringing reefs within bays, reef formations at bay entrances, and those on the submerged platform surrounding the island. On the platform, a depth gradient was detected, with generic richness, abundance and cover generally greater at deeper stations (18–20 m depth) compared with medium‐depth (10–12 m) and shallow (1–3 m) stations. A gradient was also recorded along bays, with increasing coral diversity and abundance from the bay heads to the bay entrances. The coral community at Rapa was characterized by the presence of several taxa not found in other French Polynesian archipelagos and the rarity of others that are common and abundant in the Society and the Tuamotu islands. Another distinctive feature of reef communities at Rapa is the high cover and dominance of macroalgae, particularly in the shallower parts of the surrounding platform, which probably explains the lower densities of coral colonies recorded there. These characteristics of the diversity and biogeographical composition of coral assemblages at Rapa provide considerable ecological grounds for its conservation.     

Dasymetric Spatiotemporal Interpolation

This research applies the principles of dasymetric mapping to spatiotemporal interpolation by extending the spatial concepts of zone and area to their temporal analogs of interval and duration, respectively. An example application of dasymetric spatiotemporal interpolation using crime event data is presented. Results indicate that dasymetric spatiotemporal interpolation significantly improves the accuracy of estimates over areal or duration weighting. In addition, even when dasymetric interpolation in either the spatial or temporal dimension is relatively weak, combining dasymetric estimation in both space and time dimensions simultaneously has the potential to amplify the accuracy of the overall dasymetric estimation.