Oscillation modes of microtubules.

Microtubules are long, filamentous protein complexes which play a central role in several cellular physiological processes, such as cell division transport and locomotion. Their mechanical properties are extremely important since they determine the biological function. In a recently published experiment [Phys. Rev. Lett. 89 (2002) 248101], microtubule's Young's and shear moduli were simultaneously measured, proving that they are highly anisotropic. Together with the known structure, this finding opens the way to better understand and predict their mechanical behavior under a particular set of conditions. In the present study, we modeled microtubules by using the finite elements method and analyzed their oscillation modes. The analysis revealed that oscillation modes involving a change in the diameter of the microtubules strongly depend on the shear modulus. In these modes, the correlation times of the movements are just slightly shorter than diffusion times of free molecules surrounding the microtubule. It could be therefore speculated that the matching of the two timescales could play a role in facilitating the interactions between microtubules and MT associated proteins, and between microtubules and tubulins themselves.

Orthotropic active strain models for the numerical simulation of cardiac biomechanics

An active strain formulation for orthotropic constitutive laws arising in cardiac mechanics modeling is introduced and studied. The passive mechanical properties of the tissue are described by the Holzapfel-Ogden relation. In the active strain formulation, the Euler-Lagrange equations for minimizing the total energy are written in terms of active and passive deformation factors, where the active part is assumed to depend, at the cell level, on the electrodynamics and on the specific orientation of the cardiac cells. The well-posedness of the linear system derived from a generic Newton iteration of the original problem is analyzed and different mechanical activation functions are considered. In addition, the active strain formulation is compared with the classical active stress formulation from both numerical and modeling perspectives. Taylor-Hood and MINI finite elements are employed to discretize the mechanical problem. The results of several numerical experiments show that the proposed formulation is mathematically consistent and is able to represent the main key features of the phenomenon, while allowing savings in computational costs.

Significance of serum adiponectin levels in patients with chronic liver disease.

Adiponectin, which plays a pivotal role in metabolic liver diseases, is reduced in concentration in patients with NASH (non-alcoholic steatohepatitis). The aim of the present study was to determine adiponectin concentrations in patients with different forms and stages of chronic liver diseases. Serum adiponectin concentrations were measured in 232 fasting patients with chronic liver disease: 64 with NAFLD (non-alcoholic fatty liver disease), 123 with other chronic liver disease (e.g. viral hepatitis, n=71; autoimmune disease, n=18; alcohol-induced liver disease, n=3; or elevated liver enzymes of unknown origin, n=31) and 45 with cirrhosis. Circulating adiponectin levels were significantly lower in patients with NAFLD in comparison with patients with other chronic liver disease (4.8+/-3.5 compared with 10.4+/-6.3 microg/ml respectively; P<0.0001). Circulating adiponectin levels were significantly higher in patients with cirrhosis in comparison with patients without cirrhosis (18.6+/-14.5 compared with 8.4+/-6.1 microg/ml respectively; P<0.0001). Adiponectin concentrations correlated negatively with body weight (P<0.001), serum triacylglycerols (triglycerides) (P<0.001) and, in women, with BMI (body mass index) (P<0.001). Adiponectin concentrations correlated positively with serum bile acids (P<0.001), serum hyaluronic acid (P<0.001) and elastography values (P<0.001). Adiponectin levels were decreased in patients with NAFLD. In conclusion, adiponectin levels correlate positively with surrogate markers of hepatic fibrosis (transient elastography, fasting serum bile acids and hyaluronate) and are significantly elevated in cases of cirrhosis.

Comparative Study of Chemoembolization, Loadable Beads: In vitro Drug Release and Physical Properties of DC Bead and Hepasphere, Loaded with Doxorubicin and Irinotecan

