Simulative Investigation on the Electronic, Vibrational and Optical Properties of the Ge2Sb2Te5 Chalcogenide

Chalcogenides are chemical compounds with at least one of the following three chemical elements: Sulfur (S), Selenium (Sn), and Tellurium (Te). As opposed to other materials, chalcogenide atomic arrangement can quickly and reversibly inter-change between crystalline, amorphous and liquid phases. Therefore they are also called phase change materials. As a results, chalcogenide thermal, optical, structural, electronic, electrical properties change pronouncedly and significantly with the phase they are in, leading to a host of different applications in different areas. The noticeable optical reflectivity difference between crystalline and amorphous phases has allowed optical storage devices to be made. Their very high thermal conductivity and heat fusion provided remarkable benefits in the frame of thermal energy storage for heating and cooling in residential and commercial buildings. The outstanding resistivity difference between crystalline and amorphous phases led to a significant improvement of solid state storage devices from the power consumption to the re-writability to say nothing of the shrinkability. This work focuses on a better understanding from a simulative stand point of the electronic, vibrational and optical properties for the crystalline phases (hexagonal and faced-centered cubic). The electronic properties are calculated implementing the density functional theory combined with pseudo-potentials, plane waves and the local density approximation. The phonon properties are computed using the density functional perturbation theory. The phonon dispersion and spectrum are calculated using the density functional perturbation theory. As it relates to the optical constants, the real part dielectric function is calculated through the Drude-Lorentz expression. The imaginary part results from the real part through the Kramers-Kronig transformation. The refractive index, the extinctive and absorption coefficients are analytically calculated from the dielectric function. The transmission and reflection coefficients are calculated using the Fresnel equations. All calculated optical constants compare well the experimental ones.

A new method for the estimation of the reverberation time from measured room impulse responses with application to Italian opera houses

The Schroeder's backward integration method is the most used method to extract the decay curve of an acoustic impulse response and to calculate the reverberation time from this curve. In the literature the limits and the possible improvements of this method are widely discussed. In this work a new method is proposed for the evaluation of the energy decay curve. The new method has been implemented in a Matlab toolbox. Its performance has been tested versus the most accredited literature method. The values of EDT and reverberation time extracted from the energy decay curves calculated with both methods have been compared in terms of the values themselves and in terms of their statistical representativeness. The main case study consists of nine Italian historical theatres in which acoustical measurements were performed. The comparison of the two extraction methods has also been applied to a critical case, i.e. the structural impulse responses of some building elements. The comparison underlines that both methods return a comparable value of the T30. Decreasing the range of evaluation, they reveal increasing differences; in particular, the main differences are in the first part of the decay, where the EDT is evaluated. This is a consequence of the fact that the new method returns a “locally" defined energy decay curve, whereas the Schroeder's method accumulates energy from the tail to the beginning of the impulse response. Another characteristic of the new method for the energy decay extraction curve is its independence on the background noise estimation. Finally, a statistical analysis is performed on the T30 and EDT values calculated from the impulse responses measurements in the Italian historical theatres. The aim of this evaluation is to know whether a subset of measurements could be considered representative for a complete characterization of these opera houses.

Ein Experiment zur Bestimmung des g-Faktors des gebundenen Elektrons in wasserstoff- und lithiumähnlichen mittelschweren Ionen

In dieser Arbeit werden der experimentelle Aufbau und erste Messungen für die Bestimmung des g-Faktors des Elektrons gebunden in wasserstoff- und lithiumähnlichen mittelschweren Ionen beschrieben. Mit dem hochpräzisenWert des g-Faktors können theoretische Berechnungen der Quantenelektrodynamik gebundener Zustände überprüft werden. Die Messungen werden in einem Dreifach-Penningfallen-System durchgeführt. Dort wurden im Rahmen dieser Arbeit auch erstmals hochgeladene Ionen bis 28Si13+ in einer hierfür entwickelten Elektronenstrahl-Ionenquelle/-falle erzeugt. Für die Bestimmung des g-Faktors werden die freie Zyklotronfrequenz und die Larmorfrequenz benötigt. Erstere wird aus den drei Eigenfrequenzen des in der Präzisionsfalle gespeicherten Ions berechnet. Um das Ion bei den Messungen nicht zu verlieren, werden die Eigenfrequenzen des Ions durch Kopplung an einen radiofrequenten Nachweisschwingkreis nicht-destruktiv nachgewiesen. Die freie Zyklotronfrequenz konnte dabei mit einer relativen Genauigkeit von wenigen 10E−9 bestimmt werden. Zur Bestimmung der Larmorfrequenz ist die genaue Kenntnis der Spinrichtung des Elektrons im Magnetfeld notwendig. Diese wird durch den kontinuierlichen Stern-Gerlach-Effekt in der sogenannten Analysefalle bestimmt. Hierzu muss eine hohe Stabilität der axialen Frequenz des Ions erreicht werden. Um dies sowie die Hochpräzisionsmessungen in der Präzisionsfalle zu erreichen, wurden in dieser Arbeit beide Fallen hinsichtlich ihrer elektrischen und magnetischen Eigenschaften charakterisiert.

