Sex Differences in Semantic Processing: Event-Related Brain Potentials Distinguish between Lower and Higher Order Semantic Analysis during Word Reading

Behavioral studies suggest that women and men differ in the strategic elaboration of verbally encoded information especially in the absence of external task demand. However, measuring such covert processing requires other than behavioral data. The present study used event-related potentials to compare sexes in lower and higher order semantic processing during the passive reading of semantically related and unrelated word pairs. Women and men showed the same early context effect in the P1-N1 transition period. This finding indicates that the initial lexical-semantic access is similar in men and women. In contrast, sexes differed in higher order semantic processing. Women showed an earlier and longer lasting context effect in the N400 accompanied by larger signal strength in temporal networks similarly recruited by men and women. The results suggest that women spontaneously conduct a deeper semantic analysis. This leads to faster processing of related words in the active neural networks as reflected in a shorter stability of the N400 map in women. Taken together, the findings demonstrate that there is a selective sex difference in the controlled semantic analysis during passive word reading that is not reflected in different functional organization but in the depth of processing.

How does spatial extent of fMRI datasets affect independent component analysis decomposition?

Spatial independent component analysis (sICA) of functional magnetic resonance imaging (fMRI) time series can generate meaningful activation maps and associated descriptive signals, which are useful to evaluate datasets of the entire brain or selected portions of it. Besides computational implications, variations in the input dataset combined with the multivariate nature of ICA may lead to different spatial or temporal readouts of brain activation phenomena. By reducing and increasing a volume of interest (VOI), we applied sICA to different datasets from real activation experiments with multislice acquisition and single or multiple sensory-motor task-induced blood oxygenation level-dependent (BOLD) signal sources with different spatial and temporal structure. Using receiver operating characteristics (ROC) methodology for accuracy evaluation and multiple regression analysis as benchmark, we compared sICA decompositions of reduced and increased VOI fMRI time-series containing auditory, motor and hemifield visual activation occurring separately or simultaneously in time. Both approaches yielded valid results; however, the results of the increased VOI approach were spatially more accurate compared to the results of the decreased VOI approach. This is consistent with the capability of sICA to take advantage of extended samples of statistical observations and suggests that sICA is more powerful with extended rather than reduced VOI datasets to delineate brain activity.

Semiparametric Latent Variable Regression Models for Spatio-temporal Modeling of Mobile Source Particles in the Greater Boston Area

Traffic particle concentrations show considerable spatial variability within a metropolitan area. We consider latent variable semiparametric regression models for modeling the spatial and temporal variability of black carbon and elemental carbon concentrations in the greater Boston area. Measurements of these pollutants, which are markers of traffic particles, were obtained from several individual exposure studies conducted at specific household locations as well as 15 ambient monitoring sites in the city. The models allow for both flexible, nonlinear effects of covariates and for unexplained spatial and temporal variability in exposure. In addition, the different individual exposure studies recorded different surrogates of traffic particles, with some recording only outdoor concentrations of black or elemental carbon, some recording indoor concentrations of black carbon, and others recording both indoor and outdoor concentrations of black carbon. A joint model for outdoor and indoor exposure that specifies a spatially varying latent variable provides greater spatial coverage in the area of interest. We propose a penalised spline formation of the model that relates to generalised kringing of the latent traffic pollution variable and leads to a natural Bayesian Markov Chain Monte Carlo algorithm for model fitting. We propose methods that allow us to control the degress of freedom of the smoother in a Bayesian framework. Finally, we present results from an analysis that applies the model to data from summer and winter separately

Practical Large-Scale Spatio-Temporal Modeling of Particulate Matter Concentrations

