Polar MM5 Simulations of the Winter Climate of the Laurentide Ice Sheet at the LGM

Optimized regional climate simulations are conducted using the Polar MM5, a version of the fifth-generation Pennsylvania State University-NCAR Mesoscale Model (MM5), with a 60-km horizontal resolution domain over North America during the Last Glacial Maximum (LGM, 21 000 calendar years ago), when much of the continent was covered by the Laurentide Ice Sheet (LIS). The objective is to describe the LGM annual cycle at high spatial resolution with an emphasis on the winter atmospheric circulation. Output from a tailored NCAR Community Climate Model version 3 (CCM3) simulation of the LGM climate is used to provide the initial and lateral boundary conditions for Polar MM5. LGM boundary conditions include continental ice sheets, appropriate orbital forcing, reduced CO2 concentration, paleovegetation, modified sea surface temperatures, and lowered sea level. Polar MM5 produces a substantially different atmospheric response to the LGM boundary conditions than CCM3 and other recent GCM simulations. In particular, from November to April the upper-level flow is split around a blocking anticyclone over the LIS, with a northern branch over the Canadian Arctic and a southern branch impacting southern North America. The split flow pattern is most pronounced in January and transitions into a single, consolidated jet stream that migrates northward over the LIS during summer. Sensitivity experiments indicate that the winter split flow in Polar MM5 is primarily due to mechanical forcing by LIS, although model physics and resolution also contribute to the simulated flow configuration. Polar MM5 LGM results are generally consistent with proxy climate estimates in the western United States, Alaska, and the Canadian Arctic and may help resolve some long-standing discrepancies between proxy data and previous simulations of the LGM climate.

Validation of network communicability metrics for the analysis of brain structural networks.

Computational network analysis provides new methods to analyze the brain's structural organization based on diffusion imaging tractography data. Networks are characterized by global and local metrics that have recently given promising insights into diagnosis and the further understanding of psychiatric and neurologic disorders. Most of these metrics are based on the idea that information in a network flows along the shortest paths. In contrast to this notion, communicability is a broader measure of connectivity which assumes that information could flow along all possible paths between two nodes. In our work, the features of network metrics related to communicability were explored for the first time in the healthy structural brain network. In addition, the sensitivity of such metrics was analysed using simulated lesions to specific nodes and network connections. Results showed advantages of communicability over conventional metrics in detecting densely connected nodes as well as subsets of nodes vulnerable to lesions. In addition, communicability centrality was shown to be widely affected by the lesions and the changes were negatively correlated with the distance from lesion site. In summary, our analysis suggests that communicability metrics that may provide an insight into the integrative properties of the structural brain network and that these metrics may be useful for the analysis of brain networks in the presence of lesions. Nevertheless, the interpretation of communicability is not straightforward; hence these metrics should be used as a supplement to the more standard connectivity network metrics.

A compilation of Western European terrestrial records 60–8 ka BP: towards an understanding of latitudinal climatic gradients

Terrestrial records of past climatic conditions, such as lake sediments and speleothems, provide data of great importance for understanding environmental changes. However, unlike marine and ice core records, terrestrial palaeodata are often not available in databases or in a format that is easily accessible to the non-specialist. As a consequence, many excellent terrestrial records are unknown to the broader palaeoclimate community and are not included in compilations, comparisons, or modelling exercises. Here we present a compilation of Western European terrestrial palaeo-records covering, entirely or partially, the 60–8-ka INTIMATE time period. The compilation contains 56 natural archives, including lake records, speleothems, ice cores, and terrestrial proxies in marine records. The compilation is limited to include records of high temporal resolution and/or records that provide climate proxies or quantitative reconstructions of environmental parameters, such as temperature or precipitation, and that are of relevance and interest to a broader community. We briefly review the different types of terrestrial archives, their respective proxies, their interpretation and their application for palaeoclimatic reconstructions. We also discuss the importance of independent chronologies and the issue of record synchronization. The aim of this exercise is to provide the wider palaeo-community with a consistent compilation of high-quality terrestrial records, to facilitate model-data comparisons, and to identify key areas of interest for future investigations. We use the compilation to investigate Western European latitudinal climate gradients during the deglacial period and, despite of poorly constrained chronologies for the older records, we summarize the main results obtained from NW and SW European terrestrial records before the LGM.

