An individual based model spatial distribution, production and biomass for Calanus finmarchicus in the Norwegian Sea

The copepod Calanus finmarchicus is the dominant species of the meso-zooplankton in the Norwegian Sea, and constitutes an important link between the phytoplankton and the higher trophic levels in the Norwegian Sea food chain. An individualbased model for C. finmarchicus, based on super-individuals and evolving traits for behaviour, stages, etc., is two-way coupled to the NORWegian ECOlogical Model system (NORWECOM). One year of modelled C. finmarchicus spatial distribution, production and biomass are found to represent observations reasonably well. High C. finmarchicus abundance is found along the Norwegian shelf-break in the early summer, while the overwintering population is found along the slope and in the deeper Norwegian Sea basins. The timing of the spring bloom is generally later than in the observations. Annual Norwegian Sea production is found to be 29 million tonnes of carbon and a production to biomass (P/B) ratio of 4.3 emerges. Sensitivity tests show that the modelling system is robust to initial values of behavioural traits and with regards to the number of super-individuals simulated given that this is above about 50,000 individuals. Experiments with the model system indicate that it provides a valuable tool for studies of ecosystem responses to causative forces such as prey density or overwintering population size. For example, introducing C. finmarchicus food limitations reduces the stock dramatically, but on the other hand, a reduced stock may rebuild in one year under normal conditions. The NetCDF file contains model grid coordinates and bottom topography.

Salinization MASS Model Version

This salinization simulation enables irrigation salinity, capillary rise, rainfall, leaching, salt mitigation strategies via fallowing, and other functions that are addressed in a coupled social-environment model applied to southern Mesopotamia. The simulation is applied to the Modeling Ancient Settlement Systems Project supported by the University of Chicago, Durham University, and Argonne National Laboratory. The simulation can be used for other regions.

Benthic cover map of Heron Reef derived from a high-spatial-resolution multi-spectral satellite image using object based image analysis

Providing accurate maps of coral reefs where the spatial scale and labels of the mapped features correspond to map units appropriate for examining biological and geomorphic structures and processes is a major challenge for remote sensing. The objective of this work is to assess the accuracy and relevance of the process used to derive geomorphic zone and benthic community zone maps for three western Pacific coral reefs produced from multi-scale, object-based image analysis (OBIA) of high-spatial-resolution multi-spectral images, guided by field survey data. Three Quickbird-2 multi-spectral data sets from reefs in Australia, Palau and Fiji and georeferenced field photographs were used in a multi-scale segmentation and object-based image classification to map geomorphic zones and benthic community zones. A per-pixel approach was also tested for mapping benthic community zones. Validation of the maps and comparison to past approaches indicated the multi-scale OBIA process enabled field data, operator field experience and a conceptual hierarchical model of the coral reef environment to be linked to provide output maps at geomorphic zone and benthic community scales on coral reefs. The OBIA mapping accuracies were comparable with previously published work using other methods; however, the classes mapped were matched to a predetermined set of features on the reef.

Tracer release experiments in the tropical North Atlantic

Lateral diffusivity is computed from a tracer release experiment in the northeastern tropical Atlantic thermocline. The uncertainties of the estimates are inferred from a synthetic particle release using a high-resolution ocean circulation model. The main method employed to compute zonal and meridional components of lateral diffusivity is the growth of the second moment of a cloud of tracer. The application of an areal comparison method for estimating tracer-based diffusivity in the field experiments is also discussed. The best estimate of meridional eddy diffusivity in the Guinea Upwelling region at about 300 m depth is estimated to be inline image m2 s-1. The zonal component of lateral diffusivity is estimated to be inline image m2 s-1, while areal comparison method yields areal equivalent zonal diffusivity component of inline image m2 s?1. In comparison to Ky, Kx is about twice larger, resulting from the tracer patch stretching by zonal jets. Employed conceptual jet model indicates that zonal jet velocities of about inline image m s?1 are required to explain the enhancement of the zonal eddy diffusivity component. Finally, different sampling strategies are tested on synthetic tracer release experiments. They indicate that the best sampling strategy is a sparse regular sampling grid covering most of the tracer patch.