Ancient Greek harbours used as geo-archives for palaeotsunami research : case studies from Krane (Cefalonia), Lechaion (Gulf of Corinth) and Kyllini (Peloponnese)

Since historical times, coastal areas throughout the eastern Mediterranean are exposed to tsunami hazard. For many decades the knowledge about palaeotsunamis was solely based on historical accounts. However, results from timeline analyses reveal different characteristics affecting the quality of the dataset (i.e. distribution of data, temporal thinning backward of events, local periodization phenomena) that emphasize the fragmentary character of the historical data. As an increasing number of geo-scientific studies give convincing examples of well dated tsunami signatures not reported in catalogues, the non-existing record is a major problem to palaeotsunami research. While the compilation of historical data allows a first approach in the identification of areas vulnerable to tsunamis, it must not be regarded as reliable for hazard assessment. Considering the increasing economic significance of coastal regions (e.g. for mass tourism) and the constantly growing coastal population, our knowledge on the local, regional and supraregional tsunami hazard along Mediterranean coasts has to be improved. For setting up a reliable tsunami risk assessment and developing risk mitigation strategies, it is of major importance (i) to identify areas under risk and (ii) to estimate the intensity and frequency of potential events. This approach is most promising when based on the analysis of palaeotsunami research seeking to detect areas of high palaeotsunami hazard, to calculate recurrence intervals and to document palaeotsunami destructiveness in terms of wave run-up, inundation and long-term coastal change. Within the past few years, geo-scientific studies on palaeotsunami events provided convincing evidence that throughout the Mediterranean ancient harbours were subject to strong tsunami-related disturbance or destruction. Constructed to protect ships from storm and wave activity, harbours provide especially sheltered and quiescent environments and thus turned out to be valuable geo-archives for tsunamigenic high-energy impacts on coastal areas. Directly exposed to the Hellenic Trench and extensive local fault systems, coastal areas in the Ionian Sea and the Gulf of Corinth hold a considerably high risk for tsunami events, respectively.Geo-scientific and geoarcheaological studies carried out in the environs of the ancient harbours of Krane (Cefalonia Island), Lechaion (Corinth, Gulf of Corinth) and Kyllini (western Peloponnese) comprised on-shore and near-shore vibracoring and subsequent sedimentological, geochemical and microfossil analyses of the recovered sediments. Geophysical methods like electrical resistivity tomography and ground penetrating radar were applied in order to detect subsurface structures and to verify stratigraphical patterns derived from vibracores over long distances. The overall geochronological framework of each study area is based on radiocarbon dating of biogenic material and age determination of diagnostic ceramic fragments. Results presented within this study provide distinct evidence of multiple palaeotsunami landfalls for the investigated areas. Tsunami signatures encountered in the environs of Krane, Lechaion and Kyllini include (i) coarse-grained allochthonous marine sediments intersecting silt-dominated quiescent harbour deposits and/or shallow marine environments, (ii) disturbed microfaunal assemblages and/or (iii) distinct geochemical fingerprints as well as (iv) geo-archaeological destruction layers and (v) extensive units of beachrock-type calcarenitic tsunamites. For Krane, geochronological data yielded termini ad or post quem (maximum ages) for tsunami event generations dated to 4150 ± 60 cal BC, ~ 3200 ± 110 cal BC, ~ 650 ± 110 cal BC, and ~ 930 ± 40 cal AD, respectively. Results for Lechaion suggest that the harbour was hit by strong tsunami impacts in the 8th-6th century BC, the 1st-2nd century AD and in the 6th century AD. At Kyllini, the harbour site was affected by tsunami impact in between the late 7th and early 4th cent. BC and between the 4th and 6th cent. AD. In case of Lechaion and Kyllini, the final destruction of the harbour facilities also seems to be related to the tsunami impact. Comparing the tsunami signals obtained for each study areas with geo-scientific data from palaeotsunami events from other sites indicates that the investigated harbour sites represent excellent geo-archives for supra-regional mega-tsunamis.