Local changes in neocortical circuit dynamics coincide with the spread of seizures to thalamus in a model of epilepsy

During the generalization of epileptic seizures, pathological activity in one brain area recruits distant brain structures into joint synchronous discharges. However, it remains unknown whether specific changes in local circuit activity are related to the aberrant recruitment of anatomically distant structures into epileptiform discharges. Further, it is not known whether aberrant areas recruit or entrain healthy ones into pathological activity. Here we study the dynamics of local circuit activity during the spread of epileptiform discharges in the zero-magnesium in vitro model of epilepsy. We employ high-speed multi-photon imaging in combination with dual whole-cell recordings in acute thalamocortical (TC) slices of the juvenile mouse to characterize the generalization of epileptic activity between neocortex and thalamus. We find that, although both structures are exposed to zero-magnesium, the initial onset of focal epileptiform discharge occurs in cortex. This suggests that local recurrent connectivity that is particularly prevalent in cortex is important for the initiation of seizure activity. Subsequent recruitment of thalamus into joint, generalized discharges is coincident with an increase in the coherence of local cortical circuit activity that itself does not depend on thalamus. Finally, the intensity of population discharges is positively correlated between both brain areas. This suggests that during and after seizure generalization not only the timing but also the amplitude of epileptiform discharges in thalamus is entrained by cortex. Together these results suggest a central role of neocortical activity for the onset and the structure of pathological recruitment of thalamus into joint synchronous epileptiform discharges.

Computed Ultrasound Tomography in Echo mode (CUTE) of speed of sound for diagnosis and for aberration correction in pulse-echo sonography

Sound speed as a diagnostic marker for various diseases of human tissue has been of interest for a while. Up to now, mostly transmission ultrasound computed tomography (UCT) was able to detect spatially resolved sound speed, and its promise as a diagnostic tool has been demonstrated. However, UCT is limited to acoustically transparent samples such as the breast. We present a novel technique where spatially resolved detection of sound speed can be achieved using conventional pulse-echo equipment in reflection mode. For this purpose, pulse-echo images are acquired under various transmit beam directions and a two-dimensional map of the sound speed is reconstructed from the changing phase of local echoes using a direct reconstruction method. Phantom results demonstrate that a high spatial resolution (1 mm) and contrast (0.5 % of average sound speed) can be achieved suitable for diagnostic purposes. In comparison to previous reflection-mode based methods, CUTE works also in a situation with only diffuse echoes, and its direct reconstruction algorithm enables real-time application. This makes it suitable as an addition to conventional clinical ultrasound where it has the potential to benefit diagnosis in a multimodal approach. In addition, knowledge of the spatial distribution of sound speed allows full aberration correction and thus improved spatial resolution and contrast of conventional B-mode ultrasound. © (2014) COPYRIGHT Society of Photo-Optical Instrumentation Engineers (SPIE). Downloading of the abstract is permitted for personal use only.

Reducing CT radiation dose with iterative reconstruction algorithms: the influence of scan and reconstruction parameters on image quality and CTDIvol

