Multispectral Bayesian reconstruction technique for real-time two color fluorescence microscopy

We have developed a real-time imaging method for two-color wide-field fluorescence microscopy using a combined approach that integrates multi-spectral imaging and Bayesian image reconstruction technique. To enable simultaneous observation of two dyes (primary and secondary), we exploit their spectral properties that allow parallel recording in both the channels. The key advantage of this technique is the use of a single wavelength of light to excite both the primary dye and the secondary dye. The primary and secondary dyes respectively give rise to fluorescence and bleed-through signal, which after normalization were merged to obtain two-color 3D images. To realize real-time imaging, we employed maximum likelihood (ML) and maximum a posteriori (MAP) techniques on a high-performance computing platform (GPU). The results show two-fold improvement in contrast while the signal-to-background ratio (SBR) is improved by a factor of 4. We report a speed boost of 52 and 350 for 2D and 3D images respectively. Using this system, we have studied the real-time protein aggregation in yeast cells and HeLa cells that exhibits dot-like protein distribution. The proposed technique has the ability to temporally resolve rapidly occurring biological events.

Studies on the effects of varying secondary gas properties in a low entrainment ratio supersonic ejector

A dilution cum purge ejector for application in fuel cells represents a domain of ejector operation involving low entrainment ratio with differing secondary and primary gas; which is hardly investigated and a cohesive design framework is not readily available. We comprehensively study a constant area ejector using analytical, experimental and numerical tools at low entrainment ratio (0.004-0.065) with Air, Helium and Argon as secondary gas while the primary gas is Air. For the first time, limits of operating parameters used in control volume method to design the ejector are found to be highly dependent on the secondary molecular weight. The entrainment ratio in the ejector (low for Helium and high for Argon) is affected by the molecular weight and the static pressure within the ejector (low for Air and high for Argon & Helium) by the gamma of the secondary gas. Sufficient suction pressure (0.3-0.55 bar) is generated by the ejector thereby preventing any backflow of secondary gas at all primary stagnation pressures (1.5, 2.2 and 3.1 bar). Numerical results agree well with experimental results. The ejector is shown to completely dilute and purge the secondary flow, meeting all key design requirements. (C) 2014 Elsevier Ltd. All rights reserved.

Increasing Critical Sensitivity of the Load/Unload Response Ratio before Large Earthquakes with Identified Stress Acclumulation Pattern

The Load/Unload Response Ratio (LURR) method is proposed for short-to-intermediate-term earthquake prediction [Yin, X.C., Chen, X.Z., Song, Z.P., Yin, C., 1995. A New Approach to Earthquake Prediction — The Load/Unload Response Ratio (LURR) Theory, Pure Appl. Geophys., 145, 701–715]. This method is based on measuring the ratio between Benioff strains released during the time periods of loading and unloading, corresponding to the Coulomb Failure Stress change induced by Earth tides on optimally oriented faults. According to the method, the LURR time series usually climb to an anomalously high peak prior to occurrence of a large earthquake. Previous studies have indicated that the size of critical seismogenic region selected for LURR measurements has great influence on the evaluation of LURR. In this study, we replace the circular region usually adopted in LURR practice with an area within which the tectonic stress change would mostly affect the Coulomb stress on a potential seismogenic fault of a future event. The Coulomb stress change before a hypothetical earthquake is calculated based on a simple back-slip dislocation model of the event. This new algorithm, by combining the LURR method with our choice of identified area with increased Coulomb stress, is devised to improve the sensitivity of LURR to measure criticality of stress accumulation before a large earthquake. Retrospective tests of this algorithm on four large earthquakes occurred in California over the last two decades show remarkable enhancement of the LURR precursory anomalies. For some strong events of lesser magnitudes occurred in the same neighborhoods and during the same time periods, significant anomalies are found if circular areas are used, and are not found if increased Coulomb stress areas are used for LURR data selection. The unique feature of this algorithm may provide stronger constraints on forecasts of the size and location of future large events.

Load-Unload Response Ratio and Accelerating Moment/Energy Release Critical Region Scaling and Earthquake Prediction

The main idea of the Load-Unload Response Ratio (LURR) is that when a system is stable, its response to loading corresponds to its response to unloading, whereas when the system is approaching an unstable state, the response to loading and unloading becomes quite different. High LURR values and observations of Accelerating Moment/Energy Release (AMR/AER) prior to large earthquakes have led different research groups to suggest intermediate-term earthquake prediction is possible and imply that the LURR and AMR/AER observations may have a similar physical origin. To study this possibility, we conducted a retrospective examination of several Australian and Chinese earthquakes with magnitudes ranging from 5.0 to 7.9, including Australia's deadly Newcastle earthquake and the devastating Tangshan earthquake. Both LURR values and best-fit power-law time-to-failure functions were computed using data within a range of distances from the epicenter. Like the best-fit power-law fits in AMR/AER, the LURR value was optimal using data within a certain epicentral distance implying a critical region for LURR. Furthermore, LURR critical region size scales with mainshock magnitude and is similar to the AMR/AER critical region size. These results suggest a common physical origin for both the AMR/AER and LURR observations. Further research may provide clues that yield an understanding of this mechanism and help lead to a solid foundation for intermediate-term earthquake prediction.

