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.

Model of interactive visualization of diffusion processes from different strains, and source routes and spatial diffusion of capuchin monkeys over the past 6 million years

Source routes and Spatial Diffusion of capuchin monkeys over the past 6 million years, rebuilt in the SPREAD 1.0.6 from the MCC tree. The map shows the 10 different regions to which distinctive samples were associated. The different transmission routes have been calculated from the average rate over time. Only rates with Bayes factor> 3 were considered as significantly different from zero. Significant diffusion pathways are highlighted with color varying from dark brown to red, being the dark brown less significant rates and deep red the most significant rates.

A 20 m spatial resolution seamless multisource Digital Elevation/Depth Model for Lizard Island, northern Great Barrier Reef

Hydrographers have traditionally referred to the nearshore area as the "white ribbon" area due to the challenges associated with the collection of elevation data in this highly dynamic transitional zone between terrestrial and marine environments. Accordingly, available information in this zone is typically characterised by a range of datasets from disparate sources. In this paper we propose a framework to 'fill' the white ribbon area of a coral reef system by integrating multiple elevation and bathymetric datasets acquired by a suite of remote-sensing technologies into a seamless digital elevation model (DEM). A range of datasets are integrated, including field-collected GPS elevation points, terrestrial and bathymetric LiDAR, single and multibeam bathymetry, nautical chart depths and empirically derived bathymetry estimations from optical remote sensing imagery. The proposed framework ranks data reliability internally, thereby avoiding the requirements to quantify absolute error and results in a high resolution, seamless product. Nested within this approach is an effective spatially explicit technique for improving the accuracy of bathymetry estimates derived empirically from optical satellite imagery through modelling the spatial structure of residuals. The approach was applied to data collected on and around Lizard Island in northern Australia. Collectively, the framework holds promise for filling the white ribbon zone in coastal areas characterised by similar data availability scenarios. The seamless DEM is referenced to the horizontal coordinate system MGA Zone 55 - GDA 1994, mean sea level (MSL) vertical datum and has a spatial resolution of 20 m.

Age model and color measurements of sediment cores from the eastern equatorial Pacific

High-resolution records of glacial-interglacial variations in biogenic carbonate, opal, and detritus (derived from non-destructive core log measurements of density, P-wave velocity and color; r >= 0.9) from 15 sediment sites in the eastern equatorial (sampling resolution is ~1 kyr) clear response to eccentricity and precession forcing. For the Peru Basin, we generate a high-resolution (21 kyr increment) orbitally-based chronology for the last 1.3 Ma. Spectral analysis indicates that the 100 kyr cycle became dominant at roughly 1.2 Ma, 200-300 kyr earlier than reported for other paleoclimatic records. The response to orbital forcing is weaker since the Mid-Brunhes Dissolution Event (at 400 ka). A west-east reconstruction of biogenic sedimentation in the Peru Basin (four cores; 91-85°W) distinguishes equatorial and coastal upwelling systems in the western and eastern sites, respectively. A north-south reconstruction perpendicular to the equatorial upwelling system (11 cores, 11°N-°3S) shows high carbonate contents (>= 50%) between 6°N and 4°S and highly variable opal contents between 2°N and 4°S. Carbonate cycles B-6, B-8, B-10, B-12, B-14, M-2, and M-6 are well developed with B-10 (430 ka) as the most prominent cycle. Carbonate highs during glacials and glacial-interglacial transitions extended up to 400 km north and south compared to interglacial or interglacial^glacial carbonate lows. Our reconstruction thus favors glacial-interglacial expansion and contraction of the equatorial upwelling system rather than shifting north or south. Elevated accumulation rates are documented near the equator from 6°N to 4°S and from 2°N to 4°S for carbonate and opal, respectively. Accumulation rates are higher during glacials and glacial-interglacial transitions in all cores, whereas increased dissolution is concentrated on Peru Basin sediments close to the carbonate compensation depth and occurred during interglacials or interglacial-glacial transitions.

Ocean temperature response to idealized Gleissberg and de Vries solar cycles in a comprehensive climate model

The ~90-year Gleissberg and ~200-year de Vries cycles have been identified as two distinctive quasi-periodic components of Holocene solar activity. Evidence exists for the impact of such multi-decadal to centennial-scale variability in total solar irradiance (TSI) on climate, but concerning the ocean, this evidence is mainly restricted to the surface response. Here we use a comprehensive global climate model to study the impact of idealized solar forcing, representing the Gleissberg and de Vries cycles, on global ocean potential temperature at different depth levels, after a recent proxy record indicates a signal of TSI anomalies in the northeastern Atlantic at mid-depth. Potential impacts of TSI anomalies on deeper oceanic levels are climatically relevant due to their possible effect on ocean circulation by altering water mass characteristics. Simulated solar anomalies are shown to penetrate the ocean down to at least deep-water levels. Despite the fact that the two forcing periods differ only by a factor of ~2, the spatial pattern of response is significantly distinctive between the experiments, suggesting different mechanisms for solar signal propagation. These are related to advection by North Atlantic Deep Water flow (200-year forcing), and barotropic adjustment in the South Atlantic in response to a latitudinal shift of the westerly wind belt (90-year forcing).