Purpose To characterize in vitro the loadability, physical properties, and release of irinotecan and doxorubicin from two commercially available embolization microspheres. Materials and Methods DC Bead (500-700 μm) and Hepasphere (400-600 μm) microspheres were loaded with either doxorubicin or irinotecan solutions. Drug amount was quantified with spectrophotometry, bead elasticity was measured under compression, and bead size and loading homogeneity were assessed with microscopy image analysis. Drug release was measured over 1-week periods in saline by using a pharmacopeia flow-through method. Results Almost complete drug loading was obtained for both microsphere types and drugs. Doxorubicin-loaded DC Beads maintained their spherical shape throughout the release. In contrast, Hepaspheres showed less homogeneous doxorubicin loading and, after release, some fractured microspheres. Incomplete doxorubicin release was observed in saline over 1 week (27% ± 2 for DC beads and 18% ± 7 for Hepaspheres; P = .013). About 75% of this amount was released within 2.2 hours for both beads. For irinotecan, complete release was obtained for both types of beads, in a sustained manner over 2-3 hours for DC Beads, and in a significantly faster manner as a 7-minute burst for Hepaspheres. Conclusions The two drug-eluting microspheres could be efficiently loaded with both drugs. Incomplete doxorubicin release was attributed to strong drug-bead ionic interactions. Weaker interactions were observed with irinotecan, which led to faster drug release.

Spatially-resolved eigenmode decomposition of red blood cells membrane fluctuations questions the role of ATP in flickering.

Red blood cells (RBCs) present unique reversible shape deformability, essential for both function and survival, resulting notably in cell membrane fluctuations (CMF). These CMF have been subject of many studies in order to obtain a better understanding of these remarkable biomechanical membrane properties altered in some pathological states including blood diseases. In particular the discussion over the thermal or metabolic origin of the CMF has led in the past to a large number of investigations and modeling. However, the origin of the CMF is still debated. In this article, we present an analysis of the CMF of RBCs by combining digital holographic microscopy (DHM) with an orthogonal subspace decomposition of the imaging data. These subspace components can be reliably identified and quantified as the eigenmode basis of CMF that minimizes the deformation energy of the RBC structure. By fitting the observed fluctuation modes with a theoretical dynamic model, we find that the CMF are mainly governed by the bending elasticity of the membrane and that shear and tension elasticities have only a marginal influence on the membrane fluctations of the discocyte RBC. Further, our experiments show that the role of ATP as a driving force of CMF is questionable. ATP, however, seems to be required to maintain the unique biomechanical properties of the RBC membrane that lead to thermally excited CMF.

Rules of origin : from analysis to reform

