Two-dimensional approximation for tidal dynamics in coastal aquifers: Capillarity correction

Most previous investigations on tide-induced watertable fluctuations in coastal aquifers have been based on one-dimensional models that describe the processes in the cross-shore direction alone, assuming negligible along-shore variability. A recent study proposed a two-dimensional approximation for tide-induced watertable fluctuations that took into account coastline variations. Here, we further develop this approximation in two ways, by extending the approximation to second order and by taking into account capillary effects. Our results demonstrate that both effects can markedly influence watertable fluctuations. In particular, with the first-order approximation, the local damping rate of the tidal signal could be subject to sizable errors.

Physiologic parameters that affect pulse transit time difference between the upper and lower limbs in children

Pulse wave velocity (PWV) is a known parameter that is related to arterial distensibility. However, its potential is hampered by the absence of appropriate techniques to estimate it noninvasively. PWV can be used as an assessment of increased arterial stiffness that is linked to systolic hypertension, excess cardiovascular morbidity and mortality.(1,2)

Discrete and continuous models for dry masonry columns

The dynamic response of dry masonry columns can be approximated with finite-difference equations. Continuum models follow by replacing the difference quotients of the discrete model by corresponding differential expressions. The mathematically simplest of these models is a one-dimensional Cosserat theory. Within the presented homogenization context, the Cosserat theory is obtained by making ad hoc assumptions regarding the relative importance of certain terms in the differential expansions. The quality of approximation of the various theories is tested by comparison of the dispersion relations for bending waves with the dispersion relation of the discrete theory. All theories coincide with differences of less than 1% for wave-length-block-height (L/h) ratios bigger than 2 pi. The theory based on systematic differential approximation remains accurate up to L/h = 3 and then diverges rapidly. The Cosserat model becomes increasingly inaccurate for L/h < 2 pi. However, in contrast to the systematic approximation, the wave speed remains finite. In conclusion, considering its relative simplicity, the Cosserat model appears to be the natural starting point for the development of continuum models for blocky structures.

Robust algorithms for solving stochastic partial differential equations

A robust semi-implicit central partial difference algorithm for the numerical solution of coupled stochastic parabolic partial differential equations (PDEs) is described. This can be used for calculating correlation functions of systems of interacting stochastic fields. Such field equations can arise in the description of Hamiltonian and open systems in the physics of nonlinear processes, and may include multiplicative noise sources. The algorithm can be used for studying the properties of nonlinear quantum or classical field theories. The general approach is outlined and applied to a specific example, namely the quantum statistical fluctuations of ultra-short optical pulses in chi((2)) parametric waveguides. This example uses a non-diagonal coherent state representation, and correctly predicts the sub-shot noise level spectral fluctuations observed in homodyne detection measurements. It is expected that the methods used wilt be applicable for higher-order correlation functions and other physical problems as well. A stochastic differencing technique for reducing sampling errors is also introduced. This involves solving nonlinear stochastic parabolic PDEs in combination with a reference process, which uses the Wigner representation in the example presented here. A computer implementation on MIMD parallel architectures is discussed. (C) 1997 Academic Press.

Encephalitogenicity of murine, but not bovine, DM20 in SJL mice is due to a single amino acid difference in the immunodominant encephalitogenic epitope

Myelin proteolipid protein (PLP) contains 2 immunodominant encephalitogenic epitopes in SJL mice, namely PLP residues 139-151 and 178-191. DM20, a minor isoform of PLP, lacks residues 116-150 and consequently contains only the single major encephalitogenic epitope 178-191. However, it has been found previously that bovine DM20 is not encephalitogenic in SJL mice. Since residue 188 within peptide 178-191 is phenylalanine (F) in murine DM20 and alanine (A) in bovine DM20, we tested the effect of this difference on the immune responses and induction of EAE. SJL mice were immunized with either highly purified murine or bovine DM20. Residues 178-191 were found to be immunodominant for each, but only murine and not bovine DM20 was encephalitogenic. A synthetic peptide corresponding to the murine 178-191 sequence (F188) was also encephalitogenic, whereas the peptide corresponding to the bovine sequence (A188) was not. Both F188 and A188 bind with high affinity to I-A(s) and both are recognized by the SJL T cell repertoire. A188-specific T cell lines reacted to both A188 and F188, but F188-specific T cell lines were not stimulated by A188. F188-specific T cell lines produced mRNA for the Th1 cytokines IL2 and IFN gamma and, in passive transfer experiments, were encephalitogenic upon stimulation with F188, but not A188. In contrast, A188-specific T cell lines produced mRNA for IL4, IL5 and IL10, in addition to IL2 and IFN gamma, and were not encephalitogenic after stimulation with either F188 or A188. Cotransfer of A188-specific T cell lines with F188-specific T cell lines resulted in protection from EAE. Thus, A188 induces a functionally different phenotype of T cells from that induced by F188. Taken together these data suggest that the failure of bovine DM20 to induce EAE may be attributable to induction of protective rather than pathogenic T cells by the immunodominant epitope.

