Impact of oceanic processes on the carbon cycle during the last termination

During the last termination (from ~18 000 years ago to ~9000 years ago), the climate significantly warmed and the ice sheets melted. Simultaneously, atmospheric CO2 increased from ~190 ppm to ~260 ppm. Although this CO2 rise plays an important role in the deglacial warming, the reasons for its evolution are difficult to explain. Only box models have been used to run transient simulations of this carbon cycle transition, but by forcing the model with data constrained scenarios of the evolution of temperature, sea level, sea ice, NADW formation, Southern Ocean vertical mixing and biological carbon pump. More complex models (including GCMs) have investigated some of these mechanisms but they have only been used to try and explain LGM versus present day steady-state climates. In this study we use a coupled climate-carbon model of intermediate complexity to explore the role of three oceanic processes in transient simulations: the sinking of brines, stratification-dependent diffusion and iron fertilization. Carbonate compensation is accounted for in these simulations. We show that neither iron fertilization nor the sinking of brines alone can account for the evolution of CO2, and that only the combination of the sinking of brines and interactive diffusion can simultaneously simulate the increase in deep Southern Ocean δ13C. The scenario that agrees best with the data takes into account all mechanisms and favours a rapid cessation of the sinking of brines around 18 000 years ago, when the Antarctic ice sheet extent was at its maximum. In this scenario, we make the hypothesis that sea ice formation was then shifted to the open ocean where the salty water is quickly mixed with fresher water, which prevents deep sinking of salty water and therefore breaks down the deep stratification and releases carbon from the abyss. Based on this scenario, it is possible to simulate both the amplitude and timing of the long-term CO2 increase during the last termination in agreement with ice core data. The atmospheric δ13C appears to be highly sensitive to changes in the terrestrial biosphere, underlining the need to better constrain the vegetation evolution during the termination.

A rate and state friction law for saline ice

Sea ice friction models are necessary to predict the nature of interactions between sea ice floes. These interactions are of interest on a range of scales, for example, to predict loads on engineering structures in icy waters or to understand the basin-scale motion of sea ice. Many models use Amonton's friction law due to its simplicity. More advanced models allow for hydrodynamic lubrication and refreezing of asperities; however, modeling these processes leads to greatly increased complexity. In this paper we propose, by analogy with rock physics, that a rate- and state-dependent friction law allows us to incorporate memory (and thus the effects of lubrication and bonding) into ice friction models without a great increase in complexity. We support this proposal with experimental data on both the laboratory (∼0.1 m) and ice tank (∼1 m) scale. These experiments show that the effects of static contact under normal load can be incorporated into a friction model. We find the parameters for a first-order rate and state model to be A = 0.310, B = 0.382, and μ0 = 0.872. Such a model then allows us to make predictions about the nature of memory effects in moving ice-ice contacts.

Evaluation of biospheric components in Earth system models using modern and palaeo-observations: the state-of-the-art

Earth system models (ESMs) are increasing in complexity by incorporating more processes than their predecessors, making them potentially important tools for studying the evolution of climate and associated biogeochemical cycles. However, their coupled behaviour has only recently been examined in any detail, and has yielded a very wide range of outcomes. For example, coupled climate–carbon cycle models that represent land-use change simulate total land carbon stores at 2100 that vary by as much as 600 Pg C, given the same emissions scenario. This large uncertainty is associated with differences in how key processes are simulated in different models, and illustrates the necessity of determining which models are most realistic using rigorous methods of model evaluation. Here we assess the state-of-the-art in evaluation of ESMs, with a particular emphasis on the simulation of the carbon cycle and associated biospheric processes. We examine some of the new advances and remaining uncertainties relating to (i) modern and palaeodata and (ii) metrics for evaluation. We note that the practice of averaging results from many models is unreliable and no substitute for proper evaluation of individual models. We discuss a range of strategies, such as the inclusion of pre-calibration, combined process- and system-level evaluation, and the use of emergent constraints, that can contribute to the development of more robust evaluation schemes. An increasingly data-rich environment offers more opportunities for model evaluation, but also presents a challenge. Improved knowledge of data uncertainties is still necessary to move the field of ESM evaluation away from a "beauty contest" towards the development of useful constraints on model outcomes.

Direct probing of photoinduced electron transfer in a self- assembled biomimetic [2Fe2S]-hydrogenase complex using ultrafast vibrational spectroscopy

