A climatology of the diurnal variations of stratospheric and mesospheric ozone over Bern, Switzerland

The ground-based radiometer GROMOS, stationed in Bern (47.95° N, 7.44° E), Switzerland, has a unique dataset: it obtains ozone profiles from November 1994 to present with a time resolution of 30 min and equal quality during night- and daytime. Here, we derive a monthly climatology of the daily ozone cycle from 17 yr of GROMOS observation. We present the diurnal ozone variation of the stratosphere and mesosphere. Characterizing the diurnal cycle of stratospheric ozone is important for correct trend estimates of the ozone layer derived from satellite observations. The diurnal ozone cycle from GROMOS is compared to two models: The Whole Atmosphere Community Climate Model (WACCM) and the Hamburg Model of Neutral and Ionized Atmosphere (HAMMONIA). Aura Microwave Limb Sounder (Aura/MLS) ozone data, from night- and daytime overpasses over Bern, have also been included in the comparison. Generally, observation and models show good qualitative agreement: in the lower mesosphere, daytime ozone is for both GROMOS and models around 25% less than nighttime ozone (reference is 22:30–01:30). In the stratosphere, ozone reaches its maximum in the afternoon showing values several percent larger than the midnight value. It is important that diurnal ozone variations of this order are taken into account when merging different data sets for the derivation of long-term ozone trends in the stratosphere. Further, GROMOS and models indicate a seasonal behavior of daily ozone variations in the stratosphere with a larger afternoon maximum during daytime in summer than in winter. At 0.35 hPa, observations from GROMOS and Aura/MLS show a seasonal pattern in diurnal ozone variations with larger relative amplitudes during daytime in winter (−25 ± 5%) than in summer (−18 ± 4%) (compared to mean values around midnight). For the first time, a time series of the diurnal variations in ozone is presented: 17 yr of GROMOS data show strong interannual variations in the diurnal ozone cycle for both the stratosphere and the mesosphere. There are some indications that strong temperature tides can suppress the diurnal variation of stratospheric ozone via the anticorrelation of temperature and ozone. That means the spatio-temporal variability of solar thermal tides seems to affect the diurnal cycle of stratospheric ozone.

Experimental validation of a mean field model of mineralized collagen fiber arrays at two levels of hierarchy

In the course of this study, stiffness of a fibril array of mineralized collagen fibrils modeled with a mean field method was validated experimentally at site-matched two levels of tissue hierarchy using mineralized turkey leg tendons (MTLT). The applied modeling approaches allowed to model the properties of this unidirectional tissue from nanoscale (mineralized collagen fibrils) to macroscale (mineralized tendon). At the microlevel, the indentation moduli obtained with a mean field homogenization scheme were compared to the experimental ones obtained with microindentation. At the macrolevel, the macroscopic stiffness predicted with micro finite element (μFE) models was compared to the experimental stiffness measured with uniaxial tensile tests. Elastic properties of the elements in μFE models were injected from the mean field model or two-directional microindentations. Quantitatively, the indentation moduli can be properly predicted with the mean-field models. Local stiffness trends within specific tissue morphologies are very weak, suggesting additional factors responsible for the stiffness variations. At macrolevel, the μFE models underestimate the macroscopic stiffness, as compared to tensile tests, but the correlations are strong.

Modelling Lake Kivu water level variations over the last seven decades

This study aimed at analysing the hydrological changes in the Lake Kivu Basin over the last seven decades with focus on the response of the lake water level to meteorological factors and hydropower dam construction. Historical precipitation and lake water levels were acquired from literature, local agencies and from global databases in order to compile a coherent dataset. The net lake inflow was modelled using a soil water balance model and the water levels were reconstructed using a parsimonious lake water balance model. The soil water balance shows that 370 mm yr−1 (25%) of the precipitation in the catchment contributes to the runoff and baseflow whereas 1100 mm yr−1 (75%) contributes to the evapotranspiration. A review of the lake water balance resulted in the following estimates of hydrological contributions: 55%, 25%, and 20% of the overall inputs from precipitation, surface inflows, and subaquatic groundwater discharge, respectively. The overall losses were 58% and 42% for lake surface evaporation and outflow discharge, respectively. The hydrological model used indicated a remarkable sensitivity of the lake water levels to hydrometeorological variability up to 1977, when the outflow bed was artificially widened.