General Introduction This thesis can be divided into two main parts :the first one, corresponding to the first three chapters, studies Rules of Origin (RoOs) in Preferential Trade Agreements (PTAs); the second part -the fourth chapter- is concerned with Anti-Dumping (AD) measures. Despite wide-ranging preferential access granted to developing countries by industrial ones under North-South Trade Agreements -whether reciprocal, like the Europe Agreements (EAs) or NAFTA, or not, such as the GSP, AGOA, or EBA-, it has been claimed that the benefits from improved market access keep falling short of the full potential benefits. RoOs are largely regarded as a primary cause of the under-utilization of improved market access of PTAs. RoOs are the rules that determine the eligibility of goods to preferential treatment. Their economic justification is to prevent trade deflection, i.e. to prevent non-preferred exporters from using the tariff preferences. However, they are complex, cost raising and cumbersome, and can be manipulated by organised special interest groups. As a result, RoOs can restrain trade beyond what it is needed to prevent trade deflection and hence restrict market access in a statistically significant and quantitatively large proportion. Part l In order to further our understanding of the effects of RoOs in PTAs, the first chapter, written with Pr. Olivier Cadot, Celine Carrère and Pr. Jaime de Melo, describes and evaluates the RoOs governing EU and US PTAs. It draws on utilization-rate data for Mexican exports to the US in 2001 and on similar data for ACP exports to the EU in 2002. The paper makes two contributions. First, we construct an R-index of restrictiveness of RoOs along the lines first proposed by Estevadeordal (2000) for NAFTA, modifying it and extending it for the EU's single-list (SL). This synthetic R-index is then used to compare Roos under NAFTA and PANEURO. The two main findings of the chapter are as follows. First, it shows, in the case of PANEURO, that the R-index is useful to summarize how countries are differently affected by the same set of RoOs because of their different export baskets to the EU. Second, it is shown that the Rindex is a relatively reliable statistic in the sense that, subject to caveats, after controlling for the extent of tariff preference at the tariff-line level, it accounts for differences in utilization rates at the tariff line level. Finally, together with utilization rates, the index can be used to estimate total compliance costs of RoOs. The second chapter proposes a reform of preferential Roos with the aim of making them more transparent and less discriminatory. Such a reform would make preferential blocs more "cross-compatible" and would therefore facilitate cumulation. It would also contribute to move regionalism toward more openness and hence to make it more compatible with the multilateral trading system. It focuses on NAFTA, one of the most restrictive FTAs (see Estevadeordal and Suominen 2006), and proposes a way forward that is close in spirit to what the EU Commission is considering for the PANEURO system. In a nutshell, the idea is to replace the current array of RoOs by a single instrument- Maximum Foreign Content (MFC). An MFC is a conceptually clear and transparent instrument, like a tariff. Therefore changing all instruments into an MFC would bring improved transparency pretty much like the "tariffication" of NTBs. The methodology for this exercise is as follows: In step 1, I estimate the relationship between utilization rates, tariff preferences and RoOs. In step 2, I retrieve the estimates and invert the relationship to get a simulated MFC that gives, line by line, the same utilization rate as the old array of Roos. In step 3, I calculate the trade-weighted average of the simulated MFC across all lines to get an overall equivalent of the current system and explore the possibility of setting this unique instrument at a uniform rate across lines. This would have two advantages. First, like a uniform tariff, a uniform MFC would make it difficult for lobbies to manipulate the instrument at the margin. This argument is standard in the political-economy literature and has been used time and again in support of reductions in the variance of tariffs (together with standard welfare considerations). Second, uniformity across lines is the only way to eliminate the indirect source of discrimination alluded to earlier. Only if two countries face uniform RoOs and tariff preference will they face uniform incentives irrespective of their initial export structure. The result of this exercise is striking: the average simulated MFC is 25% of good value, a very low (i.e. restrictive) level, confirming Estevadeordal and Suominen's critical assessment of NAFTA's RoOs. Adopting a uniform MFC would imply a relaxation from the benchmark level for sectors like chemicals or textiles & apparel, and a stiffening for wood products, papers and base metals. Overall, however, the changes are not drastic, suggesting perhaps only moderate resistance to change from special interests. The third chapter of the thesis considers whether Europe Agreements of the EU, with the current sets of RoOs, could be the potential model for future EU-centered PTAs. First, I have studied and coded at the six-digit level of the Harmonised System (HS) .both the old RoOs -used before 1997- and the "Single list" Roos -used since 1997. Second, using a Constant Elasticity Transformation function where CEEC exporters smoothly mix sales between the EU and the rest of the world by comparing producer prices on each market, I have estimated the trade effects of the EU RoOs. The estimates suggest that much of the market access conferred by the EAs -outside sensitive sectors- was undone by the cost-raising effects of RoOs. The chapter also contains an analysis of the evolution of the CEECs' trade with the EU from post-communism to accession. Part II The last chapter of the thesis is concerned with anti-dumping, another trade-policy instrument having the effect of reducing market access. In 1995, the Uruguay Round introduced in the Anti-Dumping Agreement (ADA) a mandatory "sunset-review" clause (Article 11.3 ADA) under which anti-dumping measures should be reviewed no later than five years from their imposition and terminated unless there was a serious risk of resumption of injurious dumping. The last chapter, written with Pr. Olivier Cadot and Pr. Jaime de Melo, uses a new database on Anti-Dumping (AD) measures worldwide to assess whether the sunset-review agreement had any effect. The question we address is whether the WTO Agreement succeeded in imposing the discipline of a five-year cycle on AD measures and, ultimately, in curbing their length. Two methods are used; count data analysis and survival analysis. First, using Poisson and Negative Binomial regressions, the count of AD measures' revocations is regressed on (inter alia) the count of "initiations" lagged five years. The analysis yields a coefficient on measures' initiations lagged five years that is larger and more precisely estimated after the agreement than before, suggesting some effect. However the coefficient estimate is nowhere near the value that would give a one-for-one relationship between initiations and revocations after five years. We also find that (i) if the agreement affected EU AD practices, the effect went the wrong way, the five-year cycle being quantitatively weaker after the agreement than before; (ii) the agreement had no visible effect on the United States except for aone-time peak in 2000, suggesting a mopping-up of old cases. Second, the survival analysis of AD measures around the world suggests a shortening of their expected lifetime after the agreement, and this shortening effect (a downward shift in the survival function postagreement) was larger and more significant for measures targeted at WTO members than for those targeted at non-members (for which WTO disciplines do not bind), suggesting that compliance was de jure. A difference-in-differences Cox regression confirms this diagnosis: controlling for the countries imposing the measures, for the investigated countries and for the products' sector, we find a larger increase in the hazard rate of AD measures covered by the Agreement than for other measures.