A pyridyl-functionalized diiron dithiolate complex, [μ-(4-pyCH2−NMI-S2)Fe2(CO)6] (3, py = pyridine(ligand), NMI = naphthalene monoimide) was synthesized and fully characterized. In the presence of zinc tetraphenylporphyrin (ZnTPP), a self-assembled 3·ZnTPP complex was readily formed in CH2Cl2 by the coordination of the pyridyl nitrogen to the porphyrin zinc center. Ultrafast photoinduced electron transfer from excited ZnTPP to complex 3 in the supramolecular assembly was observed in real time by monitoring the ν(CO) and ν(CO)NMI spectral changes with femtosecond time-resolved infrared (TRIR) spectroscopy. We have confirmed that photoinduced charge separation produced the monoreduced species by comparing the time-resolved IR spectra with the conventional IR spectra of 3•− generated by reversible electrochemical reduction. The lifetimes for the charge separation and charge recombination processes were found to be τCS = 40 ± 3 ps and τCR = 205 ± 14 ps, respectively. The charge recombination is much slower than that in an analogous covalent complex, demonstrating the potential of a supramolecular approach to extend the lifetime of the chargeseparated state in photocatalytic complexes. The observed vibrational frequency shifts provide a very sensitive probe of the delocalization of the electron-spin density over the different parts of the Fe2S2 complex. The TR and spectro-electrochemical IR spectra, electron paramagnetic resonance spectra, and density functional theory calculations all show that the spin density in 3•− is delocalized over the diiron core and the NMI bridge. This delocalization explains why the complex exhibits low catalytic dihydrogen production even though it features a very efficient photoinduced electron transfer. The ultrafast porphyrin-to-NMIS2−Fe2(CO)6 photoinduced electron transfer is the first reported example of a supramolecular Fe2S2-hydrogenase model studied by femtosecond TRIR spectroscopy. Our results show that TRIR spectroscopy is a powerful tool to investigate photoinduced electron transfer in potential dihydrogen-producing catalytic complexes, and that way to optimize their performance by rational approaches.

Multi‐method global sensitivity analysis (MMGSA) for modelling floodplain hydrological processes

When studying hydrological processes with a numerical model, global sensitivity analysis (GSA) is essential if one is to understand the impact of model parameters and model formulation on results. However, different definitions of sensitivity can lead to a difference in the ranking of importance of the different model factors. Here we combine a fuzzy performance function with different methods of calculating global sensitivity to perform a multi-method global sensitivity analysis (MMGSA). We use an application of a finite element subsurface flow model (ESTEL-2D) on a flood inundation event on a floodplain of the River Severn to illustrate this new methodology. We demonstrate the utility of the method for model understanding and show how the prediction of state variables, such as Darcian velocity vectors, can be affected by such a MMGSA. This paper is a first attempt to use GSA with a numerically intensive hydrological model.

Time-dependent flows in the coupled solar wind-magnetosphere-ionosphere system

Concepts of time-dependent flow in the coupled solar wind-magnetosphere-ionosphere system are discussed and compared with the frequently-adopted steady-state paradigm. Flows are viewed as resulting from departures of the system from equilibrium excited by dayside and nightside reconnection processes, with the flows then taking the system back towards a new equilibrium configuration. The response of the system to reconnection impulses, continuous but unbalanced reconnection and balanced steady-state reconnection are discussed in these terms. It is emphasized that in the time-dependent case the ionospheric and interplanetary electric fields are generally inductively decoupled from each other; a simple mapping of the interplanetary electric field along equipotential field lines into the ionosphere occurs only in the electrostatic steady-state case.

State-resolved gas-surface reactivity of methane in the symmetric C-H stretch vibration on Ni(100)

The state-resolved reactivity of CH4 in its totally symmetric C-H stretch vibration (�1) has been measured on a Ni(100) surface. Methane molecules were accelerated to kinetic energies of 49 and 63:5 kJ=mol in a molecular beam and vibrationally excited to �1 by stimulated Raman pumping before surface impact at normal incidence. The reactivity of the symmetric-stretch excited CH4 is about an order of magnitude higher than that of methane excited to the antisymmetric stretch (�3) reported by Juurlink et al. [Phys. Rev. Lett. 83, 868 (1999)] and is similar to that we have previously observed for the excitation of the first overtone (2�3). The difference between the state-resolved reactivity for �1 and �3 is consistent with predictions of a vibrationally adiabatic model of the methane reaction dynamics and indicates that statistical models cannot correctly describe the chemisorption of CH4 on nickel.

The slowly evolving background state of the atmosphere

