Medem Channel profiles based on stand-alone model simulation and corrections with nudging

State-of-the-art process-based models have shown to be applicable to the simulation and prediction of coastal morphodynamics. On annual to decadal temporal scales, these models may show limitations in reproducing complex natural morphological evolution patterns, such as the movement of bars and tidal channels, e.g. the observed decadal migration of the Medem Channel in the Elbe Estuary, German Bight. Here a morphodynamic model is shown to simulate the hydrodynamics and sediment budgets of the domain to some extent, but fails to adequately reproduce the pronounced channel migration, due to the insufficient implementation of bank erosion processes. In order to allow for long-term simulations of the domain, a nudging method has been introduced to update the model-predicted bathymetries with observations. The model-predicted bathymetry is nudged towards true states in annual time steps. Sensitivity analysis of a user-defined correlation length scale, for the definition of the background error covariance matrix during the nudging procedure, suggests that the optimal error correlation length is similar to the grid cell size, here 80-90 m. Additionally, spatially heterogeneous correlation lengths produce more realistic channel depths than do spatially homogeneous correlation lengths. Consecutive application of the nudging method compensates for the (stand-alone) model prediction errors and corrects the channel migration pattern, with a Brier skill score of 0.78. The proposed nudging method in this study serves as an analytical approach to update model predictions towards a predefined 'true' state for the spatiotemporal interpolation of incomplete morphological data in long-term simulations.

Planktic foraminiferal Mg/Ca and oxygen isotopes, plus estimated d18Oseawater and ice volume corrected d18Oseawater for the past 35000 years of sediment core GeoB10069-3

Evidence from geologic archives suggests that there were large changes in the tropical hydrologic cycle associated with the two prominent northern hemisphere deglacial cooling events, Heinrich Stadial 1 (HS1; ~19 to 15 kyr BP; kyr BP = 1000 yr before present) and the Younger Dryas (~12.9 to 11.7 kyr BP). These hydrologic shifts have been alternatively attributed to high and low latitude origin. Here, we present a new record of hydrologic variability based on planktic foraminifera-derived d18O of seawater (d18Osw) estimates from a sediment core from the tropical Eastern Indian Ocean, and using 12 additional d18Osw records, construct a single record of the dominant mode of tropical Eastern Equatorial Pacific and Indo-Pacific Warm Pool (IPWP) hydrologic variability. We show that deglacial hydrologic shifts parallel variations in the reconstructed interhemispheric temperature gradient, suggesting a strong response to variations in the Atlantic Meridional Overturning Circulation and the attendant heat redistribution. A transient model simulation of the last deglaciation suggests that hydrologic changes, including a southward shift in the Intertropical Convergence Zone (ITCZ) which likely occurred during these northern hemisphere cold events, coupled with oceanic advection and mixing, resulted in increased salinity in the Indonesian region of the IPWP and the eastern tropical Pacific, which is recorded by the d18Osw proxy. Based on our observations and modeling results we suggest the interhemispheric temperature gradient directly controls the tropical hydrologic cycle on these time scales, which in turn mediates poleward atmospheric heat transport.

Temperature-dependent phenology of Plutella xylostella (Lepidoptera: Plutellidae)

Temperature-dependent population growth of diamondback moth (DBM) Plutella xylostella (L.), a prolific insect pest of crucifer vegetables, was studied under six constant temperatures in the laboratory. The objective of the study was to predict the impacts of temperature changes on the population of DBM at high-resolution scales along altitudinal gradients and under climate change scenarios. Non-linear functions were fitted on the data for modeling the development, mortality, longevity and oviposition of the pest. The best-fitted functions for each life stage were compiled for estimating the life table parameters of the species by stochastic simulations. To quantify the impacts on the pest, three indices (establishment, generation and activity) were computed using the estimates of life table parameters and temperature data obtained at local scale (current scenario 2013) and downscaled climate change data (future scenario 2055) from the AFRICLIM database. To measure and represent the impacts of temperature change along the altitude on the pest; the indices were mapped along the altitudinal gradients of Kilimanjaro and Taita Hills, in Tanzania and Kenya, respectively. Potential impact of the changes between climate scenarios 2013 and 2055 was assessed. The data files included in this database were utilized for the above analysis to develop temperature dependent phenology of Plutella xylostella to assess current and future distribution along eastern African Afromontanes.