The oxidation capacity of the atmosphere - OH and HO 2 radical measurements in a boreal forest environment using laser induced fluorescence spectroscopy

Das wichtigste Oxidationsmittel für den Abbau flüchtiger Kohlenwasserstoffverbindungen (VOC, engl.: volatile organic compounds) in der Atmosphäre ist das Hydroxylradikal (OH), welches sich in einem schnellen chemischen Gleichgewicht mit dem Hydroperoxylradical (HO2) befindet. Bisherige Messungen und Modellvergleiche dieser Radikalspezies in Waldgebieten haben signifikante Lücken im Verständnis der zugrundeliegenden Prozesse aufgezeigt.rnIm Rahmen dieser Doktorarbeit wurden Messungen von OH- und HO2-Radikalen mittelsrnlaserinduzierten Fluoreszensmesstechnik (LIF, engl.: laser-induced fluorescence) in einem Nadelwald in Süd-Finnland während der Messkampagne HUMPPA–COPEC–2010 (Hyytiälä United Measurements of Photochemistry and Particles in Air – Comprehensive Organic Precursor Emission and Concentration study) im Sommer 2010 durchgeführt. Unterschiedliche Komponenten des LIF-Instruments wurden verbessert. Eine modifizierte Methode zur Bestimmung des Hintergrundsignals (engl.: InletPreInjector technique) wurde in den Messaufbaurnintegriert und erstmals zur Messung von atmosphärischem OH verwendet. Vergleichsmessungen zweier Instrumente basierend auf unterschiedlichen Methoden zur Messung von OH-Radikalen, chemische Ionisationsmassenspektrometrie (CIMS - engl.: chemical ionization mass spectrometry) und LIF-Technik, zeigten eine gute Übereinstimmung. Die Vergleichsmessungen belegen das Vermögen und die Leistungsfähigkeit des modifizierten LIF-Instruments atmosphärische OH Konzentrationen akkurat zu messen. Nachfolgend wurde das LIF-Instrument auf der obersten Plattform eines 20m hohen Turmes positioniert, um knapp oberhalb der Baumkronen die Radikal-Chemie an der Schnittstelle zwischen Ökosystem und Atmosphäre zu untersuchen. Umfangreiche Messungen - dies beinhaltet Messungen der totalen OH-Reaktivität - wurden durchgeführt und unter Verwendung von Gleichgewichtszustandsberechnungen und einem Boxmodell, in welches die gemessenen Daten als Randbedingungen eingehen, analysiert. Wenn moderate OH-Reaktivitäten(k′(OH)≤ 15 s−1) vorlagen, sind OH-Produktionsraten, die aus gemessenen Konzentrationen von OH-Vorläuferspezies berechnet wurden, konsistent mit Produktionsraten, die unter der Gleichgewichtsannahme von Messungen des totalen OH Verlustes abgeleitet wurden. Die primären photolytischen OH-Quellen tragen mit einem Anteil von bis zu einem Drittel zur Gesamt-OH-Produktion bei. Es wurde gezeigt, dass OH-Rezyklierung unter Bedingungen moderater OH-Reaktivität hauptsächlich durch die Reaktionen von HO2 mit NO oder O3 bestimmt ist. Während Zeiten hoher OH-Reaktivität (k′(OH) > 15 s−1) wurden zusätzliche Rezyklierungspfade, die nicht über die Reaktionen von HO2 mit NO oder O3, sondern direkt OH bilden, aufgezeigt.rnFür Hydroxylradikale stimmen Boxmodell-Simulationen und Messungen gut übereinrn(OHmod/OHobs=1.04±0.16), während HO2-Mischungsverhältnisse in der Simulation signifikant unterschätzt werden (HO2mod/HO2obs=0.3±0.2) und die simulierte OH-Reaktivität nicht mit der gemessenen OH-Reaktivität übereinstimmt. Die gleichzeitige Unterschätzung der HO2-Mischungsverhältnisse und der OH-Reaktivität, während OH-Konzentrationen von der Simulation gut beschrieben werden, legt nahe, dass die fehlende OH-Reaktivität in der Simulation eine noch unberücksichtigte HO2-Quelle darstellt. Zusätzliche, OH-unabhängigernRO2/HO2-Quellen, wie z.B. der thermische Zerfall von herantransportiertem peroxyacetylnitrat (PAN) und die Photolyse von Glyoxal sind indiziert.

