Algorithm for non-hydrostatic dynamics in the Regional Oceanic Modeling System

A non-hydrostatic algorithm for the Regional Oceanic Modeling System (ROMS) is proposed. It is based on a decomposition technique for hydrostatic and non-hydrostatic pressure. The algorithm has a pressure-correction scheme with split-explicit time-stepping for baroclinic and barotropic vertical modes with a free surface. The algorithm implementation requires solving a Poisson equation for a non-hydrostatic pressure that has a non-symmetric matrix in discrete form. The efficiency of a different class of solvers and preconditioners were tested. The algorithm is successfully implemented with several examples where non-hydrostatic effects are important. These include standing external gravity waves; strongly nonlinear internal wave generation and transformation; stratified shear instability and its associated mixing; and nonlinear internal tidal generation over a ridge. The corresponding changes in the pre-processing and post-processing infrastructure in the existing hydrostatic ROMS code were performed to implement parallel elliptic solvers and a new set of dynamical equations.     

Surface gravity wave interaction with elastic bottom

In the present paper, a hydroelastic model is developed to deal with surface gravity wave interaction with an elastic bed based on the small amplitude water wave theory and plate deflection in finite water depth. The elastic bottom bed is modelled as a thin elastic plate and is based on the Euler-Bernoulli beam equation. The wave characteristics in the presence of the elastic bed is analyzed in both the cases of deep and shallow water waves. Further, the linearized long wave equation is generalized to include bottom flexibility. A generalized expansion formula for the velocity potential is derived to deal with the boundary value problems associated with surface gravity waves having an elastic bed. The utility of the expansion formula is illustrated by demonstrating specific physical problems which will play significant role in the analysis of wave structure interaction problems. Behavior of the wave spectra are discussed in the case of closed basin having a free surface and an elastic bottom topography.     

Resilience assessment on coastline changes and urban settlements: A case study in Seribu Islands, Indonesia

The sustainable development in the ocean and costal areas has been an issue for the archipelago nation. Since two decades ago, some archipelago nations have attempted to implement the concept of both scientifically and politically sounding sustainability. The vulnerability assessment is one of the methods that are being used to measure the ocean and coastal sustainability in order to have better evaluation and redesign of the land development as well as policy making. Most of the vulnerability assessment has been conducted based on pressures, damages and changes that involve in the region. A common understanding of the vulnerability assessment is that there are three aspects to be considered: hazards, resilience and damages. These three aspects must be well defined at first in order to have better indicators or sub-indices for the vulnerability index. There are several issues and factors that should be considered before performing the vulnerability assessment. Firstly, each country has different coastal characteristics due to a different geologic process. Secondly, the three aspects of the vulnerability (i.e. hazards, resilience and damages) are impacting on each country at a different scale. Thirdly, the vulnerability of a small island region is different from that of a large island region. Finally, policies and regulations vary in each country.From the data analysis results, it is found that the urban settlement in Seribu Islands is one of the resilient factors in addition to the geological and geomorphological conditions. The resilience factors in Seribu Islands are classified into four categories: 1) settlements area, 2) population density, 3) hard infrastructure such as airfields, ports and roads, 4) geological process such as abrasion and erosion. Based on the island characteristics of Seribu Islands, a unique vulnerability index that fits to this locality is developed. It is shown that the vulnerability index developed in this study can measure the resilience of Seribu Islands. In addition to the aforementioned resilience factors, the unique geographical condition and the geological stability in Seribu Islands made the outer islands become a barrier from oceanographic conditions and made the inner islands protected. However, the population growth made significant changes in terms of ecology, water, sanitation and pollution within the region.     

Deep-sea nematode assemblage has not recovered 26 years after experimental mining of polymetallic nodules (Clarion-Clipperton Fracture Zone, Tropical Eastern Pacific)

We investigated nematode assemblages inhabiting the 26-year-old track created by experimental deep-sea mining of polymetallic nodules, and two adjacent, undisturbed sites, one with nodules and one without nodules. The aim was to compare density, assemblage structure, and diversity indices in order to assess the process of recovery of the nematode assemblage inhabiting the disturbed site. This experimental dredging was conducted in 1978 by the Ocean Minerals Company (USA) in the area of a French mining claim in the Clarion-Clipperton Fracture Zone (Tropical Eastern Pacific) at a depth of about 5000 m. The nematode assemblage had not returned its initial state 26 years after the experimental dredging: the total nematode density and biomass within the dredging track were significantly lower than outside the track; the biodiversity indices showed significantly lower nematode diversity within the track; and the structure of the nematode assemblage within the track differed significantly from those in the two undisturbed sites outside the track. However, there were no significant differences in the mean body volumes of adult nematodes and adult-juvenile ratios between the track and reference sites. Parameters such as the rate of sediment restoration (which depends on local hydrological conditions) and the degree and character of the disturbance appeared to be of considerable importance for the recovery rate of the deep-sea nematode assemblages and their ability to recolonize disturbed areas. The rates of recolonization and recovery may vary widely in different deep-sea regions.     

