Optimal multispecies fishery management in the Baltic Sea: A bioeconomic analysis

The aim of the thesis is to assess the fishery of Baltic cod, herring and sprat by simulation over 50 years time period. We form a bioeconomic multispecies model for the species. We include species interactions into the model because especially cod and sprat stocks have significant effects on each other. We model the development of population dynamics, catches and profits of the fishery with current fishing mortalities, as well as with the optimal profit maximizing fishing mortalities. Thus, we see how the fishery would develop with current mortalities, and how the fishery should be developed in order to yield maximal profits. Especially cod stock has been quite low recently and by optimizing the fishing mortality it could get recovered. In addition, we assess what would happen to the fisheries of the species if more favourable environmental conditions for cod recruitment dominate in the Baltic Sea. The results may yield new information for the fisheries management. According to the results the fishery of Baltic cod, herring and sprat are not at the most profitable level. The fishing mortalities of each species should be lower in order to maximize the profits. By fishing mortality optimizing the net present value would be almost three times higher in the simulation period. The lower fishing mortality of cod would result in a cod stock recovery. If the environmental conditions in the Baltic Sea improved, cod stock would recover even without a decrease in the fishing mortality. Then the increased cod stock would restrict herring and sprat stock remarkably, and harvesting of these species would not be as profitable anymore.

The roles of nitrification and nitrate reduction pathways in nitrogen cycling of Baltic Sea

The Baltic Sea is one of the largest brackish water bodies in the world. Primary production in the Baltic Sea is limited by nitrogen (N) availability with the exception of river outlets and the northernmost phosphorus limited basin. The excess human induced N load from the drainage basin has caused severe eutrophication of the sea. The excess N loads can be mitigated by microbe mediated natural N removal processes that are found in the oxic-anoxic interfaces in sediments and water column redoxclines. Such interfaces allow the close coupling between the oxic nitrification process, and anoxic denitrification and anaerobic ammonium oxidation (anammox) processes that lead to the formation of molecular nitrogen gas. These processes are governed by various environmental parameters. The effects of these parameters on N processes were investigated in the northern Baltic Sea sediments. During summer months when the sediment organic content is at its highest, nitrification and denitrification reach their maximum rates. However, nitrification had no excess potential, which was probably because of high competition for molecular oxygen (O2) between heterotrophic and nitrification microbes. Subsequently, the limited nitrate (NO3-) availability inhibited denitrification. In fall, winter and spring, nitrification was limited by ammonium availability and denitrification limited by the availability of organic carbon and occasionally by NO3-. Anaerobic ammonium oxidation (anammox) was not an important N removal process in the northern Baltic Sea. Modeling studies suggest that when hypoxia expands in the Baltic Sea, N removal intensifies. However, the results of this study suggest the opposite because bottom water hypoxia (O2< 2 ml l-1) decreased the denitrification rates in sediments. Moreover, N was recycled by the dissimilatory nitrate reduction to ammonium (DNRA) process instead of being removed from the water ecosystem. High N removal potentials were found in the anoxic water column in the deep basins of the Baltic Proper. However, the N removal in the water column appeared to be limited by low substrate availability, because the water at the depths at which the substrate producing nitrification process occurred, rarely mix with the water at the depths at which N removal processes were found. Overall, the natural N removal capacity of the northern Baltic Sea decreased compared to values measured in mid 1990s and early 2000. The reason for this appears to be increasing hypoxia.

Consistency limits behavior of bentonites exposed to sea water

