Shock-induced damage beneath impact craters is studied in this work. Two representative terrestrial rocks, San Marcos granite and Bedford limestone, are chosen as test target. Impacts into the rock targets with different combinations of projectile material, size, impact angle, and impact velocity are carried out at cm scale in the laboratory.
\r\n\r\nShock-induced damage and fracturing would cause large-scale compressional wave velocity reduction in the recovered target beneath the impact crater. The shock-induced damage is measured by mapping the compressional wave velocity reduction in the recovered target. A cm scale nondestructive tomography technique is developed for this purpose. This technique is proved to be effective in mapping the damage in San Marcos granite, and the inverted velocity profile is in very good agreement with the result from dicing method and cut open directly.
\r\n\r\nBoth compressional velocity and attenuation are measured in three orthogonal directions on cubes prepared from one granite target impacted by a lead bullet at 1200 m/s. Anisotropy is observed from both results, but the attenuation seems to be a more useful parameter than acoustic velocity in studying orientation of cracks.
\r\n\r\nOur experiments indicate that the shock-induced damage is a function of impact conditions including projectile type and size, impact velocity, and target properties. Combined with other crater phenomena such as crater diameter, depth, ejecta, etc., shock-induced damage would be used as an important yet not well recognized constraint for impact history.
\r\n\r\nThe shock-induced damage is also calculated numerically to be compared with the experiments for a few representative shots. The Johnson-Holmquist strength and failure model, initially developed for ceramics, is applied to geological materials. Strength is a complicated function of pressure, strain, strain rate, and damage. The JH model, coupled with a crack softening model, is used to describe both the inelastic response of rocks in the compressive field near the impact source and the tensile failure in the far field. The model parameters are determined either from direct static measurements, or from indirect numerical adjustment. The agreement between the simulation and experiment is very encouraging.
", "doi": "10.7907/39ZE-SY71", "publication_date": "2006", "thesis_type": "phd", "thesis_year": "2006" }, { "id": "thesis:3744", "collection": "thesis", "collection_id": "3744", "cite_using_url": "https://resolver.caltech.edu/CaltechETD:etd-09242008-095607", "type": "thesis", "title": "Collisional Processes Involving Icy Bodies in the Solar System", "author": [ { "family_name": "Stewart", "given_name": "Sarah Toby", "orcid": "0000-0001-9606-1593", "clpid": "Stewart-Sarah-Toby" } ], "thesis_advisor": [ { "family_name": "Ahrens", "given_name": "Thomas J.", "clpid": "Ahrens-T-J" }, { "family_name": "Yung", "given_name": "Yuk L.", "orcid": "0000-0002-4263-2562", "clpid": "Yung-Y-L" } ], "thesis_committee": [ { "family_name": "Ingersoll", "given_name": "Andrew P.", "orcid": "0000-0002-2035-9198", "clpid": "Ingersoll-A-P" }, { "family_name": "Ahrens", "given_name": "Thomas J.", "clpid": "Ahrens-T-J" }, { "family_name": "Yung", "given_name": "Yuk L.", "orcid": "0000-0002-4263-2562", "clpid": "Yung-Y-L" }, { "family_name": "Brown", "given_name": "Michael E.", "orcid": "0000-0002-8255-0545", "clpid": "Brown-M-E" }, { "family_name": "Goldreich", "given_name": "Peter Martin", "clpid": "Goldreich-P-M" }, { "family_name": "Kamb", "given_name": "W. Barclay", "clpid": "Kamb-W-B" } ], "local_group": [ { "literal": "div_gps" } ], "abstract": "1. The Shock Hugoniot of Solid Ice:
\r\n\r\nWe present a complete description of the solid ice Hugoniot based on new shock wave experiments conducted at an initial temperature of 100 K and previously published data obtained at 263 K. We identify five regions on the solid ice Hugoniot: (1) elastic shock waves, (2) ice Ih deformation shocks, transformation shocks to (3) ice VI, (4) ice VII, and (5) liquid water. In each region, data obtained at different initial temperatures are described by a single Us - Δup shock equation of state. The dynamic strength of ice Ih is strongly dependent on temperature. The Hugoniot Elastic Limit varies from 0.05 to 0.62 GPa, as a function of temperature and peak shock stress. We estimate the entropy and temperature along the 100 and 263 K Hugoniots and derive the critical pressures for shock-induced incipient (IM) and complete (CM) melting upon release. On the 100 K Hugoniot, the critical pressures are about 4.5 and between 5-6 GPa for IM and CM, respectively. On the 263 K Hugoniot, the critical pressures are 0.6 and 3.7 GPa for IM and CM, lower than previously suggested. Shock-induced melting of ice will be widespread in impact events.
\r\n\r\n2. Rampart Crater Formation on Mars:
\r\n\r\nWe present a model for the fluidization of Martian rampart crater ejecta blankets with liquid water based on the shock physics of cratering onto an ice-rich regolith. We conducted simulations of crater formation on Mars, explicitly accounting for the equations of state and shock-induced melting criteria for both the silicate and ice components and using strength models constrained by the observed transition diameter DTr from simple to complex craters on Mars, where DTr = 8 km corresponds to an effective yield strength of 107 Pa.
\r\n\r\nFor the observed size range of rampart craters (diameters D ≾ 30 km) and typical asteroidal impact conditions (silicate impactors, D ≾ 1 km, at 10 km s-1), we find that the hemispherical volume where subsurface ice is partially melted by the impact shock has a radius of about 15 projectile radii (rp), much larger than previous predictions of about 6 rp. The radius of the final crater is comparable to the radius of partial melting and more than half the ice within the excavated material is melted. Thus, the amount of shock-melted water incorporated into the continuous ejecta blanket is within a factor of two of the near-surface ground ice content.
\r\n\r\nWe find that fluidized ejecta blankets may form in the current climate with mean surface temperatures of 200 K. Decreasing the effective yield strength of the modeled materials, e.g., by increasing the ice content or porosity, modifies the impact-induced flow in the excavated cavity, resulting in deeper projectile penetration, steeper ejection angles, higher crater rim uplift, and reduced final crater diameter. The volume fraction of shock-melted water in the ejecta blanket increases with distance from the crater rim. The horizontal flow velocities during emplacement of fluidized ejecta (~ 10 - 1000 m s-1) is nearly constant in the continuous ejecta blanket and within the range of large terrestrial landslides. Therefore, ground-hugging debris flow conditions are achieved. The ejecta blanket properties from impacts into a Martian regolith containing 20-40%vol near-surface ice are consistent with the fraction of liquid water inferred from models of ejecta flow rheologies which produce rampart morphologies, about 10-30% liquid water by volume [Ivanov, B. A., Solar System Research, 30, 43-58, 1996].
\r\n\r\nWe present a model for the formation of different rampart ejecta morphologies which may be used in conjunction with an ejecta blanket debris flow model to map the distribution of ground ice. In addition, we find that formation of single or multiple-rampart ejecta blankets does not require pre-existing liquid water in the Martian crust. We estimate the minimum water content in observed rampart ejecta blankets to be equivalent to a global layer of water 0.6 m thick. Based on the crater sampling efficiency, the implied global Martian ice content, within the upper 2 km of the crust, is equivalent to a global layer of water 100 m deep. This result is comparable to other estimates of HM2O content in the Martian crust.
\r\n", "doi": "10.7907/0V4Q-VF61", "publication_date": "2002", "thesis_type": "phd", "thesis_year": "2002" }, { "id": "thesis:8168", "collection": "thesis", "collection_id": "8168", "cite_using_url": "https://resolver.caltech.edu/CaltechTHESIS:03252014-144322918", "primary_object_url": { "basename": "Ni_s_2001.pdf", "content": "final", "filesize": 28950625, "license": "other", "mime_type": "application/pdf", "url": "/8168/1/Ni_s_2001.pdf", "version": "v3.0.0" }, "type": "thesis", "title": "2D Modeling of Lower Mantle Structure with WKM Synthetics", "author": [ { "family_name": "Ni", "given_name": "Sidao", "orcid": "0000-0003-2988-4850", "clpid": "Ni-Sidao" } ], "thesis_advisor": [ { "family_name": "Ahrens", "given_name": "Thomas J.", "clpid": "Ahrens-T-J" } ], "thesis_committee": [ { "family_name": "Ahrens", "given_name": "Thomas J.", "clpid": "Ahrens-T-J" }, { "family_name": "Gurnis", "given_name": "Michael C.", "orcid": "0000-0003-1704-597X", "clpid": "Gurnis-M-C" }, { "family_name": "Helmberger", "given_name": "Donald V.", "clpid": "Helmberger-D-V" }, { "family_name": "Stevenson", "given_name": "David John", "orcid": "0000-0001-9432-7159", "clpid": "Stevenson-D-J" }, { "family_name": "Yung", "given_name": "Yuk L.", "orcid": "0000-0002-4263-2562", "clpid": "Yung-Y-L" } ], "local_group": [ { "literal": "div_gps" } ], "abstract": "The Earth is very heterogeneous, especially in the region close to the surface of the\r\nEarth, and in regions close to the core-mantle boundary (CMB). The lowermost\r\nmantle (bottom 300km of the mantle) is the place for fast anomaly (3% faster S\r\nvelocity than PREM, modeled from Scd), for slow anomaly (-3% slower S velocity than\r\nPREM, modeled from S,ScS), for extreme anomalous structure (ultra-low velocity\r\nzone, 30% lower inS velocity, 10% lower in P velocity). Strong anomaly with larger\r\ndimension is also observed beneath Africa and Pacific, originally modeled from travel\r\ntime of S, SKS and ScS. Given the heterogeneous nature of the earth, more accurate\r\napproach (than travel time) has to be applied to study the details of various anomalous\r\nstructures, and matching waveform with synthetic seismograms has proven effective\r\nin constraining the velocity structures. However, it is difficult to make synthetic\r\nseismograms in more than 1D cases where no exact analytical solution is possible.\r\nNumerical methods like finite difference or finite elements are too time consuming\r\nin modeling body waveforms. We developed a 2D synthetic algorithm, which is\r\nextended from 1D generalized ray theory (GRT), to make synthetic seismograms\r\nefficiently (each seismogram per minutes). This 2D algorithm is related to WKB\r\napproximation, but is based on different principles, it is thus named to be WKM, i.e.,\r\nWKB modified. WKM has been applied to study the variation of fast D\" structure\r\nbeneath the Caribbean sea, to study the plume beneath Africa. WKM is also applied\r\nto study PKP precursors which is a very important seismic phase in modeling lower\r\nmantle heterogeneity. By matching WKM synthetic seismograms with various data,\r\nwe discovered and confirmed that (a) The D\" beneath Caribbean varies laterally, and\r\nthe variation is best revealed with Scd+Sab beyond 88 degree where Sed overruns\r\nSab. (b) The low velocity structure beneath Africa is about 1500 km in height, at\r\nleast 1000km in width, and features 3% reduced S velocity. The low velocity structure\r\nis a combination of a relatively thin, low velocity layer (200 km thick or less) beneath\r\nthe Atlantic, then rising very sharply into mid mantle towards Africa. (c) At the\r\nedges of this huge Africa low velocity structures, ULVZs are found by modeling the\r\nlarge separation between S and ScS beyond 100 degree. The ULVZ to the eastern\r\nboundary was discovered with SKPdS data, and later is confirmed by PKP precursor\r\ndata. This is the first time that ULVZ is verified with distinct seismic phase.", "doi": "10.7907/J8XC-TF53", "publication_date": "2001", "thesis_type": "phd", "thesis_year": "2001" }, { "id": "thesis:775", "collection": "thesis", "collection_id": "775", "cite_using_url": "https://resolver.caltech.edu/CaltechETD:etd-02262008-153435", "type": "thesis", "title": "Impact-Induced Phase Transformations in Elastic Solids: A Continuum Study Including Numerical Simulations for GeO\u2082", "author": [ { "family_name": "Winfree", "given_name": "Nancy A.", "clpid": "Winfree-Nancy-A" } ], "thesis_advisor": [ { "family_name": "Knowles", "given_name": "James K.", "clpid": "Knowles-J-K" }, { "family_name": "Ahrens", "given_name": "Thomas J.", "clpid": "Ahrens-T-J" } ], "thesis_committee": [ { "family_name": "Knowles", "given_name": "James K.", "clpid": "Knowles-J-K" }, { "family_name": "Ahrens", "given_name": "Thomas J.", "clpid": "Ahrens-T-J" } ], "local_group": [ { "literal": "div_eng" } ], "abstract": "This thesis applies recently developed continuum theories of diffusionless phase transformations in solids to the study of impact problems involving materials which can experience such phase changes. Our objective is to compare the theoretical predictions against certain experimental results.
\r\n\r\nIn the experiments of interest, a face-to-face impact occurs between a disk of amorphous germanium dioxide and another material, either tungsten or an aluminum alloy. The GeO\u2082 is believed to transform to another phase if sufficient compressive stress is achieved.
\r\n\r\nWe model these experiments using one-dimensional finite elasticity. Phase-changing materials are represented by non-convex potential energy functions. This can produce phase boundaries that propagate subsonically or supersonically with respect to the slower longitudinal wave speed of the two phases. When a subsonic phase boundary is possible, it is not uniquely determined by the fundamental field equations and jump conditions. Uniqueness is obtained by invoking a nucleation criterion to control the initiation of the new phase, and a kinetic relation to govern its evolution.
\r\n\r\nThe experiments considered here are sufficiently long in duration (\u2248 3 \u00b5s) that several reflections and wave interactions occur, and the analysis becomes analytically intractable. Accordingly, a finite-difference method of Godunov type is employed to analyze these experiments numerically. Methods of Godunov type treat adjoining discretized spatial elements as the two sides of a Riemann problem, which is typically solved approximately by linearizing around the initial conditions on each side. Fortuitously, all constitutive models employed in this thesis are such that the required Riemann problems can be solved exactly without too much effort.
\r\n\r\nSimulations utilizing the numerical method demonstrate that the impact response of a material is sensitive to the kinetic relation that enters the model. It appears the theory may offer a plausible description of the experiments, though the restrictions placed on the constitutive models herein seem too severe to provide a good quantitative match to the experimental results.
", "doi": "10.7907/4dhf-fj83", "publication_date": "1999", "thesis_type": "phd", "thesis_year": "1999" }, { "id": "thesis:9923", "collection": "thesis", "collection_id": "9923", "cite_using_url": "https://resolver.caltech.edu/CaltechTHESIS:09122016-121754817", "type": "thesis", "title": "I. Rupture Properties of Large Subduction Earthquakes. II. Broadband Upper Mantle Structure of Western North America", "author": [ { "family_name": "Melbourne", "given_name": "Timothy Ian", "clpid": "Melbourne-Timothy-Ian" } ], "thesis_advisor": [ { "family_name": "Helmberger", "given_name": "Donald V.", "clpid": "Helmberger-D-V" }, { "family_name": "Ahrens", "given_name": "Thomas J.", "clpid": "Ahrens-T-J" } ], "thesis_committee": [ { "family_name": "Clayton", "given_name": "Robert W.", "orcid": "0000-0003-3323-3508", "clpid": "Clayton-R-W" }, { "family_name": "Farley", "given_name": "Kenneth A.", "orcid": "0000-0002-7846-7546", "clpid": "Farley-K-A" }, { "family_name": "Helmberger", "given_name": "Donald V.", "clpid": "Helmberger-D-V" }, { "family_name": "Kanamori", "given_name": "Hiroo", "orcid": "0000-0001-8219-9428", "clpid": "Kanamori-H" }, { "family_name": "Stock", "given_name": "Joann M.", "orcid": "0000-0003-4816-7865", "clpid": "Stock-J-M" } ], "local_group": [ { "literal": "div_gps" } ], "abstract": "This thesis contains two studies, one of which employs geodetic data bearing on large\r\nsubduction earthquakes to infer complexity of rupture duration, and the other is a\r\nhigh frequency seismological study of the upper mantle discontinuity structure under\r\nwestern North America and the East Pacific Rise. In the first part, we present\r\nGlobal Positioning System and tide gauge data which record the co-seismic deformation\r\nwhich accompanied the 1995 Mw8.0 Jalisco event offshore central Mexico, the\r\n1994 Mw7.5 Sanriku event offshore Northern Honshu, Japan, and the 1995 Mw8.1\r\nAntofagasta earthquake offshore Northern Chile. In two of the three cases we find\r\nthat the mainshocks were followed by significant amounts of rapid, post-seismic deformation\r\nwhich is best and most easily explained by continued slip near the co-seismic\r\nrupture patch. In Jalisco, we find that the post-seismic deformation which occurred\r\nduring the two weeks following the mainshock amounted to as much 70% of the co-seismic\r\ndeformation, from which we estimate an additional moment release of 40%,\r\nwhile in the Sanriku event an additional 30% moment release followed in the 10 days\r\nafter the mainshock. Because of the favorable geometry of the Jalisco network, we infer\r\nthat the post-seismic faulting occurred predominantly down-dip of the co-seismic\r\nrupture plane. This is the first documented case of rapid slip migration following\r\na large earthquake, and is pertinent to earthquake prediction based on precursory\r\ndeformation. Following the Antofagasta mainshock there was no rapid post-seismic\r\ndisplacement within the resolution of the GPS measurements, which equals roughly\r\n1% of the co-seismic displacement. As the three GPS data sets represent the best\r\nobservations of large subduction earthquakes to date and two of them show significant\r\namounts of aseismic energy release, they strongly suggest silent faulting may\r\nbe common in certain types of subduction zones. This, in turn, bears on estimates\r\nof global moment release, seismic coupling, and our understanding of the natural\r\nhazards associated with convergent margins.
\r\n\r\nThe second part of this dissertation utilizes high frequency body waves to infer\r\nthe upper mantle structure of western North America and the East Pacific Rise.\r\nAn uncharacteristically large Mw5.9 earthquake located in Western Texas provided a\r\nvivid topside reflection off the 410 Km velocity discontinuity (\"410\"), which we model\r\nto infer the fine details of this structure. We find that, contrary to conventional\r\nwisdom, the 410 is not sharp, and our results help reconcile seismic observations\r\nof 410 structure with laboratory predictions. By analyzing differences between our\r\nstructure and seismic 410 structure estimates under the nearby Gulf of California, we\r\nattempt to extract differences in temperature and mineralogy between subcontinental\r\nand suboceanic 410 structures.
\r\n\r\nExtending this analysis, we utilize teleseismic events from East Pacific Rise transform\r\nfaults to model multiple S upper mantle triplications. We find that for raypaths\r\ntraversing the rise crest the 1-D model TNA [Grand and Helmberger (1984)] derived\r\nfor the western US accurately predicts differential SnS-S travel times and triplication\r\nwaveform structure, implying that there is little velocity heterogeneity along the\r\nridge crest along nearly its entire length. We find that for energy traversing paths\r\nincreasingly away from the ridge axis there is no discernible change in the apparent\r\ndepth of the 410 and 670 Km discontinuities. In the shallowest mantle (uppermost 75\r\nKm), there is a strong lateral shear velocity gradient amounting to 3% over roughly\r\n150 Km. The LID, nonexistent at the ridge crest, grows slowly in thickness beyond\r\n150 Km from the axis. The compatible geodynamic model of these two results is that\r\nthe East Pacific Rise is not fed from the local lower mantle, rather, upper mantle\r\nmaterial must be transported laterally to supply the ridge axis spreading center, and\r\nthe LID reflects the source region of the East Pacific Rise magma supply.
", "doi": "10.7907/87nd-p040", "publication_date": "1999", "thesis_type": "phd", "thesis_year": "1999" }, { "id": "thesis:7626", "collection": "thesis", "collection_id": "7626", "cite_using_url": "https://resolver.caltech.edu/CaltechTHESIS:04232013-084302487", "type": "thesis", "title": "I. Rigid Body Penetration into Brittle Material. II. Phase Change Effect on Shock Wave Propagation", "author": [ { "family_name": "Liu", "given_name": "Cangli", "clpid": "Liu-Cangli" } ], "thesis_advisor": [ { "family_name": "Ahrens", "given_name": "Thomas J.", "clpid": "Ahrens-T-J" }, { "family_name": "Stock", "given_name": "Joann M.", "orcid": "0000-0003-4816-7865", "clpid": "Stock-J-M" } ], "thesis_committee": [ { "family_name": "Stock", "given_name": "Joann M.", "orcid": "0000-0003-4816-7865", "clpid": "Stock-J-M" }, { "family_name": "Ahrens", "given_name": "Thomas J.", "clpid": "Ahrens-T-J" }, { "family_name": "Clayton", "given_name": "Robert W.", "orcid": "0000-0003-3323-3508", "clpid": "Clayton-R-W" }, { "family_name": "Ravichandran", "given_name": "Guruswami", "orcid": "0000-0002-2912-0001", "clpid": "Ravichandran-G" }, { "family_name": "Ustundag", "given_name": "Ersan", "clpid": "Ustundag-E" } ], "local_group": [ { "literal": "div_gps" } ], "abstract": "Part I.
\r\n\r\nWe have developed a technique for measuring the depth time history of rigid body penetration into brittle materials (hard rocks and concretes) under a deceleration of ~ 105 g. The technique includes bar-coded projectile, sabot-projectile separation,\r\ndetection and recording systems. Because the technique can give very dense data on penetration depth time history, penetration velocity can be deduced. Error analysis shows that the technique has a small intrinsic error of ~ 3-4 % in time during penetration, and 0.3 to 0.7 mm in penetration depth. A series of 4140 steel projectile penetration into G-mixture mortar targets have been conducted using the Caltech 40 mm gas/ powder gun in the velocity range of 100 to 500 m/s.
\r\n\r\nWe report, for the first time, the whole depth-time history of rigid body penetration into brittle materials (the G-mixture mortar) under 105 g deceleration. Based on the experimental results, including penetration depth time history, damage of recovered target and projectile materials and theoretical analysis, we find:
\r\n\r\n1. Target materials are damaged via compacting in the region in front of a projectile and via brittle radial and lateral crack propagation in the region surrounding the penetration path. The results suggest that expected cracks in front of penetrators may be stopped by a comminuted region that is induced by wave propagation. Aggregate erosion on the projectile lateral surface is < 20% of the final penetration depth. This result suggests that the effect of lateral friction on the penetration process can be ignored.
\r\n\r\n2. Final penetration depth, Pmax, is linearly scaled with initial projectile energy per unit cross-section area, es , when targets are intact after impact. Based on the experimental data on the mortar targets, the relation is Pmax(mm) 1.15es (J/mm2 ) + 16.39.
\r\n\r\n3. Estimation of the energy needed to create an unit penetration volume suggests that the average pressure acting on the target material during penetration is ~ 10 to 20 times higher than the unconfined strength of target materials under quasi-static loading, and 3 to 4 times higher than the possible highest pressure due to friction and material strength and its rate dependence. In addition, the\r\nexperimental data show that the interaction between cracks and the target free surface significantly affects the penetration process.
\r\n\r\n4. Based on the fact that the penetration duration, tmax, increases slowly with es and does not depend on projectile radius approximately, the dependence of tmax on projectile length is suggested to be described by tmax(\u03bcs) = 2.08es (J/mm2 + 349.0 x m/(\u03c0R2), in which m is the projectile mass in grams and R is the projectile radius in mm. The prediction from this relation is in reasonable agreement with the experimental data for different projectile lengths.
\r\n\r\n5. Deduced penetration velocity time histories suggest that whole penetration history is divided into three stages: (1) An initial stage in which the projectile velocity change is small due to very small contact area between the projectile and target materials; (2) A steady penetration stage in which projectile velocity continues to decrease smoothly; (3) A penetration stop stage in which projectile deceleration jumps up when velocities are close to a critical value of ~ 35 m/s.
\r\n\r\n6. Deduced averaged deceleration, a, in the steady penetration stage for projectiles with same dimensions is found to be a(g) = 192.4v + 1.89 x 104, where v is initial projectile velocity in m/s. The average pressure acting on target materials during penetration is estimated to be very comparable to shock wave pressure.
\r\n\r\n7. A similarity of penetration process is found to be described by a relation between normalized penetration depth, P/Pmax, and normalized penetration time, t/tmax, as P/Pmax = f(t/tmax, where f is a function of t/tmax. After f(t/tmax is determined using experimental data for projectiles with 150 mm length, the penetration depth time history for projectiles with 100 mm length predicted by this relation is in good agreement with experimental data. This similarity also predicts that average deceleration increases with decreasing projectile length, that is verified by the experimental data.
\r\n\r\n8. Based on the penetration process analysis and the present data, a first principle model for rigid body penetration is suggested. The model incorporates the models for contact area between projectile and target materials, friction coefficient, penetration stop criterion, and normal stress on the projectile surface. The most important assumptions used in the model are: (1) The penetration process can be treated as a series of impact events, therefore, pressure normal to projectile surface is estimated using the Hugoniot relation of target material; (2) The necessary condition for penetration is that the pressure acting on target materials is not lower than the Hugoniot elastic limit; (3) The friction force on projectile lateral surface can be ignored due to cavitation during penetration. All the parameters involved in the model are determined based on independent experimental data. The penetration depth time histories predicted from the model are in good agreement with the experimental data.
\r\n\r\n9. Based on planar impact and previous quasi-static experimental data, the strain rate dependence of the mortar compressive strength is described by \u03c3f/\u03c30f = exp(0.0905(log(\u03ad/\u03ad_0) 1.14, in the strain rate range of 10-7/s to 103/s (\u03c30f and \u03ad are reference compressive strength and strain rate, respectively). The non-dispersive\r\nHugoniot elastic wave in the G-mixture has an amplitude of ~ 0.14 GPa and a velocity of ~ 4.3 km/s.
\r\n\r\nPart II.
\r\n\r\nStress wave profiles in vitreous GeO2 were measured using piezoresistance gauges in the pressure range of 5 to 18 GPa under planar plate and spherical projectile impact. Experimental data show that the response of vitreous GeO2 to planar shock loading\r\ncan be divided into three stages: (1) A ramp elastic precursor has peak amplitude of 4 GPa and peak particle velocity of 333 m/s. Wave velocity decreases from initial longitudinal elastic wave velocity of 3.5 km/s to 2.9 km/s at 4 GPa; (2) A ramp wave with amplitude of 2.11 GPa follows the precursor when peak loading pressure is 8.4 GPa. Wave velocity drops to the value below bulk wave velocity in this stage; (3) A shock wave achieving final shock state forms when peak pressure is > 6 GPa. The Hugoniot relation is D = 0.917 + 1.711u (km/s) using present data and the data of Jackson and Ahrens [1979] when shock wave pressure is between 6 and 40 GPa for \u03c10 = 3.655 gj cm3 . Based on the present data, the phase change from 4-fold to 6-fold coordination of Ge+4 with O-2 in vitreous GeO2 occurs in the pressure range of 4 to 15 \u00b1 1 GPa under planar shock loading. Comparison of the shock loading data for\r\nfused SiO2 to that on vitreous GeO2 demonstrates that transformation to the rutile structure in both media are similar. The Hugoniots of vitreous GeO2 and fused SiO2 are found to coincide approximately if pressure in fused SiO2 is scaled by the ratio of fused SiO2to vitreous GeO2 density. This result, as well as the same structure, provides the basis for considering vitreous Ge02 as an analogous material to fused SiO2 under shock loading. Experimental results from the spherical projectile impact demonstrate: (1) The supported elastic shock in fused SiO2 decays less rapidly than a linear elastic wave when elastic wave stress amplitude is higher than 4 GPa. The supported elastic shock in vitreous GeO2 decays faster than a linear elastic wave; (2) In vitreous GeO2 , unsupported shock waves decays with peak pressure in the phase transition range (4-15 GPa) with propagation distance, x, as \u03b1 1/x-3.35 , close to the prediction of Chen et al. [1998]. Based on a simple analysis on spherical wave propagation, we find that the different decay rates of a spherical elastic wave in fused SiO2 and vitreous GeO2 is predictable on the base of the compressibility variation with stress under one-dimensional strain condition in the two materials.
\r\n\r\n", "doi": "10.7907/w8gc-n615", "publication_date": "1999", "thesis_type": "phd", "thesis_year": "1999" }, { "id": "thesis:16202", "collection": "thesis", "collection_id": "16202", "cite_using_url": "https://resolver.caltech.edu/CaltechTHESIS:10062023-181634008", "primary_object_url": { "basename": "Holland_KG_1997.pdf", "content": "final", "filesize": 30615462, "license": "other", "mime_type": "application/pdf", "url": "/16202/1/Holland_KG_1997.pdf", "version": "v2.0.0" }, "type": "thesis", "title": "Phase Changes and Transport Properties of Geophysical Materials Under Shock Loading", "author": [ { "family_name": "Holland", "given_name": "Kathleen Gabrielle", "clpid": "Holland-Kathleen-Gabrielle" } ], "thesis_advisor": [ { "family_name": "Ahrens", "given_name": "Thomas J.", "clpid": "Ahrens-T-J" }, { "family_name": "Stock", "given_name": "Joann M.", "orcid": "0000-0003-4816-7865", "clpid": "Stock-J-M" } ], "thesis_committee": [ { "family_name": "Stock", "given_name": "Joann M.", "orcid": "0000-0003-4816-7865", "clpid": "Stock-J-M" }, { "family_name": "Ahrens", "given_name": "Thomas J.", "clpid": "Ahrens-T-J" }, { "family_name": "Stevenson", "given_name": "David John", "orcid": "0000-0001-9432-7159", "clpid": "Stevenson-D-J" }, { "family_name": "Wyllie", "given_name": "Peter J.", "clpid": "Wyllie-P-J" }, { "family_name": "Gurnis", "given_name": "Michael C.", "orcid": "0000-0003-1704-597X", "clpid": "Gurnis-M-C" } ], "local_group": [ { "literal": "div_gps" } ], "abstract": "The lower mantle of the Earth is believed to be largely composed of (Mg, Fe)O (magnesiowiistite) and (Mg,Fe)SiO\u2083 (perovskite); thus the high pressure phase of (Mg,Fe)\u2082SiO\u2084 (olivine), which is believed to be perovskite plus magnesiowiistite is of geophysical interest. Radiative temperatures of single-crystal olivine starting material [(Mg_(0.9), Fe_(0.1))\u2082SiO\u2084] decreased abruptly from 7040 \u00b1 315 to 4300 \u00b1 270 K upon shock compression above 80 GPa. The data indicate that an upper bound to the solidus of the magnesiowiistite and perovskite assemblage at 4300 \u00b1 270 K is 130 \u00b1 3 GPa. These conditions correspond to those for partial melting at the base of the mantle, as has been suggested to occur within the recently discovered ultra-low-velocity zone (ULVZ) beneath the Central Pacific. We construct speculative high pressure phase diagrams for the MgO - SiO\u2082 system using experimental data from our work, and other mineral physics experiments.
\r\n\r\nIn separate experiments, time dependent shock temperatures were measured for stainless steel (SS) films sandwiched between two transparent Al\u2082O\u2083 anvils. The anvil material was the same as the driver material so that there would be symmetric heat flow from the sample. Inferred Hugoniot temperatures, T_h, of 5000 - 8500\u00b1500 Kat 222- 321 GPa are consistent with previous measurements in SS. Temperatures at the film\u00ad anvil interface (T_i), which are directly measured (rather than T_h) indicate that T_i did not decrease measurably during the approximately 250 ns that the shock wave took to traverse the Al\u2082O\u2083 anvil. Thus an upper bound is obtained for the thermal diffusivity of Al\u2082O\u2083 at the metal/anvil interface of K \u2264 14 \u00b1 5 cm\u00b2/s at 208 GPa and 2110 K. This is a factor of 1.6 lower than previously calculated values, resulting in a decrease of the inferred T_h by at least 400 K. The observed shock temperatures are combined with temperatures calculated from measured Hugoniots and are used to calculate the thermal conductivity of Al\u2082O\u2083. There was no measurable radiant-intensity decrease during the time when the shock wave propagated through the anvil; we infer from this that Al\u2082O\u2083 remained transparent while in the shocked state. Thus an Al\u2082O\u2083 anvil is sufficiently transparent for shock temperature measurements for metals, to at least 240 GPa.
\r\n\r\nFinally, shock temperature experiments employing a six-channel pyrometer were conducted on 200, 500, and 1000 \u00c5 thick films of Fe sandwiched between 3 mm thick anvils of Al\u2082O\u2083 and LiF, to measure the thermal diffusivity ratios of Al\u2082O\u2083/Fe and LiF/Fe, at high temperatures and pressures. Temperature decays of 3000 \u00b1 800 K in 250 ns were observed at Fe pressures of 194 - 303 GPa, which reflect the conduction of heat from the thin metal films into the anvil material. These results were achieved in experiments employing LiF anvils at 164 - 166 GPa and 4190 - 4220 K, and Al2O3 anvils at 196 - 303 GPa and 1410 - 2750 K. Thermal modeling of interface temperature versus time yields best fit thermal diffusivity ratios ranging from 15 \u00b1 30 to 80 \u00b1 20 (Fe/anvil) over the pressure and temperature range of the experiments. Calculated thermal conductivities for Fe, using electron gas theory, of 110 - 212 W /mK are used to calculate thermal conductivities for the anvil materials ranging from 6 to 12 W/mK. Debye theory predicts higher values of 8 to 34 W/mK. Data from previous experiments on thick (\u2265 100\u00b5m) films of Fe and stainless steel are combined with our present results from experiments on thin (\u2264 1000 \u00c5) films to infer a 5860 \u00b1 390 K Hugoniot temperature for the onset of melting of iron at 243 GPa. Our results address the question of whether radiation observed in shock temperature experiments on metals originates from the metal at the metal/ anvil interface or from the shocked anvil. We conclude that the photon flux from the shocked iron/anvil sandwich recorded in all experiments originates from the metal. Within the uncertainties of the shock temperature data, the uncertainties in shock temperatures resulting from the radiation from the anvils is negligible. This is in direct disagreement with previous conclusions of Kondo.
", "doi": "10.7907/cbd3-mk54", "publication_date": "1999", "thesis_type": "phd", "thesis_year": "1999" }, { "id": "thesis:770", "collection": "thesis", "collection_id": "770", "cite_using_url": "https://resolver.caltech.edu/CaltechETD:etd-02262007-141751", "primary_object_url": { "basename": "Chen_g_1998.pdf", "content": "final", "filesize": 10120655, "license": "other", "mime_type": "application/pdf", "url": "/770/1/Chen_g_1998.pdf", "version": "v3.0.0" }, "type": "thesis", "title": "I. High pressure melting of [gamma]-iron and the thermal profile in the Earth's core. II. High pressure, high temperature equation of state of fayalite (Fe2SiO4)", "author": [ { "family_name": "Chen", "given_name": "George (Guangqing)", "clpid": "Chen-G-G" } ], "thesis_advisor": [ { "family_name": "Ahrens", "given_name": "Thomas J.", "clpid": "Ahrens-T-J" } ], "thesis_committee": [ { "family_name": "Ahrens", "given_name": "Thomas J.", "clpid": "Ahrens-T-J" }, { "family_name": "Burnett", "given_name": "Donald S.", "clpid": "Burnett-D-S" }, { "family_name": "Stolper", "given_name": "Edward M.", "clpid": "Stolper-E-M" }, { "family_name": "Knowles", "given_name": "James K.", "clpid": "Knowles-J-K" }, { "family_name": "Tombrello", "given_name": "Thomas A.", "clpid": "Tombrello-T-A" } ], "local_group": [ { "literal": "div_pma" } ], "abstract": "NOTE: Text or symbols not renderable in plain ASCII are indicated by [...]. Abstract is included in .pdf document.\n\nThe melting curve of [...]-iron in the pressure range of 100 to 300 GPa has been derived by computing Gibbs free energies at high pressures and high temperatures from thermodynamic and equations of state (EOS) data for the [...], [...] and liquid-phases. Our calculations indicate the melting curve of iron is very sensitive to the EOS of both the solid [...] and melt phase. Our best estimate of the EOS parameters for [...]-iron are: p0 = 8.775 \u00b1 0.012 Mg/m3, [...] = 205 \u00b1 4 GPa, [...] = 4.80 \u00b1 0.01 (referenced to 12 GPa and 300 K). The calculation favors the melting curve of Boehler [1993] or Saxena et al. [1993]. Shock-wave experiments on pure iron preheated to 1573 K were conducted in 17-73 GPa range. The shock-wave equation of state of [gamma]-iron at 1573 K initial temperature can be fit with [...] = 4.102(0.015) km/s + 1.610(0.014)[...] with [...] = 7.413 \u00b1 0.012 Mg/m3. [Gamma]-iron's bulk modulus and its pressure derivative are 124.7\u00b11.1 GPa and 5.44\u00b10.06 respectively.\n\nWe present new data for sound velocities in the [gamma]- and liquid-phases. In the [gamma]-phase, to a first approximation, the longitudinal sound velocity is linear with respect to density: Vp = -3.13(0.72) + 1.119(0.084) [...] (units for Vp and [...] are km/s and Mg/m3 respectively). Melting was observed in the highest pressure (about 70-73 GPa) experiments at a calculated shock temperature of about 2775 \u00b1 160 K. This result is consistent with our calculated [...]-iron melting curve which is close to those measured by Boehler [1993] and Saxena et al. [1993]. The liquid iron sound velocity data yield a Gruneisen parameter value for liquid iron of 1.63\u00b10.28 at 9.37\u00b10.02 Mg/m3 at 71.6 GPa. The quantity [...] from our data is 15.2\u00b12.6 Mg/m3, which is within the bounds of Brown and McQueen [1986] (13.3-19.6 Mg/m3). Based on upward pressure and temperature extrapolation of our melting curve of [gamma]-iron, the estimated inner core-outer core boundary temperature is 5500\u00b1400 K, the temperature at the core-mantle boundary on the outer core side is about 3930\u00b1630 K, and the thermal boundary layer at the core-mantle boundary has a temperature difference between 400 and 1400 K.\n\nThe shock-wave equation of state of initially solid (300 K) and molten (1573 K) fayalite (Fe2SiO4, Fa) are reported in the ranges 23 to 212 GPa and 5 to 47 GPa, respectively. The 300 K data appear to undergo a phase change in the 35-55 GPa range. The density of the high pressure phase (HPP) is consistent with a dense oxide mixture. Although the initially 300 K fayalite may melt along its Hugoniot, this is not explicitly detected. Fitting the HPP Hugoniot data in the shock velocity ([...])-particle velocity ([...]) plane yields:\n\n[...] = 4.375(0.027) Mg/m3, (1)\n\n[...] = 4.07(0.22) km/s + 1.43(0.06) [...], (2)\n\nwhere [...] is the initial density. The isentropic bulk modulus [...] = 72.4 \u00b1 8.0 GPa, and its pressure derivative [...] = 4.72\u00b10.24. \n\nThe 1573 K data set yields:\n \n[...] = 3.750(0.018) Mg/m3, (3)\n \n[...] = 2.63(0.02) km/s + 1.59(0.01) [...], (4)\n\nand [...] = 25.9 \u00b1 0.4 GPa, [...] = 5.36 \u00b1 0.04. The bulk modulus compares favorably with Agee [1992a]'s result (24.4 GPa), but the pressure derivative is quite different (10.1 from Agee [1992a]).\n\nAbove 50 GPa, the high pressure regime of the Hugoniot of the solid fayalite can be fit with oxide mixture models using stishovite and FeO (either LPP or HPP). The fayalitic liquid compression data above 40 GPa are well fit with ideal mixing of partial molar volumes of stishovite and FeO (LPP or HPP), in support of the hypothesis of Rigden et al. [1989].\n\nA model basalt incorporating the liquid fayalite data shows the neutral buoyancy zone of basic silicate melts of plausible terrestrial compositions is at about 250-400 km depth based on the PREM Earth model.\n", "doi": "10.7907/z7r2-xf47", "publication_date": "1998", "thesis_type": "phd", "thesis_year": "1998" }, { "id": "thesis:5932", "collection": "thesis", "collection_id": "5932", "cite_using_url": "https://resolver.caltech.edu/CaltechTHESIS:06072010-160325894", "primary_object_url": { "basename": "Zajac_bj_1998.pdf", "content": "final", "filesize": 13435545, "license": "other", "mime_type": "application/pdf", "url": "/5932/1/Zajac_bj_1998.pdf", "version": "v8.0.0" }, "type": "thesis", "title": "The State of Stress as Inferred from Deviated Boreholes: Constraints on the Tectonics of Offshore Central California and Cook Inlet, Alaska", "author": [ { "family_name": "Zajac", "given_name": "Blair J.", "clpid": "Zajac-Blair-J" } ], "thesis_advisor": [ { "family_name": "Stock", "given_name": "Joann M.", "orcid": "0000-0003-4816-7865", "clpid": "Stock-J-M" }, { "family_name": "Ahrens", "given_name": "Thomas J.", "clpid": "Ahrens-T-J" } ], "thesis_committee": [ { "family_name": "Ahrens", "given_name": "Thomas J.", "clpid": "Ahrens-T-J" }, { "family_name": "Clayton", "given_name": "Robert W.", "orcid": "0000-0003-3323-3508", "clpid": "Clayton-R-W" }, { "family_name": "Gurnis", "given_name": "Michael C.", "orcid": "0000-0003-1704-597X", "clpid": "Gurnis-M-C" }, { "family_name": "Saleeby", "given_name": "Jason B.", "clpid": "Saleeby-J-B" }, { "family_name": "Stock", "given_name": "Joann M.", "orcid": "0000-0003-4816-7865", "clpid": "Stock-J-M" } ], "local_group": [ { "literal": "div_gps" } ], "abstract": "This thesis introduces a new method of constraining the vector directions of the three principal stresses and their relative magnitudes, by using borehole breakouts in non-vertical drill holes. Unlike older stress state measurements from breakouts, this work does not presume that one of the principal stresses is Vertical. This method has important uses in complicated three-dimensional structures, such as in the Los Angeles basin, and in oil drilling applications.
\r\n\r\nChapter 1 discusses why knowledge of the three-dimensional stress tensor is relevant to today's science and examines the applications of the stress state determination technique discussed herein. The history of previous work is also described.
\r\n\r\nIn Chapter 2 I discuss the techniques of determining the stress tensor from borehole breakouts, examining the physics of borehole breakouts, the theory of the inversion technique used, and data processing issues. The theory and data processing issues are not discussed separately in this work, since data processing issues often prompted new theoretical techniques. I first examine the physics of borehole breakouts and how the orientation of breakouts on the borehole wall relates to the local stress field. A new borehole breakout selection scheme which takes into account highly non-vertical boreholes is then presented along with a discussion of the real world problems of data gathering, identification, and processing. Having selected a borehole breakout data set using the criteria, I invert for the best fitting stress state using a new technique combining genetic algorithms and non- differential function optimizers. Finally, I present a way in which 95% confidence limits can be placed on the resulting stress tensor.
\r\n\r\nWith all of the technical and theoretical pieces in place, I now examine several different data sets. Chapter 3 examines a borehole breakout data set publish by Qian and Pedersen [1991] from the Siljan Deep Drilling Project in Sweden and demonstrates that even for simple borehole breakout data sets, the stress state inversions assuming a vertical principal stress direction may fall outside of the 95% confidence limits of an inversion allowing non-vertical principal stress directions. My technique of displaying the borehole breakout data makes the data quality more obvious as compared to the way Qian and Pedersen [1991] plotted the data.
\r\n\r\nChapter 4 examines a borehole breakout data set from the offshore Santa Maria Basin, California. This analysis presents vertical borehole breakout data that represent a maximum horizontal principal stress direction of N7\u00b0E, roughly consistent with other earthquake focal mechanism, GPS, and borehole breakout studies in the area. However, the stress state inversion of breakouts identified in the vertical and a limited number of nearly horizontal boreholes suggests a stress state very different from any other stress state results. This could imply that the three dimensional stress in the Santa Maria Basin is very complicated. However, given the limited amount of borehole breakouts identified in nearly horizontal wells, the stress state results from this data set are inconclusive.
\r\n\r\nChapter 5 examines the largest data set used in this study, from a series of oil wells in Cook Inlet, Alaska. These are borehole caliper arm data from 21 different wells reaching a maximum deviation of 54\u00b0 and 3,223 m true vertical depth. Stress state inversions of 31 different subsets of the borehole breakout data were performed. Inversion of breakouts identified in the top two of three marker beds analyzed in wells drilled from the Baker platform identified nearly degenerate thrust faulting stress states with the maximum principal stress axis, S_1, oriented horizontally WNWESE, perpendicular to the NNE-trending anticlinal structures. The stress state from the deepest marker is also a nearly degenerate thrust faulting stress state with S_1 oriented NNW\u2014SSE, aligned with the regional direction of relative plate motion between the North American and Pacific plates. In between the shallow and deep stress state is an apparent normal faulting stress state with S_2 oriented subhorizontally ENE\u2014WSW. This clockwise rotation of the stress tensor as a function of depth suggests that the stress field changes with depth, from a shallow stress state responsible for the local NNE-trending structures to a deeper one from the North American and Pacific plates' collision zone. The observed normal faulting stress state between the two thrust faulting stress states is anomalous and may represent some sort of transition from the shallow to the deep stress state. Stress state profiles in 500 m true vertical depth (TVD) intervals show consistently oriented thrust faulting stress regimes with NNW\u2014SSE trending S_1 azimuths. The thrust faulting S_3 principal stress direction is consistently within 30\u00b0 of vertical, suggesting that while the assumption of a purely vertical principal stress direction is not valid, the stress tensor does not significantly rotate away from the surface conditions that require a purely vertical stress tensor. The nearly degenerate thrust faulting stress states determined from the Granite Point and the 10.8 km distant Baker platform breakouts are nearly identical, implying that the technique of using deviated borehole breakouts to invert for the regional stress is valid. The orientations of the maximum horizontal stress determined from the Cook Inlet borehole breakouts are consistent with other stress indicators in south-central Alaska and consistent with the direction of relative plate motion between the North American Plate and the Pacific plate. The S_1 axis for the Cook Inlet field trends due south plunging 3\u00b0. The 95% confidence limits allow the S_1 azimuth to vary from N156\u00b0E to N195\u00b0E and the plunge to vary from 10\u00b0 to -4\u00b0. This stress state does not appear representative of the stress field for each subset of breakouts. The Granite Point S1 axis trends N19\u00b0W plunging 3\u00b0; the 95% confidence limits allow the azimuth to vary from N42\u00b0W to N7\u00b0E and the plunge to vary from 1\u00b0 to 6\u00b0. The Baker platform S_1 axis trends N170\u00b0E plunging 8\u00b0; the 95% confidence limits on S_1 allow its azimuth to vary from N139\u00b0E to N191\u00b0E and its plunge to vary from 1\u00b0 to 15\u00b0. Finally, the Dillon platform S_1 axis trends N69\u00b0W plunging 2\u00b0; the 95% confidence limits constrain the S_1 azimuth from N268\u00b0E to N324\u00b0E and the plunge from 8\u00b0 to -4\u00b0. The more westerly orientation of S_1 at the Dillon platform may be related to the local NNE-trending anticlinal structures in the Cook Inlet Basin.
\r\n\r\nChapter 6 concludes and summarized the results and conclusions from the thesis.
\r\n\r\nThe first appendix contains in minute detail some of the mathematics describing the boreholes, breakouts, and coordinate system rotations used to perform this work. The second appendix contains the individual discussion and plots of the raw dipmeter data from all of the Cook Inlet, Alaska wells.
", "doi": "10.7907/S3V0-3M43", "publication_date": "1998", "thesis_type": "phd", "thesis_year": "1998" }, { "id": "thesis:4363", "collection": "thesis", "collection_id": "4363", "cite_using_url": "https://resolver.caltech.edu/CaltechETD:etd-11012007-084634", "primary_object_url": { "basename": "Takata_t_1995.pdf", "content": "final", "filesize": 23493649, "license": "other", "mime_type": "application/pdf", "url": "/4363/1/Takata_t_1995.pdf", "version": "v3.0.0" }, "type": "thesis", "title": "Three-dimensional analysis of impact processes on planets", "author": [ { "family_name": "Takata", "given_name": "Toshiko", "clpid": "Takata-T" } ], "thesis_advisor": [ { "family_name": "Ahrens", "given_name": "Thomas J.", "clpid": "Ahrens-T-J" } ], "thesis_committee": [ { "family_name": "Unknown", "given_name": "Unknown" } ], "local_group": [ { "literal": "div_gps" } ], "abstract": "NOTE: Text or symbols not renderable in plain ASCII are indicated by [...]. Abstract is included in .pdf document.\r\n\r\nThis thesis consists of four chapters. The first chapter describes the numerical method known as Smoothed Particle Hydrodynamics (SPH), that is used in present calculations relating to oblique impact on a planet with an atmosphere. Numerical test results are presented to demonstrate the validity of the present computer code. The SPH code is applied to Comet Shoemaker-Levy 9 impact on Jupiter, particularly the disintegration of the comet in the Jovian atmosphere and the evolution of the vapor plume. These calculations are described in chapter II. In order to interpret the phenomena observed in the impacts, numerical predictions are compared with the observational data. These comparisons are summarized in chapter III. Chapter IV considers analytical models of the formation of radar dark and/or bright halos surrounding impact craters discovered on Venus as a result of the Magellan mission. The surface features caused by the atmospheric shock waves associated with impacts are modeled and applied to the observed Venusian radar features.\r\n\r\nPaper I:\r\n\r\nThe SPH method used for a series of impact calculations is tested. From the performance tests, we infer that the tree structure method for a vectorizing computer code can be applied when the number of particles, N, [...]. Most of the calculations are performed using a CRAY-YMP at Jet Propulsion Laboratory or a CRAY-C90 at Goddard Space Flight Center. Tests for one-dimensional shock propagation show that the developed code yields an error of less than several percent in density and pressure, upon comparison with the analytical impedance-match solution. For spherical, three-dimensional impact of a projectile upon a silicate half-space, the shock wave attenuation of pressure and density calculated by the SPH code is very similar to the results obtained with a finite-difference Lagrangian continuum code (previously obtained by Ahrens and O'Keefe, [1977]). For example, the shock pressures in the near field agree within 15 % differences, whereas the shock pressure in the far field, at 10 projectile radii, differs by a factor of 2 from that calculated previously using a two-dimensional axisymmetric Lagrangian code.\r\n\r\nPaper II:\r\n\r\nThe impact of fragments of Comet Shoemaker-Levy 9 on Jupiter and the resulting vapor plume expansion are investigated by conducting three-dimensional numerical simulations using the Smoothed Particle Hydrodynamics (SPH) method. An icy body, representing the cometary fragments, with a velocity of 60 km/sec and a diameter of 2 km can penetrate to 350 km below the 1-bar pressure level in the atmosphere. Most of the initial kinetic energy of the fragment is transferred to the atmosphere between 50 km and 300 km below the 1-bar pressure level. The shock-heated atmospheric gas in the wake is totally dissociated and partially ionized. Scaling our SPH results to other sizes indicates that fragments larger than [...] 100 m in diameter can penetrate to below the visible cloud decks. The energy deposited in the atmosphere is explosively released in the upward expansion of the resulting plume. The plume preferentially expands upward rather than horizontally due to the density gradient of the ambient atmosphere. It rises ...] km in [...] seconds. Eventually the total atmospheric mass ejected to above 1 bar is [...] 40 times the initial mass of the impactor. The plume temperature at a radius ~ [...] km is > [...] K for [...] seconds, for a 2-km fragment. We predict that impact-induced plumes will be observable with the remote sensing instruments of the Galileo spacecraft. As the impact site rotates into the view of the earth some 20 minutes after the impact, the plume expansion will be observable using the Hubble Space Telescope and from visible and infrared instruments on ground-based telescopes. The rising plume reaches ~ 3000 km altitude in ~ ten minutes and will be visible from the earth.\r\n\r\nPaper III:\r\n\r\nPreliminary observational data from the impact of fragments of Comet Shoemaker-Levy 9 (SL9) are compared with smoothed particle hydrodynamic (SPH) and radiative calculations to determine the energies of individual SL9 fragments and the equivalent diameter of the SL9 progenitor. The Hubble Space Telescope (HST) images of the G impact-induced plume demonstrate that it achieved a height of ~3300 km [Hammel et al., 1994]. This is in close agreement with the SPH calculations for a 7 x [...] erg or 2 km \u00b1 0.3 km diameter solid ice impactor at 60 km/sec. Comparison of the R fragment impact-induced plume brightness in the 8 to 12 [...] band, as recorded by the NASA Infrared Telescope Facility (Mauna Kea, HI) [Orton et al. 1994], with our radiative predictive calculations yields a kinetic energy of 4 x [...] ergs or a diameter of 1.6 \u00b1 0.3 for the R fragment. Using the G and R impacts to calibrate Weaver et al.'s [1994] detailed photometric determination of the relative diameters of pre-Jovian impact projectiles in the SL9 chain from HST images yields revised values for the 11 largest fragment diameters. Adding the inferred masses from the plume and radiative calibration of Weaver et al.'s catalog yields an SL9 progenitor equivalent (ice) diameter of 4.1 \u00b1 0.6 and 3.5 \u00b1 0.5 km, respectively. This compares to a 7.7 km diameter progenitor inferred by Weaver et al. and a 2 km diameter progenitor obtained using tidal break-up modeling by Scotti and Melosh [1993] and Asphaug and Benz [1994]. We also examined the inferred position of Jupiter derived H2O impact ejecta from the hypothetical (3 to 5 bar) cloud deck and find that although Jupiter derived impact induced plumes having radius of > [...] km and achieved altitudes of 3 x [...] km, the contained Jovian water ejecta is restricted to the interior of the plume and would be masked from Earth observations by the surprisingly opaque plume. On the basis of Bjorker et al.'s [1994] observations of thirty minutes' water emissions from the plumes of the impact sites of fragments G and K Field and Ferraro's [1994] analysis which indicates that even if SL9 fragments are 50% porous, upon being subjected to entry, the ram pressures, upon reaching the level of 1 bar, squeeze-out virtually all porosity and the present, self-consistent kinemetric and radiative coupling of the SL9 impactor, we conclude: 1) there is no evidence that SL9 was anything but a comet as first suggested by the coma, 2) penetration of 2 km diameter fragments occurs to depths of ~ 300 km, even for a porous comet, as calculated by the present SPH method and by other groups using finite difference methods [Zahnle and MacLow, 1993, Boslough et al., 1994].\r\n\r\nPaper IV:\r\n\r\nA paraboloidal bow shock model was developed in order to estimate the surface distribution of gas shock-induced modifications surrounding venusian impact craters. We applied two-dimensional oblique shock dynamics to describe a three-dimensional paraboloidal - shaped bow shock impinging upon an assumed incompressible venusian surface. The effects of the hypersonic atmospheric shock acting on the venusian surface are considered in terms of maximum gas pressure, the density, the particle velocity, and the temperature, for varying angles and velocities of impact. The maximum boulder size that can be saltated by the shockwave induced gas flow and the degree of mutual collision of the surface materials are also considered. The present calculations quantitatively predict the areal extent of the gas shock perturbed surface for normal and oblique impact as a function of impact angle and velocity, and radii of impactors. For a 1 km-radius stony meteorite impacting at 20 km/s, the radius of the disturbed area extends ~ 10-17 times the 3-5 km crater radius. The perturbed surface affects the surface radar properties and the present results can provide an explanation of the wide \"dark/bright halos\" surrounding some of venusian impact craters observed via Magellan imagery. For example, a ~ 50 km radius bright halo within a ~ 20 km dark halo is observed around the 3.1 km-radius crater located at 16.5\u00b0 north latitude and 334.4\u00b0 longitude. The average value of the radar backscatter cross section of the ~ 20 km radius dark halo indicates that ~ 50 cm-thick layer of porous lithologic material is superimposed upon the undisturbed surface. The occurrence of the bright halo indicates that the surface roughness in this region is ~ 30% greater than that of the surrounding original surface. The present model can relate the observed crater halo radii to the impact parameters, such as impact velocity and angle, and the impactor radius.\r\n", "doi": "10.7907/2jtf-x438", "publication_date": "1995", "thesis_type": "phd", "thesis_year": "1995" }, { "id": "thesis:7397", "collection": "thesis", "collection_id": "7397", "cite_using_url": "https://resolver.caltech.edu/CaltechTHESIS:01142013-104747074", "primary_object_url": { "basename": "Rowan_lr_1993.pdf", "content": "final", "filesize": 38194629, "license": "other", "mime_type": "application/pdf", "url": "/7397/1/Rowan_lr_1993.pdf", "version": "v3.0.0" }, "type": "thesis", "title": "Equation of state of molten mid-ocean ridge basalt. Structure of Kilauea volcano.", "author": [ { "family_name": "Rowan", "given_name": "Linda Rose", "clpid": "Rowan-L-R" } ], "thesis_advisor": [ { "family_name": "Ahrens", "given_name": "Thomas J.", "clpid": "Ahrens-T-J" }, { "family_name": "Stolper", "given_name": "Edward M.", "clpid": "Stolper-E-M" }, { "family_name": "Clayton", "given_name": "Robert W.", "clpid": "Clayton-R-W" } ], "thesis_committee": [ { "family_name": "Ahrens", "given_name": "Thomas J.", "clpid": "Ahrens-T-J" }, { "family_name": "Stolper", "given_name": "Edward M.", "clpid": "Stolper-E-M" }, { "family_name": "Clayton", "given_name": "Robert W.", "clpid": "Clayton-R-W" }, { "family_name": "Rossman", "given_name": "George Robert", "clpid": "Rossman-G-R" }, { "family_name": "Westphal", "given_name": "James A.", "clpid": "Westphal-J-A" } ], "local_group": [ { "literal": "div_gps" } ], "abstract": "Basalts are the most ubiquitous rocks erupting at the earth's surface at the present\r\ntime and they provide an important probe of the subsurface processes occurring\r\nwithin planetary interiors. Recent advances in both mineral physics and seismic\r\nanalysis have allowed me to undertake two independent studies related to the genesis\r\nand eruption of basaltic magmas. Chapters 1 and 2 are part of an experimental study\r\nconducted in the shock wave laboratory on the equation of state of molten mid-ocean\r\nridge basalt and the chemical interactions of the shocked liquid with its Mo container.\r\nMy advisors for this project were Thomas Ahrens and Edward Stolper. Chapter 3 is a\r\ntravel time tomography study of the three-dimensional structure of Kilauea Volcano,\r\nHawaii in collaboration with Robert Clayton. Chapter 3 is currently in press in the\r\nJournal of Geophysical Research.
\r\n\r\nThe EOS of molten MORB to 20 GPa was accomplished using the innovative silicate\r\nliquid shock wave measurement technique on the 40 mm propellant gun developed\r\nby Rigden [1986] and Miller [1990). This technique has been used to determine\r\nthe EOS for four synthetic melts and this thesis applies the technique to a natural\r\nmelt, a MORB dredged from the Juan de Fuca ridge. The resulting EOS indicates\r\nthat the MORB liquid is very compressible and therefore has a low bulk modulus of\r\n11.7 GPa. These results are consistent with low pressure static compression experiments\r\non similar basalts, but are not consistent with the results of ultrasonic interferometry.\r\nThe compressible nature of the MORB liquid is related to its composition\r\nand this may be expressed best by comparing the MORB Hugoniot to the Hugoniot\r\ndetermined for An_(36)Di_(64) and komatiite. The MORB and An_(36)Di_(64) Hugoniots show\r\nsignificant increase in density at low pressure followed by a stiffening at high pressures\r\nwhere the liquid Hugoniot approaches its respective dense oxide high pressure\r\ncomposition. This may be related to gradual coordination changes from four-fold to\r\nsix-fold for the Si^(+4) and Al^(+3) which are essentially complete at the high pressure\r\nwhere the curve stiffens. The MORB is much more compressible than the komatiite\r\nand overtakes the komatiite in density at a low pressure of 2.5 GPa. This is a\r\ncompositional effect caused by the enrichment of the MORB in Al_2O_3 and SiO_2 and\r\ndepletion in MgO compared to komatiite. The compressible nature of the MORB\r\nallows it to become denser than the surrounding mantle near the base of the low\r\nvelocity zone and therefore it is unlikely that MORB can be derived from very deep\r\nin the earth's upper mantle.
\r\n\r\nFor most shock wave experiments, the sample is not recovered and nothing can\r\nbe determined about its structure or composition due to the passage of the shock\r\nwave. In a few of my EOS experiments on molten MORB, however, the shocked sample\r\nwas recovered and could be studied in detail. Observations of impact-induced\r\ninteractions between the silicate liquid and its Mo container provide insight into\r\nplanetary impact and differentiation processes involving metal-silicate partitioning.\r\nThe shocked liquids showed extreme reduction and with increasing pressure the FeO\r\ncontent of the initial melt was reduced to almost nothing by reaction with the Mo.\r\nThese reactions produced metallic particles enriched in Mo, Fe and Si. These particles\r\nshow a similar texture as those found at impact sites on the earth and moon\r\nand provide clues to the impact origin of metallic particles.
\r\n\r\nA travel time tomography study of local P wave data from Kilauea Volcano,\r\nHawaii, was undertaken to determine the lateral heterogeneities produced by its intricate\r\nmagmatic and tectonic environment. The technique proved to be a powerful\r\nprobe of the volcano's intrusive plumbing because the presence of a dense seismic\r\narray and many local earthquakes allowed for excellent coverage of complex subsurface\r\nfeatures. Analysis and interpretation of the tomographic images leads to the\r\nfollowing inferrred model. The main shallow magma reservoir is delineated by a\r\nslow anomaly centered 2 km southeast of Halemaumau caldera. There is a distinct\r\nhigh velocity region centered northwest of the summit from 0 to 2 km depth that\r\nmay represent a dense wall and/or cap of intrusive rock that acts as a barrier or\r\ncontainment structure for the northern part of the reservoir. We suggest that the\r\nshallow reservoir is a narrow, compartmentalized region of sills and dikes because of\r\nthe closely spaced high and low velocity anomalies near the summit. The rift zones\r\nof Kilauea are imaged as major, high velocity entities, widening to the south with\r\ndepth until 6 km. These fast anomalies may be related to the sheeted dike complexes\r\nalong the rifts. On a finer scale, magma pockets centered at 0-2 km depth have been\r\ninferred beneath Makaopuhi, Mauna Ulu and Puu Oo, along the east rift zone. The\r\nHilina and Kaoiki fault zones, are imaged as slow features at shallow depths (less than 6\r\nkm), related to their tensional structures that produce the open fractures and cracks\r\nin the basaltic edifice. The Koae fault system is imaged as a slightly fast shallow\r\nstructure (less than 6 km) possibly related to intrusive diking from the adjacent rift zones.\r\nContinued inversions with the immense amount of seismic data collected for Hawaiian\r\nevents will allow the detailed development of a three-dimensional velocity model\r\nfor Kilauea. Such a model will be extremely useful to seismologists and petrologists\r\nalike for understanding the tectonic growth and magmatic evolution of this dynamic\r\nshield volcano.
\r\n\r\n", "doi": "10.7907/J6D0-GV56", "publication_date": "1993", "thesis_type": "phd", "thesis_year": "1993" }, { "id": "thesis:6701", "collection": "thesis", "collection_id": "6701", "cite_using_url": "https://resolver.caltech.edu/CaltechTHESIS:10032011-142032525", "primary_object_url": { "basename": "Zhao_ls_1992.pdf", "content": "final", "filesize": 25837928, "license": "other", "mime_type": "application/pdf", "url": "/6701/1/Zhao_ls_1992.pdf", "version": "v4.0.0" }, "type": "thesis", "title": "Seismic waveform modeling of regional phases and wavefields from an off-center explosion", "author": [ { "family_name": "Zhao", "given_name": "Lian-She", "clpid": "Zhao-Lian-She" } ], "thesis_advisor": [ { "family_name": "Ahrens", "given_name": "Thomas J.", "clpid": "Ahrens-T-J" } ], "thesis_committee": [ { "family_name": "Helmberger", "given_name": "Donald V.", "clpid": "Helmberger-D-V" }, { "family_name": "Harkrider", "given_name": "David G.", "clpid": "Harkrider-D-G" } ], "local_group": [ { "literal": "div_gps" } ], "abstract": "This thesis consists of two parts. Part I (Chapter 1) is a theoretical study of wavefields radiated from an off-center explosion in an embedded solid sphere, while Part II (Chapters 2, 3 and 4) deals with seismic waveform modeling of regional phases.
\r\n\r\nChapter 1 investigates the effects of asymmetric source regions on the excitation of seismic body waves of\r\nexplosions. We give an analytic formulation for determining the wave fields from an off-center explosion in an embedded\r\nsolid sphere in an elastic whole-space. It is shown that SH wave amplitudes can be ten times larger than the P wave amplitudes if the source is properly located.
\r\n\r\nChapter 2 discusses the strategies of modeling regional Broadband recordings with the application to Harvard station recordings of November 25, 1988 Saguenay earthquake. We demonstrate that the longer period motions\r\n(WWSSN) preceding the direct S arrival can be modeled reasonably well with a crustal model consisting of a layer over half-space. While a generic model assumed to be the same everywhere can be used to model the P_nl waveshape, HeImberger and Engen (1980), some adjustments are required\r\nto fit the absolute timing of P_n and S_n for specific paths. In particular, path from the Saguenay event to Harvard is slower than expected for a shield environment and is associated with the roots of the Appalachian\r\nThrust Belt.
\r\n\r\nChapters 3 and 4 concern the velocity structures of the Tibetan Plateau and surrounding regions. A fast shield-like velocity model, TIP, is proposed for the structure beneath the Tibetan Plateau as suggested by modeling SH waveforms in Chapter 3. The crustal model is determined by modeling\r\nLove waveforms and P_nl waveforms in the time domain while the upper mantle is constrained with S and SS waveform data. Relocations of Tibetan earthquakes are discussed in Chapter 4, assuming model TIP at source\r\nregion and model JB at receiver regions. A detailed\r\ninvestigation of 100 such events yields a distinctly different picture from a random distribution of events down to a depth of about 50 kms given by ISC. Waveform modeling\r\nof depth phases such as _pP indicates that only three or four events from of this population is actually deeper than 25 km. These few events occur near the edges of the Plateau where active subduction is occurring as suggested by the thrust-like nature of their mechanisms. The events occurred earlier than indicated by the ISC by 3 seconds on average which leads to about a 1.5% and 0.5% over estimation of P_n and S_n velocities respectively applying ISC tables and standard flat-layered models.
\r\n", "doi": "10.7907/thvj-vd30", "publication_date": "1992", "thesis_type": "phd", "thesis_year": "1992" }, { "id": "thesis:1847", "collection": "thesis", "collection_id": "1847", "cite_using_url": "https://resolver.caltech.edu/CaltechETD:etd-05172007-104609", "primary_object_url": { "basename": "Duffy_ts_1992.pdf", "content": "final", "filesize": 8927678, "license": "other", "mime_type": "application/pdf", "url": "/1847/1/Duffy_ts_1992.pdf", "version": "v3.0.0" }, "type": "thesis", "title": "Elastic Properties of Metals and Minerals under Shock Compression", "author": [ { "family_name": "Duffy", "given_name": "Thomas Sheehan", "clpid": "Duffy-Thomas-Sheehan" } ], "thesis_advisor": [ { "family_name": "Ahrens", "given_name": "Thomas J.", "clpid": "Ahrens-T-J" } ], "thesis_committee": [ { "family_name": "Unknown", "given_name": "Unknown" } ], "local_group": [ { "literal": "div_gps" } ], "abstract": "NOTE: Text or symbols not renderable in plain ASCII are indicated by [...]. Abstract is included in .pdf document.\r\n\r\nComparison of laboratory elasticity data with seismic measurements of the Earth provides a means to understand the deep interior. The effect of pressure and temperature on elastic properties must be well understood for meaningful comparisons. In this work, elastic wave velocities have been measured under shock compression to 80 GPa in an Fe-Cr-Ni alloy, to 27 GPa in polycrystalline MgO, and to 81 GPa in molybdenum preheated to 1400\u00b0C. These measurements were made by recording particle velocity histories at a sample surface using the method of velocity interferometry. In addition to elastic properties, these experiments provide information on the constitutive and equation of state (EOS) properties of the sample as well as the unloading adiabats.
\r\n\r\nCompressional and bulk wave velocities in Fe-Cr-Ni alloy are consistent with third-order finite strain theory and ultrasonic data. Thermal effects on the wave velocities are less than 2% at 80 GPa. Second pressure derivatives of velocity were constrained along the Hugoniot to be: (\u22022CL/\u2202P2)H = -0.16 (0.06) GPa-1 and (\u22022KS/\u2202P2)H = -0.17 (0.08) GPa-1. The measured wave profiles can be successfully reproduced by numerical simulations utilizing elastic-plastic theory modified by a Bauschinger effect and stress relaxation. Material strength was found to increase by a factor of at least 5 up to 80 GPa and to be 2-3% of the total stress.
\r\n\r\nCompressional and bulk velocities in Fe-Cr-Ni define linear velocity-density trends and can be modeled by averaging properties of Fe, Cr, and Ni. The effect of alloying ~4 wt.% Ni with Fe would change both VP and VB by less than 1% under core conditions. Compressional velocities in Fe-Ni are compatible with inner core values when corrected for thermal effects. Shear velocities in Fe, determined from a combination of VP and VB data, are ~3.6 km/s at P=150-200 GPa. Low values are most likely caused by a weak pressure dependence of the rigidity and imply that partial melting is not required in the inner core.
\r\n\r\nWave profile and EOS measurements in polycrystalline MgO define its EOS: US = 6.77(0.08) + 1.27(0.04)up. Compressional sound velocities to 27 GPa yield the longitudinal modulus and its pressure derivative: CLo = KoS + 4/3G = 335 \u00b1 1 GPa and C'Lo = 7.4 \u00b1 0.2, which are in good agreement with ultrasonic determinations. The unloading wave profiles can be modeled using a modified elastic-plastic constitutive response originally developed for metals. Thermal expansivities in MgO have been determined to be 12 \u00b1 4 x 10-6 K-1 at P=174-200 GPa and T=3100-3600 K from shock temperature and EOS data. These results imply that the lower mantle is enriched in Si and/or Fe relative to the upper mantle.
\r\n\r\nWave profiles in molybdenum at 1400\u00b0C are the first wave profile determinations at significantly high initial temperature. The EOS determined from these measurements agrees well with previous data. The compressive yield strength of Mo is 0.79-0.94 GPa at 1400\u00b0C, and the HEL stress is 1.5-1.7 GPa. The temperature coefficient of compressional velocity, (\u2202Vp/\u2202T)p, is found to vary from -0.35(0.13) m/s/K at 12 GPa to -0.18(0.14) m/s/K at 81 GPa and compares with an ambient pressure value of -0.26 m/s/K. It is inferred that (\u2202Vp/\u2202T)p decreases with pressure, and data for Mo are shown to be consistent with trends defined by other metals.
", "doi": "10.7907/m63k-zj81", "publication_date": "1992", "thesis_type": "phd", "thesis_year": "1992" }, { "id": "thesis:6695", "collection": "thesis", "collection_id": "6695", "cite_using_url": "https://resolver.caltech.edu/CaltechTHESIS:09282011-151807180", "primary_object_url": { "basename": "Smither_cl_1992.pdf", "content": "final", "filesize": 40708385, "license": "other", "mime_type": "application/pdf", "url": "/6695/1/Smither_cl_1992.pdf", "version": "v4.0.0" }, "type": "thesis", "title": "Stress-relief displacements induced by drilling and three-dimensional modeling of planetary impacts", "author": [ { "family_name": "Smither", "given_name": "Catherine Louise", "clpid": "Smither-C-L" } ], "thesis_advisor": [ { "family_name": "Ahrens", "given_name": "Thomas J.", "clpid": "Ahrens-T-J" } ], "thesis_committee": [ { "family_name": "Clayton", "given_name": "Robert W.", "clpid": "Clayton-R-W" }, { "family_name": "Harkrider", "given_name": "David G.", "clpid": "Harkrider-D-G" } ], "local_group": [ { "literal": "div_gps" } ], "abstract": "The holographic stressmeter uses double-exposure holographic interferometry to record the displacements induced by the drilling of a sidehole into the borehole wall. The local stresses, which are the result of the far-field stresses, concentrated at the borehole, cause deformation of the surface of the borehole wall near the sidehole. To interpret the data, it is essential to understand this deformation. The initial model used a thin infinite elastic plate subjected to plane stress at infinity. Two-dimensional finite element analysis showed that the displacement depends on the depth of the\r\nsidehole. We then developed a new model for the analysis of stress-relief displacements. For holes with a depth/diameter ratio greater than unity, the simple plane\r\nstress elastic plate solution breaks down. The revised model, which gives a more exact solution of displacements near the hole, does a better job of fitting the observed\r\ndata.
\r\n\r\nA three-dimensional smoothed particle hydrodynamics code was used to model normal and oblique impacts of silicate projectiles on planetary bodies. The energy of the system is partitioned after impact into internal and kinetic energy of the both bodies. These simulations show that, unlike the case of impacts onto a half-space, up to 70% of the initial energy remains in the kinetic energy of the impactor, as parts of it travel past the main planet and escape the system. More oblique impacts retain more kinetic energy in the impactor: 6 to 75% versus 4 to 30% for the normal impacts. Higher velocity collisions also show this trend, as the kinetic energy of the impactor is 2 to 50% of the total for 5 km/s impacts, and 13 to 75% for 20 km/s\r\nimpacts. Impacts at 20 km/s with an impactor 60% the size of the target completely melted both targets. Three to 4 times more vaporization of the target material occurred on the larger targets. The amount of target material ejected at velocities greater than the escape velocity is found to be higher than that predicted by studies of impact onto a half-space, and slightly less than the amount predicted by theories of catastrophic breakup of asteroids.
\r\n", "doi": "10.7907/49fn-dg49", "publication_date": "1992", "thesis_type": "phd", "thesis_year": "1992" }, { "id": "thesis:8806", "collection": "thesis", "collection_id": "8806", "cite_using_url": "https://resolver.caltech.edu/CaltechTHESIS:03272015-161328191", "primary_object_url": { "basename": "Miller_gh_1990.pdf", "content": "final", "filesize": 7846556, "license": "other", "mime_type": "application/pdf", "url": "/8806/1/Miller_gh_1990.pdf", "version": "v3.0.0" }, "type": "thesis", "title": "The Equation of State and Petrogenesis of Komatiite", "author": [ { "family_name": "Miller", "given_name": "Gregory Hale", "clpid": "Miller-G-H" } ], "thesis_advisor": [ { "family_name": "Ahrens", "given_name": "Thomas J.", "clpid": "Ahrens-T-J" }, { "family_name": "Stolper", "given_name": "Edward M.", "clpid": "Stolper-E-M" }, { "family_name": "Rossman", "given_name": "George Robert", "clpid": "Rossman-G-R" }, { "family_name": "Goddard", "given_name": "William A., III", "clpid": "Goddard-W-A-III" } ], "thesis_committee": [ { "family_name": "Unknown", "given_name": "Unknown" } ], "local_group": [ { "literal": "div_gps" } ], "abstract": "(1) Equation of State of Komatiite
\r\n\r\nThe equation of state (EOS) of a molten komatiite (27 wt% MgO) was detennined in the 5 to\r\n36 GPa pressure range via shock wave compression from 1550\u00b0C and 0 bar. Shock wave velocity,\r\nUS, and particle velocity, UP, in km/s follow the linear relationship US = 3.13(\u00b10.03) + 1.47(\u00b10.03)\r\nUP. Based on a calculated density at 1550\u00b0C, 0 bar of 2.745\u00b10.005 glee, this US-UP relationship\r\ngives the isentropic bulk modulus KS = 27.0 \u00b1 0.6 GPa, and its first and second isentropic pressure\r\nderivatives, K'S = 4.9 \u00b1 0.1 and K\"S = -0.109 \u00b1 0.003 GPa-1.
\r\n\r\nThe calculated liquidus compression curve agrees within error with the static compression\r\nresults of Agee and Walker [1988a] to 6 GPa. We detennine that olivine (FO94) will be neutrally\r\nbuoyant in komatiitic melt of the composition we studied near 8.2 GPa. Clinopyroxene would also\r\nbe neutrally buoyant near this pressure. Liquidus garnet-majorite may be less dense than this komatiitic\r\nliquid in the 20-24 GPa interval, however pyropic-garnet and perovskite phases are denser than\r\nthis komatiitic liquid in their respective liquidus pressure intervals to 36 GPa. Liquidus perovskite\r\nmay be neutrally buoyant near 70 GPa.
\r\n\r\nAt 40 GPa, the density of shock-compressed molten komatiite would be approximately equal\r\nto the calculated density of an equivalent mixture of dense solid oxide components. This observation\r\nsupports the model of Rigden et al. [1989] for compressibilities of liquid oxide components.\r\nUsing their theoretical EOS for liquid forsterite and fayalite, we calculate the densities of a spectrum\r\nof melts from basaltic through peridotitic that are related to the experimentally studied komatiitic\r\nliquid by addition or subtraction of olivine. At low pressure, olivine fractionation lowers the density\r\nof basic magmas, but above 14 GPa this trend is reversed. All of these basic to ultrabasic liquids\r\nare predicted to have similar densities at 14 GPa, and this density is approximately equal to the bulk\r\n(PREM) mantle. This suggests that melts derived from a peridotitic mantle may be inhibited from\r\nascending from depths greater than 400 km.
\r\n\r\nThe EOS of ultrabasic magmas was used to model adiabatic melting in a peridotitic mantle.\r\nIf komatiites are formed by >15% partial melting of a peridotitic mantle, then komatiites generated\r\nby adiabatic melting come from source regions in the lower transition zone (\u2248500-670 km) or the\r\nlower mantle (>670 km). The great depth of incipient melting implied by this model, and the melt\r\ndensity constraint mentioned above, suggest that komatiitic volcanism may be gravitationally hindered.\r\nAlthough komatiitic magmas are thought to separate from their coexisting crystals at a temperature\r\n=200\u00b0C greater than that for modern MORBs, their ultimate sources are predicted to be\r\ndiapirs that, if adiabatically decompressed from initially solid mantle, were more than 700\u00b0C hotter\r\nthan the sources of MORBs and derived from great depth.
\r\n\r\nWe considered the evolution of an initially molten mantle, i.e., a magma ocean. Our model\r\nconsiders the thermal structure of the magma ocean, density constraints on crystal segregation, and\r\napproximate phase relationships for a nominally chondritic mantle. Crystallization will begin at the\r\ncore-mantle boundary. Perovskite buoyancy at > 70 GPa may lead to a compositionally stratified\r\nlower mantle with iron-enriched mangesiowiistite content increasing with depth. The upper mantle\r\nmay be depleted in perovskite components. Olivine neutral buoyancy may lead to the formation of\r\na dunite septum in the upper mantle, partitioning the ocean into upper and lower reservoirs, but this\r\nseptum must be permeable.
\r\n\r\n(2) Viscosity Measurement with Shock Waves
\r\n\r\nWe have examined in detail the analytical method for measuring shear viscosity from the\r\ndecay of perturbations on a corrugated shock front The relevance of initial conditions, finite shock\r\namplitude, bulk viscosity, and the sensitivity of the measurements to the shock boundary conditions\r\nare discussed. The validity of the viscous perturbation approach is examined by numerically solving\r\nthe second-order Navier-Stokes equations. These numerical experiments indicate that shock instabilities\r\nmay occur even when the Kontorovich-D'yakov stability criteria are satisfied. The experimental\r\nresults for water at 15 GPa are discussed, and it is suggested that the large effective viscosity\r\ndetermined by this method may reflect the existence of ice VII on the Rayleigh path of the\r\nHugoniot This interpretation reconciles the experimental results with estimates and measurements\r\nobtained by other means, and is consistent with the relationship of the Hugoniot with the phase\r\ndiagram for water. Sound waves are generated at 4.8 MHz at in the water experiments at 15 GPa.\r\nThe existence of anelastic absorption modes near this frequency would also lead to large effective\r\nviscosity estimates.
\r\n\r\n(3) Equation of State of Molybdenum at 1400\u00b0C
\r\n\r\nShock compression data to 96 GPa for pure molybdenum, initially heated to 1400\u00b0C, are\r\npresented. Finite strain analysis of the data gives a bulk modulus at 1400\u00b0C, K'S. of 244\u00b12 GPa and\r\nits pressure derivative, K'OS of 4. A fit of shock velocity to particle velocity gives the coefficients of\r\nUS = CO+S UP to be CO = 4.77\u00b10.06 km/s and S = 1.43\u00b10.05. From the zero pressure sound speed, CO, a bulk modulus of 232\u00b16 GPa is calculated that is consistent with extrapolation of ultrasonic elasticity measurements. The temperature derivative of the bulk modulus at zero pressure, \u03b8KOS\u03b8T|P, is\r\napproximately -0.012 GPa/K. A thermodynamic model is used to show that the thermodynamic\r\nGr\u00fcneisen parameter is proportional to the density and independent of temperature. The Mie-Gr\u00fcneisen\r\nequation of state adequately describes the high temperature behavior of molybdenum\r\nunder the present range of shock loading conditions.
", "doi": "10.7907/r0pt-2227", "publication_date": "1990", "thesis_type": "phd", "thesis_year": "1990" }, { "id": "thesis:8781", "collection": "thesis", "collection_id": "8781", "cite_using_url": "https://resolver.caltech.edu/CaltechTHESIS:03182015-092735941", "primary_object_url": { "basename": "Anderson_ww_1990.pdf", "content": "final", "filesize": 31311830, "license": "other", "mime_type": "application/pdf", "url": "/8781/1/Anderson_ww_1990.pdf", "version": "v3.0.0" }, "type": "thesis", "title": "High Pressure States in Condensed Matter: I. High Pressure Behavior of the Iron-Sulfur System with Applications to the Earth's Core. II. Empirical Equation of State for Organic Compounds at High Pressures", "author": [ { "family_name": "Anderson", "given_name": "William Wyatt", "clpid": "Anderson-William-Wyatt" } ], "thesis_advisor": [ { "family_name": "Ahrens", "given_name": "Thomas J.", "clpid": "Ahrens-T-J" }, { "family_name": "Stevenson", "given_name": "David John", "orcid": "0000-0001-9432-7159", "clpid": "Stevenson-D-J" } ], "thesis_committee": [ { "family_name": "Stevenson", "given_name": "David John", "orcid": "0000-0001-9432-7159", "clpid": "Stevenson-D-J" }, { "family_name": "Ahrens", "given_name": "Thomas J.", "clpid": "Ahrens-T-J" }, { "family_name": "Anderson", "given_name": "Donald L.", "clpid": "Anderson-D-L" }, { "family_name": "Burnett", "given_name": "Donald S.", "clpid": "Burnett-D-S" }, { "family_name": "Albee", "given_name": "Arden Leroy", "clpid": "Albee-A-L" } ], "local_group": [ { "literal": "div_gps" } ], "abstract": "Part I:
\r\n\r\nThe earth's core is generally accepted to be composed primarily of iron, with an admixture of other elements. Because the outer core is observed not to transmit shear waves at seismic frequencies, it is known to be liquid or primarily liquid. A new equation of state is presented for liquid iron, in the form of parameters for the 4th order Birch-Murnaghan and Mie-Gr\u00fcneisen equations of state. The parameters were constrained by a set of values for numerous properties compiled from the literature. A detailed theoretical model is used to constrain the P-T behavior of the heat capacity, based on recent advances in the understanding of the interatomic potentials for transition metals. At the reference pressure of 105 Pa and temperature of 1811 K (the normal melting point of Fe), the parameters are: \u03c1 = 7037 kg/m3, KS0 = 110 GPa, KS' = 4.53, KS\" = -.0337 GPa-1, and \u03b3 = 2.8, with \u03b3 \u221d \u03c1-1.17. Comparison of the properties predicted by this model with the earth model PREM indicates that the outer core is 8 to 10 % less dense than pure liquid Fe at the same conditions. The inner core is also found to be 3 to 5% less dense than pure liquid Fe, supporting the idea of a partially molten inner core. The density deficit of the outer core implies that the elements dissolved in the liquid Fe are predominantly of lower atomic weight than Fe. Of the candidate light elements favored by researchers, only sulfur readily dissolves into Fe at low pressure, which means that this element was almost certainly concentrated in the core at early times. New melting data are presented for FeS and FeS2 which indicate that the FeS2 is the S-hearing liquidus solid phase at inner core pressures. Consideration of the requirement that the inner core boundary be observable by seismological means and the freezing behavior of solutions leads to the possibility that the outer core may contain a significant fraction of solid material. It is found that convection in the outer core is not hindered if the solid particles are entrained in the fluid flow. This model for a core of Fe and S admits temperatures in the range 3450K to 4200K at the top of the core. An all liquid Fe-S outer core would require a temperature of about 4900 K at the top of the core.
\r\n\r\nPart II.
\r\n\r\nThe abundance of uses for organic compounds in the modern world results in many applications in which these materials are subjected to high pressures. This leads to the desire to be able to describe the behavior of these materials under such conditions. Unfortunately, the number of compounds is much greater than the number of experimental data available for many of the important properties. In the past, one approach that has worked well is the calculation of appropriate properties by summing the contributions from the organic functional groups making up molecules of the compounds in question. A new set of group contributions for the molar volume, volume thermal expansivity, heat capacity, and the Rao function is presented for functional groups containing C, H, and O. This set is, in most cases, limited in application to low molecular liquids. A new technique for the calculation of the pressure derivative of the bulk modulus is also presented. Comparison with data indicates that the presented technique works very well for most low molecular hydrocarbon liquids and somewhat less well for oxygen-bearing compounds. A similar comparison of previous results for polymers indicates that the existing tabulations of group contributions for this class of materials is in need of revision. There is also evidence that the Rao function contributions for polymers and low molecular compounds are somewhat different.
", "doi": "10.7907/7dsx-ed25", "publication_date": "1990", "thesis_type": "phd", "thesis_year": "1990" }, { "id": "thesis:7994", "collection": "thesis", "collection_id": "7994", "cite_using_url": "https://resolver.caltech.edu/CaltechTHESIS:10172013-085025235", "primary_object_url": { "basename": "Polanskey_ca_1989.pdf", "content": "final", "filesize": 19873642, "license": "other", "mime_type": "application/pdf", "url": "/7994/1/Polanskey_ca_1989.pdf", "version": "v3.0.0" }, "type": "thesis", "title": "I. Impact Spallation Experiments: Fracture Patterns and Spall Velocities. II. Craters in Carbonate Rocks: An Electron Paramagnetic Resonance Analysis of Shock Damage", "author": [ { "family_name": "Polanskey", "given_name": "Carol Ann", "clpid": "Polanskey-Carol-Ann" } ], "thesis_advisor": [ { "family_name": "Ahrens", "given_name": "Thomas J.", "clpid": "Ahrens-T-J" }, { "family_name": "Burnett", "given_name": "Donald S.", "orcid": "0000-0001-9521-8675", "clpid": "Burnett-D-S" } ], "thesis_committee": [ { "family_name": "Stevenson", "given_name": "David John", "orcid": "0000-0001-9432-7159", "clpid": "Stevenson-D-J" }, { "family_name": "Ahrens", "given_name": "Thomas J.", "clpid": "Ahrens-T-J" }, { "family_name": "Burnett", "given_name": "Donald S.", "orcid": "0000-0001-9521-8675", "clpid": "Burnett-D-S" }, { "family_name": "Ingersoll", "given_name": "Andrew P.", "orcid": "0000-0002-2035-9198", "clpid": "Ingersoll-A-P" }, { "family_name": "Westphal", "given_name": "James A.", "clpid": "Westphal-J-A" } ], "local_group": [ { "literal": "div_gps" } ], "abstract": "This work is divided into two independent papers.
\r\n\r\nPAPER 1.
\r\n\r\nSpall velocities were measured for nine experimental impacts into San Marcos gabbro targets. Impact velocities ranged from 1 to 6.5 km/sec. Projectiles were iron, aluminum, lead, and basalt of varying sizes. The projectile masses ranged from a 4 g lead bullet to a 0.04 g aluminum sphere. The velocities of fragments were measured from high-speed films taken of the events. The maximum spall velocity observed was 30 m/sec, or 0.56 percent of the 5.4 km/sec impact velocity. The measured velocities were compared to the spall velocities predicted by the spallation model of Melosh (1984). The compatibility between the spallation model for large planetary impacts and the results of these small scale experiments are considered in detail.
\r\n\r\nThe targets were also bisected to observe the pattern of internal fractures. A series of fractures were observed, whose location coincided with the boundary between rock subjected to the peak shock compression and a theoretical \"near surface zone\" predicted by the spallation model. Thus, between this boundary and the free surface, the target material should receive reduced levels of compressive stress as compared to the more highly shocked region below.
\r\n\r\nPAPER 2.
\r\n\r\nCarbonate samples from the nuclear explosion crater, OAK, and a terrestrial impact crater, Meteor Crater, were analyzed for shock damage using electron paramagnetic resonance, EPR. The first series of samples for OAK Crater were obtained from six boreholes within the crater, and the second series were ejecta samples recovered from the crater floor. The degree of shock damage in the carbonate material was assessed by comparing the sample spectra to spectra of Solenhofen limestone, which had been shocked to known pressures.
\r\n\r\nThe results of the OAK borehole analysis have identified a thin zone of highly shocked carbonate material underneath the crater floor. This zone has a maximum depth of approximately 200 ft below sea floor at the ground zero borehole and decreases in depth towards the crater rim. A layer of highly shocked material is also found on the surface in the vicinity of the reference bolehole, located outside the crater. This material could represent a fallout layer. The ejecta samples have experienced a range of shock pressures.
\r\n\r\nIt was also demonstrated that the EPR technique is feasible for the study of terrestrial impact craters formed in carbonate bedrock. The results for the Meteor Crater analysis suggest a slight degree of shock damage present in the \u03b2 member of the Kaibab Formation exposed in the crater walls.
", "doi": "10.7907/3rgr-f070", "publication_date": "1989", "thesis_type": "phd", "thesis_year": "1989" }, { "id": "thesis:16215", "collection": "thesis", "collection_id": "16215", "cite_using_url": "https://resolver.caltech.edu/CaltechTHESIS:10242023-225432851", "primary_object_url": { "basename": "Schmitt_DR_1987.pdf", "content": "final", "filesize": 52360957, "license": "other", "mime_type": "application/pdf", "url": "/16215/1/Schmitt_DR_1987.pdf", "version": "v2.0.0" }, "type": "thesis", "title": "I. Applications of Double-Exposure Holography to the Measurement of In Situ Stress and the Elastic Moduli of Rock from Boreholes. II. Shock Temperature Measurements in Fused Quartz and Crystalline NaCl to 35 GPa", "author": [ { "family_name": "Schmitt", "given_name": "Douglas Ray", "clpid": "Schmitt-Douglas-Ray" } ], "thesis_advisor": [ { "family_name": "Ahrens", "given_name": "Thomas J.", "clpid": "Ahrens-T-J" }, { "family_name": "Harkrider", "given_name": "David G.", "clpid": "Harkrider-D-G" } ], "thesis_committee": [ { "family_name": "Harkrider", "given_name": "David G.", "clpid": "Harkrider-D-G" }, { "family_name": "Ahrens", "given_name": "Thomas J.", "clpid": "Ahrens-T-J" }, { "family_name": "Clayton", "given_name": "Robert W.", "orcid": "0000-0003-3323-3508", "clpid": "Clayton-R-W" }, { "family_name": "Kanamori", "given_name": "Hiroo", "orcid": "0000-0001-8219-9428", "clpid": "Kanamori-H" }, { "family_name": "Stevenson", "given_name": "David John", "orcid": "0000-0001-9432-7159", "clpid": "Stevenson-D-J" }, { "family_name": "Stolper", "given_name": "Edward M.", "orcid": "0000-0001-8008-8804", "clpid": "Stolper-E-M" } ], "local_group": [ { "literal": "div_gps" } ], "abstract": "Part I.
\r\n\r\nThe application of a new borehole technique using holographic inter\u00ad ferometry to measure the in situ state of stress and the modulus of elasti\u00adcity of rock is discussed. The apparatus exposes two holograms which are taken both before and after micron scale displacements are induced by\r\n(1) drilling a small stress-relieving hole in the wall of a borehole, and (2) applying a normal point force to the borehole wall. Maximum induced displacements are approximately 10 microns; the holograms are sensitive to movements on the order of 0.1 micron. Raw data take the form of a series interference holograms which have dark fringes superimposed on the three dimensional holographic borehole wall image. Synthetic fringe patterns are used to forward model the observed in the present method of data analysis. Calibrations of the normal force method of measuring the elastic moduli is carried out on metals with well defined elastic properties. Typi\u00adcally each test yields elastic (Young's) moduli for brass and aluminum of 100 \u00b1 10 GPa and 70 \u00b1 5 GPa, respectively, which are in close agreement with\r\nstandard tests. Laboratory holographic measurements of the Young's modulus on a sample of keragenaceous dolomitic marlstone (taken from the same mine as which the in situ experiments were conducted) yielded 16.8 \u00b1 2.8 GPa in agreement with the predicted modulus of 17.2 \u00b1 2.0 GPa based upon published density-modulus relationships. Sonic velocity determinations of the dynamic Young's modulus on cores taken from the rock sample give values consistent with the holographic measurements of 13.5 to 19.1 GPa for assumed values of Poisson's ratio of 0.35 to 0.25. The results of field tests in a horizontal borehole in a mine pillar in the Mahogany formation of Garfield County, Colorado, are presented for both experiments. The elastic modulus was found to vary with position in the borehole from 26.9 to 36.0 GPa. The farfield stresses for a borehole station 4 m from the mine pillar free surface were found from analysis of several stress-relief holograms; the determined vertical stress within the mine pillar was -10.2 MPa (compressive) close to the predicted magnitude of -11.2 MPa.
\r\n\r\nPart II.
\r\n\r\nGreybody temperatures and emittances of fused quartz under shock compression between 10 and 30 GPa are determined. Observed radiative temperatures are higher than computed continuum temperatures for shock compressed fused quartz, however; below ~26 GPa observed emittances are < 0.02. This suggests that fused quartz deforms heterogeneously in this shock pressure range, as has been observed in other minerals. Between 10 and 16 GPa, radiative temperatures decrease from 4400 K to 3200 K, above 16 GPa to 30 GPa greybody temperatures of ~3000 K with low emit\u00adtances are observed. The emittances increase with pressure from 0.02 to 0.9. The pressure range from 10 to 16 GPa coincides with the permanent densification region while the 16 to 30 GPa range coincides with the mixed phase region along the fused quartz Hugoniot. The differing radiative behaviors relate to these modes of deformation. Based upon shock recovery experiments and a proposed model of heterogeneous deformation under shock compression, the temperatures associated with low emittances in the mixed phase region probably represents the melting temperature of the high pressure phase. Above 20 GPa to 30 GPa the melting temperature of stishovite would therefore be approximately 3000 K and almost independent of pressure. The effect of pressure on melting relations for the phase system SiO\u2082-Mg\u2082SiO\u2084 are considered together with the proposed melting curve of stishovite and suggest that maximum solidus temperatures within the mantle of ~2370 K at 12.5 GPa and ~2520 K at 20.0 GPa. Using the proposed stishovite melting temperatures (T_m) and reasonable upper mantle temperatures (T), the effective viscosity (which is a function of the homologous temperature (T_m/T)) appears to remain nearly constant from 600 to 200 km depth in the Earth.
\r\n\r\nRadiative color temperatures were measured in single crystal sodium chloride under shock compression parallel to [100] over a pressure range from 20 to 35 GPa. Color temperatures from 2500 to 4500 K and emittances from 0.003 to 0.3 were determined by fitting observed spectra (450 to 850 nm) to the Planck greybody radiation law. These data support a heterogeneous shock deformation model of shocked halite in this pressure range. A 2500 K temperature rise is observed over the Bl-B2 mixed phase region from 25 to 30 GPa. Assuming that shock deformation occurs via yielding in localized planar zones which become melt and the melting temperature at high pressure controls the temperature, we infer that the temperature of the B2 fusion curve from 30 to 35 GPa rises from 3200 to 3300 K. The Bl-B2-liquid triple point is predicted at a temperature of 2250 K and 23.5 GPa.
", "doi": "10.7907/vr7c-6j57", "publication_date": "1987", "thesis_type": "phd", "thesis_year": "1987" }, { "id": "thesis:6712", "collection": "thesis", "collection_id": "6712", "cite_using_url": "https://resolver.caltech.edu/CaltechTHESIS:10132011-090721112", "primary_object_url": { "basename": "Vidale_je_1987.pdf", "content": "final", "filesize": 24558213, "license": "other", "mime_type": "application/pdf", "url": "/6712/1/Vidale_je_1987.pdf", "version": "v4.0.0" }, "type": "thesis", "title": "Application of Two-Dimensional Finite-Difference Wave Simulation to Earthquakes, Earth Structure, and Seismic Hazard", "author": [ { "family_name": "Vidale", "given_name": "John Emilio", "clpid": "Vidale-John-Emilio" } ], "thesis_advisor": [ { "family_name": "Hager", "given_name": "Bradford H.", "clpid": "Hager-B-H" }, { "family_name": "Clayton", "given_name": "Robert W.", "orcid": "0000-0003-3323-3508", "clpid": "Clayton-R-W" }, { "family_name": "Ahrens", "given_name": "Thomas J.", "clpid": "Ahrens-T-J" } ], "thesis_committee": [ { "family_name": "Ahrens", "given_name": "Thomas J.", "clpid": "Ahrens-T-J" }, { "family_name": "Helmberger", "given_name": "Donald V.", "clpid": "Helmberger-D-V" }, { "family_name": "Kanamori", "given_name": "Hiroo", "orcid": "0000-0001-8219-9428", "clpid": "Kanamori-H" }, { "family_name": "Clayton", "given_name": "Robert W.", "orcid": "0000-0003-3323-3508", "clpid": "Clayton-R-W" }, { "family_name": "Jennings", "given_name": "Paul C.", "clpid": "Jennings-P-C" }, { "family_name": "Hager", "given_name": "Bradford H.", "clpid": "Hager-B-H" } ], "local_group": [ { "literal": "div_gps" } ], "abstract": "Although the earth is 3-dimensional (3-D), numerical simulations of wave propagation through laterally heterogeneous media are easier to formulate and more practical to use in 2-D. In this thesis, schemes to model seismic wave propagation through laterally varying structures with 2-D numerical algorithms are developed and applied to earthquake and explosion problems.
\r\n\r\nIn Chapter 1, 2-D source expressions that have the same radiation patterns as their 3-D counterparts are derived which can rectify the following three problems: the use of 2-D simulations generally results in \"line source tails\" on what would be impulsive arrivals in 3-D, 1/\u221aR rather than 1/R amplitude decay for body waves, and no decay rather than 1/\u221aR amplitude decay for surface waves. Because this technique approximately transforms waves from a cartesian 2-D grid to a cylindrically symmetric 3-D world, slightly anisotropic geometrical spreading in 2-D better approximates isotropic spreading in 3-D than simple isotropic spreading in 2-D does. In Section 1.7, a correction to the explosive source expression reduces energy traveling vertically out of the source region, but leaves unchanged the energy traveling laterally out of the source region. In some cases, this correction will significantly improve the results of using a 2-D grid to simulate elastic wave propagation from an explosive point source.
\r\n\r\nIn Chapter 2, synthetic seismograms are constructed for the strong motions of the 1968 Borrego Mountain earthquake recorded at EI Centro. A good fit to the data results from using the laterally varying model determined by a detailed refraction survey and the source parameters determined by teleseismic waveform modeling. Shallow faulting is no longer necessary to explain the long-period surface-wave development.
\r\n\r\nSynthetic seismograms calculated for the 1971 San Fernando earthquake show strong effects due to lateral variation in sediment thickness in the San Fernando valley and the Los Angeles basin. Using previously determined basin structure and teleseismically determined source parameters, two-dimensional SH and P-SV finite difference calculations can reproduce the amplitude and duration of the strong motion velocities recorded across the basins in Los Angeles in the period range from 1 to 10 seconds. The edges of basins nearest the seismic source show ground motion amplification up to a factor of three, and tend to convert direct shear waves into Love and Rayleigh waves that travel within the basins. The computed motions are sensitive to the mechanism and location of earthquakes. A strike-slip earthquake on the Newport-Inglewood fault zone, for example, would produce different patterns of peak velocity and duration of shaking across the San Fernando and Los Angeles basins.
\r\n\r\nIn Chapter 3, the effect of shallow station structure and lateral velocity variation are investigated for records of the Amchitka explosion Milrow. The differences between the Meuller-Murphy, Heimberger-Hadley, and von Seggern-Blandford reduced displacement potential (RDP) source representations are small compared to the differences between using various possible velocity structures.
\r\n\r\nCreager and Jordan (1986) propose that penetration of subducting slabs under the Kurile Islands and other subduction zones to depths of at least 1000 km is necessary to explain the t ravel time anomalies of deep earthquakes. Such penetration would also affect the amplitudes and waveforms of the body waves from these earthquakes. In Chapter 4, synthetic seismograms appropriate for a record section in a plane perpendicular to the strike of the slab are presented using a coupled finite-difference and Kirchhoff method. An inferred shear-wave version of the compressional-wave velocity structure of Creager and Jordan (1986) produces an amplitude decrease up to a factor of four and waveform broadening up to 20 seconds for SH arrivals with a take-off angle pointing straight down t he slab. Slabs that extend only 300 km below the earthquake but are half as thick and twice as anomalously fast as Creager and Jordan 's (1986) velocity model will roughly preserve the travel time variation pattern, and show less waveform broadening, but produce first arrivals that are emergent. Slabs that become thicker with depth show less waveform broadening. Reconciliation of the amplitude, waveform distort ion, and timing of body waves from deep events is necessary to understand the geometry of slabs near and below the 6.50 km discontinuity.
", "doi": "10.7907/7TQ9-X746", "publication_date": "1987", "thesis_type": "phd", "thesis_year": "1987" }, { "id": "thesis:11448", "collection": "thesis", "collection_id": "11448", "cite_using_url": "https://resolver.caltech.edu/CaltechTHESIS:04052019-172938000", "type": "thesis", "title": "Magmons: Solitary Waves Arising in the Buoyant Ascent of Magma by Porous Flow through a Viscously Deformable Matrix", "author": [ { "family_name": "Scott", "given_name": "David Russell", "clpid": "Scott-David-Russell" } ], "thesis_advisor": [ { "family_name": "Ahrens", "given_name": "Thomas J.", "clpid": "Ahrens-T-J" } ], "thesis_committee": [ { "family_name": "Ahrens", "given_name": "Thomas J.", "clpid": "Ahrens-T-J" }, { "family_name": "Anderson", "given_name": "Donald L.", "clpid": "Anderson-D-L" }, { "family_name": "Hager", "given_name": "Bradford H.", "clpid": "Hager-B-H" }, { "family_name": "Stevenson", "given_name": "David John", "clpid": "Stevenson-D-J" }, { "family_name": "Wasserburg", "given_name": "Gerald J.", "clpid": "Wasserburg-G-J" } ], "local_group": [ { "literal": "div_gps" } ], "abstract": "The dynamics of partially molten regions of the Earth's mantle are studied using a combination of theoretical, experimental, and numerical techniques. The physical model is based on experimental observations of partially molten ultramafic rocks and incorporates two elements: buoyancy-driven porous flow of magma through a viscously deformable matrix, and buoyancy-driven circulation of the whole rock.
\r\n\r\nThe first element of this model is analogous to buoyancy-driven pipe flow of a liquid through a denser and more viscous wall fluid. Laboratory experiments on this system illustrate the phenomenon of solitary waves. These are waves of larger pipe radius that ascend a uniform pipe of smaller radius. The waves are very nearly conserved in collisions. These, and the corresponding waves of higher porosity that arise in one-dimensional porous flow, are characterized further by analysis and numerical experiments.
\r\n\r\nThe full system, incorporating circulation in a multidimensional porous medium, also displays solitary waves governed by the same basic processes as the one-dimensional waves. Analysis and numerical experiments show that the multidimensional waves have a circular or spherical form.
\r\n\r\nA possible natural manifestation of this fluid dynamical phenomenon is in igneous processes. Magmons, as the waves are called in that setting, probably have wavelengths of kilometers and velocities of centimeters per year. Magma ascent in magmons may account for episodicity in igneous emplacement. Also, a magmon can collect and mobilize a small degree of partial melt without disturbing its geochemical signature. In a partially molten region the characteristic wavelength of magmons will always be superimposed on that of large scale variations in porosity.
", "doi": "10.7907/mvdj-ty73", "publication_date": "1987", "thesis_type": "phd", "thesis_year": "1987" }, { "id": "thesis:9559", "collection": "thesis", "collection_id": "9559", "cite_using_url": "https://resolver.caltech.edu/CaltechTHESIS:02122016-151411807", "type": "thesis", "title": "Dynamical Models for the Earth's Geoid", "author": [ { "family_name": "Richards", "given_name": "Mark Alan", "orcid": "0000-0002-1893-953X", "clpid": "Richards-Mark-Alan" } ], "thesis_advisor": [ { "family_name": "Hager", "given_name": "Bradford H.", "clpid": "Hager-B-H" }, { "family_name": "Ahrens", "given_name": "Thomas J.", "clpid": "Ahrens-T-J" } ], "thesis_committee": [ { "family_name": "Ahrens", "given_name": "Thomas J.", "clpid": "Ahrens-T-J" }, { "family_name": "Anderson", "given_name": "Donald L.", "clpid": "Anderson-D-L" }, { "family_name": "Clayton", "given_name": "Robert W.", "orcid": "0000-0003-3323-3508", "clpid": "Clayton-R-W" }, { "family_name": "Hager", "given_name": "Bradford H.", "clpid": "Hager-B-H" }, { "family_name": "Stevenson", "given_name": "David John", "orcid": "0000-0001-9432-7159", "clpid": "Stevenson-D-J" } ], "local_group": [ { "literal": "div_gps" } ], "abstract": "The Earth's largest geoid anomalies occur at the lowest spherical harmonic degrees, or longest wavelengths, and are primarily the result of mantle convection. Thermal density contrasts due to convection are partially compensated by boundary deformations due to viscous flow whose effects must be included in order to obtain a dynamically consistent model for the geoid. These deformations occur rapidly with respect to the timescale for convection, and we have analytically calculated geoid response kernels for steady-state, viscous, incompressible, self-gravitating, layered Earth models which include the deformation of boundaries due to internal loads. Both the sign and magnitude of geoid anomalies depend strongly upon the viscosity structure of the mantle as well as the possible presence of chemical layering.
\r\n\r\nCorrelations of various global geophysical data sets with the observed geoid can be used to construct theoretical geoid models which constrain the dynamics of mantle convection. Surface features such as topography and plate velocities are not obviously related to the low-degree geoid, with the exception of subduction zones which are characterized by geoid highs (degrees 4-9). Recent models for seismic heterogeneity in the mantle provide additional constraints, and much of the low-degree (2-3) geoid can be attributed to seismically inferred density anomalies in the lower mantle. The Earth's largest geoid highs are underlain by low density material in the lower mantle, thus requiring compensating deformations of the Earth's surface. A dynamical model for whole mantle convection with a low viscosity upper mantle can explain these observations and successfully predicts more than 80% of the observed geoid variance.
\r\n\r\nTemperature variations associated with density anomalies in the man tie cause lateral viscosity variations whose effects are not included in the analytical models. However, perturbation theory and numerical tests show that broad-scale lateral viscosity variations are much less important than radial variations; in this respect, geoid models, which depend upon steady-state surface deformations, may provide more reliable constraints on mantle structure than inferences from transient phenomena such as postglacial rebound. Stronger, smaller-scale viscosity variations associated with mantle plumes and subducting slabs may be more important. On the basis of numerical modelling of low viscosity plumes, we conclude that the global association of geoid highs (after slab effects are removed) with hotspots and, perhaps, mantle plumes, is the result of hot, upwelling material in the lower mantle; this conclusion does not depend strongly upon plume rheology. The global distribution of hotspots and the dominant, low-degree geoid highs may correspond to a dominant mode of convection stabilized by the ancient Pangean continental assemblage.
", "doi": "10.7907/ncf0-0c06", "publication_date": "1986", "thesis_type": "phd", "thesis_year": "1986" }, { "id": "thesis:7983", "collection": "thesis", "collection_id": "7983", "cite_using_url": "https://resolver.caltech.edu/CaltechTHESIS:10092013-140441093", "primary_object_url": { "basename": "Le Bras_r_1985.pdf", "content": "final", "filesize": 26322234, "license": "other", "mime_type": "application/pdf", "url": "/7983/1/Le Bras_r_1985.pdf", "version": "v3.0.0" }, "type": "thesis", "title": "Methods of Multiparameter Inversion of Seismic Data Using the Acoustic and Elastic Born Approximations", "author": [ { "family_name": "Le Bras", "given_name": "Ronan", "orcid": "0000-0003-2439-6938", "clpid": "Le-Bras-Ronan" } ], "thesis_advisor": [ { "family_name": "Ahrens", "given_name": "Thomas J.", "clpid": "Ahrens-T-J" } ], "thesis_committee": [ { "family_name": "Ahrens", "given_name": "Thomas J.", "clpid": "Ahrens-T-J" }, { "family_name": "Clayton", "given_name": "Robert W.", "orcid": "0000-0003-3323-3508", "clpid": "Clayton-R-W" }, { "family_name": "Kanamori", "given_name": "Hiroo", "orcid": "0000-0001-8219-9428", "clpid": "Kanamori-H" }, { "family_name": "Silver", "given_name": "Leon T.", "clpid": "Silver-L-T" }, { "family_name": "Helmberger", "given_name": "Donald V.", "clpid": "Helmberger-D-V" } ], "local_group": [ { "literal": "div_gps" } ], "abstract": "This thesis presents two different forms of the Born approximations for acoustic and elastic wavefields and discusses their application to the inversion of seismic data. The Born approximation is valid for small amplitude heterogeneities superimposed over a slowly varying background. The first method is related to frequency-wavenumber migration methods. It is shown to properly recover two independent acoustic parameters within the bandpass of the source time function of the experiment for contrasts of about 5 percent from data generated using an exact theory for flat interfaces. The independent determination of two parameters is shown to depend on the angle coverage of the medium. For surface data, the impedance profile is well recovered.
\r\n\r\nThe second method explored is mathematically similar to iterative tomographic methods recently introduced in the geophysical literature. Its basis is an integral relation between the scattered wavefield and the medium parameters obtained after applying a far-field approximation to the first-order Born approximation. The Davidon-Fletcher-Powell algorithm is used since it converges faster than the steepest descent method. It consists essentially of successive backprojections of the recorded wavefield, with angular and propagation weighing coefficients for density and bulk modulus. After each backprojection, the forward problem is computed and the residual evaluated. Each backprojection is similar to a before-stack Kirchhoff migration and is therefore readily applicable to seismic data. Several examples of reconstruction for simple point scatterer models are performed. Recovery of the amplitudes of the anomalies are improved with successive iterations. Iterations also improve the sharpness of the images.
\r\n\r\nThe elastic Born approximation, with the addition of a far-field approximation is shown to correspond physically to a sum of WKBJ-asymptotic scattered rays. Four types of scattered rays enter in the sum, corresponding to P-P, P-S, S-P and S-S pairs of incident-scattered rays. Incident rays propagate in the background medium, interacting only once with the scatterers. Scattered rays propagate as if in the background medium, with no interaction with the scatterers. An example of P-wave impedance inversion is performed on a VSP data set consisting of three offsets recorded in two wells.
\r\n", "doi": "10.7907/14q2-fa62", "publication_date": "1985", "thesis_type": "phd", "thesis_year": "1985" }, { "id": "thesis:7477", "collection": "thesis", "collection_id": "7477", "cite_using_url": "https://resolver.caltech.edu/CaltechTHESIS:02142013-085850235", "type": "thesis", "title": "Studies of the Crust-Mantle System Beneath Southern California", "author": [ { "family_name": "Humphreys", "given_name": "Eugene Drake", "orcid": "0000-0002-1916-8378", "clpid": "Humphreys-Eugene-Drake" } ], "thesis_advisor": [ { "family_name": "Ahrens", "given_name": "Thomas J.", "clpid": "Ahrens-T-J" } ], "thesis_committee": [ { "family_name": "Clayton", "given_name": "Robert W.", "orcid": "0000-0003-3323-3508", "clpid": "Clayton-R-W" }, { "family_name": "Hager", "given_name": "Bradford H.", "clpid": "Hager-B-H" }, { "family_name": "Silver", "given_name": "Leon T.", "clpid": "Silver-L-T" }, { "family_name": "Anderson", "given_name": "Donald L.", "clpid": "Anderson-D-L" }, { "family_name": "Ahrens", "given_name": "Thomas J.", "clpid": "Ahrens-T-J" } ], "local_group": [ { "literal": "div_gps" } ], "abstract": "A back-projection method of tomographic reconstruction is adapted to inverted seismic travel-time data. The problems encountered in inverting these data include ray set inhomogeneity and anisotropy and the three-dimensionality of the space interrogated. Jacobi iteration, deconvolution and variable ray weighting are shown to work well in augmenting the basic back-projection method to produce a well-focused image. Applications of the various focusing algorithms are shown to have a degree of success that depends on the ray geometry used. Also, the ability to reconstruct an accurate image when the data include moderate amounts of noise is shown to be good.
\r\n\r\nP-wave teleseismic travel time delays recorded by the southern California array are inverted with the tomographic method to obtain variations in the P-wave velocity structure to a depth of 750 km. Two major anomalies are imaged. A curtain-like E-W trending high velocity feature is found directly beneath the Transverse Ranges. This feature is about 50 km in thickness, extends in depth to a maximum of about 250 km on its eastern end, and attains a maximum velocity that is about 3% greater than average mantle at the same depth beneath southern California. A zone of low-velocity material is found in the uppermost 100 km beneath the region of the Salton Trough. The seismic velocities here are depressed by about 4%.
\r\n\r\nThese anomalous regions are interpreted to be related to the geologic processes that have been active recently in southern California. Scaling relations are used to estimate that the Transverse Range anomaly is about 500\u00b0C colder and 1% more dense than average southern California mantle of the same depth, while the Salton Trough anomaly is about 1/2% less dense and contains about 3% melt. The density distribution drives a flow of upper mantle material from the Salton Trough region towards the Transverse Ranges, where it sinks into the mantle to form the feature seen beneath these ranges. Mantle flow results in tractions that act on the base of the lithosphere to produce stresses within the lithosphere that are tensile in the Salton Trough and compressive in the Transverse Ranges. These stresses are thought to account for the physiography seen in these provences.
\r\n\r\nThe southern California crust is modeled using late Quaternary slip rates on major faults, and a kinematic description is determined that has: 1) only local sites of convergence in the Transverse Ranges, and 2) the occurrance of significant strain rates near to the southern California coast, including the western Transverse Ranges.
", "doi": "10.7907/j2xf-gq29", "publication_date": "1985", "thesis_type": "phd", "thesis_year": "1985" }, { "id": "thesis:3777", "collection": "thesis", "collection_id": "3777", "cite_using_url": "https://resolver.caltech.edu/CaltechETD:etd-09262002-154053", "primary_object_url": { "basename": "Boslough_mb_1984.pdf", "content": "final", "filesize": 8832552, "license": "other", "mime_type": "application/pdf", "url": "/3777/1/Boslough_mb_1984.pdf", "version": "v3.0.0" }, "type": "thesis", "title": "Shock-Wave Properties and High-Pressure Equations of State of Geophysically Important Materials", "author": [ { "family_name": "Boslough", "given_name": "Mark Bruce", "orcid": "0000-0001-6912-4608", "clpid": "Boslough-Mark-Bruce" } ], "thesis_advisor": [ { "family_name": "Ahrens", "given_name": "Thomas J.", "clpid": "Ahrens-T-J" } ], "thesis_committee": [ { "family_name": "Ahrens", "given_name": "Thomas J.", "clpid": "Ahrens-T-J" }, { "family_name": "Corngold", "given_name": "Noel Robert", "clpid": "Corngold-N-R" }, { "family_name": "Stevenson", "given_name": "David John", "orcid": "0000-0001-9432-7159", "clpid": "Stevenson-D-J" }, { "family_name": "Roshko", "given_name": "Anatol", "clpid": "Roshko-A" }, { "family_name": "Vreeland", "given_name": "Thad", "clpid": "Vreeland-T" } ], "local_group": [ { "literal": "div_gps" } ], "abstract": "Shock wave (Hugoniot), shock temperature, and release data are presented for several geophysically important, refractory materials. A sensitive multichannel optical pyrometer was developed to measure shock temperatures (2500 to 5600\u00b0K at pressures from 48 to 117 GPa) in anorthite (CaAl2Si2O8) glass. Shock temperatures of 3750 to 6000\u00b0K at pressures from 140 to 182 GPa were measured in calcium oxide (CaO). Temperature data were used to constrain the energetics of the B1-B2 phase transition at 70 GPa in CaO, and to construct a finite-strain equation of state for CaO consistent with previous Hugoniot data.
\r\n\r\nThe new CaO equation of state was used with equation of state parameters of other oxides to construct a theoretical mixed oxide Hugoniot of anorthite, which is in agreement with new Hugoniot data above about 50 GPa, determined using new experimental techniques developed in this study. The mixed oxide model, however, overestimates the shock temperatures, and does not accurately predict measured release paths. Both shock temperature and release data for anorthite indicate that several high pressure phase regions of stability exist above 50 GPa. A similar mixed oxide Hugoniot was constructed for lunar gabbroic anorthosite, and agrees with two new Hugoniot points at 120 GPa. Release data from lunar gabbroic anorthosite shocked to 120 GPa give evidence for shock vaporization.
\r\n\r\nBecause the densities and bulk properties of CaO and the high pressure phase or phases of anorthite are so close to those determined seismologically for the lower mantle, the amount of these materials present in the lower mantle is not well constrained. The possibility of significant enrichment of the lower mantle in these refractory materials, as predicted by inhomogeneous accretion models, is still open.
\r\n\r\nA simple model is developed to explain the measured time dependences of radiated light in the shock temperature experiments, and constrain the absorption coefficient of the shocked material. The absorption coefficient is found to be an increasing function of shock pressure in shocked anorthite glass.
\r\n\r\nHugoniot and release paths were determined using electromagnetic particle velocity gauges for San Gabriel anorthosite and San Marcos Gabbro shocked to peak stresses between 5 and 11 GPa. The data indicate a loss of shear strength in both rocks, and a partial phase transition of the anorthosite to a denser phase. This implies that estimates of shock wave attenuation in these materials based on elastic-plastic models are too high, and previously calculated amounts of internal energy gained by surface materials from impact or explosion events have been underestimated.
", "doi": "10.7907/C9GZ-3121", "publication_date": "1984", "thesis_type": "phd", "thesis_year": "1984" }, { "id": "thesis:1324", "collection": "thesis", "collection_id": "1324", "cite_using_url": "https://resolver.caltech.edu/CaltechETD:etd-04102003-122837", "primary_object_url": { "basename": "gillespie_ar_1982.pdf", "content": "final", "filesize": 50336898, "license": "other", "mime_type": "application/pdf", "url": "/1324/14/gillespie_ar_1982.pdf", "version": "v4.0.0" }, "type": "thesis", "title": "Quaternary Glaciation and Tectonism in the Southeastern Sierra Nevada, Inyo County, California", "author": [ { "family_name": "Gillespie", "given_name": "Alan Reed", "clpid": "Gillespie-Alan-Reed" } ], "thesis_advisor": [ { "family_name": "Ahrens", "given_name": "Thomas J.", "clpid": "Ahrens-T-J" } ], "thesis_committee": [ { "family_name": "Ahrens", "given_name": "Thomas J.", "clpid": "Ahrens-T-J" }, { "family_name": "Allen", "given_name": "Clarence R.", "clpid": "Allen-C-R" }, { "family_name": "Sieh", "given_name": "Kerry E.", "orcid": "0000-0002-7311-2447", "clpid": "Sieh-K-E" }, { "family_name": "Sharp", "given_name": "Robert P.", "clpid": "Sharp-R-P" }, { "family_name": "Wasserburg", "given_name": "Gerald J.", "orcid": "0000-0002-7957-8029", "clpid": "Wasserburg-G-J" }, { "family_name": "Clark", "given_name": "M. M.", "clpid": "Clark-M-M" } ], "local_group": [ { "literal": "div_gps" } ], "abstract": "The southeastern Sierra Nevada consists of three geographic regions. From west to east, they are: an upland region across the crest, a steep east-facing escarpment along which Owens Valley has partly subsided, and foothill blocks intermediate to the Sierra Nevada and Owens Valley. Farther east, Owens Valley is a deep graben separating the Sierra Nevada and the Inyo Range.
\r\n\r\nThe main goals of this thesis were the detailed mapping of Quaternary glacial and other deposits in these regions, dating of critical events, and geomorphic analysis of the range front. The focus was on Pleistocene moraines near the range front. The motivation of this research was to improve our understanding of the chronology of Pleistocene events, to characterize details of the tectonic history of the Sierra, to infer faulting, erosion, and deposition rates, and to provide a basis for the comparison of the Quaternary geology in the southeastern Sierra and in more intensively studied regions in the central and northern Sierra and elsewhere.
\r\n\r\nThe study area extended from the alluvial fans of Owens Valley west to the Sierra crest from latitude 36\u00b045' N to 37\u00b000' N. It included the southern part of the Big Pine volcanic field, an eruptive center for basaltic lavas for most of the Pleistocene Epoch. Elevations within the study area ranged from about 1000 m (Owens Valley) to about 4000 m (peaks along the crest).
\r\n\r\nThroughout the study area the principal rocks are granodiorite and quartz monzonite of Cretaceous age. Plutons are rather small, and individual drainages generally include more than one. In the southern part of the study area, Jurassic-Triassic metavolcanic rocks are found as roof pendants. These rocks, originally ranging in composition from basalt to rhyolite, are most common near the Sierra crest. In the northern canyons of the study area, Paleozoic metasedimentary rocks including sandy marbles and biotite schist replace the metavolcanic pendants. The foothill blocks are identical to the Sierras in composition.
\r\n\r\nBelow the foothills coalescing alluvial fans grade a few km east to the alluvium and lacustrine sediments of the Owens River and Owens Lake. These sediments have been shown in geophysical studies to mask a second escarpment as high as the one of the range front, and the total bedrock relief from the Sierra crest to the floor of the graben is as much as 6 km.
\r\n\r\nDuring the Quaternary Period the southeastern Sierra Nevada was characterized by the down-faulting of Owens Valley along two zones, one a series of normal faults along the range front (Independence Fault) and the other a series of faults along the center of the valley (Owens Valley fault zone). This same period has seen the cutting of deep canyons through the 2-km-high escarpment. During repeated glaciations these canyons were widened and deepened. Traces of at least seven glaciations were found during this study. Moraines and other deposits left during these glaciations could be distinguished based on the degree of weathering of granitic clasts, vegetative cover, and morphologic characteristics. Absolute age limits were obtained for two of the Pleistocene glaciations by radiometric dating of basalt flows interfingered with the moraines.
\r\n\r\nThree of the recognized glaciations, probably corresponding to the Matthes, Recess Peak, and Hilgard neoglaciations found by J.H. Birman in the central Sierra Nevada, occurred during the Holocene Epoch. The youngest glaciers (Matthes glaciation) left unconsolidated and unvegetated till in stagnant rock glaciers and moraines in cirques on high peaks. Some rock glaciers are still ice-cored. Extending out from the cirques and into the upper reaches of the canyons are moraines correlating to the Recess Peak glaciation. Till is generally consolidated and supports heavy lichen growth and bushes but few trees. The oldest Holocene glaciation (Hilgard) left few large moraines in the study area. Hilgard glaciers extended much farther down-canyon than the younger Holocene glaciers, sometimes within one or two km of the Tioga terminal moraines. Those Hilgard terminal moraines which were found have been barely breached by streams. Moraines tend to be heavily forested, and lakes are largely unsedimented. The Hilgard glaciation may have simply been the last stade of the Tioga glaciation from the evidence found in this study.
\r\n\r\nAt least four Pleistocene glaciations occurred in the southeastern Sierra Nevada. All four postdate most of the significant incision of streams through the escarpment. The three youngest probably correlate with the Tioga, Tenaya, and Tahoe stages (in order of increasing age) recognized throughout the Sierra. In each case, moraine morphology has been well preserved. Tioga moraines were found down to about 2200 m elevation. Nested sets of moraines were common. The terminal moraines of the youngest of these were sometimes intact. lakes were rare; one (Sawmill Meadow) was completely sedimented. Granitic boulders in the moraines were largely unweathered. Weathering of boulders in Tenaya moraines was similar, but a small fraction of granitic boulders were grusy. Boulders from Tahoe moraines were conspicuously weathered, and the moraines themselves were rounded and gullied.
\r\n\r\nThe oldest group of moraines probably significantly predates the Tahoe glaciation. It is nevertheless post-Sherwin. Moraines in this group were found in five of the eight canyons studied. While obviously eroded, these moraines still retained their original shape. All surficial boulders were heavily weathered, but some exposed in road cuts were fresh. No moraines of Sherwin age were identified, although Sherwin till is widespread only a few km to the north. However, on plateaus and ridges 200 to 300 m above the modern canyons near the Sierra crest were found ancient diamictons, some of which could be till. Remnants of U-shaped valleys preserved as high passes across the crest or as cols between canyons east of the crest may be testimonials to ancient glaciers of Sherwin age or older.
\r\n\r\nRadiometric dating (40Ar - 39Ar) of basalts interfingered with moraines in Sawmill Canyon provided a new upper limit of 0.12 my for a moraine of the Tahoe glaciation, and a range of 0.13 - 0.46 my for one pre-Tahoe glaciation. These results confirm that the Tahoe glaciation occurred during the Wisconsin stage of the continental ice age, and conclusively demonstrate the existence of pre-Wisconsin glaciers in the southern Sierra. Relative dating based on acoustic wave speeds through weathered boulders on the moraines indicates the age of the pre-Wisconsin moraine may be close to the upper limit.
\r\n\r\nAlluvial fans appear to have aggraded early in the Wisconsin glaciation (Tahoe glaciation). Subsequently, the fan heads have been incised and the locus of deposition has moved eastward down the fans. The Tenaya and Tioga glaciers during the late Wisconsin stage left outwash plains and terraces along streams cut into the older fans, but aggradation during these events was considerably less than earlier.
\r\n\r\nThree ages of fans were found. The oldest fanglomerate probably is pre-Wisconsin and is exposed in regions protected from later deposition. The heavily weathered fan deposits of this group overlie basalts which appear to be contemporaneous with dated 1.1 - 1.2 my-old basalts nearby. In the middle elevations of the fans, roughly 10 m of fanglomerate was deposited over the old fanglomerate, probably during the Tahoe glaciation. Deposition rates probably are about 0.1 mm/y for the late Pleistocene Epoch. The extent and distribution of the youngest fans (Tenaya-Tioga) are variable, but they are generally found downstream from the incised Tahoe fan heads.
\r\n\r\nFaulting along the range front appears to have been dip-slip only. The offset rate along the range-front faults was determined at several canyons where the fault crossed dated moraines or lava flows. At least during the Wisconsin glaciation faulting on this zone appears to have been erratic, with rates ranging from zero to 0.5 mm/y or more. Offset moraines and terraces at Independence Creek indicated a faulting rate of 0.1 mm/y. Only a few km to the north, Tahoe moraines of both forks of Oak Creek were not offset at all, although scarps could be seen on adjacent hillsides. At Sawmill Creek an offset lava flow gave a lower limit of 0.5 mm/y. It seems that during the late Pleistocene Epoch, offset on the range front faults has been less than on the mid-valley faults east of the study area. Geodetic studies have suggested modern strain rates of 2.2 mm/y for the Owens Valley fault zone.
\r\n\r\nBasalts found in canyons through the escarpment and on terraces and ridges in the foothills to the east document stream erosion during the Pleistocene Epoch. Ridgetop basalts, dated at 1.2 my, stand at least 125 m above the modern streams through the foothills. This indicates an erosion rate of ~ 0.1 mm/y. A comparable rate of ~ 0.15 mm/y for the last 0.46 my was found for Sawmill Creek within the Sierra Nevada. Thus at least during the late Pleistocene Epoch erosion rates in the Sierra and in the foothills have been similar.
\r\n\r\nPatches of boulders and gravels atop the basalt show that some time after 1.2 my BP the foothill block was submerged by alluvial fans. Incision may have begun in response to the inception or renewal of subsidence of the graben along the Owens Valley fault zone.
\r\n\r\nExtensive volcanism in the Big Pine Volcanic field appears to have begun at least 1.2 my ago, and has continued sporadically up to perhaps 0.05 my ago. Minor eruptions may have occurred more recently.
\r\n\r\nThe eastern escarpment of the Sierra Nevada consists of two zones of truncated ridges. Within the study area, the upper zone is about 950 m high; the lower is about 750 m high. Triangular facets of the upper zone have a gradient of only ~ 24\u00b0, lower than slopes of ~ 29\u00b0 in the lower zone. This could be explained if subsidence of Owens Valley along the range-front faults occurred in two great pulses.
\r\n\r\nBoth absolute and relative dating methods were refined for this study. Absolute dating of the K-poor basaltic lavas was done indirectly, by 40Ar - 39Ar analysis of K-rich granitic xenoliths found in the lava. These ancient xenoliths were partially degassed of their accumulated 40Ar during heating in the magma, and it proved possible to date this heating event.
\r\n\r\nIn addition to conventional relative dating methods, a new quantitative approach based on the speed of acoustic waves through individual clasts in a deposit was investigated. This method had been used only once before, on terrace deposits. The technique proved to be very useful, and was capable of discriminating moraines successfully in well-studied canyons in the central Sierra. Acoustic wave speeds may be controlled by the abundance of intergranular cracks in granitic boulders. If this is the case, then this technique relies on different processes than those exploited by conventional methods of relative dating. The successful application to moraines in this study enhances our ability to analyse glacial sequences and complements conventional semi-quantitative methods of relative dating.
", "doi": "10.7907/GNES-QH83", "publication_date": "1982", "thesis_type": "phd", "thesis_year": "1982" }, { "id": "thesis:3552", "collection": "thesis", "collection_id": "3552", "cite_using_url": "https://resolver.caltech.edu/CaltechETD:etd-09152006-102853", "type": "thesis", "title": "Seismic Source Processes and Tectonics: Observations of Four Intracontinental Earthquakes", "author": [ { "family_name": "Cipar", "given_name": "John Joseph", "clpid": "Cipar-John-Joseph" } ], "thesis_advisor": [ { "family_name": "Ahrens", "given_name": "Thomas J.", "clpid": "Ahrens-T-J" }, { "family_name": "Helmberger", "given_name": "Donald V.", "clpid": "Helmberger-D-V" }, { "family_name": "Kanamori", "given_name": "Hiroo", "orcid": "0000-0001-8219-9428", "clpid": "Kanamori-H" } ], "thesis_committee": [ { "family_name": "Ahrens", "given_name": "Thomas J.", "clpid": "Ahrens-T-J" }, { "family_name": "Harkrider", "given_name": "David G.", "clpid": "Harkrider-D-G" }, { "family_name": "Helmberger", "given_name": "Donald V.", "clpid": "Helmberger-D-V" }, { "family_name": "Kanamori", "given_name": "Hiroo", "orcid": "0000-0001-8219-9428", "clpid": "Kanamori-H" }, { "family_name": "Sieh", "given_name": "Kerry E.", "orcid": "0000-0002-7311-2447", "clpid": "Sieh-K-E" } ], "local_group": [ { "literal": "div_gps" } ], "abstract": "This thesis presents studies of the source processes of four shallow earthquakes and their relation to regional tectonics. In the first chapter, long-period teleseismic P and S waves from the Haicheng, China earthquake of February 4, 1975 are compared directly to time domain synthetic seismograms to infer source parameters. The P-wave focal mechanism indicates that faulting was dominantly left-lateral strike-slip along a northwest striking nodal plane (strike = 288\u00b0, dip = 78\u00b0N, rake = 342\u00b0). The strike of this nodal plane agrees with the trend of the aftershock distribution. Azimuthal variation of P-wave duration is attributed to fault rupture 22 km in a northwesterly direction, along strike of the aftershock zone. There is considerable discrepancy between the observed SH waves and synthetics computed using this model. These discrepancies are due to either structural complexities in the source region or change in fault mechanism as the rupture propagated along strike. Seismic moment, average dislocation and stress drop are computed to be 2.7 x 1026 dyne-cm, 2.5 meters and 48 bars, respectively.
\r\n\r\nThe remaining three chapters present a detailed examination of seismograms recorded by the 1976 Friuli, Italy earthquake (May 6, 1976, MS = 6.5) and two major aftershocks (both on September 15, 1976 at 03h 15m, MS = 6.0 and 09h 21m, MS = 5.9). Teleseismic long-period body waves and surface waves radiated by the mainshock and 09h 21m aftershock are studied in Chapter Two to determine source characteristics. Focal mechanisms along with geological evidence suggest that both events represent underthrusting of the Friuli Plain beneath the Southern Alps. The depths of both earthquakes, estimated by matching synthetic body wave seismograms to observations, are found to lie between 6 and 10 km. Synthetic seismogram calculations which include source directivity effects suggest that the fault length of the mainshock is approximately 16 to 24 km assuming a rupture velocity of 3.0 km/sec. Observations of 100 sec Rayleigh waves confirm the body wave focal mechanism, but suggest that the seismic moment of the mainshock is 5 x 1025 dyne-cm compared to 2.9 x 1025 dyne-cm estimated from body waves. The P-wave moment of the aftershock is 1 x 1025 dyne-cm.
\r\n\r\nIn Chapter Three, short-period records are modeled to obtain additional details of the source time history. Two point sources of radiation are required to adequately model the aftershock short-period records. For the 09h 21m aftershock, the model derived from short-period records also produces good fits to the long-period data. The SP model for the 03h 15m aftershock, on the other hand, predicts long-period synthetics which do not agree with the observations. In particular, the SP moment (0.37 x 1025 dyne-cm) is about 2-1/2 times smaller than the LP moment (1 x 1025 dyne-cm). Adding a long-period component to the SP model considerably improves LP waveform and moment agreement. In the case of the mainshock, a reasonable fit to the observed SP data is obtained using three point sources of radiation. However, LP synthetics computed using this model do not agree with the observations, and the SP moment (0.65 x 1025 dyne-cm) is a small fraction of the LP moment (3-5 x 1025 dyne-cm). Time function durations indicate that the individual events inferred from the SP records are radiated from patches of the fault having radii of 2 to 4 km and stress drops in the range 35 to 276 bars. In comparison, overall stress drops estimated from LP data are found to be 12 bars (mainshock) and 24 bars (09h 21m aftershock). Strong-motion accelerograms are used to put additional constraint on the source geometry of the 09h 21m aftershock.
\r\n\r\nThe 03h 15m and 09h 21m aftershocks are the culminating events of a series of large aftershocks which began on September 11, 1976. Cumulative seismic moment of the Friuli aftershock sequence was as large as the moment released by the mainshock. By comparison, aftershock moments of California earthquakes are typically 1 to 10 percent of the mainshock moment. The large size, location and focal mechanism of the aftershocks suggest that they represent failure of major stress concentrations remaining after the mainshock.
", "doi": "10.7907/bwpb-4c04", "publication_date": "1981", "thesis_type": "phd", "thesis_year": "1981" }, { "id": "thesis:3001", "collection": "thesis", "collection_id": "3001", "cite_using_url": "https://resolver.caltech.edu/CaltechETD:etd-08032004-150047", "primary_object_url": { "basename": "Lyzenga_ga_1980.pdf", "content": "final", "filesize": 9181240, "license": "other", "mime_type": "application/pdf", "url": "/3001/1/Lyzenga_ga_1980.pdf", "version": "v3.0.0" }, "type": "thesis", "title": "Shock Temperatures of Materials: Experiments and Applications to the High Pressure Equation of State", "author": [ { "family_name": "Lyzenga", "given_name": "Gregory Allen", "clpid": "Lyzenga-Gregory-Allen" } ], "thesis_advisor": [ { "family_name": "Ahrens", "given_name": "Thomas J.", "clpid": "Ahrens-T-J" } ], "thesis_committee": [ { "family_name": "Unknown", "given_name": "Unknown" } ], "local_group": [ { "literal": "div_gps" } ], "abstract": "The experimental determination of temperatures in the high-pressure shocked state of condensed matter provides a useful supplement to equation-of-state models derived from Hugoniot measurements. An optical pyrometry technique has been developed to obtain temperature measurements during impact-driven shock wave experiments with solid and liquid samples at pressures near 100 GPa. Experimental results confirm that throughout moderate ranges of shock pressure amplitude, transparent dielectrics emit thermal radiation from the region of the shock front, with a spectrum which is characteristic of the Hugoniot state temperature. Shock temperatures in sodium chloride crystals have been measured in the pressure range 70-105 GPa. The observed temperatures, between 4000 and 8000 K, are in agreement with the results of earlier determinations and with calculations assuming the occurrence of shock-induced melting. Results of experiments to measure shock temperatures in metallic silver include a successful measurement of 5950 K at a pressure of 185 GPa. This result is consistent with the melting of silver under shock, with a melting pressure dependence described by the Lindemann criterion.
\r\n\r\nShock temperature measurements in silica (SiO2) have produced anomalous results suggestive of melting occurring in the stishovite phase near 100 GPa pressure and 4700 K temperature. Experimental measurements with [alpha]-quartz and fused silica samples extend from pressures of approximately 60 GPa to 140 GPa, with shock temperatures between approximately 4500 K and 7000 K. The experimental data allow quantification of the thermodynamic relations among silica phases, including heats of transition and the Gruneisen parameter. Shock temperatures in single crystal forsterite (Mg2SiO4) between pressures of 150 GPa and 175 GPa range from 4500 K to 4950 K, a result which is consistent with occurrence of a polymorphic solid state transition accompanied by a substantial heat of transition (~1.5 MJ/kg). These results have potentially important implications for solid earth geophysics, and knowledge of the melting curves of candidate minerals of the earth's mantle provides some constraints on the geotherm.
\r\n\r\nHugoniot temperatures have been measured in liquid water between approximately 50 GPa and 80 GPa, with results ranging from 3500 K to 5400 K. The observed temperatures are well reproduced by theoretical calculations assuming a constant specific heat model, but further work is required to characterize fully the thermal variation of H2O properties at high temperature. Compression measurements in pressure-volume states other than Hugoniot shock states and, in particular, in states of compression at constant entropy can provide both independent thermal equation-of-state information and the properties of high-density condensed phases inaccessible to shock wave experiments. Numerical calculations have been carried out for hypothetical experiments on water as well as carbon dioxide and liquid molecular hydrogen. The results of these calculations indicate that various experimental impact configurations may be employed to convert shock compression into isentropic compression, with pressures on the order of 100 GPa attained via impact velocities of a few kilometers per second. In the case of water, a net entropy production of a few percent of the Hugoniot entropy at the same pressure is predicted on the basis of these calculations.
", "doi": "10.7907/VSBK-ZA43", "publication_date": "1980", "thesis_type": "phd", "thesis_year": "1980" }, { "id": "thesis:16138", "collection": "thesis", "collection_id": "16138", "cite_using_url": "https://resolver.caltech.edu/CaltechTHESIS:07072023-223914164", "type": "thesis", "title": "Physics of Mantle and Core Minerals", "author": [ { "family_name": "Jeanloz", "given_name": "Raymond", "clpid": "Jeanloz-Raymond" } ], "thesis_advisor": [ { "family_name": "Ahrens", "given_name": "Thomas J.", "clpid": "Ahrens-T-J" } ], "thesis_committee": [ { "family_name": "Unknown", "given_name": "Unknown" } ], "local_group": [ { "literal": "div_gps" } ], "abstract": "Shock-wave equation-of-state (Hugoniot) data for initially porous and nonporous samples of iron provide experimental support for theoretically calculated properties of the earth's core, and show that whereas both densities and bulk moduli in the outer core are less than those of Fe under equivalent conditions (by about 10% and 12% respectively) their gradients with pressure are consistent with gross chemical homogeneity through the outer core; simple dynamic models of the core are allowed. New Hugoniot data for w\u00fcstite show that oxygen (~11 wt. %) can be the element which reduces the density of the outer core below that of Fe. The properties of the inner core are consistent with those of iron, suggesting that the inner core-outer core boundary is both a phase and a compositional boundary. The minimum estimated temperature at the top of the outer core is ~2800K, whereas subsolidus phase equilibria of olivine indicate a temperature near 2000K in the transition zone. Hugoniot data for porous and nonporous MgO and SiO\u2082 (phases considered representative of the lower mantle) provide an experimentally\u00ad constrained (adiabatic) geothermal gradient through the lower mantle which implies the presence of one or (for a more consistent result) more thermal boundary layers in the lower mantle. These suggest that the core is a major heat source for the mantle and that a barrier to convection occurs in (or near the top of) the lower mantle: a chemical discontinuity would be a likely cause. This inference is consistent with new shock\u00ad wave data for Cao which show that calcium could be substantially enriched in the lower mantle, as suggested by inhomogeneous accretion theories. A thermal equation of state is determined for anorthite from porous and nonporous Hugoniot data which, however, show that this refractory mineral can probably not be a major Ca-bearing phase in the lower mantle, except perhaps near the core-mantle boundary. Diamond-cell and Hugoniot data show that CaO undergoes a B1/B2 transition at 70 GPa with properties well predicted theoretically. FeO undergoes a similar transition (at ~70 GPa) and these results suggest that transformation in magnesiow\u00fcstite may be important in the lowermost mantle. New Hugoniot data for bronzite are combined with previous shock-wave measurements for olivines and pyroxenes. These data are consistent with static high\u00ad pressure results, but suggest the occurrence of post-perovskite phases (density ~5% greater than perovskite) and they also provide evidence of nonequilibrium effects under shock to pressures above 100 GPa. Spectroscopic and microscopic studies of shock-compressed olivines support this evidence: the structure of olivine achieved under shock is apparently far from equilibrium, as is indicated by phase-transformation theory. Although the bulk properties measured under shock are consistent with the attainment of thermodynamic equilibrium, these properties apparently represent highly transient and nonequilibrium states.", "doi": "10.7907/dzss-k041", "publication_date": "1980", "thesis_type": "phd", "thesis_year": "1980" }, { "id": "thesis:13675", "collection": "thesis", "collection_id": "13675", "cite_using_url": "https://resolver.caltech.edu/CaltechTHESIS:04142020-164137815", "type": "thesis", "title": "I. Regional Variations in Upper Mantle Compressional Velocities beneath Southern California. II. Post-Shock Temperatures: Their Experimental Determination, Calculation, and Implications", "author": [ { "family_name": "Raikes", "given_name": "Susan Ann", "clpid": "Raikes-Susan-Ann" } ], "thesis_advisor": [ { "family_name": "Kanamori", "given_name": "Hiroo", "orcid": "0000-0001-8219-9428", "clpid": "Kanamori-H" }, { "family_name": "Ahrens", "given_name": "Thomas J.", "clpid": "Ahrens-T-J" } ], "thesis_committee": [ { "family_name": "Unknown", "given_name": "Unknown" } ], "local_group": [ { "literal": "Seismological Laboratory" }, { "literal": "div_gps" } ], "abstract": "The establishment in Southern California of a large seismographic network provides an unique opportunity for studying the seismic velocity variations within a tectonically active region that includes a major plate boundary, whose surface expression is the San Andreas Fault. In the first part of this thesis, the compressional velocity within the upper mantle beneath Southern California is investigated through observations of the dependence of teleseismic P-delays at all stations of the array on the distance and azimuth to the event. The variation of residuals with azimuth was found to be as large as 1.3 sec at a single station; the delays Here stable as a function of time, and no evidence was found for temporal velocity variations related to seismic activity in the area. These delays were used in the construction of models for the upper mantle P-velocity structure to depths of 150 km, both by ray tracing and inversion techniques. The models exhibit considerable lateral heterogeneity including a region of low velocity beneath the Imperial Valley, and regions of increased velocity beneath the Sierra Nevada and much of the Transverse \u00b7Ranges. These changes are attributed to variation in the degree of partial melting within the upper mantle; their relationship to, and implications for, regional tectonics are discussed in the final chapter of this section.
\r\n\r\nOne of the major uncertainties in the interpretation of shock wave data is the temperature reached under shock compression and subsequent release. The second half of this thesis describes the development of a technique for the experimental determination of post-shock temperatures, its application to several metals and silicates shocked to pressures in the range 5 to 30 CPa. The technique utilises an infra-red radiation detector to determine the brightness temperature of the free surface of the sample after the shock wave has passed through it, and has yielded highly reproducible results that are consistent for the wavelength ranges 4.5 to 5.75 and 7 to 14\u00b5. The comparison of these results with values calculated using conventional theories provides some insight into the thermal processes occurring in shock waves. In particular, the measured temperatures are generally higher than those calculated; this is attributed to elasto-plastic effects in metals, and is probably associated with strength effects in silicates, both of which are commonly ignored in the calculation of theoretical temperatures. The implications of these observations for the interpretation of shock-induced metamorphism and impact phenomena, and for the application of shock-wave data to the interpretation of the behaviour of silicates within the earth's mantle, are discussed in the final chapter.
", "doi": "10.7907/3323-8080", "publication_date": "1978", "thesis_type": "phd", "thesis_year": "1978" }, { "id": "thesis:4392", "collection": "thesis", "collection_id": "4392", "cite_using_url": "https://resolver.caltech.edu/CaltechETD:etd-11032005-141608", "primary_object_url": { "basename": "Gibbons_rv_1974.pdf", "content": "final", "filesize": 15332503, "license": "other", "mime_type": "application/pdf", "url": "/4392/1/Gibbons_rv_1974.pdf", "version": "v4.0.0" }, "type": "thesis", "title": "Experimental Effects of High Shock Pressure on Materials of Geological and Geophysical Interest", "author": [ { "family_name": "Gibbons", "given_name": "Rex Vincent", "clpid": "Gibbons-Rex-Vincent" } ], "thesis_advisor": [ { "family_name": "Ahrens", "given_name": "Thomas J.", "clpid": "Ahrens-T-J" } ], "thesis_committee": [ { "family_name": "Unknown", "given_name": "Unknown" } ], "local_group": [ { "literal": "div_gps" } ], "abstract": "Shock recovery studies have been performed on a number of materials of geological and geophysical significance using a shock-loading propellant gun. These materials, including silicate glasses, feldspar, orthopyroxene, and pyroxenoid, have been shock-loaded to pressures up to approximately 550 kilobars. Optical and electron microscopic, optical spectral, and x-ray diffraction studies have been carried out on the recovered samples to determine the permanent effects of the shock-loading and the applicability of such information to research on naturally shocked lunar and terrestrial rocks and meteorites. The data on the pyroxene and feldspar are discussed in terms of their usefulness in interpreting, understanding, and calibrating the pressure and temperature conditions of shock metamorphism. The most significant observations include (1) the behavior of the silicate glasses, especially the occurrence of permanent densification until high postshock temperatures cause reversion to low density glass, (2) the permanent reduction of Mn(III) to Mn(II) by shock-loading to 496 kilobars, (3) the shock vitrification of calcic plagioclase to diaplectic and shock-fused glasses at pressures above 300 kilobars, (4) the shockproduction of deformation lamellae in bronzite, and (5) the formation of shock glass in bronzite at 226 kilobars.", "doi": "10.7907/9BJQ-N519", "publication_date": "1974", "thesis_type": "phd", "thesis_year": "1974" }, { "id": "thesis:14017", "collection": "thesis", "collection_id": "14017", "cite_using_url": "https://resolver.caltech.edu/CaltechTHESIS:12072020-234014838", "primary_object_url": { "basename": "Lagus_PL_1974.pdf", "content": "final", "filesize": 30518048, "license": "other", "mime_type": "application/pdf", "url": "/14017/1/Lagus_PL_1974.pdf", "version": "v2.0.0" }, "type": "thesis", "title": "The Equations of State of Hydrogen and Argon: Applications to the Jovian Interior", "author": [ { "family_name": "Lagus", "given_name": "Peter Leonard", "clpid": "Lagus-Peter-Leonard" } ], "thesis_advisor": [ { "family_name": "Ahrens", "given_name": "Thomas J.", "clpid": "Ahrens-T-J" } ], "thesis_committee": [ { "family_name": "Archambeau", "given_name": "Charles B.", "clpid": "Archambeau-C-B" }, { "family_name": "Goldreich", "given_name": "Peter Martin", "clpid": "Goldreich-P-M" }, { "family_name": "Muhleman", "given_name": "Duane Owen", "clpid": "Muhleman-D-O" }, { "family_name": "Anderson", "given_name": "Donald L.", "clpid": "Anderson-D-L" }, { "family_name": "Ahrens", "given_name": "Thomas J.", "clpid": "Ahrens-T-J" } ], "local_group": [ { "literal": "div_gps" } ], "abstract": "Hugoniot data for solid argon (initially at 77\u00b0K and l bar) and solid hydrogen (initially at 5\u00b0K and 1 bar) have been obtained to 143 kbar and 6.4 kbar respectively utilizing a propellant gun. The argon data (at volumes of 15.28, 14.84 and 14.64 cc/mole) are in fair agreement with previous shock data, and in excellent agreement with recent theoretically predicted Hugoniots. The hydrogen data (at volumes of 17.10, 15.32, 15.27, 15.11 cc/mole) are compared with Hugoniots calculated from published isothermal compression data. For both argon and hydrogen, the present data are consistent with the assumption that \u03b3/V is constant. Furthermore, to compressions of V/V\u2080 \u2243 0.65, no gross inconsistencies exist between shock-wave and isothermal compression measurements in solid hydrogen.
\r\n\r\nA simple equation of state (E0S) for molecular hydrogen based on a spherically averaged De Boer-type repulsion potential which explicitly includes the zero point energy reproduces experimental pressure-volume data between 5 kbar and 370 kbar. This molecular equa\u00adtion of state when combined with recent metallic equations of state implies a molecular to metallic phase transition pressure of 1.9 \u00b1 0.4 Mbar at 0\u00b0K.
\r\n\r\nA thermally expanded model of Jupiter which incorporates this molecular equation of state, recent metallic hydrogen and helium equa\u00adtions of state, and a van der Waals-type atmosphere yields a model with a hydrogen abundance of x = 0.57. However, the interior temperatures are everywhere above the melting temperature of metallic hydrogen.
", "doi": "10.7907/zzpq-2876", "publication_date": "1974", "thesis_type": "phd", "thesis_year": "1974" }, { "id": "thesis:10683", "collection": "thesis", "collection_id": "10683", "cite_using_url": "https://resolver.caltech.edu/CaltechTHESIS:02072018-154915291", "primary_object_url": { "basename": "Davies_gf_1973.pdf", "content": "final", "filesize": 41100534, "license": "other", "mime_type": "application/pdf", "url": "/10683/1/Davies_gf_1973.pdf", "version": "v2.0.0" }, "type": "thesis", "title": "Elasticity of Solids at High Pressures and Temperatures: Theory, Measurement, and Geophysical Application", "author": [ { "family_name": "Davies", "given_name": "Geoffrey Frederick", "clpid": "Davies-Geoffrey-Frederick" } ], "thesis_advisor": [ { "family_name": "Ahrens", "given_name": "Thomas J.", "clpid": "Ahrens-T-J" } ], "thesis_committee": [ { "family_name": "Unknown", "given_name": "Unknown" } ], "local_group": [ { "literal": "div_gps" } ], "abstract": "A theory for describing the elasticity of solids at \r\nsimultaneous high pressures and high temperatures is developed\r\nby incorporating the fourth-order ahnarmonic theory of lattice\r\ndynamics into finite strain theory. The theory is applied\r\nto the analysis of a variety of data for MgO, SiO2 and\r\nNaCl, and the results for MgO and SiO2 used as the basis of\r\na discussion of the constitution of the lower mantle. New \r\nresults are reported of measurements of elastic properties\r\nof MgO shock-compressed to over 500 Kb.
\r\n\r\n\r\nThe condition that finite strain equations be frame-\r\nindifferent is shown to require that only strain tensors\r\nbelonging to a class of frame-indifferent strain tensors be \r\nused in finite strain expansion. It is shown that the\r\ngenerality of finite strain theory is not impaired by the \r\ninclusion of an explicit theory of thermal effects. Explicit\r\nequations for isotherms, isentropes and Hugoniots and for\r\nthe effective elastic moduli of materials of cubic symmetry\r\nunder hydrostatic stress are derived. The primary parameters\r\nof these equations are related to the elastic moduli and\r\ntheir pressure and temperature derivatives in an arbitrary\r\nreference state using thermodynamic identities, some of which\r\nare derived here.
\r\n\r\n\r\nHugoniot data corresponding to different initial sample\r\ndensities of MgO, SiO2 and NaCl and original ultrasonic data\r\nof NaCl are used to test both the compressional and thermal\r\nparts of the theory, and to refine the equations of state of\r\nthese materials. The frame-indifferent analogue, E, of the\r\nusual \"Eulerian\" strain tensor, \u03b5, is found to usually give\r\nfaster convergence of finite strain expansions than the\r\n\"Lagrangian\" strain tensor, \u03b7. The effect of usinq differ\u00adent \r\nstrain measures on the values of parameters derived from \r\ndata is demonstrated, and the adverse effects of using inappropriately \r\nderived parameters in extrapolation equations is demonstrated.\t\r\nThermal effects in Hugoniot data are reasonably well described, \r\nbut higher-order anharmonic effects appear to be required in the \r\ntheory in order to describe the high temperature ultrasonic and \r\nthermal expansion data.
\r\n\r\n\r\nMeasured velocities of rarefaction waves propagating\r\ninto shocked MgO are in accord with a two-stage longitudinal\r\n(elastic)-hydrodynamic (plastic) decompression model, and\r\nconstrain the high-pressure elastic moduli of MgO.
\r\n\r\n\r\nThe effects on the determination of the lower mantle constitution \r\nof temperature, varying composition, the presence of phases denser \r\nthan oxides mixtures, and the presence of iron in the \"low-spin\" \r\nelectronic state are estimated, and a trade-off between many \r\nof these factors demonstrated. Iron content could range between 6% and \r\n15% by weight of FeO. Silica content could range from 33% to 50% or \r\nmore by weight. Phases a few percent denser than oxides mixtures \r\nseem to be likely. The temperature is very indeterminable.
", "doi": "10.7907/SK7X-EE20", "publication_date": "1973", "thesis_type": "phd", "thesis_year": "1973" }, { "id": "thesis:11127", "collection": "thesis", "collection_id": "11127", "cite_using_url": "https://resolver.caltech.edu/CaltechTHESIS:07202018-110502938", "type": "thesis", "title": "Crystal Field Effects in Mantle Minerals", "author": [ { "family_name": "Gaffney", "given_name": "Edward Stowell", "clpid": "Gaffney-Edward-Stowell" } ], "thesis_advisor": [ { "family_name": "Ahrens", "given_name": "Thomas J.", "clpid": "Ahrens-T-J" }, { "family_name": "Rossman", "given_name": "George Robert", "orcid": "0000-0002-4571-6884", "clpid": "Rossman-G-R" } ], "thesis_committee": [ { "family_name": "Unknown", "given_name": "Unknown" } ], "local_group": [ { "literal": "div_gps" } ], "abstract": "The behavior of Fe2+ in the lower mantle will depend on the effects of crystal fields. A point charge model, scaled to fit observed spectra at low pressures is developed to predict these effects. Two of the three parameters needed to predict spin-pairing transitions can only be determined from spin-forbidden electronic transitions. The spectra of garnet, gillespite and peridot are examined and found to have such absorption features. Assignment of these spectra leads to values of the Racah parameters, B and C, as well as the crystal field parameter Dq.
\r\n\r\nA new experimental technique, which allows the measurement of optical absorption spectra of solids in the visible region during shock loading, is described. Results are discussed for periclase and ruby. The ruby data indicate that the point charge model is good to at least 15 percent (volume) compression.
\r\n\r\nThe effects of low-spin Fe2+ in the earth's lower mantle are investigated in considerable detail.\tThe existence of low-spin Fe2+ permits the formation of a separate phase since Mg2+ and low-spin Fe2+ may not form solid solutions. The bulk elastic behavior of such phases is predicted from volume-bulk modulus systematics and compared with available shock wave data. It is likely that the high pressure phases of several ferrous iron compounds involve low-spin Fe2+ Iron will be spin-paired in the mantle below 1200 km and likely at higher levels as well. The observed density and bulk modulus in the lower mantle are inconsistent with any combination of phases in a pyrolite bulk composition but can be fit quite well by a model with all Fe2+ spin-paired below 630 km and nearly olivine composition at that depth, with MgO decreasing to almost a pyroxene composition at the core.
\r\n\r\nAn origin of the upper mantle from the lower mantle by chemical fractionation is proposed. The spin-pairing of Fe2+ provides an excellent mechanism for both iron and silicon enrichment in the lower mantle by partial melting yielding a pyrolite upper mantle, and hence, a chemically inhomogeneous mantle. This removes the motivation for reducing FeO and SiO2 in the mantle to supply Fe and Si for the core.
\r\n\r\n", "doi": "10.7907/3DS0-XR25", "publication_date": "1973", "thesis_type": "phd", "thesis_year": "1973" } ]