Grazing suppression of dimethylsulphoniopropionate (DMSP) accumulation in iron-fertilised, sub-Antarctic waters

The impact of in situ iron fertilisation on the production of particulate dimethylsulphoniopropionate (DMSPp) and its breakdown product dimethyl sulphide (DMS) was monitored during the SOLAS air-sea gas exchange experiment (SAGE). The experiment was conducted in the high nitrate, low chlorophyll (HNLC) waters of the sub-Antarctic Southern Ocean (46.7°S 172.5°E) to the south-east of New Zealand, during March-April, 2004. In addition to monitoring net changes in the standing stocks of DMSPp and DMS, a series of dilution experiments were used to determine the DMSPp production and consumption rates in relation to increased iron availability. In contrast to previous experiments in the Southern Ocean, DMS concentrations decreased over the course of the 15-d iron-fertilisation experiment, from an integrated volume-specific concentration in the mixed layer on day 0 of 0.78 nM (measured values 0.65-0.91 nM) to 0.46 nM (measured values 0.42-0.47 nM) by day 15, in parallel with the surrounding waters. DMSPp, chlorophyll a and the abundance of photosynthetic picoeukaryotes exhibited indiscernible or only moderate increases in response to the raised iron availability, despite an obvious physiological response by the phytoplankton. High specific growth rates of DMSPp, equivalent to 0.8-1.2 doublings d⁻¹, occurred at the simulated 60% light level of the dilution experiments. Despite the high production rates, DMSPp accumulation was suppressed in part by microzooplankton grazers who consumed between 61% d⁻¹ and 126% d⁻¹ of the DMSPp production. Temporal trends in the rates of production and consumption illustrated a close coupling between the DMSP-producing phytoplankton and their microzooplankton grazers. Similar grazing and production rates were observed for the eukaryotic picophytoplankton that dominated the phytoplankton biomass, partial evidence that picoeukaryotes contributed a substantial proportion of the DMSP synthesis. These rates for DMSPp and picoeukaryotes were considerably higher than for chlorophyll a, indicating higher cycling rates of the DMSP-producing taxa than for the bulk phytoplankton community. When compared to the total phytoplankton community, there was no evidence of selection against the DMSP-containing phytoplankton by the microzooplankton grazers; the opposite appeared to be the case. SAGE demonstrated that increased iron availability in the HNLC waters of the Southern Ocean does not invariably lead to enhanced DMS sea-air flux. The potential suppression of DMSPp accumulation by grazers needs to be taken into account in future attempts to elevate DMS emission through in situ iron fertilisation and in understanding the hypothesised link between levels of Aeolian iron deposition in the Southern Ocean, DMS emission and global albedo.     

Impact of hydrotalcite deposition on biogeochemical processes in a shallow tropical bay

The biogeochemistry of a tropical shoal bay (Melville Bay, Australia) impacted by the effluent release, precipitation, and deposition of hydrotalcite from an alumina refinery was studied in both wet and dry seasons. Within the deposition zone, sulfate reduction dominated benthic carbon cycling accounting for ≈100% of total microbial activity, with rates greater than those measured in most other marine sediments. These rapid rates of anoxic metabolism resulted in high rates of sulfide and ammonium production and low C:S ratios, implying significant preservation of S in stable sulfide minerals. Rates of total microbial activity were significantly less in control sediments of equivalent grain size, where sulfate reduction accounted for ≈50% of total benthic metabolism. Rates of planktonic carbon cycling overlying the deposition zone were also greater than those measured in the control areas of southern Melville Bay. At the sediment surface, productive algal and cyanobacterial mats helped stabilize the sediment surface and oxidize sulfide to sulfate to maintain a fully oxygenated water-column overlying the impacted zone. The mats utilized a significant fraction of dissolved inorganic N and P released from the sea bed; some nutrients escaped to the water-column such that benthic regeneration of NH?+ and PO?3? accounted for 100% and 42% of phytoplankton requirements for N and P, respectively. These percentages are high compared to other tropical coastal environments and indicate that benthic nutrient recycling may be a significant factor driving water-column production overlying the deposition zone. With regard to remediation, it is recommended that the sea bed not be disturbed as attempts at removal may result in further environmental problems and would require specific assessment of the proposed removal process.     