Application and comparison of classification algorithms for recognition of Alzheimer's disease in electrical brain activity (EEG)

The early detection of subjects with probable Alzheimer's disease (AD) is crucial for effective appliance of treatment strategies. Here we explored the ability of a multitude of linear and non-linear classification algorithms to discriminate between the electroencephalograms (EEGs) of patients with varying degree of AD and their age-matched control subjects. Absolute and relative spectral power, distribution of spectral power, and measures of spatial synchronization were calculated from recordings of resting eyes-closed continuous EEGs of 45 healthy controls, 116 patients with mild AD and 81 patients with moderate AD, recruited in two different centers (Stockholm, New York). The applied classification algorithms were: principal component linear discriminant analysis (PC LDA), partial least squares LDA (PLS LDA), principal component logistic regression (PC LR), partial least squares logistic regression (PLS LR), bagging, random forest, support vector machines (SVM) and feed-forward neural network. Based on 10-fold cross-validation runs it could be demonstrated that even tough modern computer-intensive classification algorithms such as random forests, SVM and neural networks show a slight superiority, more classical classification algorithms performed nearly equally well. Using random forests classification a considerable sensitivity of up to 85% and a specificity of 78%, respectively for the test of even only mild AD patients has been reached, whereas for the comparison of moderate AD vs. controls, using SVM and neural networks, values of 89% and 88% for sensitivity and specificity were achieved. Such a remarkable performance proves the value of these classification algorithms for clinical diagnostics.

Reproducibility of retinal thickness measurements in healthy subjects using spectralis optical coherence tomography

PURPOSE: To test the reproducibility of retinal thickness measurements in healthy volunteers of a new Frequency-domain optical coherence tomography (OCT) device (Spectralis OCT; Heidelberg Engineering, Heidelberg, Germany). DESIGN: Prospective, observational study. METHODS: Forty-one eyes of 41 healthy subjects were included into the study. Intraobserver reproducibility was tested with 20 x 15 degree raster scans consisting of 37 high-resolution line scans that were repeated three times by one examiner (M.N.M.). Mean retinal thickness was calculated for nine areas corresponding to the Early Treatment Diabetic Retinopathy Study (ETDRS) areas. Coefficients of variation (COV) were calculated. RESULTS: Retinal thickness measurements were highly reproducible for all ETDRS areas. Mean total retinal thickness was 342 +/- 15 microm. Mean foveal thickness was 286 +/- 17 microm. COVs ranged from 0.38% to 0.86%. Lowest COV was found for the temporal outer ETDRS area (area 7; COV, 0.38%). Highest COV was found for the temporal inner ETDRS area (area 3; COV, 0.86%). Mean difference between measurement 1 and 2, measurement 1 and 3, and measurement 2 and 3 for all ETDRS areas was 1.01 microm, 0.98 microm, and 0.99 microm, respectively. CONCLUSION: Spectralis OCT retinal thickness measurements in healthy volunteers showed excellent intraobserver reproducibility with virtually identical results between retinal thickness measurements performed by one operator.

Called versus Calculated Yield Grades

An evaluation of carcass data collected over a two year period from southwest Iowa steer tests and 4-H carcass shows was conducted to compare USDA yield grades called by the Federal grader to yield grades calculated by actual carcass measurements. A regression equation was developed to predict called yield grade from carcass measurements. A comparison of the generated equation with the USDA equation used in calculating yield grades suggest that USDA graders accurately predict preliminary yield grades based on fat thickness, but may not have adequate time at line speeds to fully account for adjustments in ribeye size relative to carcass weight.

Can the carbon isotopic composition of methane be reconstructed from multi-site firn air measurements?

Methane is a strong greenhouse gas and large uncertainties exist concerning the future evolution of its atmospheric abundance. Analyzing methane atmospheric mixing and stable isotope ratios in air trapped in polar ice sheets helps in reconstructing the evolution of its sources and sinks in the past. This is important to improve predictions of atmospheric CH4 mixing ratios in the future under the influence of a changing climate. The aim of this study is to assess whether past atmospheric δ13C(CH4) variations can be reliably reconstructed from firn air measurements. Isotope reconstructions obtained with a state of the art firn model from different individual sites show unexpectedly large discrepancies and are mutually inconsistent. We show that small changes in the diffusivity profiles at individual sites lead to strong differences in the firn fractionation, which can explain a large part of these discrepancies. Using slightly modified diffusivities for some sites, and neglecting samples for which the firn fractionation signals are strongest, a combined multi-site inversion can be performed, which returns an isotope reconstruction that is consistent with firn data. However, the isotope trends are lower than what has been concluded from Southern Hemisphere (SH) archived air samples and high-accumulation ice core data. We conclude that with the current datasets and understanding of firn air transport, a high precision reconstruction of δ13C of CH4 from firn air samples is not possible, because reconstructed atmospheric trends over the last 50 yr of 0.3–1.5 ‰ are of the same magnitude as inherent uncertainties in the method, which are the firn fractionation correction (up to ~2 ‰ at individual sites), the Kr isobaric interference (up to ~0.8 ‰, system dependent), inter-laboratory calibration offsets (~0.2 ‰) and uncertainties in past CH4 levels (~0.5 ‰).