Daily ozone cycle in the stratosphere: global, regional and seasonal behaviour modelled with the Whole Atmosphere Community Climate Model

The Whole Atmosphere Community Climate Model (WACCM) is utilised to study the daily ozone cycle and underlying photochemical and dynamical processes. The analysis is focused on the daily ozone cycle in the middle stratosphere at 5 hPa where satellite-based trend estimates of stratospheric ozone are most biased by diurnal sampling effects and drifting satellite orbits. The simulated ozone cycle shows a minimum after sunrise and a maximum in the late afternoon. Further, a seasonal variation of the daily ozone cycle in the stratosphere was found. Depending on season and latitude, the peak-to-valley difference of the daily ozone cycle varies mostly between 3 and 5% (0.4 ppmv) with respect to the midnight ozone volume mixing ratio. The maximal variation of 15% (0.8 ppmv) is found at the polar circle in summer. The global pattern of the strength of the daily ozone cycle is mainly governed by the solar zenith angle and the sunshine duration. In addition, we find synoptic-scale variations in the strength of the daily ozone cycle. These variations are often anti-correlated to regional temperature anomalies and are due to the temperature dependence of the rate coefficients k2 and k3 of the Chapman cycle reactions. Further, the NOx catalytic cycle counteracts the accumulation of ozone during daytime and leads to an anti-correlation between anomalies in NOx and the strength of the daily ozone cycle. Similarly, ozone recombines with atomic oxygen which leads to an anti-correlation between anomalies in ozone abundance and the strength of the daily ozone cycle. At higher latitudes, an increase of the westerly (easterly) wind cause a decrease (increase) in the sunshine duration of an air parcel leading to a weaker (stronger) daily ozone cycle.

A climatology of the diurnal variations of stratospheric and mesospheric ozone over Bern, Switzerland

The ground-based radiometer GROMOS, stationed in Bern (47.95° N, 7.44° E), Switzerland, has a unique dataset: it obtains ozone profiles from November 1994 to present with a time resolution of 30 min and equal quality during night- and daytime. Here, we derive a monthly climatology of the daily ozone cycle from 17 yr of GROMOS observation. We present the diurnal ozone variation of the stratosphere and mesosphere. Characterizing the diurnal cycle of stratospheric ozone is important for correct trend estimates of the ozone layer derived from satellite observations. The diurnal ozone cycle from GROMOS is compared to two models: The Whole Atmosphere Community Climate Model (WACCM) and the Hamburg Model of Neutral and Ionized Atmosphere (HAMMONIA). Aura Microwave Limb Sounder (Aura/MLS) ozone data, from night- and daytime overpasses over Bern, have also been included in the comparison. Generally, observation and models show good qualitative agreement: in the lower mesosphere, daytime ozone is for both GROMOS and models around 25% less than nighttime ozone (reference is 22:30–01:30). In the stratosphere, ozone reaches its maximum in the afternoon showing values several percent larger than the midnight value. It is important that diurnal ozone variations of this order are taken into account when merging different data sets for the derivation of long-term ozone trends in the stratosphere. Further, GROMOS and models indicate a seasonal behavior of daily ozone variations in the stratosphere with a larger afternoon maximum during daytime in summer than in winter. At 0.35 hPa, observations from GROMOS and Aura/MLS show a seasonal pattern in diurnal ozone variations with larger relative amplitudes during daytime in winter (−25 ± 5%) than in summer (−18 ± 4%) (compared to mean values around midnight). For the first time, a time series of the diurnal variations in ozone is presented: 17 yr of GROMOS data show strong interannual variations in the diurnal ozone cycle for both the stratosphere and the mesosphere. There are some indications that strong temperature tides can suppress the diurnal variation of stratospheric ozone via the anticorrelation of temperature and ozone. That means the spatio-temporal variability of solar thermal tides seems to affect the diurnal cycle of stratospheric ozone.

The coupled atmosphere–chemistry–ocean model SOCOL-MPIOM

