A Nonstationary Space-Time Gaussian Process Model for Partially Converged Simulations

In the context of expensive numerical experiments, a promising solution for alleviating the computational costs consists of using partially converged simulations instead of exact solutions. The gain in computational time is at the price of precision in the response. This work addresses the issue of fitting a Gaussian process model to partially converged simulation data for further use in prediction. The main challenge consists of the adequate approximation of the error due to partial convergence, which is correlated in both design variables and time directions. Here, we propose fitting a Gaussian process in the joint space of design parameters and computational time. The model is constructed by building a nonstationary covariance kernel that reflects accurately the actual structure of the error. Practical solutions are proposed for solving parameter estimation issues associated with the proposed model. The method is applied to a computational fluid dynamics test case and shows significant improvement in prediction compared to a classical kriging model.

The effect of a myocardial infarction on the normalized time-varying elastance curve

It has been suggested that the shape of the normalized time-varying elastance curve [E(n)(t(n))] is conserved in different cardiac pathologies. We hypothesize, however, that the E(n)(t(n)) differs quantitatively after myocardial infarction (MI). Sprague-Dawley rats (n = 9) were anesthetized, and the left anterior descending coronary artery was ligated to provoke the MI. A sham-operated control group (CTRL) (n = 10) was treated without the MI. Two months later, a conductance catheter was inserted into the left ventricle (LV). The LV pressure and volume were measured and the E(n)(t(n)) derived. Slopes of E(n)(t(n)) during the preejection period (alpha(PEP)), ejection period (alpha(EP)), and their ratio (beta = alpha(EP)/alpha(PEP)) were calculated, together with the characteristic decay time during isovolumic relaxation (tau) and the normalized elastance at end diastole (E(min)(n)). MI provoked significant LV chamber dilatation, thus a loss in cardiac output (-33%), ejection fraction (-40%), and stroke volume (-30%) (P < 0.05). Also, it caused significant calcium increase (17-fold), fibrosis (2-fold), and LV hypertrophy. End-systolic elastance dropped from 0.66 +/- 0.31 mmHg/microl (CTRL) to 0.34 +/- 0.11 mmHg/microl (MI) (P < 0.05). Normalized elastance was significantly reduced in the MI group during the preejection, ejection, and diastolic periods (P < 0.05). The slope of E(n)(t(n)) during the alpha(PEP) and beta were significantly altered after MI (P < 0.05). Furthermore, tau and end-diastolic E(min)(n) were both significantly augmented in the MI group. We conclude that the E(n)(t(n)) differs quantitatively in all phases of the heart cycle, between normal and hearts post-MI. This should be considered when utilizing the single-beat concept.

Unconventional Monetary Policy and the Great Recession: Estimating the Macroeconomic Effects of a Spread Compression at the Zero Lower Bound

We explore the macroeconomic effects of a compression in the long-term bond yield spread within the context of the Great Recession of 2007–09 via a time-varying parameter structural VAR model. We identify a “pure” spread shock defined as a shock that leaves the policy rate unchanged, which allows us to characterize the macroeconomic consequences of a decline in the yield spread induced by central banks’ asset purchases within an environment in which the policy rate is constrained by the effective zero lower bound. Two key findings stand out. First, compressions in the long-term yield spread exert a powerful effect on both output growth and inflation. Second, conditional on available estimates of the impact of the Federal Reserve’s and the Bank of England’s asset purchase programs on long-term yield spreads, our counterfactual simulations suggest that U.S. and U.K. unconventional monetary policy actions have averted significant risks both of deflation and of output collapses comparable to those that took place during the Great Depression.

Time-varying signal analysis to detect high-altitude periodic breathing in climbers ascending to extreme altitude.

This work investigates the performance of cardiorespiratory analysis detecting periodic breathing (PB) in chest wall recordings in mountaineers climbing to extreme altitude. The breathing patterns of 34 mountaineers were monitored unobtrusively by inductance plethysmography, ECG and pulse oximetry using a portable recorder during climbs at altitudes between 4497 and 7546 m on Mt. Muztagh Ata. The minute ventilation (VE) and heart rate (HR) signals were studied, to identify visually scored PB, applying time-varying spectral, coherence and entropy analysis. In 411 climbing periods, 30-120 min in duration, high values of mean power (MP(VE)) and slope (MSlope(VE)) of the modulation frequency band of VE, accurately identified PB, with an area under the ROC curve of 88 and 89%, respectively. Prolonged stay at altitude was associated with an increase in PB. During PB episodes, higher peak power of ventilatory (MP(VE)) and cardiac (MP(LF)(HR) ) oscillations and cardiorespiratory coherence (MP(LF)(Coher)), but reduced ventilation entropy (SampEn(VE)), was observed. Therefore, the characterization of cardiorespiratory dynamics by the analysis of VE and HR signals accurately identifies PB and effects of altitude acclimatization, providing promising tools for investigating physiologic effects of environmental exposures and diseases.

