Dynamics of turing patterns under spatiotemporal forcing

We study, both theoretically and experimentally, the dynamical response of Turing patterns to a spatiotemporal forcing in the form of a traveling-wave modulation of a control parameter. We show that from strictly spatial resonance, it is possible to induce new, generic dynamical behaviors, including temporally modulated traveling waves and localized traveling solitonlike solutions. The latter make contact with the soliton solutions of Coullet [Phys. Rev. Lett. 56, 724 (1986)] and generalize them. The stability diagram for the different propagating modes in the Lengyel-Epstein model is determined numerically. Direct observations of the predicted solutions in experiments carried out with light modulations in the photosensitive chlorine dioxide-iodine-malonic acid reaction are also reported.

A novel approach to reducing number of sensing units for wearable gait analysis systems.

Gait analysis methods to estimate spatiotemporal measures, based on two, three or four gyroscopes attached on lower limbs have been discussed in the literature. The most common approach to reduce the number of sensing units is to simplify the underlying biomechanical gait model. In this study, we propose a novel method based on prediction of movements of thighs from movements of shanks. Datasets from three previous studies were used. Data from the first study (ten healthy subjects and ten with Parkinson's disease) were used to develop and calibrate a system with only two gyroscopes attached on shanks. Data from two other studies (36 subjects with hip replacement, seven subjects with coxarthrosis, and eight control subjects) were used for comparison with the other methods and for assessment of error compared to a motion capture system. Results show that the error of estimation of stride length compared to motion capture with the system with four gyroscopes and our new method based on two gyroscopes was close ( -0.8 ±6.6 versus 3.8 ±6.6 cm). An alternative with three sensing units did not show better results (error: -0.2 ±8.4 cm). Finally, a fourth that also used two units but with a simpler gait model had the highest bias compared to the reference (error: -25.6 ±7.6 cm). We concluded that it is feasible to estimate movements of thighs from movements of shanks to reduce number of needed sensing units from 4 to 2 in context of ambulatory gait analysis.

Visualizing the spatiotemporal dynamics of DNA damage in budding yeast.

Fluorescence microscopy has enabled the analysis of both the spatial distribution of DNA damage and its dynamics during the DNA damage response (DDR). Three microscopic techniques can be used to study the spatiotemporal dynamics of DNA damage. In the first part we describe how we determine the position of DNA double-strand breaks (DSBs) relative to the nuclear envelope. The second part describes how to quantify the co-localization of DNA DSBs with nuclear pore clusters, or other nuclear subcompartments. The final protocols describe methods for the quantification of locus mobility over time.

Spatiotemporal Features of Coastal Waters in Southwest Finland

In coastal waters, physico-chemical and biological properties and constituents vary at different time scales. In the study area of this thesis, within the Archipelago Sea in the northern Baltic Sea, seasonal cycles of light and temperature set preconditions for intra-annual variations, but developments at other temporal scales occur as well. Weather-induced runoffs and currents may alter water properties over the short term, and the consequences over time of eutrophication and global changes are to a degree unpredictable. The dynamic characteristics of northern Baltic Sea waters are further diversified at the archipelago coasts. Water properties may differ in adjacent basins, which are separated by island and underwater thresholds limiting water exchange, making the area not only a mosaic of islands but also one of water masses. Long-term monitoring and in situ observations provide an essential data reserve for coastal management and research. Since the seasonal amplitudes of water properties are so high, inter-annual comparisons of water-quality variables have to be based on observations sampled at the same time each year. In this thesis I compare areas by their temporal characteristics, using both inter-annual and seasonal data. After comparing spatial differences in seasonal cycles, I conclude that spatial comparisons and temporal generalizations have to be made with caution. In classifying areas by the state of their waters, the results may be biased even if the sampling is annually simultaneous, since the dynamics of water properties may vary according to the area. The most comprehensive view of the spatiotemporal dynamics of water properties would be achieved by means of comparisons with data consisting of multiple annual samples. For practical reasons, this cannot be achieved with conventional in situ sampling. A holistic understanding of the spatiotemporal features of the water properties of the Archipelago Sea will have to be based on the application of multiple methods, complementing each other’s spatial and temporal coverage. The integration of multi-source observational data and time-series analysis may be methodologically challenging, but it will yield new information as to the spatiotemporal regime of the Archipelago Sea.

Three-dimensional kinematic analysis of upper and lower limb motion during gait of post-stroke patients

The aim of this study was to analyze the alterations of arm and leg movements of patients during stroke gait. Joint angles of upper and lower limbs and spatiotemporal variables were evaluated in two groups: hemiparetic group (HG, 14 hemiparetic men, 53 ± 10 years) and control group (CG, 7 able-bodied men, 50 ± 4 years). The statistical analysis was based on the following comparisons (P ≤ 0.05): 1) right versus left sides of CG; 2) affected (AF) versus unaffected (UF) sides of HG; 3) CG versus both the affected and unaffected sides of HG, and 4) an intracycle comparison of the kinematic continuous angular variables between HG and CG. This study showed that the affected upper limb motion in stroke gait was characterized by a decreased range of motion of the glenohumeral (HG: 6.3 ± 4.5, CG: 20.1 ± 8.2) and elbow joints (AF: 8.4 ± 4.4, UF: 15.6 ± 7.6) on the sagittal plane and elbow joint flexion throughout the cycle (AF: 68.2 ± 0.4, CG: 46.8 ± 2.7). The glenohumeral joint presented a higher abduction angle (AF: 14.2 ± 1.6, CG: 11.5 ± 4.0) and a lower external rotation throughout the cycle (AF: 4.6 ± 1.2, CG: 22.0 ± 3.0). The lower limbs showed typical alterations of the stroke gait patterns. Thus, the changes in upper and lower limb motion of stroke gait were identified. The description of upper limb motion in stroke gait is new and complements gait analysis.

