Simulated Bromoform emission and concentration using the ocean biogeochemistry model MPIOM/HAMOCC forced by 6-hourly NCEP data

Bromoform (CHBr3) is one important precursor of atmospheric reactive bromine species that are involved in ozone depletion in the troposphere and stratosphere. In the open ocean bromoform production is linked to phytoplankton that contains the enzyme bromoperoxidase. Coastal sources of bromoform are higher than open ocean sources. However, open ocean emissions are important because the transfer of tracers into higher altitude in the air, i.e. into the ozone layer, strongly depends on the location of emissions. For example, emissions in the tropics are more rapidly transported into the upper atmosphere than emissions from higher latitudes. Global spatio-temporal features of bromoform emissions are poorly constrained. Here, a global three-dimensional ocean biogeochemistry model (MPIOM-HAMOCC) is used to simulate bromoform cycling in the ocean and emissions into the atmosphere using recently published data of global atmospheric concentrations (Ziska et al., 2013) as upper boundary conditions. Our simulated surface concentrations of CHBr3 match the observations well. Simulated global annual emissions based on monthly mean model output are lower than previous estimates, including the estimate by Ziska et al. (2013), because the gas exchange reverses when less bromoform is produced in non-blooming seasons. This is the case for higher latitudes, i.e. the polar regions and northern North Atlantic. Further model experiments show that future model studies may need to distinguish different bromoform-producing phytoplankton species and reveal that the transport of CHBr3 from the coast considerably alters open ocean bromoform concentrations, in particular in the northern sub-polar and polar regions.

(Table 2) Foehn characteristics of station in the McMurdo Dry Valleys, Antarctica

Foehn winds resulting from topographic modification of airflow in the lee of mountain barriers are frequently experienced in the McMurdo Dry Valleys (MDVs) of Antarctica. Strong foehn winds in the MDVs cause dramatic warming at onset and have significant effects on landscape forming processes; however, no detailed scientific investigation of foehn in the MDVs has been conducted. As a result, they are often misinterpreted as adiabatically warmed katabatic winds draining from the polar plateau. Herein observations from surface weather stations and numerical model output from the Antarctic Mesoscale Prediction System (AMPS) during foehn events in the MDVs are presented. Results show that foehn winds in the MDVs are caused by topographic modification of south-southwesterly airflow, which is channeled into the valleys from higher levels. Modeling of a winter foehn event identifies mountain wave activity similar to that associated with midlatitude foehn winds. These events are found to be caused by strong pressure gradients over the mountain ranges of the MDVs related to synoptic-scale cyclones positioned off the coast of Marie Byrd Land. Analysis of meteorological records for 2006 and 2007 finds an increase of 10% in the frequency of foehn events in 2007 compared to 2006, which corresponds to stronger pressure gradients in the Ross Sea region. It is postulated that the intra- and interannual frequency and intensity of foehn events in the MDVs may therefore vary in response to the position and frequency of cyclones in the Ross Sea region.

Investigation of the seasonality of N. pachyderma in the North Atlantic Ocean during Heinrich Stadial 1 using an ecosystem modeling approach

Fossil shells of planktonic foraminifera serve as the prime source of information on past changes in surface ocean conditions. Because the population size of planktonic foraminifera species changes throughout the year, the signal preserved in fossil shells is biased towards the conditions when species production was at its maximum. The amplitude of the potential seasonal bias is a function of the magnitude of the seasonal cycle in production. Here we use a planktonic foraminifera model coupled to an ecosystem model to investigate to what degree seasonal variations in production of the species Neogloboquadrina pachyderma may affect paleoceanographic reconstructions during Heinrich Stadial 1 (~18-15 cal. ka B.P.) in the North Atlantic Ocean. The model implies that during Heinrich Stadial 1 the maximum seasonal production occurred later in the year compared to the Last Glacial Maximum (~21-19 cal. ka B.P.) and the pre-industrial era north of 30 ºN. A diagnosis of the model output indicates that this change reflects the sensitivity of the species to the seasonal cycle of sea-ice cover and food supply, which collectively lead to shifts in the modeled maximum production from the Last Glacial Maximum to Heinrich Stadial 1 by up to six months. Assuming equilibrium oxygen isotopic incorporation in the shells of N. pachyderma, the modeled changes in seasonality would result in an underestimation of the actual magnitude of the meltwater isotopic signal recorded by fossil assemblages of N. pachyderma wherever calcification is likely to take place.

Modeled monthly local topographic lapse rates, Pacific Northwest, USA

An emerging approach to downscaling the projections from General Circulation Models (GCMs) to scales relevant for basin hydrology is to use output of GCMs to force higher-resolution Regional Climate Models (RCMs). With spatial resolution often in the tens of kilometers, however, even RCM output will likely fail to resolve local topography that may be climatically significant in high-relief basins. Here we develop and apply an approach for downscaling RCM output using local topographic lapse rates (empirically-estimated spatially and seasonally variable changes in climate variables with elevation). We calculate monthly local topographic lapse rates from the 800-m Parameter-elevation Regressions on Independent Slopes Model (PRISM) dataset, which is based on regressions of observed climate against topographic variables. We then use these lapse rates to elevationally correct two sources of regional climate-model output: (1) the North American Regional Reanalysis (NARR), a retrospective dataset produced from a regional forecasting model constrained by observations, and (2) a range of baseline climate scenarios from the North American Regional Climate Change Assessment Program (NARCCAP), which is produced by a series of RCMs driven by GCMs. By running a calibrated and validated hydrologic model, the Soil and Water Assessment Tool (SWAT), using observed station data and elevationally-adjusted NARR and NARCCAP output, we are able to estimate the sensitivity of hydrologic modeling to the source of the input climate data. Topographic correction of regional climate-model data is a promising method for modeling the hydrology of mountainous basins for which no weather station datasets are available or for simulating hydrology under past or future climates.