“Deborah Numbers”, Coupling Multiple Space and Time Scales and Governing Damage Evolution to Failure

Two different spatial levels are involved concerning damage accumulation to eventual failure. nucleation and growth rates of microdamage nN* and V*. It is found that the trans-scale length ratio c*/L does not directly affect the process. Instead, two independent dimensionless numbers: the trans-scale one * * ( V*)including the * **5 * N c V including mesoscopic parameters only, play the key role in the process of damage accumulation to failure. The above implies that there are three time scales involved in the process: the macroscopic imposed time scale tim = /a and two meso-scopic time scales, nucleation and growth of damage, (* *4) N N t =1 n c and tV=c*/V*. Clearly, the dimensionless number De*=tV/tim refers to the ratio of microdamage growth time scale over the macroscopically imposed time scale. So, analogous to the definition of Deborah number as the ratio of relaxation time over external one in rheology. Let De be the imposed Deborah number while De represents the competition and coupling between the microdamage growth and the macroscopically imposed wave loading. In stress-wave induced tensile failure (spallation) De* < 1, this means that microdamage has enough time to grow during the macroscopic wave loading. Thus, the microdamage growth appears to be the predominate mechanism governing the failure. Moreover, the dimensionless number D* = tV/tN characterizes the ratio of two intrinsic mesoscopic time scales: growth over nucleation. Similarly let D be the “intrinsic Deborah number”. Both time scales are relevant to intrinsic relaxation rather than imposed one. Furthermore, the intrinsic Deborah number D* implies a certain characteristic damage. In particular, it is derived that D* is a proper indicator of macroscopic critical damage to damage localization, like D* ∼ (10–3~10–2) in spallation. More importantly, we found that this small intrinsic Deborah number D* indicates the energy partition of microdamage dissipation over bulk plastic work. This explains why spallation can not be formulated by macroscopic energy criterion and must be treated by multi-scale analysis.

Second-to-Fourth Digit Ratio Has a Non-Monotonic Impact on Altruism

Gene-culture co-evolution emphasizes the joint role of culture and genes for the emergence of altruistic and cooperative behaviors and behavioral genetics provides estimates of their relative importance. However, these approaches cannot assess which biological traits determine altruism or how. We analyze the association between altruism in adults and the exposure to prenatal sex hormones, using the second-to-fourth digit ratio. We find an inverted U-shaped relation for left and right hands, which is very consistent for men and less systematic for women. Subjects with both high and low digit ratios give less than individuals with intermediate digit ratios. We repeat the exercise with the same subjects seven months later and find a similar association, even though subjects' behavior differs the second time they play the game. We then construct proxies of the median digit ratio in the population (using more than 1000 different subjects), show that subjects' altruism decreases with the distance of their ratio to these proxies. These results provide direct evidence that prenatal events contribute to the variation of altruistic behavior and that the exposure to fetal hormones is one of the relevant biological factors. In addition, the findings suggest that there might be an optimal level of exposure to these hormones from social perspective.

Space-time correlations of fluctuating velocities in turbulent shear flows

Space-time correlations or Eulerian two-point two-time correlations of fluctuating velocities are analytically and numerically investigated in turbulent shear flows. An elliptic model for the space-time correlations in the inertial range is developed from the similarity assumptions on the isocorrelation contours: they share a uniform preference direction and a constant aspect ratio. The similarity assumptions are justified using the Kolmogorov similarity hypotheses and verified using the direct numerical simulation DNS of turbulent channel flows. The model relates the space-time correlations to the space correlations via the convection and sweeping characteristic velocities. The analytical expressions for the convection and sweeping velocities are derived from the Navier-Stokes equations for homogeneous turbulent shear flows, where the convection velocity is represented by the mean velocity and the sweeping velocity is the sum of the random sweeping velocity and the shearinduced velocity. This suggests that unlike Taylor’s model where the convection velocity is dominating and Kraichnan and Tennekes’ model where the random sweeping velocity is dominating, the decorrelation time scales of the space-time correlations in turbulent shear flows are determined by the convection velocity, the random sweeping velocity, and the shear-induced velocity. This model predicts a universal form of the spacetime correlations with the two characteristic velocities. The DNS of turbulent channel flows supports the prediction: the correlation functions exhibit a fair good collapse, when plotted against the normalized space and time separations defined by the elliptic model.