High-resolution in situ oxygen microprofiles, porewater and solid phase geochemistry from the Crimean shelf (Black Sea) from Maria S. Merian cruise MSM15/1

The outer western Crimean shelf of the Black Sea is a natural laboratory to investigate effects of stable oxic versus varying hypoxic conditions on seafloor biogeochemical processes and benthic community structure. Bottom-water oxygen concentrations ranged from normoxic (175 µmol O2/L) and hypoxic (< 63 µmol O2/L) or even anoxic/sulfidic conditions within a few kilometers' distance. Variations in oxygen concentrations between 160 and 10 µmol/L even occurred within hours close to the chemocline at 134 m water depth. Total oxygen uptake, including diffusive as well as fauna-mediated oxygen consumption, decreased from 15 mmol/m**2/d on average in the oxic zone, to 7 mmol/m**2/d on average in the hypoxic zone, correlating with changes in macrobenthos composition. Benthic diffusive oxygen uptake rates, comprising respiration of microorganisms and small meiofauna, were similar in oxic and hypoxic zones (on average 4.5 mmol/m**2/d), but declined to 1.3 mmol/m**2/d in bottom waters with oxygen concentrations below 20 µmol/L. Measurements and modeling of porewater profiles indicated that reoxidation of reduced compounds played only a minor role in diffusive oxygen uptake under the different oxygen conditions, leaving the major fraction to aerobic degradation of organic carbon. Remineralization efficiency decreased from nearly 100 % in the oxic zone, to 50 % in the oxic-hypoxic zone, to 10 % in the hypoxic-anoxic zone. Overall, the faunal remineralization rate was more important, but also more influenced by fluctuating oxygen concentrations, than microbial and geochemical oxidation processes.

Chemical and Sr isotopic compositions of pore fluids and marine carbonates from ODP Hole 130-807A

The 87Sr/86Sr ratios and Sr concentrations in sediment and pore fluids are used to evaluate the rates of calcite recrystallization at ODP Site 807A on the Ontong Java Plateau, an 800-meter thick section of carbonate ooze and chalk. A numerical model is used to evaluate the pore fluid chemistry and Sr isotopes in an accumulating section. The deduced calcite recrystallization rate is 2% per million years (%/Myr) near the top of the section and decreases systematically in older parts of the section such that the rate is close to 0.1/age (in years). The deduced recrystallization rates have important implications for the interpretation of Ca and Mg concentration profiles in the pore fluids. The effect of calcite recrystallization on pore fluid chemistry is described by the reaction length, L, which varies by element, and depends on the concentration in pore fluid and solid. When L is small compared to the thickness of the sedimentary section, the pore fluid concentration is controlled by equilibrium or steady-state exchange with the solid phase, except within a distance L of the sediment-water interface. When L is large relative to the thickness of sediment, the pore fluid concentration is mostly controlled by the boundary conditions and diffusion. The values of L for Ca, Sr, and Mg are of order 15, 150, and 1500 meters, respectively. L_Sr is derived from isotopic data and modeling, and allows us to infer the values of L_Ca and L_Mg. The small value for L_Ca indicates that pore fluid Ca concentrations, which gradually increase down section, must be equilibrium values that are maintained by solution-precipitation exchange with calcite and do not reflect Ca sources within or below the sediment column. The pore fluid Ca measurements and measured alkalinity allow us to calculate the in situ pH in the pore fluids, which decreases from 7.6 near the sediment-water interface to 7.1+/-0.1 at 400-800 mbsf. While the calculated pH values are in agreement with some of the values measured during ODP Leg 130, most of the measurements are artifacts. The large value for L_Mg indicates that the pore fluid Mg concentrations at 807A are not controlled by calcite-fluid equilibrium but instead are determined by the changing Mg concentration of seawater during deposition, modified by aqueous diffusion in the pore fluids. We use the pore fluid Mg concentration profile at Site 807A to retrieve a global record for seawater Mg over the past 35 Myr, which shows that seawater Mg has increased rapidly over the past 10 Myr, rather than gradually over the past 60 Myr. This observation suggests that the Cenozoic rise in seawater Mg is controlled by continental weathering inputs rather than by exchange with oceanic crust. The relationship determined between reaction rate and age in silicates and carbonates is strikingly similar, which suggests that reaction affinity is not the primary determinant of silicate dissolution rates in nature.