The last two decades have seen intense scientific and regulatory interest in the health effects of particulate matter (PM). Influential epidemiological studies that characterize chronic exposure of individuals rely on monitoring data that are sparse in space and time, so they often assign the same exposure to participants in large geographic areas and across time. We estimate monthly PM during 1988-2002 in a large spatial domain for use in studying health effects in the Nurses' Health Study. We develop a conceptually simple spatio-temporal model that uses a rich set of covariates. The model is used to estimate concentrations of PM10 for the full time period and PM2.5 for a subset of the period. For the earlier part of the period, 1988-1998, few PM2.5 monitors were operating, so we develop a simple extension to the model that represents PM2.5 conditionally on PM10 model predictions. In the epidemiological analysis, model predictions of PM10 are more strongly associated with health effects than when using simpler approaches to estimate exposure. Our modeling approach supports the application in estimating both fine-scale and large-scale spatial heterogeneity and capturing space-time interaction through the use of monthly-varying spatial surfaces. At the same time, the model is computationally feasible, implementable with standard software, and readily understandable to the scientific audience. Despite simplifying assumptions, the model has good predictive performance and uncertainty characterization.

ESTIMATING TEMPORAL ASSOCIATIONS IN ELECTROCORTICOGRAPHIC (ECoG) TIME SERIES WITH FIRST ORDER PRUNING

Granger causality (GC) is a statistical technique used to estimate temporal associations in multivariate time series. Many applications and extensions of GC have been proposed since its formulation by Granger in 1969. Here we control for potentially mediating or confounding associations between time series in the context of event-related electrocorticographic (ECoG) time series. A pruning approach to remove spurious connections and simultaneously reduce the required number of estimations to fit the effective connectivity graph is proposed. Additionally, we consider the potential of adjusted GC applied to independent components as a method to explore temporal relationships between underlying source signals. Both approaches overcome limitations encountered when estimating many parameters in multivariate time-series data, an increasingly common predicament in today's brain mapping studies.

Multi-detector row CT of the small bowel: peak enhancement temporal window--initial experience

PURPOSE: To prospectively determine quantitatively and qualitatively the timing of maximal enhancement of the normal small-bowel wall by using contrast material-enhanced multi-detector row computed tomography (CT). MATERIALS AND METHODS: This HIPAA-compliant study was approved by the institutional review board. After information on radiation risk was given, written informed consent was obtained from 25 participants with no history of small-bowel disease (mean age, 58 years; 19 men) who had undergone single-level dynamic CT. Thirty seconds after the intravenous administration of contrast material, a serial dynamic acquisition, consisting of 10 images obtained 5 seconds apart, was performed. Enhancement measurements were obtained over time from the small-bowel wall and the aorta. Three independent readers qualitatively assessed small-bowel conspicuity. Quantitative and qualitative data were analyzed during the arterial phase, the enteric phase (which represented peak small-bowel mural enhancement), and the venous phase. Statistical analysis included paired Student t test and Wilcoxon signed rank test with Bonferroni correction. A P value less than .05 was used to indicate a significant difference. RESULTS: The mean time to peak enhancement of the small-bowel wall was 49.3 seconds +/- 7.7 (standard deviation) and 13.5 seconds +/- 7.6 after peak aortic enhancement. Enhancement values were highest during the enteric phase (P < .05). Regarding small-bowel conspicuity, images obtained during the enteric phase were most preferred qualitatively; there was a significant difference between the enteric and arterial phases (P < .001) but not between the enteric and venous phases (P = .18). CONCLUSION: At multi-detector row CT, peak mural enhancement of the normal small bowel occurs on average about 50 seconds after intravenous administration of contrast material or 14 seconds after peak aortic enhancement.

Provenance discrimination of Lower Cretaceous synrift sandstones (eastern Iberian Chain Spain) Constraints from detrital modes heavy minerals and geochemistry

