Earth, Environmental and Planetary Sciences
Permanent URI for this community
Browse
Browsing Earth, Environmental and Planetary Sciences by Title
Now showing 1 - 20 of 286
Results Per Page
Sort Options
Item 4500-year paleohurricane record from the Western Gulf of Mexico, Coastal Central TX, USA(Elsevier, 2024) Monica, Sarah B.; Wallace, Davin J.; Wallace, Elizabeth J.; Du, Xiaojing; Dee, Sylvia G.; Anderson, John B.Texas receives the second-highest number of tropical cyclone (TC) landfalls per year in the United States. At present, long-term TC projections from climate models remain uncertain due to the short and biased nature of Atlantic TC observations. Sediment archives of past storms can help extend the observational record of TC strikes over the past few millennia. When a TC makes landfall along the central Texas coast, coastal downwelling channels and storm currents transport and deposit coarse sediment to a zone of rapid accumulation along the shelf, known as the Texas Mud Blanket (TMB). This “backwash” process results in expansive storm deposits along the shelf, making this region ideal for paleotempestological reconstructions. Here, we present two sediment cores, located approximately 6 km southeast of Matagorda Island (TX), that collectively yield a ∼4500-year paleohurricane record. 210Pb and 137Cs are utilized in conjunction with radiocarbon ages to produce high-resolution Bayesian age models. One-centimeter interval grain size analyses are used to identify TC deposits. Two-centimeter interval X-Ray Fluorescence (XRF) is used as an additional measure to verify depositional mechanisms in this shelf environment. We define an intense paleohurricane event threshold through statistical analysis of mean grain size data. The sediment-derived TC record is correlated to Palmer Drought Severity Index (PDSI) data from Paleo Hydrodynamics Data Assimilation (PHYDA) to bolster our interpretation of the TC record, revealing a coupled relationship between PDSI and TCs since ∼300 yr BP. Over the ∼4500-year period, 24 intense TCs were recorded in the sediment record, yielding a long-term annual landfall probability of ∼0.53%. Additionally, comparisons between other TC records within the Atlantic establish a relationship between enhanced TC activity in the Western Gulf of Mexico (GOM) and TCs formed in the Caribbean Sea.Item A 5-km-thick reservoir with > 380,000 km3 of magma within the ancient Earth's crust(Springer Nature, 2022) Latypov, Rais; Chistyakova, Sofya; Hornsey, Richard A.; Costin, Gelu; van der Merwe, MauritzSeveral recent studies have argued that large, long-lived and molten magma chambers may not occur in the shallow Earth’s crust. Here we present, however, field-based observations from the Bushveld Complex that provide evidence to the contrary. In the eastern part of the complex, the magmatic layering continuously drapes across a ~ 4-km-high sloping step in the chamber floor. Such deposition of magmatic layering implies that the resident melt column was thicker than the stepped relief of the chamber floor. Prolonged internal differentiation within this thick magma column is further supported by evolutionary trends in crystallization sequence and mineral compositions through the sequence. The resident melt column in the Bushveld chamber during this period is estimated at > 5-km in thickness and > 380,000 km3 in volume. This volume of magma is three orders of magnitude larger than any known super-eruption in the Earth’s history and is only comparable to the extrusive volumes of some of Earth’s large igneous provinces. This suggests that super-large, entirely molten, and long-lived magma chambers occur, at least occasionally, in the geological history of our planet. Therefore, the classical view of magma chambers as ‘big magma tanks’ remains a viable research concept for some of Earth’s magmatic provinces.Item A community-based geological reconstruction of Antarctic Ice Sheet deglaciation since the Last Glacial Maximum(Elsevier, 2014) The RAISED Consortium; Bentley, Michael J.; Ó Cofaigh, Colm; Anderson, John B.; Conway, Howard; Davies, Bethan; Graham, Alastair G.C.; Hillenbrand, Claus-Dieter; Hodgson, Dominic A.; Jamieson, Stewart S.R.; Larter, Robert D.; Mackintosh, Andrew; Smith, James A.; Verleyen, Elie; Ackert, Robert P.; Bart, Philip J.; Berg, Sonja; Brunstein, Daniel; Canals, Miquel; Colhoun, Eric A.; Crosta, Xavier; Dickens, William A.; Domack, Eugene; Dowdeswell, Julian A.; Dunbar, Robert; Ehrmann, Werner; Evans, Jeffrey; Favier, Vincent; Fink, David; Fogwill, Christopher J.; Glasser, Neil F.; Gohl, Karsten; Golledge, Nicholas R.; Goodwin, Ian; Gore, Damian B.; Greenwood, Sarah L.; Hall, Brenda L.; Hall, Kevin; Hedding, David W.; Hein, Andrew S.; Hocking, Emma P.; Jakobsson, Martin; Johnson, Joanne S.; Jomelli, Vincent; Jones, R. Selwyn; Klages, Johann P.; Kristoffersen, Yngve; Kuhn, Gerhard; Leventer, Amy; Licht, Kathy; Lilly, Katherine; Lindow, Julia; Livingstone, Stephen J.; Massé, Guillaume; McGlone, Matt S.; McKay, Robert M.; Melles, Martin; Miura, Hideki; Mulvaney, Robert; Nel, Werner; Nitsche, Frank O.; O'Brien, Philip E.; Post, Alexandra L.; Roberts, Stephen J.; Saunders, Krystyna M.; Selkirk, Patricia M.; Simms, Alexander R.; Spiegel, Cornelia; Stolldorf, Travis D.; Sugden, David E.; van der Putten, Nathalie; van Ommen, Tas; Verfaillie, Deborah; Vyverman, Wim; Wagner, Bernd; White, Duanne A.; Witus, Alexandra E.; Zwartz, DanA robust understanding of Antarctic Ice Sheet deglacial history since the Last Glacial Maximum is important in order to constrain ice sheet and glacial-isostatic adjustment models, and to explore the forcing mechanisms responsible for ice sheet retreat. Such understanding can be derived from a broad range of geological and glaciological datasets and recent decades have seen an upsurge in such data gathering around the continent and Sub-Antarctic islands. Here, we report a new synthesis of those datasets, based on an accompanying series of reviews of the geological data, organised by sector. We present a series of timeslice maps for 20 ka, 15 ka, 10 ka and 5 ka, including grounding line position and ice sheet thickness changes, along with a clear assessment of levels of confidence. The reconstruction shows that the Antarctic Ice sheet did not everywhere reach the continental shelf edge at its maximum, that initial retreat was asynchronous, and that the spatial pattern of deglaciation was highly variable, particularly on the inner shelf. The deglacial reconstruction is consistent with a moderate overall excess ice volume and with a relatively small Antarctic contribution to meltwater pulse 1a. We discuss key areas of uncertainty both around the continent and by time interval, and we highlight potential priorities for future work. The synthesis is intended to be a resource for the modelling and glacial geological community.Item A Displaced Lower Mantle Source of the Hainan Plume in South China Revealed by Receiver Function Imaging of the CEArray(Wiley, 2024) Zhang, Yan; Niu, Fenglin; Ning, Jieyuan; Yu, ChunquanWe analyzed 49,592 teleseismic receiver functions (RFs) recorded by 278 CEArray stations to image the mantle transition zone (MTZ) beneath the South China Block to understand the origins of deep velocity anomalies and their potential links to subduction and intraplate volcanism. We employed a fast-marching method and a high-resolution 3-D velocity model (FWEA18) derived from full waveform inversion in computing P-to-S conversion times to better image the 410- and 660-km discontinuities. Our results indicate that the common-conversion-point stacking of RFs using 3-D conversion times yielded better migration images of the two discontinuities. The images revealed a slightly depressed 410-km with a few small uplifted patches, and showed that the 660-km beneath the western Yangtze Craton is depressed by 10–25 km, which is likely caused by the stagnant Paleo-Pacific slab. The 660-km beneath the southern Cathaysia Block has a 5–15 km high plateau with a topographic low at its central part. The lateral dimension of the topographic low is ∼150 km and is located beneath the central Pearl River Mount Basin near Hong Kong. We speculate that the topographic low occurs within the Hainan plume with a temperature excess of ∼300–400 K and is caused by the garnet phase transition. The displaced deep plume enters the MTZ and spreads nearly horizontally at the base. The plume evolves into two channels with a minor one toward the northeast and a major one toward the southwest, which keep moving upward to the 410-km. The southwest channel is likely the source that feeds the Hainan volcanoes.Item A Framework for the Systematic Selection of Biosensor Chassis for Environmental Synthetic Biology(American Chemical Society, 2022) Sridhar, Swetha; Ajo-Franklin, Caroline M.; Masiello, Caroline A.Microbial biosensors sense and report exposures to stimuli, thereby facilitating our understanding of environmental processes. Successful design and deployment of biosensors hinge on the persistence of the microbial host of the genetic circuit, termed the chassis. However, model chassis organisms may persist poorly in environmental conditions. In contrast, non-model organisms persist better in environmental conditions but are limited by other challenges, such as genetic intractability and part unavailability. Here we identify ecological, metabolic, and genetic constraints for chassis development and propose a conceptual framework for the systematic selection of environmental biosensor chassis. We identify key challenges with using current model chassis and delineate major points of conflict in choosing the most suitable organisms as chassis for environmental biosensing. This framework provides a way forward in the selection of biosensor chassis for environmental synthetic biology.Item A Framework for Understanding Whole-Earth Carbon Cycling(Cambridge University Press, 2019) Lee, Cin-Ty A.; Jiang, Hehe; Dasgupta, Rajdeep; Torres, Mark; Orcutt, Beth N.; Daniel, Isabelle; Dasgupta, RajdeepBuilding from the base of knowledge presented in the preceding chapters, this chapter explores how the cycling of carbon in subduction zones and orogenic belts varies with supercontinent cycles and mountain building. It discusses the processes that link short-term and long-term carbon cycling and the timescales of these processes, such as the response times of weathering and atmospheric drawdown at periods of enhanced volcanism. This chapter covers topics of potential fluctuations in the long-term CO2 content of Earth’s atmosphere because of mantle–climate feedback, again taking advantage of the modeling platforms available for further exploration of these topics.Item A new bubble dynamics model to study bubble growth, deformation, and coalescence(American Geophysical Union, 2014) Huber, C.; Su, Y.; Nguyen, C.T.; Parmigiani, A.; Gonnermann, H.M.; Dufek, J.We propose a new bubble dynamics model to study the evolution of a suspension of bubbles over a wide range of vesicularity, and that accounts for hydrodynamical interactions between bubbles while they grow, deform under shear flow conditions, and exchange mass by diffusion coarsening. The model is based on a lattice Boltzmann method for free surface flows. As such, it assumes an infinite viscosity contrast between the exsolved volatiles and the melt. Our model allows for coalescence when two bubbles approach each other because of growth or deformation. The parameter (disjoining pressure) that controls the coalescence efficiency, i.e., drainage time for the fluid film between the bubbles, can be set arbitrarily in our calculations. We calibrated this parameter by matching the measured time for the drainage of the melt film across a range of Bond numbers (ratio of buoyancy to surface tension stresses) with laboratory experiments of a bubble rising to a free surface. The model is then used successfully to model Ostwald ripening and bubble deformation under simple shear flow conditions. The results we obtain for the deformation of a single bubble are in excellent agreement with previous experimental and theoretical studies. For a suspension, we observe that the collective effect of bubbles is different depending on the relative magnitude of viscous and interfacial stresses (capillary number). At low capillary number, we find that bubbles deform more readily in a suspension than for the case of a single bubble, whereas the opposite is observed at high capillary number.Item A novel and efficient engine for P-/S-wave-mode vector decomposition for vertical transverse isotropic elastic reverse time migration(Society of Exploration Geophysicists, 2022) Zhang, Lele; Liu, Lu; Niu, Fenglin; Zuo, Jiahui; Shuai, Da; Jia, Wanli; Zhao, YangWave-mode decomposition plays a very important role in elastic reverse time migration (ERTM). Improved imaging quality can be achieved due to reduced wave-mode crosstalk artifacts. The current state-of-the-art methods for anisotropic wavefield separation are based on either splitting model strategy, low-rank approximation, or lower-upper (LU) factorization. Most of these involve expensive matrix computation and Fourier transforms with strong model assumptions. Based on the anisotropic-Helmholtz (ani-Helmholtz) decomposition operator and decoupled formulations, we develop a novel and efficient P-/S-wave-mode vector decomposition method in vertical transverse isotropic (VTI) media with application in ERTM. We first review the basics of classical Helmholtz decomposition and isotropic decoupled formulations. In addition, the ani-Helmholtz decomposition operator is built from the P- and S-wave polarizations of the Christoffel equation in VTI media. We then derive novel decoupled formulations of anisotropic P-/S-waves based on the obtained ani-Helmholtz operator. Moreover, we use the first-order Taylor expansion to approximate the normalization term from the derived decoupled formulations and obtain an efficient ani-Helmholtz decomposition approach, which produces vector P- and S-wavefields with correct units, phases, and amplitudes. Compared with the previous studies, our approach mitigates model assumptions, avoids intricate calculations, such as LU factorization or low-rank approximation, and only needs three fast Fourier transforms at each time step. In addition, the graphic processing unit technique is used to dramatically accelerate various functions of ERTM, such as wavefields extrapolation, decomposition, and imaging. Three synthetic examples demonstrate the effectiveness and feasibility of our proposed approach.Item A perspective on climate model hierarchies(Wiley, 2017) Jeevanjee, Nadir; Hassanzadeh, Pedram; Hill, Spencer; Sheshadri, AditiTo understand Earth's climate, climate modelers employ a hierarchy of climate models spanning a wide spectrum of complexity and comprehensiveness. This essay, inspired by the World Climate Research Programme's recent “Model Hierarchies Workshop,” attempts to survey and synthesize some of the current thinking on climate model hierarchies, especially as presented at the workshop. We give a few formal descriptions of the hierarchy and survey the various ways it is used to generate, test, and confirm hypotheses. We also discuss some of the pitfalls of contemporary climate modeling, and how the “elegance” advocated for by Held (2005) has (and has not) been used to address them. We conclude with a survey of current activity in hierarchical modeling, and offer suggestions for its continued fruitful development.Item A Pliocene Precipitation Isotope Proxy-Model Comparison Assessing the Hydrological Fingerprints of Sea Surface Temperature Gradients(Wiley, 2022) Knapp, Scott; Burls, Natalie J.; Dee, Sylvia; Feng, Ran; Feakins, Sarah J.; Bhattacharya, TriptiThe Pliocene offers insights into future climate, with near-modern atmospheric pCO2 and global mean surface temperature estimated to be 3–4°C above pre-industrial. However, the hydrological response differs between future global warming and early Pliocene climate model simulations. This discrepancy results from the use of reduced meridional and zonal sea surface temperature (SST) gradients, based on foraminiferal Mg/Ca and Alkenone proxy evidence, to force the early Pliocene simulation. Subsequent, SST reconstructions based on the organic proxy TEX86, have found warmer temperatures in the warm pool, bringing the magnitude of the gradient reductions into dispute. We design an independent test of Pliocene SST scenarios and their hydrological cycle “fingerprints.” We use an isotope-enabled General Circulation Model, iCAM5, to model the distribution of water isotopes in precipitation in response to four climatological SST and sea-ice fields representing modern, abrupt 4 × CO2, late Pliocene and early Pliocene climates. We conduct a proxy-model comparison with all the available precipitation isotope proxy data, and we identify target regions that carry precipitation isotopic fingerprints of SST gradients as priorities for additional proxy reconstructions. We identify two regions with distinct precipitation isotope (D/H) fingerprints resulting from reduced SST gradients: the Maritime Continent (D-enriched due to reduced convective rainfall) and the Sahel (wetter, more deep convection, D-depleted). The proxy-model comparison using available plant wax reconstructions, mostly from Africa, is promising but inconclusive. Additional proxy reconstructions are needed in both target regions and in much of the world for significant tests of SST scenarios and dynamical linkages to the hydrological cycle.Item A subglacial hydrologic drainage hypothesis for silt sorting and deposition during retreat in Pine Island Bay(Cambridge University Press, 2019) Schroeder, Dustin M.; MacKie, Emma J.; Creyts, Timothy T.; Anderson, John B.Late Holocene sediment deposits in Pine Island Bay, West Antarctica, are hypothesized to be linked to intensive meltwater drainage during the retreat of the paleo-Pine Island Ice Stream after the Last Glacial Maximum. The uppermost sediment units show an abrupt transition from ice-proximal debris to a draped silt during the late Holocene, which is interpreted to coincide with rapid deglaciation. The small scale and fine sorting of the upper unit could be attributed to origins in subglacial meltwater; however the thickness and deposition rate for this unit imply punctuated- rather than continuous-deposition. This, combined with the deposit's location seaward of large, bedrock basins, has led to the interpretation of this unit as the result of subglacial lake outbursts in these basins. However, the fine-scale sorting of the silt unit is problematic for this energetic interpretation, which should mobilize and deposit a wider range of sediment sizes. To resolve this discrepancy, we present an alternative mechanism in which the silt was sorted by a distributed subglacial water system, stored in bedrock basins far inland of the grounding line, and subsequently eroded at higher flow speeds during retreat. We demonstrate that this mechanism is physically plausible given the subglacial conditions during the late Holocene. We hypothesize that similar silt units observed elsewhere in Antarctica downstream of bedrock basins could be the result of the same mechanism.Item A zero-dimensional view of atmospheric degradation of levoglucosan (LEVCHEM_v1) using numerical chamber simulations(Copernicus Publications, 2021) Suciu, Loredana G.; Griffin, Robert J.; Masiello, Caroline A.Here, we developed a zero-dimensional (0-D) modeling framework (LEVCHEM_v1) to provide insights into the atmospheric degradation of a key tracer emitted during biomass burning – levoglucosan (LEV), while additionally exploring its effects on the dynamics of secondary organic aerosols (SOA) and other gases. For this, we updated existing chemical mechanisms (homogeneous gas-phase chemistry and heterogeneous chemistry) in the BOXMOXv1.7 model to include the chemical degradation of LEV and its intermediary degradation products in both phases (gas and aerosol). In addition, we added a gas-particle partitioning mechanism to the model to account for the effect of evaporation and condensation on the phase-specific concentrations of LEV and its degradation products. Comparison of simulation results with measurements from various chamber experiments (spanning summer and winter conditions) show that the degradation timescale of LEV varied by phase, with gas-phase degradation occurring over ∼1.5–5 d and aerosol-phase degradation occurring over ∼8–36 h. These relatively short timescales suggest that most of the initial LEV concentration can be lost chemically or deposited locally before being transported regionally. We varied the heterogeneous reaction rate constant in a sensitivity analysis (for summer conditions only) and found that longer degradation timescales of LEV are possible, particularly in the aerosol phase (7 d), implying that some LEV may be transported regionally. The multiphase chemical degradation of LEV has effects on SOA and other gases. Several first- or second-generation products resulted from its degradation; most of the products include one or two carbonyl groups, one product contains a nitrate group, and a few products show the cleavage of C−C bonds. The relative importance of the products varies depending on the phase and the timing of the maximum concentration achieved during the simulation. Our estimated secondary organic aerosol SOA yields (4 %–32 %) reveal that conversion of LEV to secondary products is significant and occurs rapidly in the studied scenarios. LEV degradation affected other gases by increasing the concentrations of radicals and decreasing those of reactive nitrogen species. Decreases of the mixing ratios of nitrogen oxides appear to drive a more rapid increase in ozone compared with changes in volatile organic compounds levels. An important next step to confirm longer degradation timescales will be to extend the evaluation of the modeled LEV degradation beyond 3–6 h by using more extensive data from chambers and, possibly, from fire plumes. The mechanism developed here can be used in chemical transport models applied to fire plumes to trace LEV and its degradation products from source to deposition, to assess their atmospheric implications and to answer questions relevant to fire tracing, carbon and nitrogen cycling, and climate.Item Absolute plate velocities from seismic anisotropy: Importance of correlated errors(American Geophysical Union, 2014) Zheng, Lin; Gordon, Richard G.; Kreemer, CornéThe errors in plate motion azimuths inferred from shear wave splitting beneath any one tectonic plate are shown to be correlated with the errors of other azimuths from the same plate. To account for these correlations, we adopt a two-tier analysis: First, find the pole of rotation and confidence limits for each plate individually. Second, solve for the best fit to these poles while constraining relative plate angular velocities to consistency with the MORVEL relative plate angular velocities. Our preferred set of angular velocities, SKS-MORVEL, is determined from the poles from eight plates weighted proportionally to the root-mean-square velocity of each plate. SKS-MORVEL indicates that eight plates (Amur, Antarctica, Caribbean, Eurasia, Lwandle, Somalia, Sundaland, and Yangtze) have angular velocities that differ insignificantly from zero. The net rotation of the lithosphere is 0.25 ± 0.11° Ma−1 (95% confidence limits) right handed about 57.1°S, 68.6°E. The within-plate dispersion of seismic anisotropy for oceanic lithosphere (σ = 19.2°) differs insignificantly from that for continental lithosphere (σ = 21.6°). The between-plate dispersion, however, is significantly smaller for oceanic lithosphere (σ = 7.4°) than for continental lithosphere (σ = 14.7°). Two of the slowest-moving plates, Antarctica (vRMS = 4 mm a−1, σ = 29°) and Eurasia (vRMS = 3 mm a−1, σ = 33°), have two of the largest within-plate dispersions, which may indicate that a plate must move faster than ≈ 5 mm a−1 to result in seismic anisotropy useful for estimating plate motion. The tendency of observed azimuths on the Arabia plate to be counterclockwise of plate motion may provide information about the direction and amplitude of superposed asthenospheric flow or about anisotropy in the lithospheric mantle.Item Accretion and Uneven Depletion of the Main Asteroid Belt(IOP Publishing, 2024) Deienno, Rogerio; Nesvorný, David; Clement, Matthew S.; Bottke, William F.; Izidoro, André; Walsh, Kevin J.Accretion and Uneven Depletion of the Main Asteroid Belt, Rogerio Deienno, David Nesvorný, Matthew S. Clement, William F. Bottke, André Izidoro, Kevin J. WalshItem Adjoint traveltime tomography unravels a scenario of horizontal mantle flow beneath the North China craton(Springer Nature, 2021) Dong, Xingpeng; Yang, Dinghui; Niu, Fenglin; Liu, Shaolin; Tong, PingThe North China craton (NCC) was dominated by tectonic extension from late Cretaceous to Cenozoic, yet seismic studies on the relationship between crust extension and lithospheric mantle deformation are scarce. Here we present a three dimensional radially anisotropic model of NCC derived from adjoint traveltime tomography to address this issue. We find a prominent low S-wave velocity anomaly at lithospheric mantle depths beneath the Taihang Mountains, which extends eastward with a gradually decreasing amplitude. The horizontally elongated low-velocity anomaly is also featured by a distinctive positive radial anisotropy (VSH > VSV). Combining geodetic and other seismic measurements, we speculate the presence of a horizontal mantle flow beneath central and eastern NCC, which led to the extension of the overlying crust. We suggest that the rollback of Western Pacific slab likely played a pivotal role in generating the horizontal mantle flow at lithospheric depth beneath the central and eastern NCC.Item Aged black carbon in marine sediments and sinking particles(Wiley, 2014) Coppola, Alysha I.; Ziolkowski, Lori A.; Masiello, Caroline A.; Druffel, Ellen R.M.We report measurements of oceanic black carbon (BC) to determine the sources of BC to abyssal marine sediments in the northeast Pacific Ocean. We find that the average 14C age of BC is older (by 6200 ± 2200 14C years) than that of the concurrently deposited non-BC sedimentary organic carbon. We investigate sources of aged BC to sediments by measuring a sample of sinking particulate organic carbon (POC) and find that POC may provide the main transport mechanism of BC to sediments. We suggest that aged BC is incorporated into POC from a combination of resuspended sediments and sorption of ancient dissolved organic carbon BC onto POC. Our BC flux estimate represents ~8–16% of the global burial flux of organic carbon to abyssal sediments and constitutes a minimum long-term removal estimate of 6–32% of biomass-derived BC using the present day emission flux.Item Airborne radar sounding evidence for deformable sediments and outcropping bedrock beneath Thwaites Glacier, West Antarctica(Wiley, 2014) Schroeder, Dustin M.; Blankenship, Donald D.; Young, Duncan A.; Witus, Alexandra E.; Anderson, John B.The geologic and morphologic records of prior ice sheet configurations show evidence of rapid, back-stepping, meltwater intensive retreats. However, the potential for such a retreat in a contemporary glacier depends on the lithology of the current ice sheet bed, which lies beneath kilometers of ice, making its physical properties difficult to constrain. We use radar sounding and marine bathymetry data to compare the bed configuration of Thwaites Glacier to the bed of paleo-Pine Island Glacier. Using observed and modeled radar scattering, we show that the tributaries and upper trunk of Thwaites Glacier are underlain by ice flow-aligned bedforms consistent with deformable sediment and that the lower trunk is grounded on a region of high bed roughness consistent with outcropping bedrock. This is the same configuration as paleo-Pine Island Glacier during its retreat across the inner continental shelf.Item Ammonium-Bearing Fluorapophyllite-(K) in the Magnesian Skarns from Aleului Valley, Pietroasa, Romania(MDPI, 2023) Marincea, Ştefan; Dumitraş, Delia-Georgeta; Sava Ghineţ, Cristina; Filiuţă, Andra Elena; Dal Bo, Fabrice; Hatert, Frédéric; Costin, GeluAn ammonium-bearing fluorapophyllite-(K) occurs as a late hydrothermal product in the outer endoskarn zone from Aleului Valley (N 46°37′04″, E 22°35′22″), located at the contact of the granodiorite laccolith from Pietroasa, of Upper Cretaceous age, with Anisian dolostones. Associated minerals are wollastonite, K feldspar, diopside, fluorapatite, talc, and pectolite. The chemical structural formula is [K0.985Na0.012(NH4)0.076]Σ=1.073(Ca4.009Mn0.001Fe2+0.003Mg0.002Ba0.001)Σ=4.016(Si7.953Al0.047) O20.029[F0.899(OH)0.101]·8.059H2O. The structure was successfully refined as tetragonal, space group P4/mnc, with cell parameters of a = 8.9685(1) Å and c = 15.7885(5) Å. The indices of refraction are ω = 1.534(1) and ε = 1.536(1). The calculated density is Dx = 2.381 g/cm3, in good agreement with the measured density, Dm = 2.379(4) g/cm3. The thermal analysis shows that the mineral completely dehydrates at up to 450 °C (endothermic effects at 330, 371, and 448 °C) and loses ammonium at 634 °C. In the infrared spectra, the multiplicity of the bands assumed to be silicate modes (1ν1 + 3ν3 + 2ν2 + 3ν4) agrees with the reduction in the symmetry of the SiO44− ion from Td to Cδ. Fluorapophyllite-(K) from Aleului Valley is of late hydrothermal origin and crystallized from F-rich fluids originating from the granodiorite intrusion, which mobilized K, Ca, and Si from the pre-existing feldspar.Item An alternative formula for approximate extended Born inversion(Society of Exploration Geophysicists, 2017) Hou, Jie; Symes, William W.Various modifications of reverse time migration (RTM) provide asymptotic inverses to the subsurface offset extended Born modeling operator for constant-density acoustics. These approximate inverses have the same quality (asymptotic accuracy) as do generalized Radon transform pseudoinverses, but they can be computed without any ray tracing whatsoever. We have developed an approximate inverse of this type whose additional computational cost, above that of subsurface offset extended RTM, is negligible.Item An explicit time evolution method for acoustic wave propagation(Society of Exploration Geophysicists, 2014) Liu, Huafeng; Dai, Nanxun; Niu, Fenglin; Wu, WeiCost-effective waveform modeling is the key to practical reverse time migration (RTM) and full-waveform inversion (FWI) implementations. We evaluated an explicit time evolution (ETE) method to efficiently simulate wave propagation in acoustic media with high temporal accuracy. We started from the constant-density acoustic wave equation and obtained an analytical time-marching scheme in the wave number domain. We then formulated an ETE scheme in the time-space domain by introducing a cosine function approximation. Although the ETE operator appears to be similar to the second-order temporal finite-difference (FD) operator, the exact nature of the ETE formula ensures high accuracy in time. We further introduced a set of optimum stencils and coefficients by minimizing evolution errors in a least-squares sense. Our numerical tests indicated that ETE can achieve similar waveform accuracy as FD with four times larger time steps. Meanwhile, the compact ETE operator keeps the computation efficient. The efficiency and capability to handle complex velocity field make ETE an attractive engine in acoustic RTM and FWI.