OBJECTIVES In this phantom CT study, we investigated whether images reconstructed using filtered back projection (FBP) and iterative reconstruction (IR) with reduced tube voltage and current have equivalent quality. We evaluated the effects of different acquisition and reconstruction parameter settings on image quality and radiation doses. Additionally, patient CT studies were evaluated to confirm our phantom results. METHODS Helical and axial 256 multi-slice computed tomography scans of the phantom (Catphan(®)) were performed with varying tube voltages (80-140kV) and currents (30-200mAs). 198 phantom data sets were reconstructed applying FBP and IR with increasing iterations, and soft and sharp kernels. Further, 25 chest and abdomen CT scans, performed with high and low exposure per patient, were reconstructed with IR and FBP. Two independent observers evaluated image quality and radiation doses of both phantom and patient scans. RESULTS In phantom scans, noise reduction was significantly improved using IR with increasing iterations, independent from tissue, scan-mode, tube-voltage, current, and kernel. IR did not affect high-contrast resolution. Low-contrast resolution was also not negatively affected, but improved in scans with doses <5mGy, although object detectability generally decreased with the lowering of exposure. At comparable image quality levels, CTDIvol was reduced by 26-50% using IR. In patients, applying IR vs. FBP resulted in good to excellent image quality, while tube voltage and current settings could be significantly decreased. CONCLUSIONS Our phantom experiments demonstrate that image quality levels of FBP reconstructions can also be achieved at lower tube voltages and tube currents when applying IR. Our findings could be confirmed in patients revealing the potential of IR to significantly reduce CT radiation doses.

Time variable Earth’s gravity field from SLR satellites

The time variable Earth’s gravity field contains information about the mass transport within the system Earth, i.e., the relationship between mass variations in the atmosphere, oceans, land hydrology, and ice sheets. For many years, satellite laser ranging (SLR) observations to geodetic satellites have provided valuable information of the low-degree coefficients of the Earth’s gravity field. Today, the Gravity Recovery and Climate Experiment (GRACE) mission is the major source of information for the time variable field of a high spatial resolution. We recover the low-degree coefficients of the time variable Earth’s gravity field using SLR observations up to nine geodetic satellites: LAGEOS-1, LAGEOS-2, Starlette, Stella, AJISAI, LARES, Larets, BLITS, and Beacon-C. We estimate monthly gravity field coefficients up to degree and order 10/10 for the time span 2003–2013 and we compare the results with the GRACE-derived gravity field coefficients. We show that not only degree-2 gravity field coefficients can be well determined from SLR, but also other coefficients up to degree 10 using the combination of short 1-day arcs for low orbiting satellites and 10-day arcs for LAGEOS-1/2. In this way, LAGEOS-1/2 allow recovering zonal terms, which are associated with long-term satellite orbit perturbations, whereas the tesseral and sectorial terms benefit most from low orbiting satellites, whose orbit modeling deficiencies are minimized due to short 1-day arcs. The amplitudes of the annual signal in the low-degree gravity field coefficients derived from SLR agree with GRACE K-band results at a level of 77 %. This implies that SLR has a great potential to fill the gap between the current GRACE and the future GRACE Follow-On mission for recovering of the seasonal variations and secular trends of the longest wavelengths in gravity field, which are associated with the large-scale mass transport in the system Earth.

Understanding measured water rotational temperatures and column densities in the very innermost coma of Comet 73P/Schwassmann–Wachmann 3 B

Direct sublimation of a comet nucleus surface is usually considered to be the main source of gas in the coma of a comet. However, evidence from a number of comets including the recent spectacular images of Comet 103P/Hartley 2 by the EPOXI mission indicates that the nucleus alone may not be responsible for all, or possibly at times even most, of the total amount of gas seen in the coma. Indeed, the sublimation of icy grains, which have been injected into the coma, appears to constitute an important source. We use the fully-kinetic Direct Simulation Monte Carlo model of Tenishev et al. (Tenishev, V.M., Combi, M.R., Davidsson, B. [2008]. Astrophys. J., 685, 659−677; Tenishev, V.M., Combi, M.R., Rubin, M. [2011]. Astrophys. J., 732) to reproduce the measurements of column density and rotational temperature of water in Comet 73P-B/Schwassmann–Wachmann 3 obtained with a very high spatial resolution of ∼30 km using IRCS/Subaru in May 2006 (Bonev, B.P., Mumma, M.J., Kawakita, H., Kobayashi, H., Villanueva, G.L. [2008]. Icarus, 196, 241−248). For gas released solely from the cometary nucleus at a heliocentric distance of 1 AU, modeled rotational temperatures start at 110 K close to the surface and decrease to only several tens of degrees by 10–20 nucleus radii. However, the measured decay of both rotational temperature and column density with distance from the nucleus is much slower than predicted by this simple model. The addition of a substantial (distributed) source of gas from icy grains in the model slows the decay in rotational temperature and provides a more gradual drop in column density profiles. Together with a contribution of rotational heating of water molecules by electrons, the combined effects allow a much better match to the IRCS/Subaru observations. From the spatial distributions of water abundance and temperature measured in 73P/SW3-B, we have identified and quantified multiple mechanisms of release. The application of this tool to other comets may permit such studies over a range of heliocentric and geocentric distances.

