Modelling and mass balance assessments of nutrient retention in a seasonally-flowing estuary (Swan River Estuary,Western Australia)

An integrated mass balance and modelling approach for analysis of estuarine nutrient fluxes is demonstrated in the Swan River Estuary, a microtidal system with strong hydrological dependence on seasonal river inflows. Mass balance components included estimation of gauged and ungauged inputs to the estuary and losses to the ocean (outflow and tidal exchange). Modelling components included estimation of atmospheric (N fixation, denitrification) and sediment–water column nutrient exchanges. Gross and net denitrification derived using two independent methods were significantly correlated (r² = 0.49, p < 0.01) with net rates averaging 40% of gross. Annual nitrogen (N) and phosphorus (P) loads from major tributaries were linearly correlated with annual freshwater discharge and were 3-fold higher in wet years than in dry years. Urban drains and groundwater contributed, on average, 26% of N inputs and 19% of P inputs, with higher relative contributions in years of low river discharge. Overall, ungauged inputs accounted for almost 35% of total nitrogen loads. For N, elevated loading in wet years was accompanied by large increases in outflow (7x) and tidal flushing (2x) losses and resulted in overall lower retention efficiency (31%) relative to dry years (70%). For P, tidal flushing losses were similar in wet and dry years, while outflow losses (4-fold higher) were comparable in magnitude to increases in loading. As a result, P retention within the estuary was not substantially affected by inter-annual variation in water and P loading (ca. 50% in all years). Sediment nutrient stores increased in most years (remineralisation efficiency ca. 50%), but sediment nutrient releases were significant and in some circumstances were a net source of nutrients to the water column.     

Floating and fixed artificial habitats: Spatial and temporal patterns of benthic communities in a coral reef environment

While natural marine habitats with motion capabilities, e.g., kelps and seaweeds, have been studied alongside their associated fouling communities, little is known of the effect of motion on the communities of floating artificial habitats such as buoys, rafts, and pontoons, particularly in tropical systems. Hydrodynamic features greatly differ between floating and fixed artificial substrata, which in turn affect the structure of their associated communities. This study tested the hypothesis that floating and fixed artificial installations in a tropical reef system (Eilat, Red Sea) would support different benthic communities throughout space and time. Specifically, we examined differences in communities recruited onto settlement plates between floating and fixed installations deployed at three different sites, along a two-year monitoring period. The three sites exhibited distinct differences in species assemblages between the monitoring dates (6, 12, 18 and 24 months post deployment), mainly between the first and the last two dates. The average level of dissimilarity between floating and fixed installations increased over time at all sites. Over 50% of the dissimilarity between the floating and fixed installations resulted from five taxonomic groups i.e., bryozoans, bivalves, barnacles, sponges, including the amount of bare space on the settlement plates. The contribution of these groups to the dissimilarity changed both temporally within each site, and spatially among sites. The observed differences were related to the hydrodynamic characteristics of floating and fixed habitats, interacting with biotic features such as predation, successional processes and seasonality; and abiotic features including small-scale spatial changes, light, and position in the water column.     

Simulation of nonlinear wave run-up with a high-order Boussinesq model

This paper considers the numerical simulation of nonlinear wave run-up within a highly accurate Boussinesq-type model. Moving wet–dry boundary algorithms based on so-called extrapolating boundary techniques are utilized, and a new variant of this approach is proposed in two horizontal dimensions. As validation, computed results involving the nonlinear run-up of periodic as well as transient waves on a sloping beach are considered in a single horizontal dimension, demonstrating excellent agreement with analytical solutions for both the free surface and horizontal velocity. In two horizontal dimensions cases involving long wave resonance in a parabolic basin, solitary wave evolution in a triangular channel, and solitary wave run-up on a circular conical island are considered. In each case the computed results compare well against available analytical solutions or experimental measurements. The ability to accurately simulate a moving wet–dry boundary is of considerable practical importance within coastal engineering, and the extension described in this work significantly improves the nearshore versatility of the present high-order Boussinesq approach.     