From EEG to BOLD: Brain mapping and estimating transfer functions in simultaneous EEG-fMRI acquisitions

Simultaneous acquisition of electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) aims to disentangle the description of brain processes by exploiting the advantages of each technique. Most studies in this field focus on exploring the relationships between fMRI signals and the power spectrum at some specific frequency bands (alpha, beta, etc.). On the other hand, brain mapping of EEG signals (e.g., interictal spikes in epileptic patients) usually assumes an haemodynamic response function for a parametric analysis applying the GLM, as a rough approximation. The integration of the information provided by the high spatial resolution of MR images and the high temporal resolution of EEG may be improved by referencing them by transfer functions, which allows the identification of neural driven areas without strong assumptions about haemodynamic response shapes or brain haemodynamic`s homogeneity. The difference on sampling rate is the first obstacle for a full integration of EEG and fMRI information. Moreover, a parametric specification of a function representing the commonalities of both signals is not established. In this study, we introduce a new data-driven method for estimating the transfer function from EEG signal to fMRI signal at EEG sampling rate. This approach avoids EEG subsampling to fMRI time resolution and naturally provides a test for EEG predictive power over BOLD signal fluctuations, in a well-established statistical framework. We illustrate this concept in resting state (eyes closed) and visual simultaneous fMRI-EEG experiments. The results point out that it is possible to predict the BOLD fluctuations in occipital cortex by using EEG measurements. (C) 2010 Elsevier Inc. All rights reserved.

On the validity of the dispersion model of hepatic drug elimination when intravascular transit time densities are long-tailed

The dispersion model with mixed boundary conditions uses a single parameter, the dispersion number, to describe the hepatic elimination of xenobiotics and endogenous substances. An implicit a priori assumption of the model is that the transit time density of intravascular indicators is approximated by an inverse Gaussian distribution. This approximation is limited in that the model poorly describes the tail part of the hepatic outflow curves of vascular indicators. A sum of two inverse Gaussian functions is proposed as ail alternative, more flexible empirical model for transit time densities of vascular references. This model suggests that a more accurate description of the tail portion of vascular reference curves yields an elimination rate constant (or intrinsic clearance) which is 40% less than predicted by the dispersion model with mixed boundary conditions. The results emphasize the need to accurately describe outflow curves in using them as a basis for determining pharmacokinetic parameters using hepatic elimination models. (C) 1997 Society for Mathematical Biology.

The subcycled Newmark algorithm

The popular Newmark algorithm, used for implicit direct integration of structural dynamics, is extended by means of a nodal partition to permit use of different timesteps in different regions of a structural model. The algorithm developed has as a special case an explicit-explicit subcycling algorithm previously reported by Belytschko, Yen and Mullen. That algorithm has been shown, in the absence of damping or other energy dissipation, to exhibit instability over narrow timestep ranges that become narrower as the number of degrees of freedom increases, making them unlikely to be encountered in practice. The present algorithm avoids such instabilities in the case of a one to two timestep ratio (two subcycles), achieving unconditional stability in an exponential sense for a linear problem. However, with three or more subcycles, the trapezoidal rule exhibits stability that becomes conditional, falling towards that of the central difference method as the number of subcycles increases. Instabilities over narrow timestep ranges, that become narrower as the model size increases, also appear with three or more subcycles. However by moving the partition between timesteps one row of elements into the region suitable for integration with the larger timestep these the unstable timestep ranges become extremely narrow, even in simple systems with a few degrees of freedom. As well, accuracy is improved. Use of a version of the Newmark algorithm that dissipates high frequencies minimises or eliminates these narrow bands of instability. Viscous damping is also shown to remove these instabilities, at the expense of having more effect on the low frequency response.

The use of sodium-chloride difference and chloride-sodium ratio as strong ion difference surrogates in the evaluation of metabolic acidosis in critically ill patients

Purpose: Inorganic apparent strong ion difference (SIDai) improves chloride-associated acidosis recognition in dysnatremic patients. We investigated whether the difference between sodium and chloride (Na+-C1-) or the ratio between chloride and sodium (Cl-/Na+) could be used as SIDai surrogates in mixed and dysnatremic patients. Patients and Methods: Two arterial blood samples were collected from 128 patients. Physicochemical analytical approach was used. Correlation, agreement, accuracy, sensitivity, and specificity were measured to examine whether Na(+)-C1(-) and CI(-)/Na(+) could be used instead of SIDai in the diagnosis of acidosis. Results: Na(+)-C1(-) and CF/Na+ were well correlated with SIDai (R = 0.987, P < 0.001 and R = 0.959, P < 0.001, respectively). Bias between Na(+)-C1(-) and SIDai was high (6.384 with a limit of agreement of 4.4638.305 mEq/L). Accuracy values for the identification of SIDai acidosis (<38.9 mEq/L) were 0.989 (95% confidence interval [CI], 0.980-0.998) for Na+-C1- and 0.974 (95% CI, 0.959-0.989) for Cr/Na+. Receiver operator characteristic curve showed that values revealing SIDai acidosis were less than 32.5 mEq/L for Nata- and more than 0.764 for C17Na+ with sensitivities of 94.0% and 92.0% and specificities of 97.0% and 90.0%, respectively. Nata- was a reliable S IDai surrogate in dysnatremic patients. Conclusions: Nata- and CI-/Na+ are good tools to disclose S IDai acidosis. In patients with dysnatremia, Nata- is an accurate tool to diagnose SIDai acidosis. (C) 2010 Elsevier Inc. All rights reserved.

Comparison of renal transplantation outcomes in children with and without bladder dysfunction. A customized approach equals the difference

Purpose: We examined the development of urological abnormalities in a group of pediatric renal transplant recipients. Materials and Methods: We reviewed 211 patients younger than 19 years who underwent 226 renal transplants. Three groups of patients were studied-136 children with end stage renal disease due to a nonurological cause (group 1), 56 children with a urological disorder but with an adequate bladder (group 2a) and 19 children with lower urinary tract dysfunction and/or inadequate bladder drainage (group 2b). A total of 15 children in group 2b underwent bladder augmentation (ureterocystoplasty in 6, enterocystoplasty in 9), 2 underwent continent urinary diversion, 1 underwent autoaugmentation and 1 underwent a Mitrofanoff procedure at the bladder for easier drainage. Kidney transplantation was performed in the classic manner by extraperitoneal access, and whenever possible the ureter was reimplanted using an antireflux procedure. Results: At a mean followup of 75 months 13 children had died, 59 grafts were lost and 15 children had received a second transplant. Two patients in group 2a required a complementary urological procedure to preserve renal function (1 enterocystoplasty, 1 vesicostomy). A total of 12 major surgical complications occurred in 226 kidney transplants (5.3%), with a similar incidence in all groups. The overall graft survival at 5 years was 75%, 74% and 84%, respectively, in groups 1, 2a and 2b. Conclusions: With individualized treatment children with severely inferior lower urinary tract function may undergo renal transplantation with a safe and adequate outcome.

Representational Difference Analysis (RDA) reveals differential expression of conserved as well as novel genes during caste-specific development of the honey bee (Apis mellifera L.) ovary

In highly eusocial insects, such as the honey bee, Apis mellifera, the reproductive bias has become embedded in morphological caste differences. These are most expressively denoted in ovary size, with adult queens having large ovaries consisting of 150-200 ovarioles each, while workers typically have only 1-20 ovarioles per ovary. This morphological differentiation is a result of hormonal signals triggered by the diet change in the third larval instar, which eventually generate caste-specific gene expression patterns. To reveal these we produced differential gene expression libraries by Representational Difference Analysis (RDA) for queen and worker ovaries in a developmental stage when cell death is a prominent feature in the ovarioles of workers, whereas all ovarioles are maintained and extend in length in queens. In the queen library, 48% of the gene set represented homologs of known Drosophila genes, whereas in the worker ovary, the largest set (59%) were ESTs evidencing novel genes, not even computationally predicted in the honey bee genome. Differential expression was confirmed by quantitative RT-PCR for a selected gene set, denoting major differences for two queen and two worker library genes. These included two unpredicted genes located in chromosome 11 (Group11.35 and Group11.31, respectively) possibly representing long non-coding RNAs. Being candidates as modulators of ovary development, their expression and functional analysis should be a focal point for future studies. (C) 2011 Elsevier Ltd. All rights reserved.