LPJ-LMfire simulations for Europe for the Last Glacial Maximum and preindustrial control, link to netCDF files

This data sets contains LPJ-LMfire dynamic global vegetation model output covering Europe and the Mediterranean for the Last Glacial Maximum (LGM; 21 ka) and for a preindustrial control simulation (20th century detrended climate). The netCDF data files are time averages of the final 30 years of the model simulation. Each netCDF file contains four or five variables: fractional cover of 9 plant functional types (PFTs; cover), total fractional coverage of trees (treecover), population density of hunter-gatherers (foragerPD; only for the "people" simulations), fraction of the gridcell burned on 30-year average (burnedf), and vegetation net primary productivity (NPP). The model spatial resolution is 0.5-degrees For the LGM simulations, LPJ-LMfire was driven by the PMIP3 suite of eight GCMs for which LGM climate simulations were available. Also provided in this archive is the result of an LPJ-LMfire run that was forced by the average climate of all GCMs (the "GCM-mean" files), and the average of each of the individual LPJ-LMfire runs over the eight LGM scenarios individually (the "LPJ-mean" files). The model simulations are provided that include the influence of human presence on the landscape (the "people" files), and in a "world without humans" scenario (the "natural" files). Finally this archive contains the preindustrial reference simulation with and without human influence ("PI_reference_people" and "PI_reference_nat", respectively). There are therefore 22 netCDF files in this archive: 8 each of LGM simulations with and without people (total 16) and the "GCM mean" simulation (2 files) and the "LPJ mean" aggregate (2 files), and finally the two preindustrial "control" simulations ("PI"), with and without humans (2 files). In addition to the LPJ-LMfire model output (netCDF files), this archive also contains a table of arboreal pollen percent calculated from pollen samples dated to the LGM at sites throughout (lgmAP.txt), and a table containing the location of archaeological sites dated to the LGM (LGM_archaeological_site_locations.txt).

Particle flux in the NE Atlantic as recorded by sediment traps

Downward particle flux was measured using sediment traps at various depths over the Porcupine Abyssal Plain (water depth ab. 4850 m) for prolonged periods from 1989 to 1999. A strong seasonal pattern of flux was evident reaching a maximum in mid-summer. The composition of the material changed with depth, reflecting the processes of remineralisation and dissolution as the material sank through the water column. However, there was surprisingly little seasonal variation in its composition to reflect changes in the biology of the euphotic zone. Currents at the site have a strong tidal component with speeds almost always less than 15 cm/sec. In the deeper part of the water column they tend to be northerly in direction, when averaged over periods of several months. A model of upper ocean biogeochemistry forced by meteorology was run for the decade in order to provide an estimate of flux at 3000 m depth. Agreement with measured organic carbon flux is good, both in terms of the timings of the annual peaks and in the integrated annual flux. Interannual variations in the integrated flux are of similar magnitude for both the model output and sediment trap measurements, but there is no significant relationship between these two sets of estimates. No long-term trend in flux is evident, either from the model, or from the measurements. During two spring/summer periods, the marine snow concentration in the water column was assessed by time-lapse photography and showed a strong peak at the start of the downward pulse of material at 3000 m. This emphasises the importance of large particles during periods of maximum flux and at the start of flux peaks. Time lapse photographs of the seabed show a seasonal cycle of coverage of phytodetrital material, in agreement with the model output both in terms of timing and magnitude of coverage prior to 1996. However, after a change in the structure of the benthic community in 1996 no phytodetritus was evident on the seabed. The model output shows only a single peak in flux each year, whereas the measured data usually indicated a double peak. It is concluded that the observed double peak may be a reflection of lowered sediment trap efficiency when flux is very high and is dominated by large marine snow particles. Resuspension into the trap 100 m above the seabed, when compared to the primary flux at 3000 m depth (1800 mab) was lower during periods of high primary flux probably because of a reduction in the height of resuspension when the material is fresh. At 2 mab, the picture is more complex with resuspension being enhanced during the periods of higher flux in 1997, which is consistent with this hypothesis. However there was rather little relationship to flux at 3000 m in 1998. At 3000 m depth, the Flux Stability Index (FSI), which provides a measure of the constancy of the seasonal cycle of flux, exhibited an inverse relationship with flux, such that the highest flux of organic carbon was recorded during the year with the greatest seasonal variation.

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.

Model simulation output data (CCSM3) related to uplift of African topography and the intensification of upwelling in the Benguela region

The Benguela Current, located off the west coast of southern Africa, is tied to a highly productive upwelling system**1. Over the past 12 million years, the current has cooled, and upwelling has intensified**2, 3, 4. These changes have been variously linked to atmospheric and oceanic changes associated with the glaciation of Antarctica and global cooling**5, the closure of the Central American Seaway**1, 6 or the further restriction of the Indonesian Seaway**3. The upwelling intensification also occurred during a period of substantial uplift of the African continent**7, 8. Here we use a coupled ocean-atmosphere general circulation model to test the effect of African uplift on Benguela upwelling. In our simulations, uplift in the East African Rift system and in southern and southwestern Africa induces an intensification of coastal low-level winds, which leads to increased oceanic upwelling of cool subsurface waters. We compare the effect of African uplift with the simulated impact of the Central American Seaway closure9, Indonesian Throughflow restriction10 and Antarctic glaciation**11, and find that African uplift has at least an equally strong influence as each of the three other factors. We therefore conclude that African uplift was an important factor in driving the cooling and strengthening of the Benguela Current and coastal upwelling during the late Miocene and Pliocene epochs.