Complex Transition of Double-Diffusive Convection in a Rectangular Enclosure with Height-to-Length Ratio Equal to 4: Part I

This is the first part of direct numerical simulation (DNS) of double-diffusive convection in a slim rectangular enclosure with horizontal temperature and concentration gradients. We consider the case with the thermal Rayleigh number of 10^5, the Pradtle number of 1, the Lewis number of 2, the buoyancy ratio of composition to temperature being in the range of [0,1], and height-to-width aspect ration of 4. A new 7th order upwind compact scheme was developed for approximation of convective terms, and a three-stage third-order Runge-Kutta method was employed for time advancement. Our DNS suggests that with the buoyancy ratio increasing form 0 to 1, the flow of transition is a complex series changing fromthe steady to periodic, chaotic, periodic, quasi-periodic, and finally back to periodic. There are two types of periodic flow, one is simple periodic flow with single fundamental frequency (FF), and another is complex periodic flow with multiple FFs. This process is illustrated by using time-velocity histories, Fourier frequency spectrum analysis and the phase-space rajectories.

Hunting spectro-temporal information in unevenly spaced paleoclimate time series

18 p.

Real-Time Bayesian Analysis of Ground Motion Envelopes for Earthquake Early Warning

Current earthquake early warning systems usually make magnitude and location predictions and send out a warning to the users based on those predictions. We describe an algorithm that assesses the validity of the predictions in real-time. Our algorithm monitors the envelopes of horizontal and vertical acceleration, velocity, and displacement. We compare the observed envelopes with the ones predicted by Cua & Heaton's envelope ground motion prediction equations (Cua 2005). We define a "test function" as the logarithm of the ratio between observed and predicted envelopes at every second in real-time. Once the envelopes deviate beyond an acceptable threshold, we declare a misfit. Kurtosis and skewness of a time evolving test function are used to rapidly identify a misfit. Real-time kurtosis and skewness calculations are also inputs to both probabilistic (Logistic Regression and Bayesian Logistic Regression) and nonprobabilistic (Least Squares and Linear Discriminant Analysis) models that ultimately decide if there is an unacceptable level of misfit. This algorithm is designed to work at a wide range of amplitude scales. When tested with synthetic and actual seismic signals from past events, it works for both small and large events.

The potential use of time-area closures to reduce catches of bigeye tuna (Thunnus obesus) in the purse-seine fishery of the eastern Pacific Ocean

Skipjack (Katsuwonus pelamis), yellowfin (Thunnus albacares), and bigeye (Thunnus obesus) tunas are caught by purse-seine vessels in the eastern Pacific Ocean (EPO). Although there is no evidence to indicate that current levels of fishing-induced mortality will affect the sustainability of skipjack or yellowfin tunas, fishing mortality on juvenile (younger than 5 years of age) bigeye tuna has increased, and overall fishing mortality is greater than that necessary to produce the maximum sustainable yield of this species. We investigated whether time-area closures have the potential to reduce purse-seine bigeye catches without significantly reducing skipjack catches. Using catch and effort data for 1995–2002, we identified regions where the ratio of bigeye to skipjack tuna catches was high and applied simple closed-area models to investigate the possible benefits of time-area closures. We estimated that the most optimistic and operationally feasible 3-month closures, covering the equatorial region of the EPO during the third quarter of the year, could reduce bigeye catches by 11.5%, while reducing skipjack tuna catches by 4.3%. Because this level of bigeye tuna catch reduction is insufficient to address sustainability concerns, and larger and longer closures would reduce catches of this species signficantly, we recommend that future research be directed toward gear technology solutions because these have been successful in many other fisheries. In particular, because over 50% of purse-seine catches of bigeye tuna are taken in sets in which bigeye tuna are the dominant species, methods to allow the determination of the species composition of aggregations around floating objects may be important.

A general framework for integrating environmental time series into stock assessment models: model description, simulation testing, and example

We present a method to integrate environmental time series into stock assessment models and to test the significance of correlations between population processes and the environmental time series. Parameters that relate the environmental time series to population processes are included in the stock assessment model, and likelihood ratio tests are used to determine if the parameters improve the fit to the data significantly. Two approaches are considered to integrate the environmental relationship. In the environmental model, the population dynamics process (e.g. recruitment) is proportional to the environmental variable, whereas in the environmental model with process error it is proportional to the environmental variable, but the model allows an additional temporal variation (process error) constrained by a log-normal distribution. The methods are tested by using simulation analysis and compared to the traditional method of correlating model estimates with environmental variables outside the estimation procedure. In the traditional method, the estimates of recruitment were provided by a model that allowed the recruitment only to have a temporal variation constrained by a log-normal distribution. We illustrate the methods by applying them to test the statistical significance of the correlation between sea-surface temperature (SST) and recruitment to the snapper (Pagrus auratus) stock in the Hauraki Gulf–Bay of Plenty, New Zealand. Simulation analyses indicated that the integrated approach with additional process error is superior to the traditional method of correlating model estimates with environmental variables outside the estimation procedure. The results suggest that, for the snapper stock, recruitment is positively correlated with SST at the time of spawning.