This article examines the changes in interparticle forces brought about on prolonged contact (1 year period) of a bentonite clay with artificial seawater. The study is undertaken with the purpose of identifying the physico-chemical factors that impart a nonswelling character to smectite clays deposited in marine environments. Results show that equilibration of the bentonite clay with artificial seawater (total pore salinity approximately 42 gL-1) for a 1 year period does not lead to any mineralogical changes in the clay specimens; however, their exchangeable cation positions become prominently dominated by magnesium ions. The consistency limits of the seawater-equilibrated bentonite was determined on stepwise leaching to lower salinities. The predominance of diffuse double-layer repulsion forces in the pore salt concentration range of 42 gL-1 to 1.1 gL-1 caused an increase in the liquid limits of the seawater-equilibrated bentonite specimens on reducing the salinity in the corresponding range (42 gL-1 to 1.1 gL-1). The attraction forces, however, prevail over the repulsion forces at salt concentrations <1.1 gL-1 and cause a decrease in liquid limit of the clay specimens with reduction in pore salinity, which is typical of nonswelling clays. The attraction forces cause aggregation of the clay unit layers into domains that break down on sodium saturation of the clay specimens. It is inferred that the physico-chemical factors responsible for the nonswelling character of the seawater-equilibrated bentonite specimens at pore salt concentrations below 1.1 gL-1 are inadequate to explain the nonswelling character of smectite-rich Ariake marine clays. The lower consistency limits of the Ariake marine clays in comparison to the nonswelling character, seawater-equilibrated bentonite specimens is attributed to a relative deficiency of interparticle forces in the Ariake marine clay.

Denitrification in the River Estuaries of the Northern Baltic Sea

Estuaries have been suggested to have an important role in reducing the nitrogen load transported to the sea. We measured denitrification rates in six estuaries of the northern Baltic Sea. Four of them were river mouths in the Bothnian Bay (northern Gulf of Bothnia), and two were estuary bays, one in the Archipelago Sea (southern Gulf of Bothnia) and the other in the Gulf of Finland. Denitrification rates in the four river mouths varied between 330 and 905 mu mol N m(-2) d(-1). The estuary bays at the Archipelago Sea and the Gulf of Bothnia had denitrification rates from 90 mu mol N m(-2) d(-1) to 910 mu mol N m(-2) d(-1) and from 230 mu mol N m(-2) d(-1) to 320 mu mol N m(-2) d(-1), respectively. Denitrification removed 3.6-9.0% of the total nitrogen loading in the river mouths and in the estuary bay in the Gulf of Finland, where the residence times were short. In the estuary bay with a long residence time, in the Archipelago Sea, up to 4.5% of nitrate loading and 19% of nitrogen loading were removed before entering the sea. According to our results, the sediments of the fast-flowing rivers and them estuary areas with short residence times have a limited capacity to reduce the nitrogen load to the Baltic Sea.

Seasonal and short-term variation in denitrification and anammox at a coastal station on the Gulf of Finland, Baltic Sea

Benthic processes were measured at a coastal deposition area in the northern Baltic Sea, covering all seasons. The N-2 production rates, 90-400 mu mol N m(-2) d(-1), were highest in autumn-early winter and lowest in spring. Heterotrophic bacterial production peaked unexpectedly late in the year, indicating that in addition to the temperature, the availability of carbon compounds suitable for the heterotrophic bacteria also plays a major role in regulating the denitrification rate. Anaerobic ammonium oxidation (anammox) was measured in spring and autumn and contributed 10% and 15%, respectively, to the total N-2 production. The low percentage did, however, result in a significant error in the total N-2 production rate estimate, calculated using the isotope pairing technique. Anammox must be taken into account in the Gulf of Finland in future sediment nitrogen cycling research.

Tropical deep convection, convective available potential energy and sea surface temperature

The convective available potential energy (CAFE) based on monthly mean sounding has been shown to be relevant to deep convection in the tropics. The variation of CAFE with SST has been found to be similar to the variation of the frequency of deep convection at one station each in the tropical Atlantic and W. Pacific oceans. This suggests a strong link between the frequency of tropical convection and CAFE. It has been shown that CAFE so derived can be interpreted as the work potential of the atmosphere above the boundary layer with ascent in the convective region and subsidence in the surrounding cloud-free region.

Simulation of ENSO-Related Surface Winds in the Tropical Pacific by an Atmospheric General Circulation Model Forced by Observed Sea Surface Temperatures