The theory of wave–mean flow interaction requires a partition of the atmospheric flow into a notional background state and perturbations to it. Here, a background state, known as the Modified Lagrangian Mean (MLM), is defined as the zonally symmetric state obtained by requiring that every potential vorticity (PV) contour lying within an isentropic layer encloses the same mass and circulation as in the full flow. For adiabatic and frictionless flow, these two integral properties are time-invariant and the MLM state is a steady solution of the primitive equations. The time dependence in the adiabatic flow is put into the perturbations, which can be described by a wave-activity conservation law that is exact even at large amplitude. Furthermore, the effects of non-conservative processes on wave activity can be calculated from the conservation law. A new method to calculate the MLM state is introduced, where the position of the lower boundary is obtained as part of the solution. The results are illustrated using Northern Hemisphere ERA-Interim data. The MLM state evolves slowly, implying that the net non-conservative effects are weak. Although ‘adiabatic eddy fluxes’ cannot affect the MLM state, the effects of Rossby-wave breaking, PV filamentation and subsequent dissipation result in sharpening of the polar vortex edge and meridional shifts in the MLM zonal flow, both at tropopause level and on the winter stratospheric vortex. The rate of downward migration of wave activity during stratospheric sudden warmings is shown to be given by the vertical scale associated with polar vortex tilt divided by the time-scale for wave dissipation estimated from the wave-activity conservation law. Aspects of troposphere–stratosphere interaction are discussed. The new framework is suitable to examine the climate and its interactions with disturbances, such as midlatitude storm tracks, and makes a clean partition between adiabatic and non-conservative processes.

Multi-method global sensitivity analysis (MMGSA) for modelling floodplain hydrological processes

When studying hydrological processes with a numerical model, global sensitivity analysis (GSA) is essential if one is to understand the impact of model parameters and model formulation on results. However, different definitions of sensitivity can lead to a difference in the ranking of importance of the different model factors. Here we combine a fuzzy performance function with different methods of calculating global sensitivity to perform a multi-method global sensitivity analysis (MMGSA). We use an application of a finite element subsurface flow model (ESTEL-2D) on a flood inundation event on a floodplain of the River Severn to illustrate this new methodology. We demonstrate the utility of the method for model understanding and show how the prediction of state variables, such as Darcian velocity vectors, can be affected by such a MMGSA. This paper is a first attempt to use GSA with a numerically intensive hydrological model

Atmosphere-ocean coupled processes in the Madden-Julian oscillation

The Madden-Julian oscillation (MJO) is a convectively coupled 30-70 day (intraseasonal) tropical atmospheric mode that drives variations in global weather, but which is poorly simulated in most atmospheric general circulation models. Over the past two decades, field campaigns and modeling experiments have suggested that tropical atmosphere-ocean interactions may sustain or amplify the pattern of enhanced and suppressed atmospheric convection that defines the MJO, and encourage its eastward propagation through the Indian and Pacific Oceans. New observations collected during the past decade have advanced our understand of the ocean response to atmospheric MJO forcing and the resulting intraseasonal sea surface temperature (SST) fluctuations. Numerous modeling studies have revealed a considerable impact of the mean state on MJO ocean-atmosphere coupled processes, as well as the importance of resolving the diurnal cycle of atmosphere--upper-ocean interactions. New diagnostic methods provide insight to atmospheric variability and physical processes associated with the MJO, but offer limited insight on the role of ocean feedbacks. Consequently, uncertainty remains concerning the role of the ocean in MJO theory. Our understanding of how atmosphere-ocean coupled processes affect the MJO can be improved by collecting observations in poorly sampled regions of MJO activity, assessing oceanic and atmospheric drivers of surface fluxes, improving the representation of upper-ocean mixing in coupled-model simulations, designing model experiments that minimize mean-state differences, and developing diagnostic tools to evaluate the nature and role of coupled ocean-atmosphere processes over the MJO cycle.

The impact of westerly wind bursts and ocean initial state on the development, and diversity of El Niño events

Westerly wind bursts (WWBs) that occur in the western tropical Pacific are believed to play an important role in the development of El Niño events. Here, following the study of Lengaigne et al. (Clim Dyn 23(6):601–620, 2004), we conduct numerical simulations in which we reexamine the response of the climate system to an observed wind burst added to a coupled general circulation model. Two sets of twin ensemble experiments are conducted (each set has control and perturbed experiments). In the first set, the initial ocean heat content of the system is higher than the model climatology (recharged), while in the second set it is nearly normal (neutral). For the recharged state, in the absence of WWBs, a moderate El Niño with a maximum warming in the central Pacific (CP) develops in about a year. In contrast, for the neutral state, there develops a weak La Niña. However, when the WWB is imposed, the situation dramatically changes: the recharged state slides into an El Niño with a maximum warming in the eastern Pacific, while the neutral set produces a weak CP El Niño instead of previous La Niña conditions. The different response of the system to the exact same perturbations is controlled by the initial state of the ocean and the subsequent ocean–atmosphere interactions involving the interplay between the eastward shift of the warm pool and the warming of the eastern equatorial Pacific. Consequently, the observed diversity of El Niño, including the occurrence of extreme events, may depend on stochastic atmospheric processes, modulating El Niño properties within a broad continuum.

The sensitivity of convective aggregation to diabatic processes in idealized radiative-convective equilibrium simulations