Petrography, geochemical whole-rock composition, and chemical analyses of tourmaline were performed in order to determine the source areas of Lower Cretaceous Mora, El Castellar, and uppermost Camarillas Formation sandstones from the Iberian Chain, Spain. Sandstones were deposited in intraplate subbasins, which are bound by plutonic and volcanic rocks of Permian, Triassic, and Jurassic age, Paleozoic metamorphic rocks, and Triassic sedimentary rocks. Modal analyses together with petrographic and cathodoluminescence observations allowed us to define three quartz-feldspathic petrofacies and recognize diagenetic processes that modified the original framework composition. Results from average restored petrofacies are: Mora petrofacies = P/F >1 and Q(r)70 F(r)22 R(r)9; El Castellar petrofacies = P/F >1 and Q(r)57 F(r)25 R(r)18; and Camarillas petrofacies = P/F ∼ zero and Q(r)64 F(r)28 R(r)7 (P—plagioclase; F—feldspar; Q—quartz; R—rock fragments; r—restored composition). Trace-element and rare earth element abundances of whole-rock analyses discriminate well between the three petrofacies based on: (1) the Rb concentration, which is indicative of the K content and reflects the amount of K-feldspar modal abundance, and (2) the relative modal abundance of heavy minerals (tourmaline, zircon, titanite, and apatite), which is reproduced by the elements hosted in the observed heavy mineral assemblage (i.e., B and Li for tourmaline; Zr, Hf, and Ta for zircon; Ti, Ta, Nb, and their rare earth elements for titanite; and P, Y, and their rare earth elements for apatite). Tourmaline chemical composition for the three petrofacies ranges from Fe-tourmaline of granitic to Mg-tourmaline of metamorphic origin. The three defined petrofacies suggest a mixed provenance from plutonic and metamorphic source rocks. However, a progressively major influence of granitic source rocks was detected from the lowermost Mora petrofacies toward the uppermost Camarillas petrofacies. This provenance trend is consistent with the uplift and erosion of the Iberian Massif, which coincided with the development of the latest Berriasian synrift regional unconformity and affected all of the Iberian intraplate basins. The uplifting stage of Iberian Massif pluton caused a significant dilution of Paleozoic metamorphic source areas, which were dominant during the sedimentation of the lowermost Mora and El Castellar petrofacies. The association of petrographic data with whole-rock geochemical compositions and tourmaline chemical analysis has proved to be useful for determining source area characteristics, their predominance, and the evolution of source rock types during the deposition of quartz-feldspathic sandstones in intraplate basins. This approach ensures that provenance interpretation is consistent with the geological context.

BOLD correlates of continuously fluctuating epileptic activity isolated by independent component analysis

Combined EEG/fMRI recordings offer a promising opportunity to detect brain areas with altered BOLD signal during interictal epileptic discharges (IEDs). These areas are likely to represent the irritative zone, which is itself a reflection of the epileptogenic zone. This paper reports on the imaging findings using independent component analysis (ICA) to continuously quantify epileptiform activity in simultaneously acquired EEG and fMRI. Using ICA derived factors coding for the epileptic activity takes into account that epileptic activity is continuously fluctuating with each spike differing in amplitude, duration and maybe topography, including subthreshold epileptic activity besides clear IEDs and may thus increase the sensitivity and statistical power of combined EEG/fMRI in epilepsy. Twenty patients with different types of focal and generalized epilepsy syndromes were investigated. ICA separated epileptiform activity from normal physiological brain activity and artifacts. In 16/20 patients, BOLD correlates of epileptic activity matched the EEG sources, the clinical semiology, and, if present, the structural lesions. In clinically equivocal cases, the BOLD correlates aided to attribute proper diagnosis of the underlying epilepsy syndrome. Furthermore, in one patient with temporal lobe epilepsy, BOLD correlates of rhythmic delta activity could be employed to delineate the affected hippocampus. Compared to BOLD correlates of manually identified IEDs, the sensitivity was improved from 50% (10/20) to 80%. The ICA EEG/fMRI approach is a safe, non-invasive and easily applicable technique, which can be used to identify regions with altered hemodynamic effects related to IEDs as well as intermittent rhythmic discharges in different types of epilepsy.

Estimation of scots pine defoliation by the common pine sawfly (Diprion pini L.) using multi-temporal radar data