Chemical concentrations and fluxes from EPICA ice cores EDML and EDC

Continuous sea salt and mineral dust aerosol records have been studied on the two EPICA (European Project for Ice Coring in Antarctica) deep ice cores. The joint use of these records from opposite sides of the East Antarctic plateau allows for an estimate of changes in dust transport and emission intensity as well as for the identification of regional differences in the sea salt aerosol source. The mineral dust flux records at both sites show a strong coherency over the last 150 kyr related to dust emission changes in the glacial Patagonian dust source with three times higher dust fluxes in the Atlantic compared to the Indian Ocean sector of the Southern Ocean (SO). Using a simple conceptual transport model this indicates that transport can explain only 40% of the atmospheric dust concentration changes in Antarctica, while factor 5-10 changes occurred. Accordingly, the main cause for the strong glacial dust flux changes in Antarctica must lie in environmental changes in Patagonia. Dust emissions, hence environmental conditions in Patagonia, were very similar during the last two glacials and interglacials, respectively, despite 2-4 °C warmer temperatures recorded in Antarctica during the penultimate interglacial than today. 2-3 times higher sea salt fluxes found in both ice cores in the glacial compared to the Holocene are difficult to reconcile with a largely unchanged transport intensity and the distant open ocean source. The substantial glacial enhancements in sea salt aerosol fluxes can be readily explained assuming sea ice formation as the main sea salt aerosol source with a significantly larger expansion of (summer) sea ice in the Weddell Sea than in the Indian Ocean sector. During the penultimate interglacial, our sea salt records point to a 50% reduction of winter sea ice coverage compared to the Holocene both in the Indian and Atlantic Ocean sector of the SO. However, from 20 to 80 ka before present sea salt fluxes show only very subdued millennial changes despite pronounced temperature fluctuations, likely due to the large distance of the sea ice salt source to our drill sites.

Results of hydrological and geological investigations in Benin

Focus of this study is the analysis of a local hydrogeological system in the subhumid outer tropics in the western African country of Benin. The aim was to characterize, qualify and quantify the hydrogeological and hydrological properties of the approx. 30 km2 big study area and to develop a conceptual hydrogeological model. This model should provide the basis for further studies on a regional scale. The main goal was to obtain the process knowledge of the hydrogeological system and to determine the process and the quantity of the groundwater recharge in the working area. According to the objectives, a broad hydrogeological approach was chosen. In a spacious network on the local scale TDR probes, suction cups and groundwater observation bores were installed. Also in a multidisciplinary cooperation with hydrology, geography, soil science, biology, meteorology and plant nutrition sciences, instruments like discharge gauging stations, tensiometers, lysimeter, climate stations, runoff plots and erosion pins were installed in the test site for the investigation of the relevant parameters of the hydrological cycle.

(Table T1) Coral taxa and their ecological significance in ODP Site 194-1196 and Hole 194-1199A sediments

The Southern Marion Plateau (SMP) represents a vertical stacking of Miocene carbonate platform deposits. Two sites (1196 and 1199) were drilled on top of this plateau, penetrating a 663-m carbonate succession of bioclastic and reefal sedimentary bodies. The study focuses on the least dolomitized 410-m-thick upper part of the succession, which is middle to late Miocene in age. Sedimentological and paleontological studies were conducted at both sites in order to propose a paleoenvironmental model and its evolution through the Miocene age. Six main microfacies of possible environmental significance were defined using statistical multivariate analyses, based on the recognition and point counting of 24 biogenic components. Depositional environment reconstructions are proposed as well as the biosedimentary and the environmental evolution regarding seismic architectures, stratigraphy, biosedimentology, and microfacies analysis. The SMP platform mainly results from a vertical stacking of lens-shaped bodies in homoclinal to distally steepened ramp settings.

