Quantitative 1H-magnetic resonance spectroscopy of human brain: Influence of composition and parameterization of the basis set in linear combination model-fitting

Localized short-echo-time (1)H-MR spectra of human brain contain contributions of many low-molecular-weight metabolites and baseline contributions of macromolecules. Two approaches to model such spectra are compared and the data acquisition sequence, optimized for reproducibility, is presented. Modeling relies on prior knowledge constraints and linear combination of metabolite spectra. Investigated was what can be gained by basis parameterization, i.e., description of basis spectra as sums of parametric lineshapes. Effects of basis composition and addition of experimentally measured macromolecular baselines were investigated also. Both fitting methods yielded quantitatively similar values, model deviations, error estimates, and reproducibility in the evaluation of 64 spectra of human gray and white matter from 40 subjects. Major advantages of parameterized basis functions are the possibilities to evaluate fitting parameters separately, to treat subgroup spectra as independent moieties, and to incorporate deviations from straightforward metabolite models. It was found that most of the 22 basis metabolites used may provide meaningful data when comparing patient cohorts. In individual spectra, sums of closely related metabolites are often more meaningful. Inclusion of a macromolecular basis component leads to relatively small, but significantly different tissue content for most metabolites. It provides a means to quantitate baseline contributions that may contain crucial clinical information.

Reactive mass transport modelling of a three-dimensional vertical fault zone with a finger-like convective flow regime

For a three-dimensional vertically-oriented fault zone, we consider the coupled effects of fluid flow, heat transfer and reactive mass transport, to investigate the patterns of fluid flow, temperature distribution, mineral alteration and chemically induced porosity changes. We show, analytically and numerically, that finger-like convection patterns can arise in a vertically-oriented fault zone. The onset and patterns of convective fluid flow are controlled by the Rayleigh number which is a function of the thermal properties of the fluid and the rock, the vertical temperature gradient, and the height and the permeability of the fault zone. Vigorous fluid flow causes low temperature gradients over a large region of the fault zone. In such a case, flow across lithological interfaces becomes the most important mechanism for the formation of sharp chemical reaction fronts. The degree of rock buffering, the extent and intensity of alteration, the alteration mineralogy and in some cases the formation of ore deposits are controlled by the magnitude of the flow velocity across these compositional interfaces in the rock. This indicates that alteration patterns along compositional boundaries in the rock may provide some insights into the convection pattern. The advective mass and heat exchanges between the fault zone and the wallrock depend on the permeability contrast between the fault zone and the wallrock. A high permeability contrast promotes focussed convective flow within the fault zone and diffusive exchange of heat and chemical reactants between the fault zone and the wallrock. However, a more gradual permeability change may lead to a regional-scale convective flow system where the flow pattern in the fault affects large-scale fluid flow, mass transport and chemical alteration in the wallrocks

Parameterization of temperature sensitivity of spring phenology and its application in explaining diverse phenological responses to temperature change

Existing evidence of plant phenological change to temperature increase demonstrates that the phenological responsiveness is greater at warmer locations and in early-season plant species. Explanations of these findings are scarce and not settled. Some studies suggest considering phenology as one functional trait within a plant's life history strategy. In this study, we adapt an existing phenological model to derive a generalized sensitivity in space (SpaceSens) model for calculating temperature sensitivity of spring plant phenophases across species and locations. The SpaceSens model have three parameters, including the temperature at the onset date of phenophases (Tp), base temperature threshold (Tb) and the length of period (L) used to calculate the mean temperature when performing regression analysis between phenology and temperature. A case study on first leaf date of 20 plant species from eastern China shows that the change of Tp and Tb among different species accounts for interspecific difference in temperature sensitivity. Moreover, lower Tp at lower latitude is the main reason why spring phenological responsiveness is greater there. These results suggest that spring phenophases of more responsive, early-season plants (especially in low latitude) will probably continue to diverge from the other late-season plants with temperatures warming in the future.

A new, fast and flexible radiative transfer method for Venus general circulation models