The newly developed atmosphere–ocean-chemistry-climate model SOCOL-MPIOM is presented by demonstrating the influence of the interactive chemistry module on the climate state and the variability. Therefore, we compare pre-industrial control simulations with (CHEM) and without (NOCHEM) interactive chemistry. In general, the influence of the chemistry on the mean state and the variability is small and mainly restricted to the stratosphere and mesosphere. The largest differences are found for the atmospheric dynamics in the polar regions, with slightly stronger northern and southern winter polar vortices in CHEM. The strengthening of the vortex is related to larger stratospheric temperature gradients, which are attributed to a parametrization of the absorption of ozone and oxygen in the Lyman-alpha, Schumann–Runge, Hartley, and Higgins bands. This effect is parametrized in the version with interactive chemistry only. A second reason for the temperature differences between CHEM and NOCHEM is related to diurnal variations in the ozone concentrations in the higher atmosphere, which are missing in NOCHEM. Furthermore, stratospheric water vapour concentrations differ substantially between the two experiments, but their effect on the temperatures is small. In both setups, the simulated intensity and variability of the northern polar vortex is inside the range of present day observations. Sudden stratospheric warming events are well reproduced in terms of their frequency, but the distribution amongst the winter months is too uniform. Additionally, the performance of SOCOL-MPIOM under changing external forcings is assessed for the period 1600–2000 using an ensemble of simulations driven by a spectral solar forcing reconstruction. The amplitude of the reconstruction is large in comparison to other state-of-the-art reconstructions, providing an upper limit for the importance of the solar signal. In the pre-industrial period (1600–1850) the simulated surface temperature trends are in reasonable agreement with temperature reconstructions, although the multi-decadal variability is more pronounced. This enhanced variability can be attributed to the variability in the solar forcing. The simulated temperature reductions during the Maunder Minimum are in the lowest probability range of the proxy records. During the Dalton Minimum, when also volcanic forcing is an important driver of temperature variations, the agreement is better. In the industrial period from 1850 onward SOCOL-MPIOM overestimates the temperature increase in comparison to observational data sets. Sensitivity simulations show that this overestimation can be attributed to the increasing trend in the solar forcing reconstruction that is used in this study and an additional warming induced by the simulated ozone changes.

Impact of epinephrine and norepinephrine on two dynamic indices in a porcine hemorrhagic shock model

Abstract BACKGROUND: Pulse pressure variations (PPVs) and stroke volume variations (SVVs) are dynamic indices for predicting fluid responsiveness in intensive care unit patients. These hemodynamic markers underscore Frank-Starling law by which volume expansion increases cardiac output (CO). The aim of the present study was to evaluate the impact of the administration of catecholamines on PPV, SVV, and inferior vena cava flow (IVCF). METHODS: In this prospective, physiologic, animal study, hemodynamic parameters were measured in deeply sedated and mechanically ventilated pigs. Systemic hemodynamic and pressure-volume loops obtained by inferior vena cava occlusion were recorded. Measurements were collected during two conditions, that is, normovolemia and hypovolemia, generated by blood removal to obtain a mean arterial pressure value lower than 60 mm Hg. At each condition, CO, IVCF, SVV, and PPV were assessed by catheters and flow meters. Data were compared between the conditions normovolemia and hypovolemia before and after intravenous administrations of norepinephrine and epinephrine using a nonparametric Wilcoxon test. RESULTS: Eight pigs were anesthetized, mechanically ventilated, and equipped. Both norepinephrine and epinephrine significantly increased IVCF and decreased PPV and SVV, regardless of volemic conditions (p < 0.05). However, epinephrine was also able to significantly increase CO regardless of volemic conditions. CONCLUSION: The present study demonstrates that intravenous administrations of norepinephrine and epinephrine increase IVCF, whatever the volemic conditions are. The concomitant decreases in PPV and SVV corroborate the fact that catecholamine administration recruits unstressed blood volume. In this regard, understanding a decrease in PPV and SVV values, after catecholamine administration, as an obvious indication of a restored volemia could be an outright misinterpretation.

A multi-model comparison of Atlantic multidecadal variability