Idealized explicit convection simulations of the Met Office Unified Model exhibit spontaneous self-aggregation in radiative-convective equilibrium, as seen in other models in previous studies. This self-aggregation is linked to feedbacks between radiation, surface fluxes, and convection, and the organization is intimately related to the evolution of the column water vapor field. Analysis of the budget of the spatial variance of column-integrated frozen moist static energy (MSE), following Wing and Emanuel [2014], reveals that the direct radiative feedback (including significant cloud longwave effects) is dominant in both the initial development of self-aggregation and the maintenance of an aggregated state. A low-level circulation at intermediate stages of aggregation does appear to transport MSE from drier to moister regions, but this circulation is mostly balanced by other advective effects of opposite sign and is forced by horizontal anomalies of convective heating (not radiation). Sensitivity studies with either fixed prescribed radiative cooling, fixed prescribed surface fluxes, or both do not show full self-aggregation from homogeneous initial conditions, though fixed surface fluxes do not disaggregate an initialized aggregated state. A sensitivity study in which rain evaporation is turned off shows more rapid self-aggregation, while a run with this change plus fixed radiative cooling still shows strong self-aggregation, supporting a “moisture memory” effect found in Muller and Bony [2015]. Interestingly, self-aggregation occurs even in simulations with sea surface temperatures (SSTs) of 295 K and 290 K, with direct radiative feedbacks dominating the budget of MSE variance, in contrast to results in some previous studies.

Utilisation of by-products from sunflower-based biodiesel production processes for the production of fermentation feedstock

By-products streams from a sunflower-based biodiesel plant were utilised for the production of fermentation media that can be used for the production of polyhydroxyalkanoates (PHA). Sunflower meal was utilised as substrate for the production of crude enzyme consortia through solid state fermentation (SSF) with the fungal strain Aspergillus oryzae. Fermented solids were subsequently mixed with unprocessed sunflower meal aiming at the production of a nutrient-rich fermentation feedstock. The highest free amino nitrogen (FAN) and inorganic phosphorus concentrations achieved were 1.5 g L-1 and 246 mg L-1, respectively, when an initial proteolytic activity of 6.4 U mL-1 was used. The FANconcentrationwas increased to 2.3 g L-1 when the initial proteolytic activity was increased to 16 U mL-1. Sunflower meal hydrolysates were mixed with crude glycerol to provide fermentationmedia that were evaluated for the production of poly(3-hydroxybutyrateco- 3-hydroxyvalerate) (P(3HB-co-3HV)) using Cupriavidus necator DSM545. The P(3HB-co-3HV) (9.9 g l-1) produced contained 3HB and 3HV units with 97 and 3 mol %, respectively. Integrating PHA production in existing 1st generation biodiesel production plants through valorisation of by-product streams could improve their sustainability.

Two case studies of sulfate scavenging processes in the Amazon region (Rondonia)

The scavenging processes of chemical species have been previously studied with numerical modeling, in order to understand the gas and particulate matter intra-reservoir transferences. In this study, the atmospheric (RAMS) and scavenging (B.V.2) models were used, in order to simulate sulfate concentrations in rainwater using scavenging processes as well as the local atmospheric conditions obtained within the LBA Project in the State of Rondonia, during a dry-to-wet transition season. Two case studies were conducted. The RAMS atmospheric simulation of these events presented satisfactory results, showing the detailed microphysical processes of clouds in the Amazonian region. On the other hand, with cloud entrainments, observed values have been overestimated. Modeled sulfate rainwater concentration, using exponential decay and cloud heights of 16 km and no entrainments, presented the best results, reaching 97% of the observed value. The results, using shape parameter 5, are the best, improving the overall result. (C) 2008 Elsevier Ltd. All rights reserved.

Phosphoproteomics Profiling Suggests a Role for Nuclear beta IPKC in Transcription Processes of Undifferentiated Murine Embryonic Stem Cells

Protein kinase C (PKC) plays a key role in embryonic stem cell (ESC) proliferation, self-renewal and differentiation However, the function of specific PKC Isoenzymes have yet to be determined Of the PKCs expressed in undifferentiated ESCs, beta IPKC was the only isoenzyme abundantly expressed in the nuclei To investigate the role of beta IPKC in these cells, we employed a phosphoproteomics strategy and used two classical (cPKC) peptide modulators and one beta IPKC-specific inhibitor peptide We identified 13 nuclear proteins that are direct or indirect beta IPKC substrates in undifferentiated ESCs These proteins are known to be involved in regulating transcription, splicing, and chromatin remodeling during proliferation and differentiation Inhibiting beta IPKC had no effect on DNA synthesis in undifferentiated ESCs However, upon differentiation many cells seized to express beta IPKC and beta IPKC was frequently found in the cytoplasm Taken together, our results suggest that beta IPKC takes part in the processes that maintain ESCs in their undifferentiated state