Spatiotemporal behavior of the TIGGE medium-range ensemble forecasts

Using the recently-developed mean–variance of logarithms (MVL) diagram, together with the TIGGE archive of medium-range ensemble forecasts from nine different centres, an analysis is presented of the spatiotemporal dynamics of their perturbations, showing how the differences between models and perturbation techniques can explain the shape of their characteristic MVL curves. In particular, a divide is seen between ensembles based on singular vectors or empirical orthogonal functions, and those based on bred vector, Ensemble Transform with Rescaling or Ensemble Kalman Filter techniques. Consideration is also given to the use of the MVL diagram to compare the growth of perturbations within the ensemble with the growth of the forecast error, showing that there is a much closer correspondence for some models than others. Finally, the use of the MVL technique to assist in selecting models for inclusion in a multi-model ensemble is discussed, and an experiment suggested to test its potential in this context.

Automatic spiral analysis for objective assessment of motor symptoms in Parkinson's disease

Objective: To develop a method for objective quantification of PD motor symptoms related to Off episodes and peak dose dyskinesias, using spiral data gathered by using a touch screen telemetry device. The aim was to objectively characterize predominant motor phenotypes (bradykinesia and dyskinesia), to help in automating the process of visual interpretation of movement anomalies in spirals as rated by movement disorder specialists. Background: A retrospective analysis was conducted on recordings from 65 patients with advanced idiopathic PD from nine different clinics in Sweden, recruited from January 2006 until August 2010. In addition to the patient group, 10 healthy elderly subjects were recruited. Upper limb movement data were collected using a touch screen telemetry device from home environments of the subjects. Measurements with the device were performed four times per day during week-long test periods. On each test occasion, the subjects were asked to trace pre-drawn Archimedean spirals, using the dominant hand. The pre-drawn spiral was shown on the screen of the device. The spiral test was repeated three times per test occasion and they were instructed to complete it within 10 seconds. The device had a sampling rate of 10Hz and measured both position and time-stamps (in milliseconds) of the pen tip. Methods: Four independent raters (FB, DH, AJ and DN) used a web interface that animated the spiral drawings and allowed them to observe different kinematic features during the drawing process and to rate task performance. Initially, a number of kinematic features were assessed including ‘impairment’, ‘speed’, ‘irregularity’ and ‘hesitation’ followed by marking the predominant motor phenotype on a 3-category scale: tremor, bradykinesia and/or choreatic dyskinesia. There were only 2 test occasions for which all the four raters either classified them as tremor or could not identify the motor phenotype. Therefore, the two main motor phenotype categories were bradykinesia and dyskinesia. ‘Impairment’ was rated on a scale from 0 (no impairment) to 10 (extremely severe) whereas ‘speed’, ‘irregularity’ and ‘hesitation’ were rated on a scale from 0 (normal) to 4 (extremely severe). The proposed data-driven method consisted of the following steps. Initially, 28 spatiotemporal features were extracted from the time series signals before being presented to a Multilayer Perceptron (MLP) classifier. The features were based on different kinematic quantities of spirals including radius, angle, speed and velocity with the aim of measuring the severity of involuntary symptoms and discriminate between PD-specific (bradykinesia) and/or treatment-induced symptoms (dyskinesia). A Principal Component Analysis was applied on the features to reduce their dimensions where 4 relevant principal components (PCs) were retained and used as inputs to the MLP classifier. Finally, the MLP classifier mapped these components to the corresponding visually assessed motor phenotype scores for automating the process of scoring the bradykinesia and dyskinesia in PD patients whilst they draw spirals using the touch screen device. For motor phenotype (bradykinesia vs. dyskinesia) classification, the stratified 10-fold cross validation technique was employed. Results: There were good agreements between the four raters when rating the individual kinematic features with intra-class correlation coefficient (ICC) of 0.88 for ‘impairment’, 0.74 for ‘speed’, 0.70 for ‘irregularity’, and moderate agreements when rating ‘hesitation’ with an ICC of 0.49. When assessing the two main motor phenotype categories (bradykinesia or dyskinesia) in animated spirals the agreements between the four raters ranged from fair to moderate. There were good correlations between mean ratings of the four raters on individual kinematic features and computed scores. The MLP classifier classified the motor phenotype that is bradykinesia or dyskinesia with an accuracy of 85% in relation to visual classifications of the four movement disorder specialists. The test-retest reliability of the four PCs across the three spiral test trials was good with Cronbach’s Alpha coefficients of 0.80, 0.82, 0.54 and 0.49, respectively. These results indicate that the computed scores are stable and consistent over time. Significant differences were found between the two groups (patients and healthy elderly subjects) in all the PCs, except for the PC3. Conclusions: The proposed method automatically assessed the severity of unwanted symptoms and could reasonably well discriminate between PD-specific and/or treatment-induced motor symptoms, in relation to visual assessments of movement disorder specialists. The objective assessments could provide a time-effect summary score that could be useful for improving decision-making during symptom evaluation of individualized treatment when the goal is to maximize functional On time for patients while minimizing their Off episodes and troublesome dyskinesias.