Performance of community-based coastal resource management (CBCRM) programs in the Philippines: A meta-analysis

Community-based coastal resource management (CBCRM) is a major conservation and fisheries management strategy in the tropics. In this study, the performance of 16 CBCRM programs in the Philippines was assessed using a meta-analysis of eight indicators that represented the perceptions of local resource users. Overall, the CBCRM programs in the Philippines were perceived to have a significant positive impact. However, the performance of each of the indicators was mixed. Although the CBCRM programs were perceived to be effective in empowering the local fishing communities, their perceived impact on improving the state of the local fisheries resources remained limited. This highlights the importance of incorporating ecological and socio-economic considerations in setting fisheries management regimes.     

Application of a vanishing, quasi-sigma, vertical coordinate for simulation of high-speed, deep currents over the Sigsbee Escarpment in the Gulf of Mexico

Recent observations over the Sigsbee Escarpment in the Gulf of Mexico have revealed extremely energetic deep currents (near 1 m s⁻¹), which are trapped along the escarpment. Both scientific interest and engineering needs demand dynamical understanding of these extreme events, and can benefit from a numerical model designed to complement observational and theoretical investigations in this region of complicated topography. The primary objective of this study is to develop a modeling methodology capable of simulating these physical processes and apply the model to the Sigsbee Escarpment region. The very steep slope of the Sigsbee Escarpment (0.05–0.1) limits the application of ocean models with traditional terrain-following (sigma) vertical coordinates, which may represent the very complicated topography in the region adequately, can result in large truncation errors during calculation of the horizontal pressure gradient. A new vertical coordinate system, termed a vanishing quasi-sigma coordinate, is implemented in the Navy Coastal Ocean Model for application to the Sigsbee Escarpment region. Vertical coordinate surfaces for this grid have noticeably gentler slopes than a traditional sigma grid, while still following the terrain near the ocean bottom. The new vertical grid is tested with a suite of numerical experiments and compared to a classical sigma-layer model. The numerical error is substantially reduced in the model with the new vertical grid. A one-year, realistic, numerical simulation is performed to simulate strong, deep currents over the Escarpment using a very-high-resolution nested modeling approach. The model results are analyzed to demonstrate that the deep-ocean currents in the simulation replicate the prominent dynamical features of the observed intense currents in the region.     

Different approaches to model error formulation in 4D-Var: a study with high-resolution advection schemes

All numerical models are imperfect. Weak constraint variational data assimilation ( VDA ), which provides a treatment of the modelling errors, is studied; building on the approach of Vidard et al. (Tellus, 56 A, pp. 177–188, 2004). The evolution of model error ( ME ) is modelled using ordinary differential equations, which involve a scalar parameter. These approaches were tested using different high-resolution advection schemes. The first set of experiments were constructed to see if it is possible to account for (numerical) discretization error within such a framework. In other set of experiments, a systematic source of modelling error was introduced by deliberately specifying an incorrect value for the Coriolis parameter in the model. Results with observational state at half of the model state resolution, are also presented. We also discuss a method of estimating the scalar parameter in the ME through VDA . In all cases, the inclusion of ME provides reduction in forecasting errors. Also, our experiments indicate that different settings of the model (e.g. using different high-resolution advection schemes) would need different ME formulation. Results presented in this paper could be used to formulate sophisticated ME forms to account for systematic errors in higher dimensional models with complex advection schemes.     

Interaction of linear waves with an array of infinitely long horizontal circular cylinders in a two-layer fluid of infinite depth