Estimation of changes in characteristics of the climate and carbon cycle in the 21st century accounting for the uncertainty of terrestrial biota parameter values

ensemble simulations with the A.M. Obukhov Institute of Atmospheric Physics, Russian Academy of Sciences (IAP RAS) climate model (CM) for the 21st century are analyzed taking into account anthropogenic forcings in accordance with the Special Report on Emission Scenarios (SRES) A2, A1B, and B1, whereas agricultural land areas were assumed to change in accordance with the Land Use Harmonization project scenarios. Different realizations within these ensemble experiments were constructed by varying two governing parameters of the terrestrial carbon cycle. The ensemble simulations were analyzed with the use of Bayesian statistics, which makes it possible to suppress the influence of unrealistic members of these experiments on their results. It is established that, for global values of the main characteristics of the terrestrial carbon cycle, the SRES scenarios used do not differ statistically from each other, so within the framework of the model, the primary productivity of terrestrial vegetation will increase in the 21st century from 74 ± 1 to 102 ± 13 PgC yr⁻¹ and the carbon storage in terrestrial vegetation will increase from 511 ± 8 to 611 ± 8 PgC (here and below, we indicate the mean ± standard deviations). The mutual compensation of changes in the soil carbon stock in different regions will make global changes in the soil carbon storage in the 21st century statistically insignificant. The global CO₂ uptake by terrestrial ecosystems will increase in the first half of the 21st century, whereupon it will decrease. The uncertainty interval of this variable in the middle (end) of the 21st century will be from 1.3 to 3.4 PgC yr⁻¹ (from 0.3 to 3.1 PgC yr⁻¹). In most regions, an increase in the net productivity of terrestrial vegetation (especially outside the tropics), the accumulation of carbon in this vegetation, and changes in the amount of soil carbon stock (with the total carbon accumulation in soils of the tropics and subtropics and the regions of both accumulation and loss of soil carbon at higher latitudes) will be robust within the ensemble in the 21st century, as will the CO₂ uptake from the atmosphere only by terrestrial ecosystems located at extratropical latitudes of Eurasia, first and foremost by the Siberian taiga. However, substantial differences in anthropogenic emissions between the SRES scenarios in the 21st century lead to statistically significant differences between these scenarios in the carbon dioxide uptake by the ocean, the carbon dioxide content in the atmosphere, and changes in the surface air temperature. In particular, according to the SRES A2 (A1B, B1) scenario, in 2071–2100 the carbon flux from the atmosphere to the ocean will be 10.6 ± 0.6 PgC yr⁻¹ (8.3 ± 0.5, 5.6 ± 0.3 PgC yr⁻¹), and the carbon dioxide concentration in the atmosphere will reach 773 ± 28 ppmv (662 ± 24, 534 ± 16 ppmv) by 2100. The annual mean warming in 2071–2100 relatively to 1961–1990 will be 3.19 ± 0.09 K (2.52 ± 0.08, 1.84 ± 0.06 K).     

Development of the summer breeze circulation in the west part of the Black Sea

On the basis of the results of regional reanalysis of the atmospheric circulation presented with a resolution of 9 km, the process of formation of breezes is studied for the case of weak synoptic background activity in the rectilinear part of the west coast of the Black Sea for the period 01–04.07.2007. It is shown that the gravitational currents, breeze fronts, and intense internal waves are formed under these conditions in the troposphere during the daytime. The Hovm?ller diagrams of the wind velocity and the maps of the vertical sections of potential temperature, vertical velocity of the air masses, and other parameters are presented. On the basis of these diagrams and maps, we obtain quantitative estimates of the wind velocity in the breeze and of the velocity of propagation of the breeze front and compare these estimates with the available literature data.     

Analysis of volcanogenic perturbations of the subarctic ozonosphere on the basis of space monitoring data

The response of the total ozone (TO) at subarctic latitudes to volcanic eruptions, products of which were injected into the stratosphere, is analyzed. It is established that the behavior of the series of average annual TO values according to the TOMS, SCIAMACHY, and GOME space equipment data averaged for 55°–65° N latitudes agrees with the activity of explosive volcanic eruptions. The series of the TO satellite monitoring instrumental data are extended by 200 years into the past using a reconstruction from the dendro-chronologic data. An analysis of the series of TO reconstructed values indicates that volcanogenic perturbations of the subarctic ozonosphere initiate long-term negative TO deviations. In this case, the TO negative deviation depth depends on the frequency of the ozonosphere volcanogenic perturbations and the phase of quasiperiodic oscillation cycles rather than on the strength of a single volcanic explosion.