In 1998-2001 Finland suffered the most severe insect outbreak ever recorded, over 500,000 hectares. The outbreak was caused by the common pine sawfly (Diprion pini L.). The outbreak has continued in the study area, Palokangas, ever since. To find a good method to monitor this type of outbreaks, the purpose of this study was to examine the efficacy of multi-temporal ERS-2 and ENVISAT SAR imagery for estimating Scots pine (Pinus sylvestris L.) defoliation. Three methods were tested: unsupervised k-means clustering, supervised linear discriminant analysis (LDA) and logistic regression. In addition, I assessed if harvested areas could be differentiated from the defoliated forest using the same methods. Two different speckle filters were used to determine the effect of filtering on the SAR imagery and subsequent results. The logistic regression performed best, producing a classification accuracy of 81.6% (kappa 0.62) with two classes (no defoliation, >20% defoliation). LDA accuracy was with two classes at best 77.7% (kappa 0.54) and k-means 72.8 (0.46). In general, the largest speckle filter, 5 x 5 image window, performed best. When additional classes were added the accuracy was usually degraded on a step-by-step basis. The results were good, but because of the restrictions in the study they should be confirmed with independent data, before full conclusions can be made that results are reliable. The restrictions include the small size field data and, thus, the problems with accuracy assessment (no separate testing data) as well as the lack of meteorological data from the imaging dates.

Automation of shear-wave splitting parameter determination of local earthquakes at Yellowstone : application as indicator of crustal stress and temporal variation

Shear-wave splitting can be a useful technique for determining crustal stress fields in volcanic settings and temporal variations associated with activity. Splitting parameters were determined for a subset of local earthquakes recorded from 2000-2010 at Yellowstone. Analysis was automated using an unsupervised cluster analysis technique to determine optimum splitting parameters from 270 analysis windows for each event. Six stations clearly exhibit preferential fast polarization values sub-orthogonal to the direction of minimum horizontal compression. Yellowstone deformation results in a local crustal stress field differing from the regional field dominated by NE-SW extension, and fast directions reflect this difference rotating around the caldera maintaining perpendicularity to the rim. One station exhibits temporal variations concordant with identified periods of caldera subsidence and uplift. From splitting measurements, we calculated a crustal anisotropy of ~17-23% and crack density ~0.12-0.17 possibly resulting from stress-aligned fluid filled microcracks in the upper crust and an active hydrothermal system.

Triggering of volcanic activity by large earthquakes

Statistical analyses of temporal relationships between large earthquakes and volcanic eruptions suggest seismic waves may trigger eruptions even over great (>1000 km) distances, although the causative mechanism is not well constrained. In this study the relationship between large earthquakes and subtle changes in volcanic activity was investigated in order to gain greater insight into the relationship between dynamic stresses propagated by surface waves and volcanic response. Daily measurements from the Ozone Monitoring Instrument (OMI), onboard the Aura satellite, provide constraints on volcanic sulfur-dioxide (SO2) emission rates as a measure of subtle changes in activity. Time series of SO2 emission rates were produced from OMI data for thirteen persistently active volcanoes from 1 October 2004 to 30 September 2010. In order to quantify the affect of earthquakes at teleseismic distances, we modeled surface-wave amplitudes from the source mechanisms of moment magnitude (Mw) ≥7 earthquakes, and calculated the Peak Dynamic Stress (PDS). We assessed the influence of earthquakes on volcanic activity in two ways: 1) by identifying increases in the SO2 time series data and looking for causative earthquakes and 2) by examining the average emission rate before and after each earthquake. In the first, the SO2 time series for each volcano was used to calculate a baseline threshold for comparison with post-earthquake emission. Next, we generated a catalog of responses based on sustained SO2 emission increases above this baseline. Delay times between each SO2 response and each prior earthquake were analyzed using both the actual earthquake catalog, and a randomly generated catalog of earthquakes. This process was repeated for each volcano. Despite varying multiple parameters, this analysis did not demonstrate a clear relationship between earthquake-generated PDS and SO2 emission. However, the second analysis, which was based on the occurrence of large earthquakes indicated a response at most volcanoes. Using the PDS calculations as a filtering criterion for the earthquake catalog, the SO2 mass for each volcano was analyzed in 28-day windows centered on the earthquake origin time. If the average SO2 mass after the earthquake was greater than an arbitrary percentage of pre-earthquake mass, we identified the volcano as having a response to the event. This window analysis provided insight on what type of volcanic activity is more susceptible to triggering by dynamic stress. The volcanoes with very open systems included in this study, Ambrym, Gaua, Villarrica, Erta Ale and, Turrialba, showed a clear response to dynamic stress while the volcanoes with more closed systems, Merapi, Semeru, Fuego, Pacaya, and Bagana, showed no response.