Effects of the transition time between muscle-tendon stretch and shortening on mechanical efficiency.

The net mechanical efficiency of positive work (eta(pos)) has been shown to increase if it is immediately preceded by negative work. This phenomenon is explained by the storage of elastic energy during the negative phase and its release during the subsequent positive phase. If a transition time (T) takes place, the elastic energy is dissipated into heat. The aim of the present study was to investigate the relationship between eta(pos) and T, and to determine the minimal T required so that eta(pos) reached its minimal value. Seven healthy male subjects were tested during four series of lowering-raising of the body mass. In the first series (S (0)), the negative and positive phases were executed without any transition time. In the three other series, T was varied by a timer (0.12, 0.24 and 0.56 s for series S (1), S (2) and S (3), respectively). These exercises were performed on a force platform sensitive to vertical forces to measure the mechanical work and a gas analyser was used to determine the energy expenditure. The results indicated that eta(pos) was the highest (31.1%) for the series without any transition time (S (0)). The efficiencies observed with transition times (S (1), S (2) and S (3)) were 27.7, 26.0 and 23.8%, respectively, demonstrating that T plays an important role for mechanical efficiency. The investigation of the relationship between eta(pos) and T revealed that the minimal T required so that eta(pos) reached its minimal value is 0.59 s.

Chronic rejection of allogenic dermal grafts in 73 % TBSA burn patient

Rationale: Allogenic grafts are an excellent way to temporarily cover a wound. It prevents the loss of electrolytes and water, reduces the risk of infection and diminishes pain. Another advantage of the allograft is in circumventing problems such as the morbidity of skin graft donor sites. We present here the case of a patient grafted in 1991 with cultured epidermal autografts (CEA) and allogenic skin transplants on his legs, outlining the risks and potential long-term complications. Methods: The 40-year-old male patient was treated with allogenic Split Thickness Skin Graft (STSG) transplantations, CEA and Cyclosporine-A therapy. Allogenic STSG for lower extremities were harvested from a female HIV-negative organ donor. They were transplanted, de-epithelialized and subsequently covered with CEAs. Cyclosporine-A was administered systemically from the first day following transplantation until three weeks after the last CEAs were placed on the allogenic dermis. Results: Immediate results showed a 90% successful grafting under cyclosporine therapy. However, some lesions were still present 16 months later. The skin was hard with little or no elasticity. Five years after the transplantation there were no more lesions. However, a 10-year follow-up showed new ulcers on both lower extremities. All the skin of the right leg was removed and replaced by STSG from the patient's back. Postoperative results were excellent with a 100% graft take. The anatomopathology showed dermo-hypodermic tissue with fibrosis of the dermis, vasculopathy and chronic ulcers compatible with chronic rejection. Conclusion: While early functional results of the allografts may seem encouraging, their long-term evolution remains uncertain and, in this case, presents complications. The apparent antigenic effect of the dermal tissue may be controlled with long-term immunosuppression which may cause important secondary effects. Even with such treatments, 15 years after organ transplantation, about 35% of a transplant is no longer functional. It is therefore important to take these long-term observations into consideration when treating sensitive areas such as hands or a face.

Measurements of seismic attenuation and transient fluid pressure in partially saturated Berea sandstone: Evidence of fluid flow on the mesoscopic scale

A novel laboratory technique is proposed to investigate wave-induced fluid flow on the mesoscopic scale as a mechanism for seismic attenuation in partially saturated rocks. This technique combines measurements of seismic attenuation in the frequency range from 1 to 100?Hz with measurements of transient fluid pressure as a response of a step stress applied on top of the sample. We used a Berea sandstone sample partially saturated with water. The laboratory results suggest that wave-induced fluid flow on the mesoscopic scale is dominant in partially saturated samples. A 3-D numerical model representing the sample was used to verify the experimental results. Biot's equations of consolidation were solved with the finite-element method. Wave-induced fluid flow on the mesoscopic scale was the only attenuation mechanism accounted for in the numerical solution. The numerically calculated transient fluid pressure reproduced the laboratory data. Moreover, the numerically calculated attenuation, superposed to the frequency-independent matrix anelasticity, reproduced the attenuation measured in the laboratory in the partially saturated sample. This experimental?numerical fit demonstrates that wave-induced fluid flow on the mesoscopic scale and matrix anelasticity are the dominant mechanisms for seismic attenuation in partially saturated Berea sandstone.

