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).

Gravity potential spherical harmonic series

Time variable gravity fields, reflecting variations of mass distribution in the system Earth is one of the key parameters to understand the changing Earth. Mass variations are caused either by redistribution of mass in, on or above the Earth's surface or by geophysical processes in the Earth's interior. The first set of observations of monthly variations of the Earth gravity field was provided by the US/German GRACE satellite mission beginning in 2002. This mission is still providing valuable information to the science community. However, as GRACE has outlived its expected lifetime, the geoscience community is currently seeking successor missions in order to maintain the long time series of climate change that was begun by GRACE. Several studies on science requirements and technical feasibility have been conducted in the recent years. These studies required a realistic model of the time variable gravity field in order to perform simulation studies on sensitivity of satellites and their instrumentation. This was the primary reason for the European Space Agency (ESA) to initiate a study on ''Monitoring and Modelling individual Sources of Mass Distribution and Transport in the Earth System by Means of Satellites''. The goal of this interdisciplinary study was to create as realistic as possible simulated time variable gravity fields based on coupled geophysical models, which could be used in the simulation processes in a controlled environment. For this purpose global atmosphere, ocean, continental hydrology and ice models were used. The coupling was performed by using consistent forcing throughout the models and by including water flow between the different domains of the Earth system. In addition gravity field changes due to solid Earth processes like continuous glacial isostatic adjustment (GIA) and a sudden earthquake with co-seismic and post-seismic signals were modelled. All individual model results were combined and converted to gravity field spherical harmonic series, which is the quantity commonly used to describe the Earth's global gravity field. The result of this study is a twelve-year time-series of 6-hourly time variable gravity field spherical harmonics up to degree and order 180 corresponding to a global spatial resolution of 1 degree in latitude and longitude. In this paper, we outline the input data sets and the process of combining these data sets into a coherent model of temporal gravity field changes. The resulting time series was used in some follow-on studies and is available to anybody interested.

Atmospheric, oceanic and hydrological polar motion excitation mechanisms based on a combination of geometric and gravimetric space observations

The goal of our study is to determine accurate time series of geophysical Earth rotation excitations to learn more about global dynamic processes in the Earth system. For this purpose, we developed an adjustment model which allows to combine precise observations from space geodetic observation systems, such as Satellite Laser Ranging (SLR), Global Navigation Satellite Systems (GNSS), Very Long Baseline Interferometry (VLBI), Doppler Orbit determination and Radiopositioning Integrated on Satellite (DORIS), satellite altimetry and satellite gravimetry in order to separate geophysical excitation mechanisms of Earth rotation. Three polar motion time series are applied to derive the polar motion excitation functions (integral effect). Furthermore we use five time variable gravity field solutions from Gravity Recovery and Climate Experiment (GRACE) to determine not only the integral mass effect but also the oceanic and hydrological mass effects by applying suitable filter techniques and a land-ocean mask. For comparison the integral mass effect is also derived from degree 2 potential coefficients that are estimated from SLR observations. The oceanic mass effect is also determined from sea level anomalies observed by satellite altimetry by reducing the steric sea level anomalies derived from temperature and salinity fields of the oceans. Due to the combination of all geodetic estimated excitations the weaknesses of the individual processing strategies can be reduced and the technique-specific strengths can be accounted for. The formal errors of the adjusted geodetic solutions are smaller than the RMS differences of the geophysical model solutions. The improved excitation time series can be used to improve the geophysical modeling.

Barium measured on water bottle samples during POLARSTERN cruise ARK-XXII/2 (SPACE)

Dissolved barium has been shown to have the potential to distinguish Eurasian from North American (NA) river runoff. As part of the ARK-XXII/2 Polarstern expedition in summer 2007, Ba was analyzed in the Barents, Kara, Laptev seas, and the Eurasian Basins as well as the Makarov Basin up to the Alpha and Mendeleyev Ridges. By combining salinity, d18O and initial phosphate corrected for mineralization with oxygen (PO4*) or N/P ratios we identified the water mass fractions of meteoric water, sea ice meltwater, and marine waters of Atlantic as well as Pacific origin in the upper water column. In all basins inside the lower halocline layer and the Arctic intermediate waters we find Ba concentrations close to those of the Fram Strait branch of the lower halocline (41-45 nM), reflecting the composition of the incoming Atlantic water. A layer of upper halocline water (UHW) with higher Ba concentrations (45-55 nM) is identified in the Makarov Basin. Atop of the UHW, the Surface Mixed Layer (SML), including the summer and winter mixed layers, has high concentrations of Ba (58-67 nM). In the SML of the investigated area of the central Arctic the meteoric fraction can be identified by assuming a conservative behavior of Ba to be primarily of Eurasian river origin. However, in productive coastal regions biological removal compromises the use of Ba to distinguish between Eurasian and NA rivers. As a consequence, the NA river water fraction is underestimated in productive surface waters or waters that have passed a productive region, whereas this fraction is overestimated in subsurface waters containing remineralised Ba, particularly when these waters have passed productive shelf regions. Especially in the Laptev Sea and small regions in the Barents Sea, Ba concentrations are low in surface waters. In the Laptev Sea exceptionally high Ba concentrations in shelf bottom waters indicate that Ba is removed from surface waters to deep waters by biological activity enhanced by increasing ice-free conditions as well as by scavenging by organic matter of terrestrial origin. We interpret high Ba concentrations in the UHW of the Makarov Basin to result from enrichment by remineralisation in bottom waters on the shelf of the Chukchi Sea and therefore the calculated NA runoff is an artefact. We conclude that no NA runoff can be demonstrated unequivocally anywhere during our expedition with the set of tracers considered here. Small contributions of NA runoff may have been masked by Ba depletion and could only be resolved by supportive tracers on the uptake history. We thus suggest that Ba has to be used with care as it can put limits but not yield quantitative water mass distributions. Only if the extra Ba inputs exceed the cumulative biological uptake the signal can be unequivocally attributed to NA runoff.