Automatic spiral analysis for objective assessment of motor symptoms in Parkinson's disease

A challenge for the clinical management of advanced Parkinson’s disease (PD) patients is the emergence of fluctuations in motor performance, which represents a significant source of disability during activities of daily living of the patients. There is a lack of objective measurement of treatment effects for in-clinic and at-home use that can provide an overview of the treatment response. The objective of this paper was to develop a method for objective quantification of advanced PD motor symptoms related to off episodes and peak dose dyskinesia, using spiral data gathered by a touch screen telemetry device. More specifically, the aim was to objectively characterize motor symptoms (bradykinesia and dyskinesia), to help in automating the process of visual interpretation of movement anomalies in spirals as rated by movement disorder specialists. Digitized upper limb movement data of 65 advanced PD patients and 10 healthy (HE) subjects were recorded as they performed spiral drawing tasks on a touch screen device in their home environment settings. Several spatiotemporal features were extracted from the time series and used as inputs to machine learning methods. The methods were validated against ratings on animated spirals scored by four movement disorder specialists who visually assessed a set of kinematic features and the motor symptom. The ability of the method to discriminate between PD patients and HE subjects and the test-retest reliability of the computed scores were also evaluated. Computed scores correlated well with mean visual ratings of individual kinematic features. The best performing classifier (Multilayer Perceptron) classified the motor symptom (bradykinesia or dyskinesia) with an accuracy of 84% and area under the receiver operating characteristics curve of 0.86 in relation to visual classifications of the raters. In addition, the method provided high discriminating power when distinguishing between PD patients and HE subjects as well as had good test-retest reliability. This study demonstrated the potential of using digital spiral analysis for objective quantification of PD-specific and/or treatment-induced motor symptoms.

Using Spatiotemporal Methods to Fill Gaps In Energy Usage Interval Data

Researchers analyzing spatiotemporal or panel data, which varies both in location and over time, often find that their data has holes or gaps. This thesis explores alternative methods for filling those gaps and also suggests a set of techniques for evaluating those gap-filling methods to determine which works best.

Balance and ankle muscle strength predict spatiotemporal gait parameters in individuals with diabetic peripheral neuropathy