Euphausiid respiration model revamped, link to model results

Euphausiids constitute major biomass component in shelf ecosystems and play a fundamental role in the rapid vertical transport of carbon from the ocean surface to the deeper layers during their daily vertical migration (DVM). DVM depth and migration patterns depend on oceanographic conditions with respect to temperature, light and oxygen availability at depth, factors that are highly dependent on season in most marine regions. Changes in the abiotic conditions also shape Euphausiid metabolism including aerobic and anaerobic energy production. Here we introduce a global krill respiration model which includes the effect of latitude (LAT), the day of the year of interest (DoY), and the number of daylight hours on the day of interest (DLh), in addition to the basal variables that determine ectothermal oxygen consumption (temperature, body mass and depth) in the ANN model (Artificial Neural Networks). The newly implemented parameters link space and time in terms of season and photoperiod to krill respiration. The ANN model showed a better fit (r**2=0.780) when DLh and LAT were included, indicating a decrease in respiration with increasing LAT and decreasing DLh. We therefore propose DLh as a potential variable to consider when building physiological models for both hemispheres. We also tested for seasonality the standard respiration rate of the most common species that were investigated until now in a large range of DLh and DoY with Multiple Linear Regression (MLR) or General Additive model (GAM). GAM successfully integrated DLh (r**2= 0.563) and DoY (r**2= 0.572) effects on respiration rates of the Antarctic krill, Euphausia superba, yielding the minimum metabolic activity in mid-June and the maximum at the end of December. Neither the MLR nor the GAM approach worked for the North Pacific krill Euphausia pacifica, and MLR for the North Atlantic krill Meganyctiphanes norvegica remained inconclusive because of insufficient seasonal data coverage. We strongly encourage comparative respiration measurements of worldwide Euphausiid key species at different seasons to improve accuracy in ecosystem modelling.

Model experiment on optimal size of a fish stock with environmental uncertainties

We analyze the effect of environmental uncertainties on optimal fishery management in a bio-economic fishery model. Unlike most of the literature on resource economics, but in line with ecological models, we allow the different biological processes of survival and recruitment to be affected differently by environmental uncertainties. We show that the overall effect of uncertainty on the optimal size of a fish stock is ambiguous, depending on the prudence of the value function. For the case of a risk-neutral fishery manager, the overall effect depends on the relative magnitude of two opposing effects, the 'convex-cost effect' and the 'gambling effect'. We apply the analysis to the Baltic cod and the North Sea herring fisheries, concluding that for risk neutral agents the net effect of environmental uncertainties on the optimal size of these fish stocks is negative, albeit small in absolute value. Under risk aversion, the effect on optimal stock size is positive for sufficiently high coefficients of constant relative risk aversion.

Lake Baikal amphipods under climate change: Thermal constraints and ecological consequences, links to supplementary material