The linear water wave scattering and radiation by an array of infinitely long horizontal circular cylinders in a two-layer fluid of infinite depth is investigated by use of the multipole expansion method. The diffracted and radiated potentials are expressed as a linear combination of infinite multipoles placed at the centre of each cylinder with unknown coefficients to be determined by the cylinder boundary conditions. Analytical expressions for wave forces, hydrodynamic coefficients, reflection and transmission coefficients and energies are derived. Comparisons are made between the present analytical results and those obtained by the boundary element method, and some examples are presented to illustrate the hydrodynamic behavior of multiple horizontal circular cylinders in a two-layer fluid. It is found that for two submerged circular cylinders the influence of the fluid density ratio on internal-mode wave forces is more appreciable than surface-mode wave forces, and the periodic oscillations of hydrodynamic results occur with the increase of the distance between two cylinders; for four submerged circular cylinders the influence of adding two cylinders on the wave forces of the former cylinders is small in low and high wave frequencies, but the influence is appreciable in intermediate wave frequencies.     

Soft computing approach for real-time estimation of missing wave heights

This paper presents soft computing approach for estimation of missing wave heights at a particular location on a real-time basis using wave heights at other locations. Six such buoy networks are developed in Eastern Gulf of Mexico using soft computing techniques of Artificial Neural Networks (ANN) and Genetic Programming (GP). Wave heights at five stations are used to estimate wave height at the sixth station. Though ANN is now an established tool in time series analysis, use of GP in the field of time series forecasting/analysis particularly in the area of Ocean Engineering is relatively new and needs to be explored further. Both ANN and GP approach perform well in terms of accuracy of estimation as evident from values of various statistical parameters employed. The GP models work better in case of extreme events. Results of both approaches are also compared with the performance of large-scale continuous wave modeling/forecasting system WAVEWATCH III. The models are also applied on real time basis for 3 months in the year 2007. A software is developed using evolved GP codes (C++) as back end with Visual Basic as the Front End tool for real-time application of wave estimation model.     

Mesozooplankton biomass and grazing responses to Cyclone Opal, a subtropical mesoscale eddy

As part of E-Flux III cruise studies in March 2005, plankton net collections were made to assess the effects of a cyclonic cold-core eddy (Cyclone Opal) on the biomass and grazing of mesozooplankton. Mesozooplankton biomass in the central region of Cyclone Opal, an area of uplifted nutricline and a subsurface diatom bloom, averaged 0.80±0.24 and 1.51±0.59 g DW m⁻², for day and night tows, respectively. These biomass estimates were about 80% higher than control (OUT) stations, with increases more or less proportionately distributed among size classes from 0.2 to >5 mm. Though elevated relative to surrounding waters south of the Hawaiian Islands (Hawai’i lee), total biomass and size distribution in Cyclone Opal were almost exactly the same as contemporary measurements made at Stn. ALOHA, 100 km north of the islands, by the HOT (Hawaii Ocean Time-series) Program. Mesozooplankton biomass and community composition at the OUT stations were also similar to ALOHA values from 1994 to 1996, preceding a recent decadal increase. These comparisons may therefore provide insight into production characteristics or biomass gradients associated with decadal changes at Stn. ALOHA. Gut fluorescence estimates were higher in Opal than in ambient waters, translating to grazing impacts of 0.11±0.02 d⁻¹ (IN) versus 0.03±0.01 d⁻¹ (OUT). Over the depth-integrated euphotic zone, mesozooplankton accounted for 30% of the combined grazing losses of phytoplankton to micro- and meso-herbivores in Opal, as compared to 13% at control stations. Estimates of active export flux by migrating zooplankton averaged 0.81 mmol C m⁻² d⁻¹ in Cyclone Opal and 0.37 mmol C m⁻² d⁻¹ at OUT stations, 53% and 24%, respectively, of the carbon export measured by passive sediment traps. Migrants also exported 0.18 mmol N m⁻² d⁻¹ (117% of trap N flux) in Cyclone Opal compared to 0.08 mmol N m⁻² d⁻¹ (51% of trap flux) at control stations. Overall, the food-web importance of mesozooplankton increased in Cyclone Opal both in absolute and relative terms. Diel migrants provided evidence for enhanced export flux in the eddy that was missed by sediment trap and ²³⁴Th techniques, and migrant-mediated flux was the major export term in the observed bloom-perturbation response and N mass balance of the eddy.