The aims of this study were to evaluate aspects of balance, ankle muscle strength and spatiotemporal gait parameters in individuals with diabetic peripheral neuropathy (DPN) and verify whether deficits in spatiotemporal gait parameters were associated with ankle muscle strength and balance performance. Thirty individuals with DPN and 30 control individuals have participated. Spatiotemporal gait parameters were evaluated by measuring the time to walk a set distance during self-selected and maximal walking speeds. Functional mobility and balance performance were assessed using the Functional Reach and the Time Up and Go tests. Ankle isometric muscle strength was assessed with a handheld digital dynamometer. Analyses of variance were employed to verify possible differences between groups and conditions. Multiple linear regression analysis was employed to uncover possible predictors of gait deficits. Gait spatiotemporal, functional mobility, balance performance and ankle muscle strength were affected in individuals with DPN. The Time Up and Go test performance and ankle muscle isometric strength were associated to spatiotemporal gait changes, especially during maximal walking speed condition. Functional mobility and balance performance are damaged in DPN and balance performance and ankle muscle strength can be used to predict spatiotemporal gait parameters in individuals with DPN.

Dynamic texture analysis and segmentation using deterministic partially self-avoiding walks

Dynamic texture is a recent field of investigation that has received growing attention from computer vision community in the last years. These patterns are moving texture in which the concept of selfsimilarity for static textures is extended to the spatiotemporal domain. In this paper, we propose a novel approach for dynamic texture representation, that can be used for both texture analysis and segmentation. In this method, deterministic partially self-avoiding walks are performed in three orthogonal planes of the video in order to combine appearance and motion features. We validate our method on three applications of dynamic texture that present interesting challenges: recognition, clustering and segmentation. Experimental results on these applications indicate that the proposed method improves the dynamic texture representation compared to the state of the art.

The spatiotemporal pattern of auditory cortical responses during verbal hallucinations

Functional magnetic resonance imaging (fMRI) studies can provide insight into the neural correlates of hallucinations. Commonly, such studies require self-reports about the timing of the hallucination events. While many studies have found activity in higher-order sensory cortical areas, only a few have demonstrated activity of the primary auditory cortex during auditory verbal hallucinations. In this case, using self-reports as a model of brain activity may not be sensitive enough to capture all neurophysiological signals related to hallucinations. We used spatial independent component analysis (sICA) to extract the activity patterns associated with auditory verbal hallucinations in six schizophrenia patients. SICA decomposes the functional data set into a set of spatial maps without the use of any input function. The resulting activity patterns from auditory and sensorimotor components were further analyzed in a single-subject fashion using a visualization tool that allows for easy inspection of the variability of regional brain responses. We found bilateral auditory cortex activity, including Heschl's gyrus, during hallucinations of one patient, and unilateral auditory cortex activity in two more patients. The associated time courses showed a large variability in the shape, amplitude, and time of onset relative to the self-reports. However, the average of the time courses during hallucinations showed a clear association with this clinical phenomenon. We suggest that detection of this activity may be facilitated by examining hallucination epochs of sufficient length, in combination with a data-driven approach.

ANALYSIS AND VISUALIZATION OF FLOW FIELDS USING INFORMATION-THEORETIC TECHNIQUES AND GRAPH-BASED REPRESENTATIONS