A multi-model analysis of Atlantic multidecadal variability is performed with the following aims: to investigate the similarities to observations; to assess the strength and relative importance of the different elements of the mechanism proposed by Delworth et al. (J Clim 6:1993–2011, 1993) (hereafter D93) among coupled general circulation models (CGCMs); and to relate model differences to mean systematic error. The analysis is performed with long control simulations from ten CGCMs, with lengths ranging between 500 and 3600 years. In most models the variations of sea surface temperature (SST) averaged over North Atlantic show considerable power on multidecadal time scales, but with different periodicity. The SST variations are largest in the mid-latitude region, consistent with the short instrumental record. Despite large differences in model configurations, we find quite some consistency among the models in terms of processes. In eight of the ten models the mid-latitude SST variations are significantly correlated with fluctuations in the Atlantic meridional overturning circulation (AMOC), suggesting a link to northward heat transport changes. Consistent with this link, the three models with the weakest AMOC have the largest cold SST bias in the North Atlantic. There is no linear relationship on decadal timescales between AMOC and North Atlantic Oscillation in the models. Analysis of the key elements of the D93 mechanisms revealed the following: Most models present strong evidence that high-latitude winter mixing precede AMOC changes. However, the regions of wintertime convection differ among models. In most models salinity-induced density anomalies in the convective region tend to lead AMOC, while temperature-induced density anomalies lead AMOC only in one model. However, analysis shows that salinity may play an overly important role in most models, because of cold temperature biases in their relevant convective regions. In most models subpolar gyre variations tend to lead AMOC changes, and this relation is strong in more than half of the models.

Reconstructing the subsurface ocean decadal variability using surface nudging in a perfect model framework

Initialising the ocean internal variability for decadal predictability studies is a new area of research and a variety of ad hoc methods are currently proposed. In this study, we explore how nudging with sea surface temperature (SST) and salinity (SSS) can reconstruct the threedimensional variability of the ocean in a perfect model framework. This approach builds on the hypothesis that oceanic processes themselves will transport the surface information into the ocean interior as seen in ocean-only simulations. Five nudged simulations are designed to reconstruct a 150 years ‘‘target’’ simulation, defined as a portion of a long control simulation. The nudged simulations differ by the variables restored to, SST or SST + SSS, and by the area where the nudging is applied. The strength of the heat flux feedback is diagnosed from observations and the restoring coefficients for SSS use the same time-scale. We observed that this choice prevents spurious convection at high latitudes and near sea-ice border when nudging both SST and SSS. In the tropics, nudging the SST is enough to reconstruct the tropical atmosphere circulation and the associated dynamical and thermodynamical impacts on the underlying ocean. In the tropical Pacific Ocean, the profiles for temperature show a significant correlation from the surface down to 2,000 m, due to dynamical adjustment of the isopycnals. At mid-tohigh latitudes, SSS nudging is required to reconstruct both the temperature and the salinity below the seasonal thermocline. This is particularly true in the North Atlantic where adding SSS nudging enables to reconstruct the deep convection regions of the target. By initiating a previously documented 20-year cycle of the model, the SST + SSS nudging is also able to reproduce most of the AMOC variations, a key source of decadal predictability. Reconstruction at depth does not significantly improve with amount of time spent nudging and the efficiency of the surface nudging rather depends on the period/events considered. The joint SST + SSS nudging applied verywhere is the most efficient approach. It ensures that the right water masses are formed at the right surface density, the subsequent circulation, subduction and deep convection further transporting them at depth. The results of this study underline the potential key role of SSS for decadal predictability and further make the case for sustained largescale observations of this field.

Modelling Nd-isotopes with a coarse resolution ocean circulation model: Sensitivities to model parameters and source/sink distributions