Lake Baikal, the world's most voluminous freshwater lake, has experienced unprecedented warming during the last decades. A uniquely diverse amphipod fauna inhabits the littoral zone and can serve as a model system to identify the role of thermal tolerance under climate change. This study aimed to identify sublethal thermal constraints in two of the most abundant endemic Baikal amphipods, Eulimnogammarus verrucosus and Eulimnogammarus cyaneus, and Gammarus lacustris, a ubiquitous gammarid of the Holarctic. As the latter is only found in some shallow isolated bays of the lake, we further addressed the question whether rising temperatures could promote the widespread invasion of this non-endemic species into the littoral zone. Animals were exposed to gradual temperature increases (4 week, 0.8 °C/d; 24 h, 1 °C/h) starting from the reported annual mean temperature of the Baikal littoral (6 °C). Within the framework of oxygen- and capacity-limited thermal tolerance (OCLTT), we used a nonlinear regression approach to determine the points at which the changing temperature-dependence of relevant physiological processes indicates the onset of limitation. Limitations in ventilation representing the first limits of thermal tolerance (pejus (= "getting worse") temperatures (Tp)) were recorded at 10.6 (95% confidence interval; 9.5, 11.7), 19.1 (17.9, 20.2), and 21.1 (19.8, 22.4) °C in E. verrucosus, E. cyaneus, and G. lacustris, respectively. Field observations revealed that E. verrucosus retreated from the upper littoral to deeper and cooler waters once its Tp was surpassed, identifying Tp as the ecological thermal boundary. Constraints in oxygen consumption at higher than critical temperatures (Tc) led to an exponential increase in mortality in all species. Exposure to short-term warming resulted in higher threshold values, consistent with a time dependence of thermal tolerance. In conclusion, species-specific limits to oxygen supply capacity are likely key in the onset of constraining (beyond pejus) and then life-threatening (beyond critical) conditions. Ecological consequences of these limits are mediated through behavioral plasticity in E. verrucosus. However, similar upper thermal limits in E. cyaneus (endemic, Baikal) and G. lacustris (ubiquitous, Holarctic) indicate that the potential invader G. lacustris would not necessarily benefit from rising temperatures. Secondary effects of increasing temperatures remain to be investigated.

Simulation Model of Benthic Antarctic Assemblages (SIMBAA), link to PC-executable

SIMBAA is a spatially explicit, individual-based simulation model. It was developed to analyse the response of populations of Antarctic benthic species and their diversity to iceberg scouring. This disturbance is causing a high local mortality providing potential space for new colonisation. Traits can be attributed to model species, e.g. in terms of reproduction, dispersal, and life span. Physical disturbances can be designed in space and time, e.g. in terms of size, shape, and frequency. Environmental heterogeneity can be considered by cell-specific capacities to host a certain number of individuals. When grid cells become empty (after a disturbance event or due to natural mortality of of an individual), a lottery decides which individual from which species stored in a pool of candidates (for this cell) will recruit in that cell. After a defined period the individuals become mature and their offspring are dispersed and stored in the pool of candidates. The biological parameters and disturbance regimes decide on how long an individual lives. Temporal development of single populations of species as well as Shannon diversity are depicted in the main window graphically and primary values are listed. Examples for simulations can be loaded and saved as sgf-files. The results are also shown in an additional window in a dimensionless area with 50 x 50 cells, which contain single individuals depicted as circles; their colour indicates the assignment to the self-designed model species and the size represents their age. Dominant species per cell and disturbed areas can also be depicted. Output of simulation runs can be saved as images, which can be assembled to video-clips by standard computer programs (see GIF-examples of which "Demo 1" represents the response of the Antarctic benthos to iceberg scouring and "Demo 2" represents a simulation of a deep-sea benthic habitat).

Terrain model of the Håkon Mosby Mud Volcano (10 m grid size)

The Håkon Mosby Mud Volcano is a natural laboratory to study geological, geochemical, and ecological processes related to deep-water mud volcanism. High resolution bathymetry of the Håkon Mosby Mud Volcano was recorded during RV Polarstern expedition ARK-XIX/3 utilizing the multibeam system Hydrosweep DS-2. Dense spacing of the survey lines and slow ship speed (5 knots) provided necessary point density to generate a regular 10 m grid. Generalization was applied to preserve and represent morphological structures appropriately. Contour lines were derived showing detailed topography at the centre of the Håkon Mosby Mud Volcano and generalized contours in the vicinity. We provide a brief introduction to the Håkon Mosby Mud Volcano area and describe in detail data recording and processing methods, as well as the morphology of the area. Accuracy assessment was made to evaluate the reliability of a 10 m resolution terrain model. Multibeam sidescan data were recorded along with depth measurements and show reflectivity variations from light grey values at the centre of the Håkon Mosby Mud Volcano to dark grey values (less reflective) at the surrounding moat.