Three-dimensional flow visualization plays an essential role in many areas of science and engineering, such as aero- and hydro-dynamical systems which dominate various physical and natural phenomena. For popular methods such as the streamline visualization to be effective, they should capture the underlying flow features while facilitating user observation and understanding of the flow field in a clear manner. My research mainly focuses on the analysis and visualization of flow fields using various techniques, e.g. information-theoretic techniques and graph-based representations. Since the streamline visualization is a popular technique in flow field visualization, how to select good streamlines to capture flow patterns and how to pick good viewpoints to observe flow fields become critical. We treat streamline selection and viewpoint selection as symmetric problems and solve them simultaneously using the dual information channel [81]. To the best of my knowledge, this is the first attempt in flow visualization to combine these two selection problems in a unified approach. This work selects streamline in a view-independent manner and the selected streamlines will not change for all viewpoints. My another work [56] uses an information-theoretic approach to evaluate the importance of each streamline under various sample viewpoints and presents a solution for view-dependent streamline selection that guarantees coherent streamline update when the view changes gradually. When projecting 3D streamlines to 2D images for viewing, occlusion and clutter become inevitable. To address this challenge, we design FlowGraph [57, 58], a novel compound graph representation that organizes field line clusters and spatiotemporal regions hierarchically for occlusion-free and controllable visual exploration. We enable observation and exploration of the relationships among field line clusters, spatiotemporal regions and their interconnection in the transformed space. Most viewpoint selection methods only consider the external viewpoints outside of the flow field. This will not convey a clear observation when the flow field is clutter on the boundary side. Therefore, we propose a new way to explore flow fields by selecting several internal viewpoints around the flow features inside of the flow field and then generating a B-Spline curve path traversing these viewpoints to provide users with closeup views of the flow field for detailed observation of hidden or occluded internal flow features [54]. This work is also extended to deal with unsteady flow fields. Besides flow field visualization, some other topics relevant to visualization also attract my attention. In iGraph [31], we leverage a distributed system along with a tiled display wall to provide users with high-resolution visual analytics of big image and text collections in real time. Developing pedagogical visualization tools forms my other research focus. Since most cryptography algorithms use sophisticated mathematics, it is difficult for beginners to understand both what the algorithm does and how the algorithm does that. Therefore, we develop a set of visualization tools to provide users with an intuitive way to learn and understand these algorithms.

A community-based assessment of the spatiotemporal distribution of influenza in a U.S.-Mexico border county during the spring 2009 novel H1N1 swine-origin influenza A pandemic

This study retrospectively evaluated the spatial and temporal disease patterns associated with influenza-like illness (ILI), positive rapid influenza antigen detection tests (RIDT), and confirmed H1N1 S-OIV cases reported to the Cameron County Department of Health and Human Services between April 26 and May 13, 2009 using the space-time permutation scan statistic software SaTScan in conjunction with geographical information system (GIS) software ArcGIS 9.3. The rate and age-adjusted relative risk of each influenza measure was calculated and a cluster analysis was conducted to determine the geographic regions with statistically higher incidence of disease. A Poisson distribution model was developed to identify the effect that socioeconomic status, population density, and certain population attributes of a census block-group had on that area's frequency of S-OIV confirmed cases over the entire outbreak. Predominant among the spatiotemporal analyses of ILI, RIDT and S-OIV cases in Cameron County is the consistent pattern of a high concentration of cases along the southern border with Mexico. These findings in conjunction with the slight northward space-time shifts of ILI and RIDT cluster centers highlight the southern border as the primary site for public health interventions. Finally, the community-based multiple regression model revealed that three factors—percentage of the population under age 15, average household size, and the number of high school graduates over age 25—were significantly associated with laboratory-confirmed S-OIV in the Lower Rio Grande Valley. Together, these findings underscore the need for community-based surveillance, improve our understanding of the distribution of the burden of influenza within the community, and have implications for vaccination and community outreach initiatives.^

Analysis of an ancient global warming event, the Paleocene-Eocene thermal maximum, under examination from ODP Holes 113-689B and 113-690B

Spatiotemporal patterns of carbonate dissolution provide a critical constraint on carbon input during an ancient (~55.5 Ma) global warming event known as the Paleocene-Eocene thermal maximum (PETM), yet the magnitude of lysocline shoaling in the Southern Ocean is poorly constrained due to limited spatial coverage in the circum-Antarctic region. This shortcoming is partially addressed by comparing patterns of carbonate sedimentation at the Site 690 PETM reference section to those herein reconstructed for nearby Site 689. Biochemostratigraphic correlation of the two records reveals that the first ~36 ka of the carbon isotope excursion (CIE) signaling PETM conditions is captured by the Site 689 section, while the remainder of the CIE interval and nearly all of the CIE recovery are missing due to a coring gap. A relatively expanded stratigraphy and higher carbonate content at mid-bathyal Site 689 indicate that dissolution was less severe than at Site 690. Thus, the bathymetric transect delimited by these two PETM records indicates that the lysocline shoaled above Site 689 (~1,100 m) while the calcite compensation depth remained below Site 690 (~1,900 m) in the Weddell Sea region. The ensuing recovery of carbonate sedimentation conforms to a bathymetric trend best explained by gradual lysocline deepening as negative feedback mechanisms neutralized ocean acidification. Further, biochemostratigraphic evidence indicates the tail end of the CIE recovery interval at both sites has been truncated by a hiatus most likely related to vigorous production and advection of intermediate waters.