New formulation of vasoactive intestinal peptide using liposomes in hyaluronic acid gel for uveitis.

We evaluated the benefits of a novel formulation of vasoactive intestinal peptide (VIP) based on the incorporation of VIP-loaded rhodamine-conjugated liposomes (VIP-Rh-Lip) within hyaluronic acid (HA) gel (Gel-VIP-Rh-Lip) for the treatment of endotoxin-induced uveitis (EIU) in comparison with VIP-Rh-Lip alone. In vitro release study and rheological analysis showed that interactions between HA chains and liposomes resulted in increased viscosity and reinforced elasticity of the gel. In vivo a single intravitreal injection of Gel-VIP-Rh-Lip was performed in rats 7 days prior to uveitis induction by subcutaneous lipopolysaccharide injection. The maximal ocular inflammation occurs within 16-24 h in controls (VIP-Rh-Lip, unloaded-Rh-Lip). Whereas intraocular injection of VIP-Rh-Lip had no effect on EIU severity compared with controls, Gel-VIP-Rh-Lip reduced significantly the clinical score and number of inflammatory cells infiltrating the eye. The fate of liposomes, VIP and HA in the eyes, regional and inguinal lymph nodes and spleen was analyzed by immunostaining and fluorescence microscopy. Retention of liposomes by HA gel was observed in vitro and in vivo. Inflammation severity seemed to impact on system stability resulting in the delayed release of VIP. Thus, HA gel containing VIP-Rh-Lip is an efficient strategy to obtain a sustained delivery of VIP in ocular and lymph node tissues.

Looking for orders of magnitude : five essays in trade and the environment and economic geography