Sable isotope record and sediment composition of Cretaceous samples of Deremara Rise

We estimate tropical Atlantic upper ocean temperatures using oxygen isotope and Mg/Ca ratios in well-preserved planktonic foraminifera extracted from Albian through Santonian black shales recovered during Ocean Drilling Program Leg 207 (North Atlantic Demerara Rise). On the basis of a range of plausible assumptions regarding seawater composition at the time the data support temperatures between 33° and 42°C. In our low-resolution data set spanning ~84-100 Ma a local temperature maximum occurs in the late Turonian, and a possible minimum occurs in the mid to early late Cenomanian. The relation between single species foraminiferal d18O and Mg/Ca suggests that the ratio of magnesium to calcium in the Turonian-Coniacian ocean may have been lower than in the Albian-Cenomanian ocean, perhaps coincident with an ocean 87Sr/86Sr minimum. The carbon isotopic compositions of distinct marine algal biomarkers were measured in the same sediment samples. The d13C values of phytane, combined with foraminiferal d13C and inferred temperatures, were used to estimate atmospheric carbon dioxide concentrations through this interval. Estimates of atmospheric CO2 concentrations range between 600 and 2400 ppmv. Within the uncertainty in the various proxies, there is only a weak overall correspondence between higher (lower) tropical temperatures and more (less) atmospheric CO2. The GENESIS climate model underpredicts tropical Atlantic temperatures inferred from ODP Leg 207 foraminiferal d18O and Mg/Ca when we specify approximate CO2 concentrations estimated from the biomarker isotopes in the same samples. Possible errors in the temperature and CO2 estimates and possible deficiencies in the model are discussed. The potential for and effects of substantially higher atmospheric methane during Cretaceous anoxic events, perhaps derived from high fluxes from the oxygen minimum zone, are considered in light of recent work that shows a quadratic relation between increased methane flux and atmospheric CH4 concentrations. With 50 ppm CH4, GENESIS sea surface temperatures approximate the minimum upper ocean temperatures inferred from proxy data when CO2 concentrations specified to the model are near those inferred using the phytane d13C proxy. However, atmospheric CO2 concentrations of 3500 ppm or more are still required in the model in order to reproduce inferred maximum temperatures.

Temperature, water level and bathymetry of thermokarst lakes in the continuous permafrost zone of northern Siberia - Lena River Delta, Siberia

Thermokarst lakes are typical features of the northern permafrost ecosystems, and play an important role in the thermal exchange between atmosphere and subsurface. The objective of this study is to describe the main thermal processes of the lakes and to quantify the heat exchange with the underlying sediments. The thermal regimes of five lakes located within the continuous permafrost zone of northern Siberia (Lena River Delta) were investigated using hourly water temperature and water level records covering a 3-year period (2009-2012), together with bathymetric survey data. The lakes included thermokarst lakes located on Holocene river terraces that may be connected to Lena River water during spring flooding, and a thermokarst lake located on deposits of the Pleistocene Ice Complex. Lakes were covered by ice up to 2 m thick that persisted for more than 7 months of the year, from October until about mid-June. Lake-bottom temperatures increased at the start of the ice-covered period due to upward-directed heat flux from the underlying thawed sediment. Prior to ice break-up, solar radiation effectively warmed the water beneath the ice cover and induced convective mixing. Ice break-up started at the beginning of June and lasted until the middle or end of June. Mixing occurred within the entire water column from the start of ice break-up and continued during the ice-free periods, as confirmed by the Wedderburn numbers, a quantitative measure of the balance between wind mixing and stratification that is important for describing the biogeochemical cycles of lakes. The lake thermal regime was modeled numerically using the FLake model. The model demonstrated good agreement with observations with regard to the mean lake temperature, with a good reproduction of the summer stratification during the ice-free period, but poor agreement during the ice-covered period. Modeled sensitivity to lake depth demonstrated that lakes in this climatic zone with mean depths > 5 m develop continuous stratification in summer for at least 1 month. The modeled vertical heat flux across the bottom sediment tends towards an annual mean of zero, with maximum downward fluxes of about 5 W/m**2 in summer and with heat released back into the water column at a rate of less than 1 W/m**2 during the ice-covered period. The lakes are shown to be efficient heat absorbers and effectively distribute the heat through mixing. Monthly bottom water temperatures during the ice-free period range up to 15 °C and are therefore higher than the associated monthly air or ground temperatures in the surrounding frozen permafrost landscape. The investigated lakes remain unfrozen at depth, with mean annual lake-bottom temperatures of between 2.7 and 4 °C.

