Development and extension of an aggregated scale model: Part 2 — Extensions to ASMITA

Abstract

Whilst much attention has been given to models that describe wave, tide and sediment transport processes in sufficient detail to determine the local changes in bed level over a relatively detailed representation of the bathymetry, far less attention has been given to models that consider the problem at a much larger scale (e.g. that of geomorphological elements such as a tidal flat and tidal channel). Such aggregated or lumped models tend not to represent the processes in detail but rather capture the behaviour at the scale of interest. One such model developed using the concept of an equilibrium concentration is the Aggregated Scale Morphological Interaction between Tidal basin and Adjacent coast (ASMITA). In a companion paper (Part 1), we detail the original model and provide some new insights into the concepts of equilibrium, and horizontal and vertical exchange that are key components of this modelling approach. In this paper, we summarise a range of developments that have been undertaken to extend the original model concept, to illustrate the flexibility and power of the conceptual framework. However, adding detail progressively moves the model in the direction of the more detailed process-based models and we give some consideration to the boundary between the two.

Highlights

• A range of extensions to the original ASMITA model are presented that allow additional processes or features to be represented in the model.
• The merits of ever increasing complexity within an aggregated model, versus the use of a suitable local-scale and more detailed process-based model are discussed.

     

From Jurassic rifting to Cretaceous subduction in NW Iranian Azerbaijan: geochronological and geochemical signals from granitoids

Previous interpretations of a Jurassic subduction in Iran were based on trace element classification diagrams for granitoids, but their reliability is questionable, underscored by modern examples of continental break-up zones such as the Baja California. We present new field observations, bulk rock geochemistry, Sr and Nd isotope analyses and U–Pb zircon geochronology to assess the age and tectonic setting of previously undated intermediate to felsic magmatic rocks cropping out in the Precambrian basement of NW Iranian Azerbaijan. The geochronology revealed an uneven distribution in space and time: Late Jurassic (159–154 Ma) intrusions and dikes are alkaline to calc-alkaline. Their melt source is mantle dominated with a distinct continental contribution disclosed by radiogenic isotopes and abundant inherited zircon cores. Mid-Cretaceous (112–96 Ma) plutonic bodies and associated volcanic rocks occur only to the east of the major Siah Cheshmeh–Khoy Fault. They have geochemical signatures typical of a metasomatized mantle. In consistence with the sedimentation history of the area, our new interpretation attributes the Late Jurassic magmatism to thinning of a continental lithosphere in a rift-related setting. Mid-Cretaceous magmatism was produced by oceanic subduction beneath the Central Iran continent. We interpret the 40-Ma age gap between the two magmatic episodes as the time of opening of the oceanic basin witnessed by the Khoy ophiolite in the study area.     

The geoarchaeological and paleoenvironmental context of Paleoindian sites in western Middle Park,Colorado, USA

Geoarchaeological investigations in western Middle Park provide important information for understanding the soil‐stratigraphic context of Paleoindian components, as well as the latest Quaternary environmental change and landscape evolution in a Southern Rocky Mountain intermontane basin. Paleoindian components are associated with the oldest two of four latest Quaternary stratigraphic units (1–4) recognized in co‐alluvial mantles (combined slopewash and colluvium) in uplands and in alluvial valley fills. Limited data suggest accumulation of unit 1 as early as ∼12,500 ¹⁴C yr B.P. in alluvial valleys and by at least ∼11,000 ¹⁴C yr B.P. in uplands was followed by brief stability and soil formation. A relatively widespread disconformity marks earliest Holocene erosion and substantial removal of latest Pleistocene deposits in upland and alluvial settings followed by unit 2 deposition ∼10,000–9000 ¹⁴C yr B.P., perhaps signaling the abrupt onset of an intensified summer monsoon. In situ Paleoindian components in uplands are found in a moderately developed buried soil (the Kremmling soil) formed in units 1 and 2 in thin (≤1m) hillslope co‐alluvial mantles. The Kremmling soil reflects geomorphic stability in upland and alluvial settings ∼9000–4500 ¹⁴C yr BP, and represents a buried landscape with the potential to contain additional Paleoindian components, although elsewhere in western Middle Park Early Archaic components are documented in morphologically similar soils. Kremmling soil morphology, the relative abundance of charcoal in unit 2 relative to younger units, and charcoal morphology indicate the expansion of forest cover, including Pinus, and grass cover during the early and middle Holocene, suggesting conditions moister than present. © 2010 Wiley Periodicals, Inc.     