We present a new radiation scheme for the Oxford Planetary Unified Model System for Venus, suitable for the solar and thermal bands. This new and fast radiative parameterization uses a different approach in the two main radiative wavelength bands: solar radiation (0.1-5.5 mu m) and thermal radiation (1.7-260 mu m). The solar radiation calculation is based on the delta-Eddington approximation (two-stream-type) with an adding layer method. For the thermal radiation case, a code based on an absorptivity/emissivity formulation is used. The new radiative transfer formulation implemented is intended to be computationally light, to allow its incorporation in 3D global circulation models, but still allowing for the calculation of the effect of atmospheric conditions on radiative fluxes. This will allow us to investigate the dynamical-radiative-microphysical feedbacks. The model flexibility can be also used to explore the uncertainties in the Venus atmosphere such as the optical properties in the deep atmosphere or cloud amount. The results of radiative cooling and heating rates and the global-mean radiative-convective equilibrium temperature profiles for different atmospheric conditions are presented and discussed. This new scheme works in an atmospheric column and can be easily implemented in 3D Venus global circulation models. (C) 2014 Elsevier Ltd. All rights reserved.

PM(10) remote sensing from geostationary SEVIRI and polar-orbiting MODIS sensors over the complex terrain of the European Alpine region

The subject of this study is to investigate the capability of spaceborne remote sensing data to predict ground concentrations of PM10 over the European Alpine region using satellite derived Aerosol Optical Depth (AOD) from the geostationary Spinning Enhanced Visible and InfraRed Imager (SEVIRI) and the polar-orbiting MODerate resolution Imaging Spectroradiometer (MODIS). The spatial and temporal resolutions of these aerosol products (10 km and 2 measurements per day for MODIS, ∼ 25 km and observation intervals of 15 min for SEVIRI) permit an evaluation of PM estimation from space at different spatial and temporal scales. Different empirical linear relationships between coincident AOD and PM10 observations are evaluated at 13 ground-based PM measurement sites, with the assumption that aerosols are vertically homogeneously distributed below the planetary Boundary Layer Height (BLH). The BLH and Relative Humidity (RH) variability are assessed, as well as their impact on the parameterization. The BLH has a strong influence on the correlation of daily and hourly time series, whilst RH effects are less clear and smaller in magnitude. Despite its lower spatial resolution and AOD accuracy, SEVIRI shows higher correlations than MODIS (rSEV∼ 0.7, rMOD∼ 0.6) with regard to daily averaged PM10. Advantages from MODIS arise only at hourly time scales in mountainous locations but lower correlations were found for both sensors at this time scale (r∼ 0.45). Moreover, the fraction of days in 2008 with at least one satellite observation was 27% for SEVIRI and 17% for MODIS. These results suggest that the frequency of observations plays an important role in PM monitoring, while higher spatial resolution does not generally improve the PM estimation. Ground-based Sun Photometer (SP) measurements are used to validate the satellite-based AOD in the study region and to discuss the impact of aerosols' micro-physical properties in the empirical models. A lower error limit of 30 to 60% in the PM10 assessment from space is estimated in the study area as a result of AOD uncertainties, variability of aerosols properties and the heterogeneity of ground measurement sites. It is concluded that SEVIRI has a similar capacity to map PM as sensors on board polar-orbiting platforms, with the advantage of a higher number of observations. However, the accuracy represents a serious limitation to the applicability of satellites for ground PM mapping, especially in mountainous areas.

The simultaneous role of an alveolus as flow mixer and flow feeder for the deposition of inhaled submicron particles