Optical recording of calcium currents during impulse conduction in cardiac tissue

We explore the feasibility of obtaining a spatially resolved picture of Ca2+Ca2+ inward currents (ICaICa) in multicellular cardiac tissue by differentiating optically recorded Ca2+Ca2+ transients that accompany propagating action potentials. Patterned growth strands of neonatal rat ventricular cardiomyocytes were stained with the Ca2+Ca2+ indicators Fluo-4 or Fluo-4FF. Preparations were stimulated at 1 Hz, and Ca2+Ca2+ transients were recorded with high spatiotemporal resolution (50  μm50  μm, 2 kHz analog bandwidth) with a photodiode array. Signals were differentiated after appropriate digital filtering. Differentiation of Ca2+Ca2+ transients resulted in optically recorded calcium currents (ORCCs) that carried the temporal and pharmacological signatures of L-type Ca2+Ca2+ inward currents: the time to peak amounted to ∼2.1  ms∼2.1  ms (Fluo-4FF) and ∼2.4  ms∼2.4  ms (Fluo-4), full-width at half-maximum was ∼8  ms∼8  ms, and ORCCs were completely suppressed by 50  μmol/L50  μmol/LCdCl2CdCl2. Also, and as reported before from patch-clamp studies, caffeine reversibly depressed the amplitude of ORCCs. The results demonstrate that the differentiation of Ca2+Ca2+ transients can be used to obtain a spatially resolved picture of the initial phase of ICaICa in cardiac tissue and to assess relative changes of activation/fast inactivation of ICaICa following pharmacological interventions.

Photometric monitoring of non-resolved space debris and databases of optical light curves.

The population of space debris increased drastically during the last years. Collisions involving massive objects may produce large number of fragments leading to significantly growth of the space debris population. An effective remediation measure in order to stabilize the population in LEO, is therefore the removal of large, massive space debris. To remove these objects, not only precise orbits, but also more detailed information about their attitude states will be required. One important property of an object targeted for removal is its spin period and spin axis orientation. If we observe a rotating object, the observer sees different surface areas of the object which leads to changes in the measured intensity. Rotating objects will produce periodic brightness vari ations with frequencies which are related to the spin periods. Photometric monitoring is the real tool for remote diagnostics of the satellite rotation around its center of mass. This information is also useful, for example, in case of contingency. Moreover, it is also important to take into account the orientation of non-spherical body (e.g. space debris) in the numerical integration of its motion when a close approach with the another spacecr aft is predicted. We introduce the two databases of light curves: the AIUB data base, which contains about a thousand light curves of LEO, MEO and high-altitude debris objects (including a few functional objects) obtained over more than seven years, and the data base of the Astronomical Observatory of Odessa University (Ukraine), which contains the results of more than 10 years of photometric monitoring of functioning satellites and large space debris objects in low Earth orbit. AIUB used its 1m ZIMLAT telescope for all light curves. For tracking low-orbit satellites, the Astronomical Observatory of Odessa used the KT-50 telescope, which has an alt-azimuth mount and allows tracking objects moving at a high angular velocity. The diameter of the KT-50 main mirror is 0.5 m, and the focal length is 3 m. The Odessa's Atlas of light curves includes almost 5,5 thousand light curves for ~500 correlated objects from a time period of 2005-2014. The processing of light curves and the determination of the rotation period in the inertial frame is challenging. Extracted frequencies and reconstructed phases for some interesting targets, e.g. GLONASS satellites, for which also SLR data were available for confirmation, will be presented. The rotation of the Envisat satellite after its sudden failure will be analyzed. The deceleration of its rotation rate within 3 years is studied together with the attempt to determine the orientation of the rotation axis.