COMPARISON OF MEASUREMENTS OF ABSORBED DOSE TO WATER USING A WATER CALORIMETER AND IONIZATION CHAMBERS FOR CLINICAL RADIOTHERAPY PHOTON AND ELECTRON BEAMS

With the development of the water calorimeter direct measurement of absorbed dose in water becomes possible. This could lead to the establishment of an absorbed dose rather than an exposure related standard for ionization chambers for high energy electrons and photons. In changing to an absorbed dose standard it is necessary to investigate the effect of different parameters, among which are the energy dependence, the air volume, wall thickness and material of the chamber. The effect of these parameters is experimentally studied and presented for several commercially available chambers and one experimental chamber, for photons up to 25 MV and electrons up to 20 MeV, using a water calorimeter as the absorbed dose standard and the most recent formalism to calculate the absorbed dose with ion chambers.^ For electron beams, the dose measured with the calorimeter was 1% lower than the dose calculated with the chambers, independent of beam energy and chamber.^ For photon beams, the absorbed dose measured with the calorimeter was 3.8% higher than the absorbed dose calculated from the chamber readings. Such differences were found to be chamber and energy independent.^ The results for the photons were found to be statistically different from the results with the electron beams. Such difference could not be attributed to a difference in the calorimeter response. ^

Zircons, isotope ratios and element analyses from tonalite and oxide gabbro of the mid-ocean ridge from ODP Hole 176-735B

A limiting factor in the accuracy and precision of U/Pb zircon dates is accurate correction for initial disequilibrium in the 238U and 235U decay chains. The longest-lived-and therefore most abundant-intermediate daughter product in the 235U isotopic decay chain is 231Pa (T1/2 = 32.71 ka), and the partitioning behavior of Pa in zircon is not well constrained. Here we report high-precision thermal ionization mass spectrometry (TIMS) U-Pb zircon data from two samples from Ocean Drilling Program (ODP) Hole 735B, which show evidence for incorporation of excess 231Pa during zircon crystallization. The most precise analyses from the two samples have consistent Th-corrected 206Pb/238U dates with weighted means of 11.9325 ± 0.0039 Ma (n = 9) and 11.920 ± 0.011 Ma (n = 4), but distinctly older 207Pb/235U dates that vary from 12.330 ± 0.048 Ma to 12.140 ± 0.044 Ma and 12.03 ± 0.24 to 12.40 ± 0.27 Ma, respectively. If the excess 207Pb is due to variable initial excess 231Pa, calculated initial (231Pa)/(235U) activity ratios for the two samples range from 5.6 ± 1.0 to 9.6 ± 1.1 and 3.5 ± 5.2 to 11.4 ± 5.8. The data from the more precisely dated sample yields estimated DPazircon/DUzircon from 2.2-3.8 and 5.6-9.6, assuming (231Pa)/(235U) of the melt equal to the global average of recently erupted mid-ocean ridge basaltic glasses or secular equilibrium, respectively. High precision ID-TIMS analyses from nine additional samples from Hole 735B and nearby Hole 1105A suggest similar partitioning. The lower range of DPazircon/DUzircon is consistent with ion microprobe measurements of 231Pa in zircons from Holocene and Pleistocene rhyolitic eruptions (Schmitt (2007; doi:10.2138/am.2007.2449) and Schmitt (2011; doi:10.1146/annurev-earth-040610-133330)). The data suggest that 231Pa is preferentially incorporated during zircon crystallization over a range of magmatic compositions, and excess initial 231Pa may be more common in zircons than acknowledged. The degree of initial disequilibrium in the 235U decay chain suggested by the data from this study, and other recent high precision datasets, leads to resolvable discordance in high precision dates of Cenozoic to Mesozoic zircons. Minor discordance in zircons of this age may therefore reflect initial excess 231Pa and does not require either inheritance or Pb loss.