In an effort to understand the fate of inhaled submicron particles in the small sacs, or alveoli, comprising the gas-exchange region of the lung, we calculated the flow in three-dimensional (3D) rhythmically expanding models of alveolated ducts. Since convection toward the alveolar walls is a precursor to particle deposition, it was the goal of this paper to investigate the streamline maps' dependence upon alveoli location along the acinar tree. On the alveolar midplane, the recirculating flow pattern exhibited closed streamlines with a stagnation saddle point. Off the midplane we found no closed streamlines but nested, funnel-like, spiral, structures (reminiscent of Russian nesting dolls) that were directed towards the expanding walls in inspiration, and away from the contracting walls in expiration. These nested, funnel-like, structures were surrounded by air that flowed into the cavity from the central channel over inspiration and flowed from the cavity to the central channel over expiration. We also found that fluid particle tracks exhibited similar nested funnel-like spiral structures. We conclude that these unique alveolar flow structures may be of importance in enhancing deposition. In addition, due to inertia, the nested, funnel-like, structures change shape and position slightly during a breathing cycle, resulting in flow mixing. Also, each inspiration feeds a fresh supply of particle-laden air from the central channel to the region surrounding the mixing region. Thus, this combination of flow mixer and flow feeder makes each individual alveolus an effective mixing unit, which is likely to play an important role in determining the overall efficiency of convective mixing in the acinus.

Correspondence establishment in statistical modeling of shapes with arbitrary topology

Correspondence establishment is a key step in statistical shape model building. There are several automated methods for solving this problem in 3D, but they usually can only handle objects with simple topology, like that of a sphere or a disc. We propose an extension to correspondence establishment over a population based on the optimization of the minimal description length function, allowing considering objects with arbitrary topology. Instead of using a fixed structure of kernel placement on a sphere for the systematic manipulation of point landmark positions, we rely on an adaptive, hierarchical organization of surface patches. This hierarchy can be built on surfaces of arbitrary topology and the resulting patches are used as a basis for a consistent, multi-scale modification of the surfaces' parameterization, based on point distribution models. The feasibility of the approach is demonstrated on synthetic models with different topologies.

Zielinstruktionen, räumliche Quiet-Eye-Verankerung und Bewegungsparametrisierung: Hinweise auf einen Wirkmechanismus

Die motorikwissenschaftliche Befundlage zum sogenannten „Quiet Eye“ weist darauf hin, dass hohe sportmotorische Leistungen, insbesondere in Präzisionsaufgaben, mit einer langen finalen Fixation vor der Bewegungsentfaltung einhergehen. Ein Mechanismus, der diesen Zusammenhang aus einer kognitionspsychologischen Perspektive erklären könnte, ist die Optimierung von Informationsverarbeitungsprozessen der Bewegungsparametrisierung. Diese Annahme wurde durch eine experimentelle Manipulation von Zielinstruktionen in einer Ballwurfaufgabe untersucht. Zum einen zeigen die Ergebnisse, dass sich die räumliche Verankerung des Quiet Eye in Abhängigkeit der variierten Aufgabenziele verändert; zum anderen deuten die Befunde darauf hin, dass sich Veränderungen der Verankerung im Bewegungsresultat niederschlagen. Damit wird ein kognitiver Wirkmechanismus plausibilisiert, nach dem die Bewegungsgenauigkeit durch Zielinstruktion via räumliche Quiet-Eye-Verankerung bestimmt wird.

SLR-derived terrestrial reference frame using the observations to LAGEOS-1/2, Starlette, Stella, and AJISAI

Currently, the contributions of Starlette, Stella, and AJISAI are not taken into account when defining the International Terrestrial Reference Frame (ITRF), despite the large amount of data collected in a long time-span. Consequently, the SLR-derived parameters and the SLR part of the ITRF are almost exclusively defined by LAGEOS-1 and LAGEOS-2. We investigate the potential of combining the observations to several SLR satellites with different orbital characteristics. Ten years of SLR data are homogeneously processed using the development version 5.3 of the Bernese GNSS Software. Special emphasis is put on orbit parameterization and the impact of LEO data on the estimation of the geocenter coordinates, Earth rotation parameters, Earth gravity field coefficients, and the station coordinates in one common adjustment procedure. We find that the parameters derived from the multi-satellite solutions are of better quality than those obtained in single satellite solutions or solutions based on the two LAGEOS satellites. A spectral analysis of the SLR network scale w.r.t. SLRF2008 shows that artifacts related to orbit perturbations in the LAGEOS-1/2 solutions, i.e., periods related to the draconitic years of the LAGEOS satellites, are greatly reduced in the combined solutions.