Experimental constraints on Pt, Pd and Au partitioning and fractionation in silicate melt–sulfide–oxide–aqueous fluid systems at 800 °C, 150 MPa and variable sulfur fugacity

We have performed experiments to constrain the effect of sulfur fugacity (fS₂) and sulfide saturation on the fractionation and partitioning behavior of Pt, Pd and Au in a silicate melt–sulfide crystal/melt–oxide–supercritical aqueous fluid phase–Pt–Pd–Au system. Experiments were performed at 800 °C, 150 MPa, with oxygen fugacity (fO₂) fixed at approximately the nickel–nickel oxide buffer (NNO). Sulfur fugacity in the experiments was varied five orders of magnitude from approximately log fS₂ = 0 to log fS₂ = −5 by using two different sulfide phase assemblages. Assemblage one consisted initially of chalcopyrite plus pyrrhotite and assemblage two was loaded with chalcopyrite plus bornite. At run conditions pyrrhotite transformed compositionally to monosulfide solid solution (mss), chalcopyrite to intermediate solid solution (iss), and in assemblage two chalcopyrite and bornite formed a sulfide melt. Run-product silicate glass (i.e., quenched silicate melt) and crystalline materials were analyzed by using both electron probe microanalysis and laser ablation inductively coupled plasma mass spectrometry. The measured concentrations of Pt, Pd and Au in quenched silicate melt in runs with log fS₂ values ranging from approximately 0.0 to −5.0 do not exhibit any apparent dependence on fS₂. The measured Pt, Pd and Au concentrations in mss do vary as a function of fS₂. The measured Pt, Pd and Au concentrations in iss do not appear dependent on fS₂. The data suggest that fS₂, working in concert with fO₂, via the determinant role that these variables play in controlling the magmatic sulfide phase assemblage and the solubility of Pt, Pd and Au as lattice bound components in magmatic sulfide phases, is a controlling factor on the budgets of Pt, Pd and Au during the evolution of magmatic systems.     

On the photodissociation of water vapour in the mesosphere

The photodissociation of water vapour in the mesosphere depends on the absorption of solar radiation in the region (175–200 nm) of the O₂ Schumann-Runge band system and also at H-Lyman alpha. The photodissociation products are OH + H, OH + H, O + 2H and H₂ + O at Lyman alpha; the percentages for these four channels are 70, 8, 12 and 10%, respectively, but OH + H is the only channel between 175 and 200 nm. Such proportions lead to a production of H atoms corresponding to practically the total photodissociation of H₂O, while the production of H₂ molecules is only 10% of the H₂O photodissociation by Lyman alpha.The photodissociation frequency (s^?1) at Lyman alpha can be expressed by a simple formula

$\text{J}{}_{\mathrm{<mtext>Ly</mtext><mtext>\alpha </mtext>}}\text{H}{}_2\text{O=4.5 ×10}{}^{\mathrm{?6}}\text{1+0.2}\text{F}{}_{10.7}\text{?65}\text{100}\text{exp}\text{[?4.4 ×10}{}^{\mathrm{?19}}\text{N}{}^{0.917}\text{]}$

where F_{10.7 cm} is the solar radioflux at 10.7 cm and N the total number of O₂ molecules (cm^?2), and when the following conventional value is accepted for the Lyman alpha solar irradiance at the top of the Earth's atmosphere ($\text{Δλ = 3.5}\text{A}\text{?}$) q_Lyα,∞ = 3 × 10¹¹ photons cm^?2 s^1?.The photodissociation frequency for the Schumann-Runge band region is also given for mesospheric conditions by a simple formula

$\text{J}{}_{\mathrm{<mtext>SRB</mtext>}}\text{(}\text{H}{}_2\text{O}\text{) = J}{}_{\mathrm{<mtext>SRB</mtext><mtext>,\infty </mtext>}}\text{(}\text{H}{}_2\text{O}\text{) exp [?10}{}^{\mathrm{?7}}\text{N}{}^{0.35}\text{]}$

where J_SRB,∞(H₂O) = 1.2 × 10^?6 and 1.4 × 10^?6 s^?1 for quiet and active sun conditions, respectively.The precision of both formulae is good, with an uncertainty less than 10%, but their accuracy depends on the accuracy of observational and experimental parameters such as the absolute solar irradiances, the variable transmittance of O₂ and the H₂O effective absorption cross sections. The various uncertainties are discussed. As an example, the absolute values deduced from the above formulae could be decreased by about 25-20% if the possible minimum values of the solar irradiances were used.     