Coastal climate reflected in carbon-13/carbon-12 ratio of organic carbon in varved sediment from Santa Barbara basin

A 1844-1987 time-series of carbon stable isotope ratios from dated sedimentary total organic carbon from the center of the Santa Barbara basin is compared with historical climate and oceanographic records. Carbon derived from carbon-13-depleted phytoplankton and carbon-13-enriched kelp appear responsible for a large part of the isotopic variance in sedimentary total organic carbon. El Niño/Southern Oscillation events are recorded by the isotopic response of marine organic carbon in sediments.

The effect of an upstream turbine on a low-aspect ratio vane

The interaction between a high-pressure rotor and a downstream vane is dominated by vortex-blade interaction. Each rotor blade passing period two co-rotating vortex pairs, the tip-leakage and upper passage vortex and the lower passage and trailing shed vortex, impinge on, and are cut by, the vane leading edge. In addition to the streamwise vortex the tip-leakage flow also contains a large velocity deficit. This causes the interaction of the tip-leakage flow with a downstream vane to differ from typical vortex blade interaction. This paper investigates the effect these interaction mechanisms have on a downstream vane. The test geometry considered was a low aspect ratio second stage vane located within a S-shaped diffuser with large radius change mounted downstream of a shroudless high-pressure turbine stage. Experimental measurements were conducted at engine-representative Mach and Reynolds numbers, and data was acquired using a fast-response aerodynamic probe upstream and downstream of the vane. Time-resolved numerical simulations were undertaken with and without a rotor tip gap in order to investigate the relative magnitude of the interaction mechanisms. The presence of the upstream stage is shown to significantly change the structure of the secondary flow in the vane and to cause a small drop in its performance.

Sperm cryopreservation of Indian major carp, Labeo rohita: cryodiluents, sperm: cryodiluent dilution ratio and cryoprotectan t concentration

Cryogenic preservation trials of spermatozoa of Labeo rohita were carried out. Twenty four cryodiluents (extender + cryoprotectant), with the combination of six extenders such as egg-yolk citrate, urea-egg-yolk, 0.9% NaCl, Kurokura-2, Ma and Mb and four cryoprotectants viz. DMSO, glycerol, methanol and ethanol, were used to screen out the suitable cryodiluents. Sperm was preserved in 0.25ml plastic straw in programmable freezer. Two step freezing method was followed. Sperm preserved with egg-yolk citrate and urea-egg-yolk containing 10% DMSO showed best post-thaw motility (80%) followed by 0.9% NaCl (60%) and Kurokura-2(30%) solutions. Sperm with the extenders M" and Mb clotted at the time of equilibration and also after few days of preservation. Egg-yolk citrate mixed with ethanol and methanol also showed good percentage of motility (80%) but egg-yolk citrate with glycerol showed less sperm motility (>60%). To determine suitable dilution ratio of milt and cryodiluent two best extender eggyolk citrate and urea-egg-yolk with four cryoprotectants such as DMSO, glycerol, methanol and ethanol at different ratio viz 1:2,1:4,1:7,1:10,1:15 and 1:20 were used. Highest post-thaw motility (>80%) was observed when milt was preserved with egg-yolk citrate containing 10% DMSO at 1:2, 1:4, 1:7 and 1:10 dilutions. Meanwhile using glycerol as cryoprotectants provided less post thaw motility at lower dilution ratio but with the increase of its dilution showed good sperm motility compared with other cryoprotectants. Finally, evaluation on the effect of cryoprotectant concentration on post-thaw sperm motility was conducted. Egg-yolk citrate and four cryoprotectant i.e. DMSO, glycerol, methanol and ethanol with six different concentrations namely 5%,7%, 10%, 15%, 20% and 30%.were evaluated. Among the cryoprotectants DMSO, methanol and ethanol showed highest post-thaw motility (about 80%) at 7% and 10% concentrations. Although glycerol was not suitable at low concentration but its 20% and 30% concentration levels provided best post-thaw motility. No post-thaw motility was obtained with DMSO at 30% concentration. The overall analysis on cryoprotectant concentration indicated that below 5% and above 20% cryoprotectant concentrations could not be suitable for effective cryopreservation of spermatozoa.