Modeling the sensitivity of precipitation oxygen isotopes to the Last Glacial Maximum boundary conditions: A study using Community Atmosphere Model (CAM3.0)

To understand the validity of d18O proxy records as indicators of past temperature change, a series of experiments was conducted using an atmospheric general circulation model fitted with water isotope tracers (Community Atmosphere Model version 3.0, IsoCAM). A pre-industrial simulation was performed as the control experiment, as well as a simulation with all the boundary conditions set to Last Glacial Maximum (LGM) values. Results from the pre-industrial and LGM simulations were compared to experiments in which the influence of individual boundary conditions (greenhouse gases, ice sheet albedo and topography, sea surface temperature (SST), and orbital parameters) were changed each at a time to assess their individual impact. The experiments were designed in order to analyze the spatial variations of the oxygen isotopic composition of precipitation (d18Oprecip) in response to individual climate factors. The change in topography (due to the change in land ice cover) played a significant role in reducing the surface temperature and d18Oprecip over North America. Exposed shelf areas and the ice sheet albedo reduced the Northern Hemisphere surface temperature and d18Oprecip further. A global mean cooling of 4.1 °C was simulated with combined LGM boundary conditions compared to the control simulation, which was in agreement with previous experiments using the fully coupled Community Climate System Model (CCSM3). Large reductions in d18Oprecip over the LGM ice sheets were strongly linked to the temperature decrease over them. The SST and ice sheet topography changes were responsible for most of the changes in the climate and hence the d18Oprecip distribution among the simulations.

Digital surface model, hillshade and Lambertian reflectance model of the rock glaciers Oelgrube and Aeusseres Hochebenkar (Oetztal Alps, Tyrol, Austria)

With full-waveform (FWF) lidar systems becoming increasingly available from different commercial manufacturers, the possibility for extracting physical parameters of the scanned surfaces in an area-wide sense, as addendum to their geometric representation, has risen as well. The mentioned FWF systems digitize the temporal profiles of the transmitted laser pulse and of its backscattered echoes, allowing for a reliable determination of the target distance to the instrument and of physical target quantities by means of radiometric calibration, one of such quantities being the diffuse Lambertian reflectance. The delineation of glaciers is a time-consuming task, commonly performed manually by experts and involving field trips as well as image interpretation of orthophotos, digital terrain models and shaded reliefs. In this study, the diffuse Lambertian reflectance was compared to the glacier outlines mapped by experts. We start the presentation with the workflow for analysis of FWF data, their direct georeferencing and the calculation of the diffuse Lambertian reflectance by radiometric calibration; this workflow is illustrated for a large FWF lidar campaign in the Ötztal Alps (Tyrol, Austria), operated with an Optech ALTM 3100 system. The geometric performance of the presented procedure was evaluated by means of a relative and an absolute accuracy assessment using strip differences and orthophotos, resp. The diffuse Lambertian reflectance was evaluated at two rock glaciers within the mentioned lidar campaign. This feature showed good performance for the delineation of the rock glacier boundaries, especially at their lower parts.

Digital elevation model (DEM) of the ice sheet in western Dronning Maud Land, Antarctica

In this paper, a new digital elevation model (DEM) is derived for the ice sheet in western Dronning Maud Land, Antarctica. It is based on differential interferometric synthetic aperture radar (SAR) from the European Remote Sensing 1/2 (ERS-1/2) satellites, in combination with ICESat's Geoscience Laser Altimeter System (GLAS). A DEM mosaic is compiled out of 116 scenes from the ERS-1 ice phase in 1994 and the ERS-1/2 tandem mission between 1996 and 1997 with the GLAS data acquired in 2003 that served as ground control. Using three different SAR processors, uncertainties in phase stability and baseline model, resulting in height errors of up to 20 m, are exemplified. Atmospheric influences at the same order of magnitude are demonstrated, and corresponding scenes are excluded. For validation of the DEM mosaic, covering an area of about 130,000 km**2 on a 50-m grid, independent ICESat heights (2004-2007), ground-based kinematic GPS (2005), and airborne laser scanner data (ALS, 2007) are used. Excluding small areas with low phase coherence, the DEM differs in mean and standard deviation by 0.5 +/- 10.1, 1.1 +/- 6.4, and 3.1 +/- 4.0 m from ICESat, GPS, and ALS, respectively. The excluded data points may deviate by more than 50 m. In order to suppress the spatially variable noise below a 5-m threshold, 18% of the DEM area is selectively averaged to a final product at varying horizontal spatial resolution. Apart from mountainous areas, the new DEM outperforms other currently available DEMs and may serve as a benchmark for future elevation models such as from the TanDEM-X mission to spatially monitor ice sheet elevation.