GROMOS-C, a novel ground-based microwave radiometer for ozone measurement campaigns

Stratospheric ozone is of major interest as it absorbs most harmful UV radiation from the sun, allowing life on Earth. Ground-based microwave remote sensing is the only method that allows for the measurement of ozone profiles up to the mesopause, over 24 hours and under different weather conditions with high time resolution. In this paper a novel ground-based microwave radiometer is presented. It is called GROMOS-C (GRound based Ozone MOnitoring System for Campaigns), and it has been designed to measure the vertical profile of ozone distribution in the middle atmosphere by observing ozone emission spectra at a frequency of 110.836 GHz. The instrument is designed in a compact way which makes it transportable and suitable for outdoor use in campaigns, an advantageous feature that is lacking in present day ozone radiometers. It is operated through remote control. GROMOS-C is a total power radiometer which uses a pre-amplified heterodyne receiver, and a digital fast Fourier transform spectrometer for the spectral analysis. Among its main new features, the incorporation of different calibration loads stands out; this includes a noise diode and a new type of blackbody target specifically designed for this instrument, based on Peltier elements. The calibration scheme does not depend on the use of liquid nitrogen; therefore GROMOS-C can be operated at remote places with no maintenance requirements. In addition, the instrument can be switched in frequency to observe the CO line at 115 GHz. A description of the main characteristics of GROMOS-C is included in this paper, as well as the results of a first campaign at the High Altitude Research Station at Jungfraujoch (HFSJ), Switzerland. The validation is performed by comparison of the retrieved profiles against equivalent profiles from MLS (Microwave Limb Sounding) satellite data, ECMWF (European Centre for Medium-Range Weather Forecast) model data, as well as our nearby NDACC (Network for the Detection of Atmospheric Composition Change) ozone radiometer measuring at Bern.

An automated method for mapping physical soil and water conservation structures on cultivated land using GIS and remote sensing techniques

An efficient and reliable automated model that can map physical Soil and Water Conservation (SWC) structures on cultivated land was developed using very high spatial resolution imagery obtained from Google Earth and ArcGIS, ERDAS IMAGINE, and SDC Morphology Toolbox for MATLAB and statistical techniques. The model was developed using the following procedures: (1) a high-pass spatial filter algorithm was applied to detect linear features, (2) morphological processing was used to remove unwanted linear features, (3) the raster format was vectorized, (4) the vectorized linear features were split per hectare (ha) and each line was then classified according to its compass direction, and (5) the sum of all vector lengths per class of direction per ha was calculated. Finally, the direction class with the greatest length was selected from each ha to predict the physical SWC structures. The model was calibrated and validated on the Ethiopian Highlands. The model correctly mapped 80% of the existing structures. The developed model was then tested at different sites with different topography. The results show that the developed model is feasible for automated mapping of physical SWC structures. Therefore, the model is useful for predicting and mapping physical SWC structures areas across diverse areas.

Palynostratigraphy of the last centuries in Switzerland based on 23 lake and mire deposits: chronostratigraphic pollen markers, regional patterns, and local histories