Orbital-resolution Pliocene-Pleistocene simulations of CO2, sea level, global temperature and benthic d18O, and d11B-based proxy-CO2 data

During the past five million yrs, benthic d18O records indicate a large range of climates, from warmer than today during the Pliocene Warm Period to considerably colder during glacials. Antarctic ice cores have revealed Pleistocene glacial-interglacial CO2 variability of 60-100 ppm, while sea level fluctuations of typically 125 m are documented by proxy data. However, in the pre-ice core period, CO2 and sea level proxy data are scarce and there is disagreement between different proxies and different records of the same proxy. This hampers comprehensive understanding of the long-term relations between CO2, sea level and climate. Here, we drive a coupled climate-ice sheet model over the past five million years, inversely forced by a stacked benthic d18O record. We obtain continuous simulations of benthic d18O, sea level and CO2 that are mutually consistent. Our model shows CO2 concentrations of 300 to 470 ppm during the Early Pliocene. Furthermore, we simulate strong CO2 variability during the Pliocene and Early Pleistocene. These features are broadly supported by existing and new d11B-based proxy CO2 data, but less by alkenone-based records. The simulated concentrations and variations therein are larger than expected from global mean temperature changes. Our findings thus suggest a smaller Earth System Sensitivity than previously thought. This is explained by a more restricted role of land ice variability in the Pliocene. The largest uncertainty in our simulation arises from the mass balance formulation of East Antarctica, which governs the variability in sea level, but only modestly affects the modeled CO2 concentrations.

Field measurements of the atmosphere, ocean, sea ice and sub-ice platelet layer at Atka Bay in 2013

Ice shelves strongly interact with coastal Antarctic sea ice and the associated ecosystem by creating conditions favourable to the formation of a sub-ice platelet layer. The close investigation of this phenomenon and its seasonal evolution remain a challenge due to logistical constraints and a lack of suitable methodology. In this study, we characterize the seasonal cycle of Antarctic fast ice adjacent to the Ekström Ice Shelf in the eastern Weddell Sea. We used a thermistor chain with the additional ability to record the temperature response induced by cyclic heating of resistors embedded in the chain. Vertical sea-ice temperature and heating profiles obtained daily between November 2012 and February 2014 were analyzed to determine sea-ice and snow evolution, and to calculate the basal energy budget. The residual heat flux translated into an ice-volume fraction in the platelet layer of 0.18 ± 0.09, which we reproduced by a independent model simulation and agrees with earlier results. Manual drillings revealed an average annual platelet-layer thickness increase of at least 4m, and an annual maximum thickness of 10m beneath second-year sea ice. The oceanic contribution dominated the total sea-ice production during the study, effectively accounting for up to 70% of second-year sea-ice growth. In summer, an oceanic heat flux of 21 W/m**2 led to a partial thinning of the platelet layer. Our results further show that the active heating method, in contrast to the acoustic sounding approach, is well suited to derive the fast-ice mass balance in regions influenced by ocean/ice-shelf interaction, as it allows sub-diurnal monitoring of the platelet-layer thickness.