Lagrangian back-tracing results for Kohnen station; output of time dependent numerical experiment with a 3D nested Antarctic ice sheet model

A nested ice flow model was developed for eastern Dronning Maud Land to assist with the dating and interpretation of the EDML deep ice core. The model consists of a high-resolution higher-order ice dynamic flow model that was nested into a comprehensive 3-D thermomechanical model of the whole Antarctic ice sheet. As the drill site is on a flank position the calculations specifically take into account the effects of horizontal advection as deeper ice in the core originated from higher inland. First the regional velocity field and ice sheet geometry is obtained from a forward experiment over the last 8 glacial cycles. The result is subsequently employed in a Lagrangian backtracing algorithm to provide particle paths back to their time and place of deposition. The procedure directly yields the depth-age distribution, surface conditions at particle origin, and a suite of relevant parameters such as initial annual layer thickness. This paper discusses the method and the main results of the experiment, including the ice core chronology, the non-climatic corrections needed to extract the climatic part of the signal, and the thinning function. The focus is on the upper 89% of the ice core (appr. 170 kyears) as the dating below that is increasingly less robust owing to the unknown value of the geothermal heat flux. It is found that the temperature biases resulting from variations of surface elevation are up to half of the magnitude of the climatic changes themselves.

Planimetric ice-flow convergence and flow-orthonormal strain rate across the Antarctic Ice Sheet, with links to model result files

Fast-flowing ice streams discharge most of the ice from the interior of the Antarctic Ice Sheet coastward. Understanding how their tributary organisation is governed and evolves is essential for developing reliable models of the ice sheet's response to climate change. Despite much research on ice-stream mechanics, this problem is unsolved, because the complexity of flow within and across the tributary networks has hardly been interrogated. Here I present the first map of planimetric flow convergence across the ice sheet, calculated from satellite measurements of ice surface velocity, and use it to explore this complexity. The convergence map of Antarctica elucidates how ice-stream tributaries draw ice from the interior. It also reveals curvilinear zones of convergence along lateral shear margins of streaming, and abundant convergence ripples associated with nonlinear ice rheology and changes in bed topography and friction. Flow convergence on ice-stream tributaries and their feeding zones is markedly uneven, and interspersed with divergence at distances of the order of kilometres. For individual drainage basins as well as the ice sheet as a whole, the range of convergence and divergence decreases systematically with flow speed, implying that fast flow cannot converge or diverge as much as slow flow. I therefore deduce that flow in ice-stream networks is subject to mechanical regulation that limits flow-orthonormal strain rates. These properties and the gridded data of convergence and flow-orthonormal strain rate in this archive provide targets for ice- sheet simulations and motivate more research into the origin and dynamics of tributarization.

Global LAI uncertainty and relative uncertainty datasets derived from a Triple Collocation Error Model (TCEM)

Uncertainty information for global leaf area index (LAI) products is important for global modeling studies but usually difficult to systematically obtain at a global scale. Here, we present a new method that cross-validates existing global LAI products and produces consistent uncertainty information. The method is based on a triple collocation error model (TCEM) that assumes errors among LAI products are not correlated. Global monthly absolute and relative uncertainties, in 0.05° spatial resolutions, were generated for MODIS, CYCLOPES, and GLOBCARBON LAI products, with reasonable agreement in terms of spatial patterns and biome types. CYCLOPES shows the lowest absolute and relative uncertainties, followed by GLOBCARBON and MODIS. Grasses, crops, shrubs, and savannas usually have lower uncertainties than forests in association with the relatively larger forest LAI. With their densely vegetated canopies, tropical regions exhibit the highest absolute uncertainties but the lowest relative uncertainties, the latter of which tend to increase with higher latitudes. The estimated uncertainties of CYCLOPES generally meet the quality requirements (± 0.5) proposed by the Global Climate Observing System (GCOS), whereas for MODIS and GLOBCARBON only non-forest biome types have met the requirement. Nevertheless, none of the products seems to be within a relative uncertainty requirements of 20%. Further independent validation and comparative studies are expected to provide a fair assessment of uncertainties derived from TCEM. Overall, the proposed TCEM is straightforward and could be automated for the systematic processing of real time remote sensing observations to provide theoretical uncertainty information for a wider range of land products.