Glacial–interglacial dynamics of Antarctic firn columns: comparison between simulations and ice core air-δ15N measurements

Correct estimation of the firn lock-in depth is essential for correctly linking gas and ice chronologies in ice core studies. Here, two approaches to constrain the firn depth evolution in Antarctica are presented over the last deglaciation: outputs of a firn densification model, and measurements of δ15N of N2 in air trapped in ice core, assuming that δ15N is only affected by gravitational fractionation in the firn column. Since the firn densification process is largely governed by surface temperature and accumulation rate, we have investigated four ice cores drilled in coastal (Berkner Island, BI, and James Ross Island, JRI) and semi-coastal (TALDICE and EPICA Dronning Maud Land, EDML) Antarctic regions. Combined with available ice core air-δ15N measurements from the EPICA Dome C (EDC) site, the studied regions encompass a large range of surface accumulation rates and temperature conditions. Our δ15N profiles reveal a heterogeneous response of the firn structure to glacial–interglacial climatic changes. While firn densification simulations correctly predict TALDICE δ15N variations, they systematically fail to capture the large millennial-scale δ15N variations measured at BI and the δ15N glacial levels measured at JRI and EDML – a mismatch previously reported for central East Antarctic ice cores. New constraints of the EDML gas–ice depth offset during the Laschamp event (~41 ka) and the last deglaciation do not favour the hypothesis of a large convective zone within the firn as the explanation of the glacial firn model–δ15N data mismatch for this site. While we could not conduct an in-depth study of the influence of impurities in snow for firnification from the existing datasets, our detailed comparison between the δ15N profiles and firn model simulations under different temperature and accumulation rate scenarios suggests that the role of accumulation rate may have been underestimated in the current description of firnification models.

Two Stable Equilibria of the Atlantic Subpolar Gyre

The cyclonic circulation of the Atlantic subpolar gyre is a key mechanism for North Atlantic climate variability on a wide range of time scales. It is generally accepted that it is driven by both cyclonic winds and buoyancy forcing, yet the individual importance and dynamical interactions of the two contributions remain unclear. The authors propose a simplified four-box model representing the convective basin of the Labrador Sea and its shallow and deep boundary current system, the western subpolar gyre. Convective heat loss drives a baroclinic flow of relatively light water around the dense center. Eddy salt flux from the boundary current to the center increases with a stronger circulation, favors the formation of dense waters, and thereby sustains a strong baroclinic flow, approximately 10%–25% of the total. In contrast, when the baroclinic flow is not active, surface waters may be too fresh to convect, and a buoyancy-driven circulation cannot develop. This situation corresponds to a second stable circulation mode. A hysteresis is found for variations in surface freshwater flux and the salinity of the near-surface boundary current. An analytical solution is presented and analyzed.

Transconvolution and the virtual positron emission tomograph--a new method for cross calibration in quantitative PET∕CT imaging