Quantitative OCT analysis of idiopathic perifoveal telangiectasia

PURPOSE: To identify and quantitate specific changes in optical coherence tomography (OCT) images of patients with type 2 idiopathic perifoveal telangiectasia (IPT). METHODS: In a prospectively designed, observational, case-control study, 28 eyes of 14 consecutive patients with IPT were examined with OCT and compared with eyes of 14 unaffected control subjects. Light reflectivity profiles of raw scan data of OCT images were quantitatively analyzed for differences in distance between different retinal reflectivity layers and their respective reflectivities. Maculae were examined in four separate regions: (1) central fovea, (2) nasal perifovea, (3) temporal perifovea, and (4) outside the fovea. RESULTS: Retinal thinning, shortening of the photoreceptor outer segments and loss of reflectivity of the photoreceptor ellipsoid region were found in the central foveal region as well as the nasal and temporal perifoveal regions in eyes with IPT. In addition, increased reflectivity of the outer nuclear layer was found in a sharply demarcated area of the inferotemporal perifoveal region in all affected eyes. Retinal tissue located more than 2000 mum away from the foveola was indistinguishable from that in normal eyes. CONCLUSIONS: Quantitative OCT analysis shows unique and specific changes in the photoreceptors of the central macula in IPT which can be detected from first clinical presentation. These changes may be of use as an additional diagnostic tool. Correlation of the findings in the outer nuclear layer with histologic studies may help identify the nature of the reflectivity increase and define more clearly the type of damage sustained by the photoreceptors in this condition.

HIGH PERFORMANCE, LOW COST SUBSPACE DECOMPOSITION AND POLYNOMIAL ROOTING FOR REAL TIME DIRECTION OF ARRIVAL ESTIMATION: ANALYSIS AND IMPLEMENTATION

This thesis develops high performance real-time signal processing modules for direction of arrival (DOA) estimation for localization systems. It proposes highly parallel algorithms for performing subspace decomposition and polynomial rooting, which are otherwise traditionally implemented using sequential algorithms. The proposed algorithms address the emerging need for real-time localization for a wide range of applications. As the antenna array size increases, the complexity of signal processing algorithms increases, making it increasingly difficult to satisfy the real-time constraints. This thesis addresses real-time implementation by proposing parallel algorithms, that maintain considerable improvement over traditional algorithms, especially for systems with larger number of antenna array elements. Singular value decomposition (SVD) and polynomial rooting are two computationally complex steps and act as the bottleneck to achieving real-time performance. The proposed algorithms are suitable for implementation on field programmable gated arrays (FPGAs), single instruction multiple data (SIMD) hardware or application specific integrated chips (ASICs), which offer large number of processing elements that can be exploited for parallel processing. The designs proposed in this thesis are modular, easily expandable and easy to implement. Firstly, this thesis proposes a fast converging SVD algorithm. The proposed method reduces the number of iterations it takes to converge to correct singular values, thus achieving closer to real-time performance. A general algorithm and a modular system design are provided making it easy for designers to replicate and extend the design to larger matrix sizes. Moreover, the method is highly parallel, which can be exploited in various hardware platforms mentioned earlier. A fixed point implementation of proposed SVD algorithm is presented. The FPGA design is pipelined to the maximum extent to increase the maximum achievable frequency of operation. The system was developed with the objective of achieving high throughput. Various modern cores available in FPGAs were used to maximize the performance and details of these modules are presented in detail. Finally, a parallel polynomial rooting technique based on Newton’s method applicable exclusively to root-MUSIC polynomials is proposed. Unique characteristics of root-MUSIC polynomial’s complex dynamics were exploited to derive this polynomial rooting method. The technique exhibits parallelism and converges to the desired root within fixed number of iterations, making this suitable for polynomial rooting of large degree polynomials. We believe this is the first time that complex dynamics of root-MUSIC polynomial were analyzed to propose an algorithm. In all, the thesis addresses two major bottlenecks in a direction of arrival estimation system, by providing simple, high throughput, parallel algorithms.

Comparison of voxel-based 3-D MRI analysis and subtraction ictal SPECT coregistered to MRI in focal epilepsy