Radiocarbon ages and 230Th data of five sediment cores from the Björn and Gardar Drift, North Atlantic

In locations of rapid sediment accumulation receiving substantial amounts of laterally transported material the timescales of transport and accurate quantification of the transported material are at the focus of intense research. Here we present radiocarbon data obtained on co-occurring planktic foraminifera, marine haptophyte biomarkers (alkenones) and total organic carbon (TOC) coupled with excess Thorium-230 (230Thxs) measurements on four sediment cores retrieved in 1649-2879 m water depth from two such high accumulation drift deposits in the Northeast Atlantic, Björn and Gardar Drifts. While 230Thxs inventories imply strong sediment focussing, no age offsets are observed between planktic foraminifera and alkenones, suggesting that redistribution of sediments is rapid and occurs soon after formation of marine organic matter, or that transported material contains negligible amounts of alkenones. An isotopic mass balance calculation based on radiocarbon concentrations of co-occurring sediment components leads us to estimate that transported sediment components contain up to 12% of fossil organic matter that is free of or very poor in alkenones, but nevertheless appears to consist of a mixture of fresh and eroded fossil material. Considering all available constraints to characterize transported material, our results show that although focussing factors calculated from bulk sediment 230Thxs inventories may allow useful approximations of bulk redeposition, they do not provide a unique estimate of the amount of each laterally transported sediment component. Furthermore, our findings provide evidence that the occurrence of lateral sediment redistribution alone does not always hinder the use of multiple proxies but that individual sediment fractions are affected to variable extents by sediment focussing.

Age models and stabile isotope ratios from SO35 cores

Stable isotope measurements on the planktonic foraminifer Globigerinoides ruber (white) have been carried out on a number of selected deep-seas sediment cores from the South Lau and Norlh Fiji Basins. The d18O-curves show good correlation with the inter-ocean oraphic correlation composite d18O-record of the standard reference section (Prell et al. 1986), which, in combination with the chronostratigraphic classifications of Herterich & Sarnthein (1984, modified) and Imbrie et al. 1984), allows a detailed dating of the sedimentary sequences. The deepest layers in core no. 119 (southern Lau Basin) could be assigned to Isotope Stage 24. Measurements made on bulk carbonate in two cores show a much higher glacial-interglacial amplitude, allowing the general identification of the conventional oxygen isotope stages. The d13C-values of the benthic foraminifer Cibicidoides wuellerstorfi show progressively lighter values northwards reflecting an increasing contribution of the isotopically lighter CO2 from the remineralisation of organic matter during the general northward movement of the deep water masses. Cyclicities in the sedimentation rates were observed in core nos. 117 and 119 (both southern Lau Basin) where the interglacials exhibit higher levels than the glacials. Calculated new or export paleoproductivity show that the glacials had higher productivity in the euphotic zone. From the oxygen isotope stratigraphy, the five ash layers in core nos. 117 and 119 could be dated as about 530 ka B.P. in Stage 14, 695 ka B.P. in Stage 18, 775 ka B.P. in Stage 21, 790 ka B.P. and 825 ka B.P. in Stage 22. Carbonate dissolution occurred during stages 5, 8 and 10 to 12.

(Table 1) Physical properties of basalts from DSDP Holes 69-504B and 70-504B

A detailed study of physical properties was made on core samples from Deep Sea Drilling Project Hole 504B. The measured properties are density, porosity, sonic velocity, electrical resistivity, and fluid permeability. Basalts from this young oceanic crust have higher density and sonic velocity than the average DSDP basalts. Porosity (and temperature) dependences of physical properties are given by V = Vo - a-phi; roo = roo-0 exp(E*/RT)phi**-q; k = k0' phi**2q-1; where V is the sonic velocity (km/s), Vo = 6.45 (km/s), a = 0.111 (km/s %), phi is the porosity (%), roo is the electrical resistivity (ohm m), roo-0 = 0.002 (ohm m), E* = 2.7 (Kcal/mol) for fresh basalts, RT has its usual meaning, q = 1.67 ± 0.27, k is the permeability, k0' = (1 to about 10) x 10**-12 (cm**2). Porosity distribution in the crust in this area is estimated by combining the seismic velocity distribution and velocity-porosity relation. Because of the rapid decrease in porosity with depth, resistivity increases and permeability decreases rapidly with depth. The decreasing rate of permeability with increasing depth is approximately given by k(cm**2) = 2 x 10**-10 exp(-z (km)/0.3).