The authors present the simulation of the tropical Pacific surface wind variability by a low-resolution (R15 horizontal resolution and 18 vertical levels) version of the Center for Ocean-Land-Atmosphere Interactions, Maryland, general circulation model (GCM) when forced by observed global sea surface temperature. The authors have examined the monthly mean surface winds acid precipitation simulated by the model that was integrated from January 1979 to March 1992. Analyses of the climatological annual cycle and interannual variability over the Pacific are presented. The annual means of the simulated zonal and meridional winds agree well with observations. The only appreciable difference is in the region of strong trade winds where the simulated zonal winds are about 15%-20% weaker than observed, The amplitude of the annual harmonics are weaker than observed over the intertropical convergence zone and the South Pacific convergence zone regions. The amplitudes of the interannual variation of the simulated zonal and meridional winds are close to those of the observed variation. The first few dominant empirical orthogonal functions (EOF) of the simulated, as well as the observed, monthly mean winds are found to contain a targe amount of high-frequency intraseasonal variations, While the statistical properties of the high-frequency modes, such as their amplitude and geographical locations, agree with observations, their detailed time evolution does not. When the data are subjected to a 5-month running-mean filter, the first two dominant EOFs of the simulated winds representing the low-frequency EI Nino-Southern Oscillation fluctuations compare quite well with observations. However, the location of the center of the westerly anomalies associated with the warm episodes is simulated about 15 degrees west of the observed locations. The model simulates well the progress of the westerly anomalies toward the eastern Pacific during the evolution of a warm event. The simulated equatorial wind anomalies are comparable in magnitude to the observed anomalies. An intercomparison of the simulation of the interannual variability by a few other GCMs with comparable resolution is also presented. The success in simulation of the large-scale low-frequency part of the tropical surface winds by the atmospheric GCM seems to be related to the model's ability to simulate the large-scale low-frequency part of the precipitation. Good correspondence between the simulated precipitation and the highly reflective cloud anomalies is seen in the first two EOFs of the 5-month running means. Moreover, the strong correlation found between the simulated precipitation and the simulated winds in the first two principal components indicates the primary role of model precipitation in driving the surface winds. The surface winds simulated by a linear model forced by the GCM-simulated precipitation show good resemblance to the GCM-simulated winds in the equatorial region. This result supports the recent findings that the large-scale part of the tropical surface winds is primarily linear.

Enhanced aerosol loading over Arabian Sea during the pre-monsoon season: Natural or anthropogenic?

[1] Recent experiments conducted over the oceanic regions adjacent to the Indian sub continent have revealed the presence of anthropogenic aerosol haze during January to March. It has been suggested that the major source of this aerosol is South and Southeast Asia. Here we show from long term, multi-station and ship borne observations that aerosols transported from regions northwest of Indian subcontinent especially Arabian and Saharan regions (mostly natural dust) along with the locally produced sea-salt aerosols by sea-surface winds constitute a more significant source of aerosols during April-May period. The radiative forcing due to Arabian/Saharan aerosols (mostly natural) during April May period is comparable and often exceed (as much as 1.5 times) the forcing due to anthropogenic aerosols during January to March period. The presence of dust load over the Arabian Sea can influence the temperature profile and radiative balance in this region.

Spring Warming of the Eastern Arabian Sea and Bay of Bengal from Buoy Data

Observations from moored buoys during spring of 1998-2000 suggest that the warming of the mixed layer (similar to20 m deep) of the north Indian Ocean warm pool is a response to net surface heat flux Q(net) (similar to100 W m(-2)) minus penetrative solar radiation Q(pen) (similar to45 W m(-2)). A residual cooling due to vertical mixing and advection is indirectly estimated to be about 25 W m(-2). The rate of warming due to typical values of Q(net) minus Q(pen) is not very sensitive to the depth of the mixed layer if it lies between 10 m and 30 m.

Modeling Indian Ocean circulation: Bay of Bengal fresh water plume and the Arabian Sea Mini Warm Pool