The neodymium (Nd) isotopic composition (Nd) of seawater is a quasi-conservative tracer of water mass mixing and is assumed to hold great potential for paleoceanographic studies. Here we present a comprehensive approach for the simulation of the two neodymium isotopes 143Nd, and 144Nd using the Bern3D model, a low resolution ocean model. The high computational efficiency of the Bern3D model in conjunction with our comprehensive approach allows us to systematically and extensively explore the sensitivity of Nd concentrations and Nd to the parametrisation of sources and sinks. Previous studies have been restricted in doing so either by the chosen approach or by computational costs. Our study thus presents the most comprehensive survey of the marine Nd cycle to date. Our model simulates both Nd concentrations as well as Nd in good agreement with observations. Nd covaries with salinity, thus underlining its potential as a water mass proxy. Results confirm that the continental margins are required as a Nd source to simulate Nd concentrations and Nd consistent with observations. We estimate this source to be slightly smaller than reported in previous studies and find that above a certain magnitude its magnitude affects Nd only to a small extent. On the other hand, the parametrisation of the reversible scavenging considerably affects the ability of the model to simulate both, Nd concentrations and Nd. Furthermore, despite their small contribution, we find dust and rivers to be important components of the Nd cycle. In additional experiments, we systematically varied the diapycnal diffusivity as well as the Atlantic-to-Pacific freshwater flux to explore the sensitivity of Nd concentrations and its isotopic signature to the strength and geometry of the overturning circulation. These experiments reveal that Nd concentrations and Nd are comparatively little affected by variations in diapycnal diffusivity and the Atlantic-to-Pacific freshwater flux. In contrast, an adequate representation of Nd sources and sinks is crucial to simulate Nd concentrations and Nd consistent with observations. The good agreement of our results with observations paves the way for the evaluation of the paleoceanographic potential of Nd in further model studies.

Effect of surface restoring on subsurface variability in a climate model during 1949–2005

Initializing the ocean for decadal predictability studies is a challenge, as it requires reconstructing the little observed subsurface trajectory of ocean variability. In this study we explore to what extent surface nudging using well-observed sea surface temperature (SST) can reconstruct the deeper ocean variations for the 1949–2005 period. An ensemble made with a nudged version of the IPSLCM5A model and compared to ocean reanalyses and reconstructed datasets. The SST is restored to observations using a physically-based relaxation coefficient, in contrast to earlier studies, which use a much larger value. The assessment is restricted to the regions where the ocean reanalyses agree, i.e. in the upper 500 m of the ocean, although this can be latitude and basin dependent. Significant reconstruction of the subsurface is achieved in specific regions, namely region of subduction in the subtropical Atlantic, below the thermocline in the equatorial Pacific and, in some cases, in the North Atlantic deep convection regions. Beyond the mean correlations, ocean integrals are used to explore the time evolution of the correlation over 20-year windows. Classical fixed depth heat content diagnostics do not exhibit any significant reconstruction between the different existing bservation-based references and can therefore not be used to assess global average time-varying correlations in the nudged simulations. Using the physically based average temperature above an isotherm (14°C) alleviates this issue in the tropics and subtropics and shows significant reconstruction of these quantities in the nudged simulations for several decades. This skill is attributed to the wind stress reconstruction in the tropics, as already demonstrated in a perfect model study using the same model. Thus, we also show here the robustness of this result in an historical and observational context.

On the seasonal variations of salinity of the tropical Atlantic mixed layer

The physical processes controlling the mixed layer salinity (MLS) seasonal budget in the tropical Atlantic Ocean are investigated using a regional configuration of an ocean general circulation model. The analysis reveals that the MLS cycle is generally weak in comparison of individual physical processes entering in the budget because of strong compensation. In evaporative regions, around the surface salinity maxima, the ocean acts to freshen the mixed layer against the action of evaporation. Poleward of the southern SSS maxima, the freshening is ensured by geostrophic advection, the vertical salinity diffusion and, during winter, a dominant contribution of the convective entrainment. On the equatorward flanks of the SSS maxima, Ekman transport mainly contributes to supply freshwater from ITCZ regions while vertical salinity diffusion adds on the effect of evaporation. All these terms are phase locked through the effect of the wind. Under the seasonal march of the ITCZ and in coastal areas affected by river (7°S:15°N), the upper ocean freshening by precipitations and/or runoff is attenuated by vertical salinity diffusion. In the eastern equatorial regions, seasonal cycle of wind forced surface currents advect freshwaters, which are mixed with subsurface saline water because of the strong vertical turbulent diffusion. In all these regions, the vertical diffusion presents an important contribution to the MLS budget by providing, in general, an upwelling flux of salinity. It is generally due to vertical salinity gradient and mixing due to winds. Furthermore, in the equator where the vertical shear, associated to surface horizontal currents, is developed, the diffusion depends also on the sheared flow stability.