PURPOSE Positron emission tomography (PET)∕computed tomography (CT) measurements on small lesions are impaired by the partial volume effect, which is intrinsically tied to the point spread function of the actual imaging system, including the reconstruction algorithms. The variability resulting from different point spread functions hinders the assessment of quantitative measurements in clinical routine and especially degrades comparability within multicenter trials. To improve quantitative comparability there is a need for methods to match different PET∕CT systems through elimination of this systemic variability. Consequently, a new method was developed and tested that transforms the image of an object as produced by one tomograph to another image of the same object as it would have been seen by a different tomograph. The proposed new method, termed Transconvolution, compensates for differing imaging properties of different tomographs and particularly aims at quantitative comparability of PET∕CT in the context of multicenter trials. METHODS To solve the problem of image normalization, the theory of Transconvolution was mathematically established together with new methods to handle point spread functions of different PET∕CT systems. Knowing the point spread functions of two different imaging systems allows determining a Transconvolution function to convert one image into the other. This function is calculated by convolving one point spread function with the inverse of the other point spread function which, when adhering to certain boundary conditions such as the use of linear acquisition and image reconstruction methods, is a numerically accessible operation. For reliable measurement of such point spread functions characterizing different PET∕CT systems, a dedicated solid-state phantom incorporating (68)Ge∕(68)Ga filled spheres was developed. To iteratively determine and represent such point spread functions, exponential density functions in combination with a Gaussian distribution were introduced. Furthermore, simulation of a virtual PET system provided a standard imaging system with clearly defined properties to which the real PET systems were to be matched. A Hann window served as the modulation transfer function for the virtual PET. The Hann's apodization properties suppressed high spatial frequencies above a certain critical frequency, thereby fulfilling the above-mentioned boundary conditions. The determined point spread functions were subsequently used by the novel Transconvolution algorithm to match different PET∕CT systems onto the virtual PET system. Finally, the theoretically elaborated Transconvolution method was validated transforming phantom images acquired on two different PET systems to nearly identical data sets, as they would be imaged by the virtual PET system. RESULTS The proposed Transconvolution method matched different PET∕CT-systems for an improved and reproducible determination of a normalized activity concentration. The highest difference in measured activity concentration between the two different PET systems of 18.2% was found in spheres of 2 ml volume. Transconvolution reduced this difference down to 1.6%. In addition to reestablishing comparability the new method with its parameterization of point spread functions allowed a full characterization of imaging properties of the examined tomographs. CONCLUSIONS By matching different tomographs to a virtual standardized imaging system, Transconvolution opens a new comprehensive method for cross calibration in quantitative PET imaging. The use of a virtual PET system restores comparability between data sets from different PET systems by exerting a common, reproducible, and defined partial volume effect.

Influence of the background wind on the local soil moisture-precipitation feedback

The importance of soil moisture anomalies on airmass convection over semiarid regions has been recognized in several studies. The underlying mechanisms remain partly unclear. An open question is why wetter soils can result in either an increase or a decrease of precipitation (positive or negative soil moisture–precipitation feedback, respectively). Here an idealized cloud-resolving modeling framework is used to explore the local soil moisture–precipitation feedback. The approach is able to replicate both positive and negative feedback loops, depending on the environmental parameters. The mechanism relies on horizontal soil moisture variations, which may develop and intensify spontaneously. The positive expression of the feedback is associated with the initiation of convection over dry soil patches, but the convective cells then propagate over wet patches where they strengthen and preferentially precipitate. The negative feedback may occur when the wind profile is too weak to support the propagation of convective features from dry to wet areas. Precipitation is then generally weaker and falls preferentially over dry patches. The results highlight the role of the midtropospheric flow in determining the sign of the feedback. A key element of the positive feedback is the exploitation of both low convective inhibition (CIN) over dry patches (for the initiation of convection) and high CAPE over wet patches (for the generation of precipitation).

Influence of microphysical schemes on atmospheric water in the Weather Research and Forecasting model

This study examines how different microphysical parameterization schemes influence orographically induced precipitation and the distributions of hydrometeors and water vapour for midlatitude summer conditions in the Weather Research and Forecasting (WRF) model. A high-resolution two-dimensional idealized simulation is used to assess the differences between the schemes in which a moist air flow is interacting with a bell-shaped 2 km high mountain. Periodic lateral boundary conditions are chosen to recirculate atmospheric water in the domain. It is found that the 13 selected microphysical schemes conserve the water in the model domain. The gain or loss of water is less than 0.81% over a simulation time interval of 61 days. The differences of the microphysical schemes in terms of the distributions of water vapour, hydrometeors and accumulated precipitation are presented and discussed. The Kessler scheme, the only scheme without ice-phase processes, shows final values of cloud liquid water 14 times greater than the other schemes. The differences among the other schemes are not as extreme, but still they differ up to 79% in water vapour, up to 10 times in hydrometeors and up to 64% in accumulated precipitation at the end of the simulation. The microphysical schemes also differ in the surface evaporation rate. The WRF single-moment 3-class scheme has the highest surface evaporation rate compensated by the highest precipitation rate. The different distributions of hydrometeors and water vapour of the microphysical schemes induce differences up to 49 W m−2 in the downwelling shortwave radiation and up to 33 W m−2 in the downwelling longwave radiation.