While voxel-based 3-D MRI analysis methods as well as assessment of subtracted ictal versus interictal perfusion studies (SISCOM) have proven their potential in the detection of lesions in focal epilepsy, a combined approach has not yet been reported. The present study investigates if individual automated voxel-based 3-D MRI analyses combined with SISCOM studies contribute to an enhanced detection of mesiotemporal epileptogenic foci. Seven consecutive patients with refractory complex partial epilepsy were prospectively evaluated by SISCOM and voxel-based 3-D MRI analysis. The functional perfusion maps and voxel-based statistical maps were coregistered in 3-D space. In five patients with temporal lobe epilepsy (TLE), the area of ictal hyperperfusion and corresponding structural abnormalities detected by 3-D MRI analysis were identified within the same temporal lobe. In two patients, additional structural and functional abnormalities were detected beyond the mesial temporal lobe. Five patients with TLE underwent epileptic surgery with favourable postoperative outcome (Engel class Ia and Ib) after 3-5 years of follow-up, while two patients remained on conservative treatment. In summary, multimodal assessment of structural abnormalities by voxel-based analysis and SISCOM may contribute to advanced observer-independent preoperative assessment of seizure origin.

Effects of Chemical Deicers on Amphibian Communites

Roads and highways present a unique challenge to wildlife as they exhibit substantial impacts on the surrounding ecosystem through the interruption of a number of ecological processes. With new roads added to the national highway system every year, an understanding of these impacts is required for effective mitigation of potential environmental impacts. A major contributor to these negative effects is the deposition of chemicals used in winter deicing activities to nearby surface waters. These chemicals often vary in composition and may affect freshwater species differently. The negative impacts of widespread deposition of sodium chloride (NaCl) have prompted a search for an `environmentally friendly' alternative. However, little research has investigated the potential environmental effects of widespread use of these alternatives. Herein, I detail the results of laboratory tests and field surveys designed to determine the impacts of road salt (NaCl) and other chemical deicers on amphibian communities in Michigan's Upper Peninsula. Using larval amphibians I demonstrate the lethal impacts of a suite of chemical deicers on this sensitive, freshwater species. Larval wood frogs (Lithobates sylvatica) were tolerant of short-term (96 hours) exposure to urea (CH4N2O), sodium chloride (NaCl), and magnesium chloride (MgCl2). However, these larvae were very sensitive to acetate products (C8H12CaMgO8, CH3COOK) and calcium chloride (CaCl2). These differences in tolerance suggest that certain deicers may be more harmful to amphibians than others. Secondly, I expanded this analysis to include an experiment designed to determine the sublethal effects of chronic exposure to environmentally realistic concentrations of NaCl on two unique amphibian species, L. sylvatica and green frogs (L. clamitans). L. sylvatica tend to breed in small, ephemeral wetlands and metamorphose within a single season. However, L. clamitans breed primarily in more permanent wetlands and often remain as tadpoles for one year or more. These species employ different life history strategies in this region which may influence their response to chronic NaCl exposure. Both species demonstrated potentially harmful effects on individual fitness. L. sylvatica larvae had a high incidence of edema suggesting the NaCl exposure was a significant physiologic stressor to these larvae. L. clamitans larvae reduced tail length during their exposure which may affect adult fitness of these individuals. In order to determine the risk local amphibians face when using these roadside pools, I conducted a survey of the spatial distribution of chloride in the three northernmost counties of Michigan. This area receives a relatively low amount of NaCl which is confined to state and federal highways. The chloride concentrations in this region were much lower than those in urban systems; however, amphibians breeding in the local area may encounter harmful chloride levels arising from temporal variations in hydroperiods. Spatial variation of chloride levels suggests the road-effect zone for amphibians may be as large as 1000 m from a salt-treated highway. Lastly, I performed an analysis of the use of specific conductance to predict chloride concentrations in natural surface water bodies. A number of studies have used this regression to predict chloride concentrations from measurements of specific conductance. This method is often chosen in the place of ion chromatography due to budget and time constraints. However, using a regression method to characterize this relationship does not result in accurate chloride ion concentration estimates.