Model results of sensitivity experiments for marine nitrogen cycle in four NetCDF files

The marine nitrogen (N) inventory is thought to be stabilized by negative feedback mechanisms that reduce N inventory excursions relative to the more slowly overturning phosphorus inventory. Using a global biogeochemical ocean circulation model we show that negative feedbacks stabilizing the N inventory cannot persist if a close spatial association of N2 fixation and denitrification occurs. In our idealized model experiments, nitrogen deficient waters, generated by denitrification, stimulate local N2 fixation activity. But, because of stoichiometric constraints, the denitrification of newly fixed nitrogen leads to a net loss of N. This can enhance the N deficit, thereby triggering additional fixation in a vicious cycle, ultimately leading to a runaway N loss. To break this vicious cycle, and allow for stabilizing negative feedbacks to occur, inputs of new N need to be spatially decoupled from denitrification. Our idealized model experiments suggest that factors such as iron limitation or dissolved organic matter cycling can promote such decoupling and allow for negative feedbacks that stabilize the N inventory. Conversely, close spatial co-location of N2 fixation and denitrification could lead to net N loss.

Transient Holocene experiments under orbital forcing using the comprehensive global climate model CCSM3 (Community Climate System Model 3)

The Southern Hemisphere Westerly Winds (SWW) have been suggested to exert a critical influence on global climate through wind-driven upwelling of deep water in the Southern Ocean and the potentially resulting atmospheric CO2 variations. The investigation of the temporal and spatial evolution of the SWW along with forcings and feedbacks remains a significant challenge in climate research. In this study, the evolution of the SWW under orbital forcing from the early Holocene (9 kyr BP) to pre-industrial modern times is examined with transient experiments using the comprehensive coupled global climate model CCSM3. Analyses of the model results suggest that the annual and seasonal mean SWW were subject to an overall strengthening and poleward shifting trend during the course of the early-to-late Holocene under the influence of orbital forcing, except for the austral spring season, where the SWW exhibited an opposite trend of shifting towards the equator.

Age model of firn core DML05C98_32 (B32)

In the framework of the European Project for Ice Coring in Antarctica (EPICA), a comprehensive glaciological pre-site survey has been carried out on Amundsenisen, Dronning Maud Land, East Antarctica, in the past decade. Within this survey, four intermediate-depth ice cores and 13 snow pits were analyzed for their ionic composition and interpreted with respect to the spatial and temporal variability of volcanic sulphate deposition. The comparison of the non-sea-salt (nss)-sulphate peaks that are related to the well-known eruptions of Pinatubo and Cerro Hudson in AD 1991 revealed sulphate depositions of comparable size (15.8 ± 3.4 kg/km**2) in 11 snow pits. There is a tendency to higher annual concentrations for smaller snow-accumulation rates. The combination of seasonal sodium and annually resolved nss-sulphate records allowed the establishment of a time-scale derived by annual-layer counting over the last 2000 years and thus a detailed chronology of annual volcanic sulphate deposition. Using a robust outlier detection algorithm, 49 volcanic eruptions were identified between AD 165 and 1997. The dating uncertainty is ±3 years between AD 1997 and 1601, around ±5 years between AD 1601 and 1257, and increasing to ±24 years at AD 165, improving the accuracy of the volcanic chronology during the penultimate millennium considerably.

Temporal and spatial LGM variability of biomarker records in the Arctic Ocean

Using the sea ice proxy IP25 and phytoplankton-derived biomarkers (brassicasterol and dinosterol) Arctic sea-ice conditions were reconstructed for Marine Isotope Stage (MIS) 3 to 1 in sediment cores from the north of Barents Sea continental margin across the Central Arctic to the Southern Mendeleev Ridge. Our results suggest more extensive sea-ice cover than present-day during MIS 3, increasing sea-ice growth during MIS 2 and decreased sea-ice cover during the last deglacial. The summer ice edge sustained north of the Barents Sea even during extremely cold (i.e., Last Glacial Maximum (LGM)) as well as warm periods (i.e., Bølling-Allerød). During the LGM, the western Svalbard margin and the northern Barents Sea margin areas were characterized by high concentrations of both IP25 and phytoplankton biomarkers, interpreted as a productive ice-edge situation, caused by the inflow of warm Atlantic Water. In contrast, the LGM high Arctic proper (north of 84°N) was covered by thick permanent sea ice throughout the year with rare break up, indicated by zero or near-zero biomarker concentrations. The spring/summer sea-ice margin significantly extended southwards to the southern Lomonosov Ridge and Mendeleev Ridge during the LGM. Our proxy reconstructions are very consistent with published model results based on the North Atlantic/Arctic Ocean Sea Ice Model (NAOSIM).