Atmospheric input of nitrogen into the North Sea: ANICE project overview

The aim of the atmospheric nitrogen inputs into the coastal ecosystem (ANICE) project is to improve transport–chemistry models that estimate nitrogen deposition to the sea. To achieve this, experimental and modelling work is being conducted which aims to improve understanding of the processes involved in the chemical transformation, transport and deposition of atmospheric nitrogen compounds. Of particular emphasis within ANICE is the influence of coastal zone processes. Both short episodes with high deposition and chronic nitrogen inputs are considered in the project. The improved transport–chemistry models will be used to assess the atmospheric inputs of nitrogen compounds into the European regional seas (the North Sea is studied as a prototype) and evaluate the impact of various emission reduction strategies on the atmospheric nitrogen loads. Assessment of the impact of atmospheric nitrogen on coastal ecosystems will be based on comparisons of phytoplankton nitrogen requirements, other external nitrogen inputs to the ANICE area of interest and the direct nitrogen fluxes provided by ANICE. Selected results from both the experimental and modelling components are presented here. The experimental results show the large spatial and temporal variability in the concentrations of gaseous nitrogen compounds, and their influences on fluxes. Model calculations show the strong variation of both concentrations and gradients of nitric acid at fetches of up to 25 km. Aerosol concentrations also show high temporal variability and experimental evidence for the reaction between nitric acid and sea salt aerosol is provided by size-segregated aerosol composition measured at both sides of the North Sea. In several occasions throughout the experimental period, air mass back trajectory analysis showed connected flow between the two sampling sites (the Weybourne Atmospheric Observatory on the North Norfolk coast of the UK and Meetpost Noordwijk, a research tower at 9 km off the Dutch coast). Results from the METRAS/SEMA mesoscale chemistry transport model system for one of these cases are presented. Measurements of aerosol and rain chemical composition, using equipment mounted on a commercial ferry, show variations in composition across the North Sea. These measurements have been compared to results obtained with the transport–chemistry model ACDEP which calculates the atmospheric inputs into the whole North Sea area. Finally, the results will be made available for the assessment of the impact of atmospheric nitrogen on coastal ecosystems.     

Cross-shore mean flow in the surf zone

The seaward mean returnflow or undertow in the surf zone, which compensates for the shoreward mass flux above the wave trough level, is found to be driven by the force imbalance between the wave momentum flux on the one hand and the set-up on the other hand, in qualitative agreement with the model of Dyhr-Nielsen and Sørensen (1970). A model, which quantifies this imbalance, is shown to yield theoretical results in good accordance with experiments when the proper boundary conditions are accounted for.     

Nonlinear saturation-based whitecapping dissipation in SWAN for deep and shallow water

This study investigates the effectiveness of a revised whitecapping source term in the spectral wind wave model SWAN (Simulating WAves Nearshore) that is local in frequency space, nonlinear with respect to the variance density and weakly dependent on the wave age. It is investigated whether this alternative whitecapping expression is able to correct the tendency towards underprediction of period measures that has been identified in the default SWAN model. This whitecapping expression is combined with an alternative wind input source term that is more accurate for young waves than the default expression. The shallow water source terms of bottom friction, depth-induced breaking and triad interaction are left unaltered. It is demonstrated that this alternative source term combination yields improved agreement with fetch- and depth-limited growth curves. Moreover, it is shown, by means of a field case over a shelf sea, that the investigated model corrects the erroneous overprediction of wind-sea energy displayed by the default model under combined swell-sea conditions. For a selection of field cases recorded at two shallow lakes, the investigated model generally improves the agreement with observed spectra and integral parameters. The improvement is most notable in the prediction of period measures.     

Neural network modelling of wave overtopping at coastal structures

A method has been developed to estimate wave overtopping discharges for a wide range of coastal structures. The prediction method is based on Neural Network modelling. For this purpose use is made of a data set obtained from a large number of physical model tests (collected within the framework of the European project CLASH, see e.g. [Steendam, G.J., Van der Meer, J.W., Verhaeghe, H., Besley, P., Franco, L. and Van Gent, M.R.A. (2004). The international database on wave overtopping. World Scientific, Proc. 29th ICCE, vol. 4, pp. 4301–4313, Lisbon, Portugal.]). Moreover, a method was developed to obtain confidence intervals for the overtopping predictions of the neural network.