A total of 23 pollen diagrams [stored in the Alpine Palynological Data-Base (ALPADABA), Geobotanical Institute, Bern] cover the last 100 to over 1000 years. The sites include 15 lakes, seven mires, and one soil profile distributed in the Jura Mts (three sites), Swiss Plateau (two sites), northern Pre-Alps and Alps (six sites), central Alps (five sites), southern Alps (three sites), and southern Pre-Alps (four sites) in the western and southern part of Switzerland or just outside the national borders. The pollen diagrams have both a high taxonomic resolution and a high temporal resolution, with sampling distances of 0.5–3 cm, equivalent to 1 to 11 years for the last 100 years and 8 to 130 years for earlier periods. The chronology is based on absolute dating (14 sites: 210Pb 11 sites; 14C six sites; varve counting two sites) or on biostratigraphic correlation among pollen diagrams. The latter relies mainly on trends in Cannabis sativa, Ambrosia, Mercurialis annua, and Ostrya-type pollen. Individual pollen stratigraphies are discussed and sites are compared within each region. The principle of designating local, extra-local, and regional pollen signals and vegetation is exemplified by two pairs of sites lying close together. Trends in biostratigraphies shared by a major part of the pollen diagrams allow the following generalisations. Forest declined in phases since medieval times up to the late 19th century. Abies and Fagus declined consistently, whereas the behaviour of short-lived trees and trees of moist habitats differed among sites (Alnus glutinosa-type, Alnus viridis, Betula, Corylus avellana). In the present century, however, Picea and Pinus increased, followed by Fraxinus excelsior in the second half of this century. Grassland (traced by Gramineae and Plantago lanceolata-type pollen) increased, replacing much of the forest, and declined again in the second half of this century. Nitrate enrichment of the vegetation (traced by Urtica) took place in the first half of this century. These trends reflect the intensification of forest use and the expansion of grassland from medieval times up to the end of the last century, whereas subsequently parts of the grassland became used more intensively and the marginal parts were abandoned for forest regrowth. In most pollen diagrams human impact is the dominant factor in explaining inferred changes in vegetation, but climatic change plays a role at three sites.

Geomorphology of the Imhotep region on comet 67P/Churyumov-Gerasimenko from OSIRIS observations

Context. Since August 2014, the OSIRIS Narrow Angle Camera (NAC) onboard the Rosetta spacecraft has acquired high spatial resolution images of the nucleus of comet 67P/Churyumov-Gerasimenko, down to the decimeter scale. This paper focuses on the Imhotep region, located on the largest lobe of the nucleus, near the equator. Aims. We map, inventory, and describe the geomorphology of the Imhotep region. We propose and discuss some processes to explain the formation and ongoing evolution of this region. Methods. We used OSIRIS NAC images, gravitational heights and slopes, and digital terrain models to map and measure the morphologies of Imhotep. Results. The Imhotep region presents a wide variety of terrains and morphologies: smooth and rocky terrains, bright areas, linear features, roundish features, and boulders. Gravity processes such as mass wasting and collapse play a significant role in the geomorphological evolution of this region. Cometary processes initiate erosion and are responsible for the formation of degassing conduits that are revealed by elevated roundish features on the surface. We also propose a scenario for the formation and evolution of the Imhotep region; this implies the presence of large primordial voids inside the nucleus, resulting from its formation process.

Detection efficiency of microchannel plates for e(-) and pi(-) in the momentum range from 17.5 to 345 MeV/c

High-energy e(-) and pi(-) were measured by the multichannel plate (MCP) detector at the PiM1 beam line of the High Intensity Proton Accelerator Facilities located at the Paul Scherrer Institute, Villigen, Switzerland. The measurements provide the absolute detection efficiencies for these particles: 5.8% +/- 0.5% for electrons in the beam momenta range 17.5-300 MeV/c and 6.0% +/- 1.3% for pions in the beam momenta range 172-345 MeV/c. The pulse height distribution determined from the measurements is close to an exponential function with negative exponent, indicating that the particles penetrated the MCP material before producing the signal somewhere inside the channel. Low charge extraction and nominal gains of the MCP detector observed in this study are consistent with the proposed mechanism of the signal formation by penetrating radiation. A very similar MCP ion detector will be used in the Neutral Ion Mass (NIM) spectrometer designed for the JUICE mission of European Space Agency (ESA) to the Jupiter system, to perform measurements of the chemical composition of the Galilean moon exospheres. The detection efficiency for penetrating radiation determined in the present studies is important for the optimisation of the radiation shielding of the NIM detector against the high-rate and high-energy electrons trapped in Jupiter's magnetic field. Furthermore, the current studies indicate that MCP detectors can be useful to measure high-energy particle beams at high temporal resolution. (C) 2015 AIP Publishing LLC.