General Summary Although the chapters of this thesis address a variety of issues, the principal aim is common: test economic ideas in an international economic context. The intention has been to supply empirical findings using the largest suitable data sets and making use of the most appropriate empirical techniques. This thesis can roughly be divided into two parts: the first one, corresponding to the first two chapters, investigates the link between trade and the environment, the second one, the last three chapters, is related to economic geography issues. Environmental problems are omnipresent in the daily press nowadays and one of the arguments put forward is that globalisation causes severe environmental problems through the reallocation of investments and production to countries with less stringent environmental regulations. A measure of the amplitude of this undesirable effect is provided in the first part. The third and the fourth chapters explore the productivity effects of agglomeration. The computed spillover effects between different sectors indicate how cluster-formation might be productivity enhancing. The last chapter is not about how to better understand the world but how to measure it and it was just a great pleasure to work on it. "The Economist" writes every week about the impressive population and economic growth observed in China and India, and everybody agrees that the world's center of gravity has shifted. But by how much and how fast did it shift? An answer is given in the last part, which proposes a global measure for the location of world production and allows to visualize our results in Google Earth. A short summary of each of the five chapters is provided below. The first chapter, entitled "Unraveling the World-Wide Pollution-Haven Effect" investigates the relative strength of the pollution haven effect (PH, comparative advantage in dirty products due to differences in environmental regulation) and the factor endowment effect (FE, comparative advantage in dirty, capital intensive products due to differences in endowments). We compute the pollution content of imports using the IPPS coefficients (for three pollutants, namely biological oxygen demand, sulphur dioxide and toxic pollution intensity for all manufacturing sectors) provided by the World Bank and use a gravity-type framework to isolate the two above mentioned effects. Our study covers 48 countries that can be classified into 29 Southern and 19 Northern countries and uses the lead content of gasoline as proxy for environmental stringency. For North-South trade we find significant PH and FE effects going in the expected, opposite directions and being of similar magnitude. However, when looking at world trade, the effects become very small because of the high North-North trade share, where we have no a priori expectations about the signs of these effects. Therefore popular fears about the trade effects of differences in environmental regulations might by exaggerated. The second chapter is entitled "Is trade bad for the Environment? Decomposing worldwide SO2 emissions, 1990-2000". First we construct a novel and large database containing reasonable estimates of SO2 emission intensities per unit labor that vary across countries, periods and manufacturing sectors. Then we use these original data (covering 31 developed and 31 developing countries) to decompose the worldwide SO2 emissions into the three well known dynamic effects (scale, technique and composition effect). We find that the positive scale (+9,5%) and the negative technique (-12.5%) effect are the main driving forces of emission changes. Composition effects between countries and sectors are smaller, both negative and of similar magnitude (-3.5% each). Given that trade matters via the composition effects this means that trade reduces total emissions. We next construct, in a first experiment, a hypothetical world where no trade happens, i.e. each country produces its imports at home and does no longer produce its exports. The difference between the actual and this no-trade world allows us (under the omission of price effects) to compute a static first-order trade effect. The latter now increases total world emissions because it allows, on average, dirty countries to specialize in dirty products. However, this effect is smaller (3.5%) in 2000 than in 1990 (10%), in line with the negative dynamic composition effect identified in the previous exercise. We then propose a second experiment, comparing effective emissions with the maximum or minimum possible level of SO2 emissions. These hypothetical levels of emissions are obtained by reallocating labour accordingly across sectors within each country (under the country-employment and the world industry-production constraints). Using linear programming techniques, we show that emissions are reduced by 90% with respect to the worst case, but that they could still be reduced further by another 80% if emissions were to be minimized. The findings from this chapter go together with those from chapter one in the sense that trade-induced composition effect do not seem to be the main source of pollution, at least in the recent past. Going now to the economic geography part of this thesis, the third chapter, entitled "A Dynamic Model with Sectoral Agglomeration Effects" consists of a short note that derives the theoretical model estimated in the fourth chapter. The derivation is directly based on the multi-regional framework by Ciccone (2002) but extends it in order to include sectoral disaggregation and a temporal dimension. This allows us formally to write present productivity as a function of past productivity and other contemporaneous and past control variables. The fourth chapter entitled "Sectoral Agglomeration Effects in a Panel of European Regions" takes the final equation derived in chapter three to the data. We investigate the empirical link between density and labour productivity based on regional data (245 NUTS-2 regions over the period 1980-2003). Using dynamic panel techniques allows us to control for the possible endogeneity of density and for region specific effects. We find a positive long run elasticity of density with respect to labour productivity of about 13%. When using data at the sectoral level it seems that positive cross-sector and negative own-sector externalities are present in manufacturing while financial services display strong positive own-sector effects. The fifth and last chapter entitled "Is the World's Economic Center of Gravity Already in Asia?" computes the world economic, demographic and geographic center of gravity for 1975-2004 and compares them. Based on data for the largest cities in the world and using the physical concept of center of mass, we find that the world's economic center of gravity is still located in Europe, even though there is a clear shift towards Asia. To sum up, this thesis makes three main contributions. First, it provides new estimates of orders of magnitudes for the role of trade in the globalisation and environment debate. Second, it computes reliable and disaggregated elasticities for the effect of density on labour productivity in European regions. Third, it allows us, in a geometrically rigorous way, to track the path of the world's economic center of gravity.

Frequency-dependent attenuation as a potential indicator of oil saturation

At seismic frequencies, wave-induced fluid flow is a major cause of P-wave attenuation in partially saturated porous rocks. Attenuation is of great importance for the oil industry in the interpretation of seismic field data. Here, the effects on P-wave attenuation resulting from changes in oil saturation are studied for media with coexisting water, oil, and gas. For that, creep experiments are numerically simulated by solving Biot's equations for consolidation of poroelastic media with the finite-element method. The experiments yield time-dependent stress?strain relations that are used to calculate the complex P-wave modulus from which frequency-dependent P-wave attenuation is determined. The models are layered media with periodically alternating triplets of layers. Models consisting of triplets of layers having randomly varying layer thicknesses are also considered. The layers in each triplet are fully saturated with water, oil, and gas. The layer saturated with water has lower porosity and permeability than the layers saturated with oil and gas. These models represent hydrocarbon reservoirs in which water is the wetting fluid preferentially saturating regions of lower porosity. The results from the numerical experiments showed that increasing oil saturation, connected to a decrease in gas saturation, resulted in a significant increase of attenuation at low frequencies (lower than 2 Hz). Furthermore, replacing the oil with water resulted in a distinguishable behavior of the frequency-dependent attenuation. These results imply that, according to the physical mechanism of wave-induced fluid flow, frequency-dependent attenuation in media saturated with water, oil, and gas is a potential indicator of oil saturation.