Racing against Time: Varying Patterns of Social Mobilization in Serbia and Georgia

After the collapse of the Soviet Union and Yugoslavia, a number of actors started to engage in the power struggle for the opportunities to shape the new order in successive nation-states. In Serbia and Georgia historically hegemonic Orthodox Christian churches were among the firsts in the frontlines for political and economic power. More than a decade has passed since the so-called Coloured Revolutions in Georgia and Serbia, and the Orthodox churches still remain participants of an ongoing socio-political transition of these states. The revival of public role of religion appeared temporary in Serbia followed by a gradual decline of an influence of the Orthodox Church over political life and legal process. However, in Georgia the public and political role of religion increased rather than declined albeit changed shape. Examining the degree to which the two Orthodox churches can influence the political agenda in Serbia and Georgia, the paper attempts to understand how church-State relations work in practice. By bringing rich empirical data from the field (70 interviews with (arch)bishops, priests and religious clerics in Georgia and Serbia added to field observations), the paper reflects on the themes under which the two Orthodox churches mobilize public protest in Serbia and Georgia. The paper further looks at varying State responses and their broader implication for church-state problematique.

A numerical model for inert gas transport in the lung based on fractal airway morphology

Patients suffering from cystic fibrosis (CF) show thick secretions, mucus plugging and bronchiectasis in bronchial and alveolar ducts. This results in substantial structural changes of the airway morphology and heterogeneous ventilation. Disease progression and treatment effects are monitored by so-called gas washout tests, where the change in concentration of an inert gas is measured over a single or multiple breaths. The result of the tests based on the profile of the measured concentration is a marker for the severity of the ventilation inhomogeneity strongly affected by the airway morphology. However, it is hard to localize underlying obstructions to specific parts of the airways, especially if occurring in the lung periphery. In order to support the analysis of lung function tests (e.g. multi-breath washout), we developed a numerical model of the entire airway tree, coupling a lumped parameter model for the lung ventilation with a 4th-order accurate finite difference model of a 1D advection-diffusion equation for the transport of an inert gas. The boundary conditions for the flow problem comprise the pressure and flow profile at the mouth, which is typically known from clinical washout tests. The natural asymmetry of the lung morphology is approximated by a generic, fractal, asymmetric branching scheme which we applied for the conducting airways. A conducting airway ends when its dimension falls below a predefined limit. A model acinus is then connected to each terminal airway. The morphology of an acinus unit comprises a network of expandable cells. A regional, linear constitutive law describes the pressure-volume relation between the pleural gap and the acinus. The cyclic expansion (breathing) of each acinus unit depends on the resistance of the feeding airway and on the flow resistance and stiffness of the cells themselves. Special care was taken in the development of a conservative numerical scheme for the gas transport across bifurcations, handling spatially and temporally varying advective and diffusive fluxes over a wide range of scales. Implicit time integration was applied to account for the numerical stiffness resulting from the discretized transport equation. Local or regional modification of the airway dimension, resistance or tissue stiffness are introduced to mimic pathological airway restrictions typical for CF. This leads to a more heterogeneous ventilation of the model lung. As a result the concentration in some distal parts of the lung model remains increased for a longer duration. The inert gas concentration at the mouth towards the end of the expirations is composed of gas from regions with very different washout efficiency. This results in a steeper slope of the corresponding part of the washout profile.

Accounting for baseline differences and measurement error in the analysis of change over time

If change over time is compared in several groups, it is important to take into account baseline values so that the comparison is carried out under the same preconditions. As the observed baseline measurements are distorted by measurement error, it may not be sufficient to include them as covariate. By fitting a longitudinal mixed-effects model to all data including the baseline observations and subsequently calculating the expected change conditional on the underlying baseline value, a solution to this problem has been provided recently so that groups with the same baseline characteristics can be compared. In this article, we present an extended approach where a broader set of models can be used. Specifically, it is possible to include any desired set of interactions between the time variable and the other covariates, and also, time-dependent covariates can be included. Additionally, we extend the method to adjust for baseline measurement error of other time-varying covariates. We apply the methodology to data from the Swiss HIV Cohort Study to address the question if a joint infection with HIV-1 and hepatitis C virus leads to a slower increase of CD4 lymphocyte counts over time after the start of antiretroviral therapy.

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.