Cosmogenic 180W variations in meteorites and re-assessment of a possible 184Os–180W decay system

We measured tungsten (W) isotopes in 23 iron meteorites and the metal phase of the CB chondrite Gujba in order to ascertain if there is evidence for a large-scale nucleosynthetic heterogeneity in the p-process isotope 180W in the solar nebula as recently suggested by Schulz et al. (2013). We observed large excesses in 180W (up to ≈ 6 ε) in some irons. However, significant within-group variations in magmatic IIAB and IVB irons are not consistent with a nucleosynthetic origin, and the collateral effects on 180W from an s-deficit in IVB irons cannot explain the total variation. We present a new model for the combined effects of spallation and neutron capture reactions on 180W in iron meteorites and show that at least some of the observed within-group variability is explained by cosmic ray effects. Neutron capture causes burnout of 180W, whereas spallation reactions lead to positive shifts in 180W. These effects depend on the target composition and cosmic-ray exposure duration; spallation effects increase with Re/W and Os/W ratios in the target and with exposure age. The correlation of 180W/184W with Os/W ratios in iron meteorites results in part from spallogenic production of 180W rather than from 184Os decay, contrary to a recent study by Peters et al. (2014). Residual ε180W excesses after correction for an s-deficit and for cosmic ray effects may be due to ingrowth of 180W from 184Os decay, but the magnitude of this ingrowth is at least a factor of ≈2 smaller than previously suggested. These much smaller effects strongly limit the applicability of the putative 184Os-180W system to investigate geological problems.

Atmospheric water parameters measured by a ground-based microwave radiometer and compared with the WRF model

The microwave radiometer TROWARA measures integrated water vapour (IWV) and integrated cloud liquid water (ILW) at Bern since 1994 with a time resolution of 7 s. In this study, we compare TROWARA measurements with a simulation of summer 2012 in Switzerland performed with the Weather Research and Forecasting (WRF) model. It is found that the WRF model agrees very well with TROWARA’s IWV variations with a mean bias of only 0.7 mm. The ILW distribution of the WRF model, although similar in shape to TROWARA’s distribution, overestimates the fraction of clear sky periods (83% compared to 60%).

Correlated cosmogenic W and Os isotopic variations in Carbo and implications for Hf-W chronology

An obstacle for establishing the chronology of iron meteorite formation using 182Hf-182W systematics (t1/2 = 8.9 Myr) is to find proper neutron fluence monitors to correct for cosmic ray modification of W isotopic composition. Recent studies showed that siderophile elements such as Pt and Os could serve such a purpose. To test and calibrate these neutron dosimeters, the isotopic compositions of W and Os were measured in a slab of the IID iron meteorite Carbo. This slab has a well-characterized noble gas depth profile reflecting different degrees of shielding to cosmic rays. The results show that W and Os isotopic ratios correlate with distance from the pre-atmospheric center. Negative correlations, barely resolved within error, were found between epsilo190Os-epsilo189Os and epsilo186Os-epsilo189Os with slopes of -0.64 ± 0.45 and -1.8(+1.9/-2.1), respectively. These Os isotope correlations broadly agree with model predictions for capture of secondary neutrons produced by cosmic ray irradiation and results reported previously for other groups of iron meteorites. Correlations were also found between epsilo182W-epsilo189Os (slope = 1.02 ± 0.37) and epsilo182W-epsilo190Os (slope = -1.38 ± 0.58). Intercepts of these two correlations yield pre-exposure epsilo182W values of -3.32 ± 0.51 and -3.62 ± 0.23, respectively (weighted average epsilo182W = -3.57 ± 0.21). This value relies on a large extrapolation leading to a large uncertainty but gives a metal-silicate segregation age of -0.5 ± 2.4 Myr after formation of the solar system. Combining the iron meteorite measurements with simulations of cosmogenic effects in iron meteorites, equations are presented to calculate and correct for cosmogenic effects on 182W using Os isotopes.