Constraints on the permeability structure of alluvial aquifers from the poro-elastic inversion of multi-frequency P-wave velocity logs

We have explored the possibility of obtaining first-order permeability estimates for saturated alluvial sediments based on the poro-elastic interpretation of the P-wave velocity dispersion inferred from sonic logs. Modern sonic logging tools designed for environmental and engineering applications allow one for P-wave velocity measurements at multiple emitter frequencies over a bandwidth covering 5 to 10 octaves. Methodological considerations indicate that, for saturated unconsolidated sediments in the silt to sand range and typical emitter frequencies ranging from approximately 1 to 30 kHz, the observable velocity dispersion should be sufficiently pronounced to allow one for reliable first-order estimations of the permeability structure. The corresponding predictions have been tested on and verified for a borehole penetrating a typical surficial alluvial aquifer. In addition to multifrequency sonic logs, a comprehensive suite of nuclear and electrical logs, an S-wave log, a litholog, and a limited number laboratory measurements of the permeability from retrieved core material were also available. This complementary information was found to be essential for parameterizing the poro-elastic inversion procedure and for assessing the uncertainty and internal consistency of corresponding permeability estimates. Our results indicate that the thus obtained permeability estimates are largely consistent with those expected based on the corresponding granulometric characteristics, as well as with the available evidence form laboratory measurements. These findings are also consistent with evidence from ocean acoustics, which indicate that, over a frequency range of several orders-of-magnitude, the classical theory of poro-elasticity is generally capable of explaining the observed P-wave velocity dispersion in medium- to fine-grained seabed sediments

Direct three-dimensional myocardial strain tensor quantification and tracking using zHARP.

Images of myocardial strain can be used to diagnose heart disease, plan and monitor treatment, and to learn about cardiac structure and function. Three-dimensional (3D) strain is typically quantified using many magnetic resonance (MR) images obtained in two or three orthogonal planes. Problems with this approach include long scan times, image misregistration, and through-plane motion. This article presents a novel method for calculating cardiac 3D strain using a stack of two or more images acquired in only one orientation. The zHARP pulse sequence encodes in-plane motion using MR tagging and out-of-plane motion using phase encoding, and has been previously shown to be capable of computing 3D displacement within a single image plane. Here, data from two adjacent image planes are combined to yield a 3D strain tensor at each pixel; stacks of zHARP images can be used to derive stacked arrays of 3D strain tensors without imaging multiple orientations and without numerical interpolation. The performance and accuracy of the method is demonstrated in vitro on a phantom and in vivo in four healthy adult human subjects.

Simulation of surface waves in porous media

We present a novel numerical algorithm for the simulation of seismic wave propagation in porous media, which is particularly suitable for the accurate modelling of surface wave-type phenomena. The differential equations of motion are based on Biot's theory of poro-elasticity and solved with a pseudospectral approach using Fourier and Chebyshev methods to compute the spatial derivatives along the horizontal and vertical directions, respectively. The time solver is a splitting algorithm that accounts for the stiffness of the differential equations. Due to the Chebyshev operator the grid spacing in the vertical direction is non-uniform and characterized by a denser spatial sampling in the vicinity of interfaces, which allows for a numerically stable and accurate evaluation of higher order surface wave modes. We stretch the grid in the vertical direction to increase the minimum grid spacing and reduce the computational cost. The free-surface boundary conditions are implemented with a characteristics approach, where the characteristic variables are evaluated at zero viscosity. The same procedure is used to model seismic wave propagation at the interface between a fluid and porous medium. In this case, each medium is represented by a different grid and the two grids are combined through a domain-decomposition method. This wavefield decomposition method accounts for the discontinuity of variables and is crucial for an accurate interface treatment. We simulate seismic wave propagation with open-pore and sealed-pore boundary conditions and verify the validity and accuracy of the algorithm by comparing the numerical simulations to analytical solutions based on zero viscosity obtained with the Cagniard-de Hoop method. Finally, we illustrate the suitability of our algorithm for more complex models of porous media involving viscous pore fluids and strongly heterogeneous distributions of the elastic and hydraulic material properties.