The Indian subcontinent divides the north Indian Ocean into two tropical basins, namely the Arabian Sea and the Bay of Bengal. The Arabian Sea has high salinity whereas the salinity of the Bay of Bengal is much lower due to the contrast in freshwater forcing of the two basins. The freshwater received by the Bay in large amounts during the summer monsoon through river discharge is flushed out annually by ocean circulation. After the withdrawal of the summer monsoon, the Ganga – Brahmaputra river plume flows first along the Indian coast and then around Sri Lanka into the Arabian Sea creating a low salinity pool in the southeastern Arabian Sea (SEAS). In the same region, during the pre-monsoon months of February – April, a warm pool, known as the Arabian Sea Mini Warm Pool (ASMWP), which is distinctly warmer than the rest of the Indian Ocean, takes shape. In fact, this is the warmest region in the world oceans during this period. Simulation of the river plume and its movement as well as its implications to thermodynamics has been a challenging problem for models of Indian Ocean. Here we address these issues using an ocean general circulation model – first we show that the model is capable of reproducing fresh plumes in the Bay of Bengal as well as its movement and then we use the model to determine the processes that lead to formation of the ASMWP. Hydrographic observations from the western Bay of Bengal have shown the presence of a fresh plume along the northern part of the Indian coast during summer monsoon. The Indian Ocean model when forced by realistic winds and climatological river discharge reproduces the fresh plume with reasonable accuracy. The fresh plume does not advect along the Indian coast until the end of summer monsoon. The North Bay Monsoon Current, which flows eastward in the northern Bay, separates the low salinity water from the more saline southern parts of the bay and thus plays an important role in the fresh water budget of the Bay of Bengal. The model also reproduces the surge of the fresh-plume along the Indian coast, into the Arabian Sea during northeast monsoon. Mechanisms that lead to the formation of the Arabian Sea Mini Warm Pool are investigated using several numerical experiments. Contrary to the existing theories, we find that salinity effects are not necessary for the formation of the ASMWP. The orographic effects of the Sahyadris (Western Ghats) and resulting reduction in wind speed leads to the formation of the ASMWP. During November – April, the SEAS behave as a low-wind heatdominated regime where the evolution of sea surface temperature is solely determined by atmospheric forcing. In such regions the evolution of surface layer temperature is not dependent on the characteristics of the subsurface ocean such as the barrier layer and temperature inversion.

Mechanisms of summer intraseasonal sea surface temperature oscillations in the Bay of Bengal

Intraseasonal variations (ISV) of sea surface temperature (SST) in the Bay of Bengal (BoB) is highest in its northwestern part. An Indian Ocean model forced by QuikSCAT winds and climatological river discharge (QR run) reproduces ISV of SST, albeit with weaker magnitude. Air-sea fluxes, in the presence of a shallow mixed layer, efficiently effect intraseasonal SST fluctuations. Warming during intraseasonal events is smaller (<1°C) for June - July period and larger (1.5° to 2°C) during September, the latter due to a thinner mixed layer. To examine the effect of salinity on ISV, the model was run by artificially increasing the salinity (NORR run) and by decreasing it (MAHA10 run). In NORR, both rainfall and river discharge were switched off and in MAHA10 the discharge by river Mahanadi was increased tenfold. The spatial pattern of ISV as well as its periodicity was similar in QR, NORR and MAHA10. The ISV was stronger in NORR and weaker in MAHA10, compared to QR. In NORR, both intraseasonal warming and cooling were higher than in QR, the former due to reduced air-sea heat loss as the mean SST was lower, and the latter due to enhanced subsurface processes resulting from weaker stratification. In MAHA10, both warming and cooling were lower than in QR, the former due to higher air-sea heat loss owing to higher mean SST, and the latter due to weak subsurface processes resulting from stronger stratification. These model experiments suggest that salinity effects are crucial in determining amplitudes of intraseasonal SST variations in the BoB.