Showing 1–50 of 84 results for author: McClure, M

JWST study of the DG Tau B disk wind candidate: I -- Overview and Nested H$_2$/CO outflows

Authors: Valentin Delabrosse, Catherine Dougados, Sylvie Cabrit, Benoit Tabone, Lukasz Tychoniec, Tom Ray, Linda Podio, Melissa McClure

Abstract: The origin and impact of outflows on proto-planetary disks and planet formation are key open questions. DG Tau B, a Class I protostar with a structured disk and a striking rotating conical CO outflow, recently identified with ALMA as one of the best MHD disk wind candidate, is an ideal target for studying these phenomena. Our aim is to analyse the outflow components intermediate between the fast a… ▽ More The origin and impact of outflows on proto-planetary disks and planet formation are key open questions. DG Tau B, a Class I protostar with a structured disk and a striking rotating conical CO outflow, recently identified with ALMA as one of the best MHD disk wind candidate, is an ideal target for studying these phenomena. Our aim is to analyse the outflow components intermediate between the fast axial jet and the wider molecular CO outflow to discriminate between the different scenarios at their origin (irradiated/shocked disk wind or swept-up material). Using observations from JWST NIRSpec-IFU, NIRCam and SINFONI/VLT, we investigate the morphology, kinematics and excitation conditions of H$_2$ emission lines of the red-shifted outflow lobe. We find an onion-like structure of the outflows with increasing temperature, velocity and collimation towards the flow axis. The red-shifted H$_2$ emission reveals a narrow conical cavity nested inside the CO outflow and originating from the inner disk regions (< 6 au). The H$_2$ shell exhibits a constant vertical velocity ($\simeq$22 km/s), twice faster that of the CO flow and an average mass flux of $\dot{M}$(H$_2$) = 3e-11 M$_\odot$/yr significantly lower than the jet and CO values, suggesting low H$_2$ abundance. The global layered structure of the H$_2$/CO outflows is consistent with an MHD disk wind scenario, with the hot H$_2$ possibly tracing an inner dense photodissociation layer of the wind coming from a launching radius in the disk of 0.2-0.4 au. Further analysis, including MIRI observations will provide additional insights into the H$_2$ excitation mechanisms and the origin of the layered outflows observed in DG Tau B. △ Less

Submitted 28 March, 2024; originally announced March 2024.

Comments: 23 pages, 21 figures, 2 tables, submitted and reviewed to A&A

JWST MIRI MRS Images Disk Winds, Water, and CO in an Edge-On Protoplanetary Disk

Authors: Nicole Arulanantham, M. K. McClure, Klaus Pontoppidan, Tracy L. Beck, J. A. Sturm, D. Harsono, A. C. A. Boogert, M. Cordiner, E. Dartois, M. N. Drozdovskaya, C. Espaillat, G. J. Melnick, J. A. Noble, M. E. Palumbo, Y. J. Pendleton, H. Terada, E. F. van Dishoeck

Abstract: We present JWST MIRI MRS observations of the edge-on protoplanetary disk around the young sub-solar mass star Tau 042021, acquired as part of the Cycle 1 GO program "Mapping Inclined Disk Astrochemical Signatures (MIDAS)." These data resolve the mid-IR spatial distributions of H$_2$, revealing X-shaped emission extending to ~200 au above the disk midplane with a semi-opening angle of $35 \pm 5$ de… ▽ More We present JWST MIRI MRS observations of the edge-on protoplanetary disk around the young sub-solar mass star Tau 042021, acquired as part of the Cycle 1 GO program "Mapping Inclined Disk Astrochemical Signatures (MIDAS)." These data resolve the mid-IR spatial distributions of H$_2$, revealing X-shaped emission extending to ~200 au above the disk midplane with a semi-opening angle of $35 \pm 5$ degrees. We do not velocity-resolve the gas in the spectral images, but the measured semi-opening angle of the H$_2$ is consistent with an MHD wind origin. A collimated, bipolar jet is seen in forbidden emission lines from [Ne II], [Ne III], [Ni II], [Fe II], [Ar II], and [S III]. Extended H$_2$O and CO emission lines are also detected, reaching diameters between ~90 and 190 au, respectively. Hot molecular emission is not expected at such radii, and we interpret its extended spatial distribution as scattering of inner disk molecular emission by dust grains in the outer disk surface. H I recombination lines, characteristic of inner disk accretion shocks, are similarly extended, and are likely also scattered light from the innermost star-disk interface. Finally, we detect extended PAH emission at 11.3 microns co-spatial with the scattered light continuum, making this the first low-mass T Tauri star around which extended PAHs have been confirmed, to our knowledge. MIRI MRS line images of edge-on disks provide an unprecedented window into the outflow, accretion, and scattering processes within protoplanetary disks, allowing us to constrain the disk lifetimes and accretion and mass loss mechanisms. △ Less

Submitted 20 March, 2024; v1 submitted 19 February, 2024; originally announced February 2024.

Comments: Accepted for publication in ApJL on March 13th, 2024

JWST Detects Neon Line Variability in a Protoplanetary Disk

Authors: C. C. Espaillat, T. Thanathibodee, C. V. Pittman, J. A. Sturm, M. K. McClure, N. Calvet, F. M. Walter, R. Franco-Hernandez, J. Muzerolle Page

Abstract: We report the first detection of variability in the mid-infrared neon line emission of a protoplanetary disk by comparing a JWST MIRI MRS spectrum of SZ Cha taken in 2023 with a Spitzer IRS SH spectrum of this object from 2008. We measure the [Ne III]-to-[Ne II] line flux ratio, which is a diagnostic of the high-energy radiation field, to distinguish between the dominance of EUV- or X-ray-driven d… ▽ More We report the first detection of variability in the mid-infrared neon line emission of a protoplanetary disk by comparing a JWST MIRI MRS spectrum of SZ Cha taken in 2023 with a Spitzer IRS SH spectrum of this object from 2008. We measure the [Ne III]-to-[Ne II] line flux ratio, which is a diagnostic of the high-energy radiation field, to distinguish between the dominance of EUV- or X-ray-driven disk photoevaporation. We find that the [Ne III]-to-[Ne II] line flux ratio changes significantly from $\sim1.4$ in 2008 to $\sim0.2$ in 2023. This points to a switch from EUV-dominated to X-ray-dominated photoevaporation of the disk. We present contemporaneous ground-based optical spectra of the Halpha emission line that show the presence of a strong wind in 2023. We propose that this strong wind prevents EUV radiation from reaching the disk surface while the X-rays permeate the wind and irradiate the disk. We speculate that at the time of the Spitzer observations, the wind was suppressed and EUV radiation reached the disk. These observations confirm that the MIR neon emission lines are sensitive to changes in high-energy radiation reaching the disk surface. This highlights the [Ne III]-to-[Ne II] line flux ratio as a tool to gauge the efficiency of disk photoevaporation in order to provide constraints on the planet-formation timescale. However, multiwavelength observations are crucial to interpret the observations and properly consider the star-disk connection. △ Less

Submitted 13 November, 2023; originally announced November 2023.

Comments: Accepted to ApJL

StageOpt technical write-up

Authors: Egor Dontsov, Mark McClure

Abstract: This document provides technical background for StageOpt. This tool helps to better understand fluid and proppant distribution between perforations and clusters within a fracturing stage. Limited entry technique is often used to achieve uniform fluid distribution between clusters. This, however, does not automatically guarantee uniform proppant distribution. In view of this observation, StageOpt a… ▽ More This document provides technical background for StageOpt. This tool helps to better understand fluid and proppant distribution between perforations and clusters within a fracturing stage. Limited entry technique is often used to achieve uniform fluid distribution between clusters. This, however, does not automatically guarantee uniform proppant distribution. In view of this observation, StageOpt allows users to overcome this obstacle and to better design perforation design by focusing on both fluid and proppant distribution between clusters. The algorithm in StageOpt is based primarily on the paper Dontsov (2023), but it also incorporates practically important features, such as the effect of stress shadow from previous stage, perforation erosion, and other phenomena. In this write-up, flow distribution between perforations is described first, and the effects of perforation friction, near-wellbore pressure drop, as well as perforation break down and stress shadow are included. Then, proppant distribution between perforations is discussed as a function of perforation orientation. Finally, perforation erosion is added to the model. The erosion changes perforation diameter as a function of several parameters, including the flow rate through the perforation as well as the amount of proppant that flows through it. The change of perforation size in turn changes magnitude of perforation friction, which then changes the distribution of slurry, that consequently changes proppant flow distribution. Such a coupled system often leads to complex results that are hard to deduce without systematic physical modeling, which is exactly what StageOpt provides. △ Less

Submitted 9 April, 2024; v1 submitted 2 November, 2023; originally announced November 2023.

A JWST inventory of protoplanetary disk ices: The edge-on protoplanetary disk HH 48 NE, seen with the Ice Age ERS program

Authors: J. A. Sturm, M. K. McClure, T. L. Beck, D. Harsono, J. B. Bergner, E. Dartois, A. C. A. Boogert, J. E. Chiar, M. A. Cordiner, M. N. Drozdovskaya, S. Ioppolo, C. J. Law, H. Linnartz, D. C. Lis, G. J. Melnick, B. A. McGuire, J. A. Noble, K. I. Öberg, M. E. Palumbo, Y. J. Pendleton, G. Perotti, K. M. Pontoppidan, D. Qasim, W. R. M. Rocha, H. Terada , et al. (2 additional authors not shown)

Abstract: Ices are the main carriers of volatiles in protoplanetary disks and are crucial to our understanding of the chemistry that ultimately sets the organic composition of planets. The ERS program Ice Age on the JWST follows the ice evolution through all stages of star and planet formation. JWST/NIRSpec observations of the edge-on Class II protoplanetary disk HH~48~NE reveal spatially resolved absorptio… ▽ More Ices are the main carriers of volatiles in protoplanetary disks and are crucial to our understanding of the chemistry that ultimately sets the organic composition of planets. The ERS program Ice Age on the JWST follows the ice evolution through all stages of star and planet formation. JWST/NIRSpec observations of the edge-on Class II protoplanetary disk HH~48~NE reveal spatially resolved absorption features of the major ice components H$_2$O, CO$_2$, CO, and multiple weaker signatures from less abundant ices NH$_3$, OCN$^-$, and OCS. Isotopologue $^{13}$CO$_2$ ice has been detected for the first time in a protoplanetary disk. Since multiple complex light paths contribute to the observed flux, the ice absorption features are filled in by ice-free scattered light. The $^{12}$CO$_2$/$^{13}$CO$_2$ ratio of 14 implies that the $^{12}$CO$_2$ feature is saturated, without the flux approaching 0, indicative of a very high CO$_2$ column density on the line of sight, and a corresponding abundance with respect to hydrogen that is higher than ISM values by a factor of at least a few. Observations of rare isotopologues are crucial, as we show that the $^{13}$CO$_2$ observation allows us to determine the column density of CO$_2$ to be at an order of magnitude higher than the lower limit directly inferred from the observed optical depth. Radial variations in ice abundance, e.g., snowlines, are significantly modified since all observed photons have passed through the full radial extent of the disk. CO ice is observed at perplexing heights in the disk, extending to the top of the CO-emitting gas layer. We argue that the most likely interpretation is that we observe some CO ice at high temperatures, trapped in less volatile ices like H$_2$O and CO$_2$. Future radiative transfer models will be required to constrain the implications on our current understanding of disk physics and chemistry. △ Less

Submitted 14 September, 2023; originally announced September 2023.

Comments: 16 pages, 8 figures, accepted for publication in A&A

Water ice: temperature-dependent refractive indexes and their astrophysical implications

Authors: W. R. M. Rocha, M. G. Rachid, M. K. McClure, J. He, H. Linnartz

Abstract: Interstellar ices are largely composed of frozen water. It is important to derive fundamental parameters for H$_2$O ice such as absorption and scattering opacities for which accurate complex refractive indexes are needed. The primary goal of this work is to derive ice-grain opacities based on accurate H$_2$O ice complex refractive indexes and to assess their impact on the derivation of ice column… ▽ More Interstellar ices are largely composed of frozen water. It is important to derive fundamental parameters for H$_2$O ice such as absorption and scattering opacities for which accurate complex refractive indexes are needed. The primary goal of this work is to derive ice-grain opacities based on accurate H$_2$O ice complex refractive indexes and to assess their impact on the derivation of ice column densities and porosity in space. We use the \texttt{optool} code to derive ice-grain opacities values based on new mid-IR complex refractive index measurements of H$_2$O ice. Next, we use those opacities in the \texttt{RADMC-3D} code to run a radiative transfer simulation of a protostellar envelope containing H$_2$O ice. This is used to calculate water ice column densities. We find that the real refractive index in the mid-IR of H$_2$O ice at 30~K is $\sim$14\% lower than previously reported in the literature. This has a direct impact on the ice column densities derived from the simulations of embedded protostars. We find that ice porosity plays a significant role in the opacity of icy grains and that the H$_2$O libration mode can be used as a diagnostic tool to constrain the porosity level. Finally, the refractive indexes presented here allow us to estimate a grain size detection limit of 18~$μみゅー$m based on the 3~$μみゅー$m band whereas the 6~$μみゅー$m band allows tracing grain sizes larger than 20~$μみゅー$m. Based on radiative transfer simulations using new mid-IR refractive indexes, we conclude that H$_2$O ice leads to more absorption of infrared light than previously estimated. This implies that the 3 and 6~$μみゅー$m bands remain detectable in icy grains with sizes larger than 10~$μみゅー$m. Finally, we propose that also the H$_2$O ice libration band can be a diagnostic tool to constrain the porosity level of the interstellar ice, in addition to the OH dangling bond, which is routinely used for this purpose. △ Less

Submitted 23 August, 2023; originally announced August 2023.

Comments: Accepted for publication in A&A, 12 pages, 15 figures

Journal ref: A&A 681, A9 (2024)

Thermoelastic fracturing and buoyancy-driven convection: Surprising sources of longevity for EGS circulation

Authors: Mark McClure

Abstract: To maximize value from an EGS system, engineers need to optimize flow rate, well spacing, and well configuration. Numerical simulations yield interesting and surprising results regarding the effect of coupled processes on long term thermal performance. Thermoelastic fracture opening and propagation can have a significant negative effect on the uniformity of flow. On the other hand, interactions be… ▽ More To maximize value from an EGS system, engineers need to optimize flow rate, well spacing, and well configuration. Numerical simulations yield interesting and surprising results regarding the effect of coupled processes on long term thermal performance. Thermoelastic fracture opening and propagation can have a significant negative effect on the uniformity of flow. On the other hand, interactions between fracture opening and buoyancy driven fluid circulation cause downward fracture propagation during long-term circulation that greatly improves the thermal longevity of the system. Passive inflow control design can significantly mitigate the negative effect of thermoelastic fracture opening on flow uniformity, while maintaining the positive effects of thermoelastic fracture opening and propagation on flow rate and thermal longevity. Overall, simulations suggest that an EGS doublet with 8000 ft laterals at 475' F, using inflow control at the production well, could sustain electricity generation rates of 8 to 10 MWe for more than 30 years. Without inflow control, 6 to 8 MWe over 30 years is possible; however, there is greater risk of uncontrolled thermal breakthrough. △ Less

Submitted 4 August, 2023; originally announced August 2023.

Comments: Originally posted as a ResFrac Blog post on July 24 2023

The edge-on protoplanetary disk HH 48 NE II. Modeling ices and silicates

Authors: J. A. Sturm, M. K. McClure, J. B. Bergner, D. Harsono, E. Dartois, M. N. Drozdovskaya, S. Ioppolo, K. I. Öberg, C. J. Law, M. E. Palumbo, Y. J. Pendleton, W. R. M. Rocha, H. Terada, R. G. Urso

Abstract: The abundance and distribution of ice in protoplanetary disks (PPD) is critical to understand the linkage between the composition of circumstellar matter and the composition of exoplanets. Edge-on PPDs are a useful tool to constrain such ice composition and its location in the disk, as ice spectral signatures can be observed in absorption against the continuum emission arising from the warmer cent… ▽ More The abundance and distribution of ice in protoplanetary disks (PPD) is critical to understand the linkage between the composition of circumstellar matter and the composition of exoplanets. Edge-on PPDs are a useful tool to constrain such ice composition and its location in the disk, as ice spectral signatures can be observed in absorption against the continuum emission arising from the warmer central disk regions. The aim of this work is to model ice absorption features in PPDs and determine how well the abundance of the main ice species across the disk can be determined within the uncertainty of the physical parameter space. The edge-on PPD around HH 48 NE, a target of the JWST ERS program IceAge, is used as a reference system. We use RADMC-3D to raytrace the mid-infrared continuum. Using a constant parameterized ice abundance, ice opacities are added to the dust opacity in regions wherever the disk is cold enough for the main carbon, oxygen and nitrogen carriers to freeze out. The global abundance of the main ice carriers in HH 48 NE can be determined within a factor of 3, when taking the uncertainty of the physical parameters into account. Ice features in PPDs can be saturated at an optical depth <1, due to local saturation. Spatially observed ice optical depths cannot be directly related to column densities due to radiative transfer effects. Vertical snowlines will not be a clear transition due to the radially increasing height of the snowsurface, but their location may be constrained from observations using radiative transfer modeling. Radial snowlines are not really accesible. Not only the ice abundance, but also inclination, settling, grain size distribution and disk mass have strong impact on the observed ice absorption features in disks. Relative changes in ice abundance can be inferred from observations only if the source structure is well constrained △ Less

Submitted 8 May, 2023; v1 submitted 3 May, 2023; originally announced May 2023.

Comments: 18 pages, 12 figures, accepted for publication in Astronomy & Astrophysics

The edge-on protoplanetary disk HH 48 NE I. Modeling the geometry and stellar parameters

Authors: J. A. Sturm, M. K. McClure, C. J. Law, D. Harsono, J. B. Bergner, E. Dartois, M. N. Drozdovskaya, S. Ioppolo, K. I. Öberg, M. E. Palumbo, Y. J. Pendleton, W. R. M. Rocha, H. Terada, R. G. Urso

Abstract: Context. Observations of edge-on disks are an important tool for constraining general protoplanetary disk properties that cannot be determined in any other way. However, most radiative transfer models cannot simultaneously reproduce the spectral energy distributions (SEDs) and resolved scattered light and submillimeter observations of these systems, due to the differences in geometry and dust prop… ▽ More Context. Observations of edge-on disks are an important tool for constraining general protoplanetary disk properties that cannot be determined in any other way. However, most radiative transfer models cannot simultaneously reproduce the spectral energy distributions (SEDs) and resolved scattered light and submillimeter observations of these systems, due to the differences in geometry and dust properties at different wavelengths. Aims. We simultaneously constrain the geometry of the edge-on protoplanetary disk HH 48 NE and the characteristics of the host star. HH 48 NE is part of the JWST early release science program Ice Age. This work serves as a stepping stone towards a better understanding of the disk physical structure and icy chemistry in this particular source. This kind of modeling lays the groundwork for studying other edge-on sources to be observed with the JWST. Methods. We fit a parameterized dust model to HH 48 NE by coupling the radiative transfer code RADMC-3D and an MCMC framework. The dust structure was fitted independently to a compiled SED, a scattered light image at 0.8 $μみゅー$m and an ALMA dust continuum observation at 890 $μみゅー$m. Results. We find that 90% of the dust mass in HH 48 NE is settled to the disk midplane, less than in average disks, and that the atmospheric layers of the disk contain exclusively large grains (0.3-10 $μみゅー$m). The exclusion of small grains in the upper atmosphere likely has important consequences for the chemistry due to the deep penetration of high-energy photons. The addition of a relatively large cavity (ca. 50 au in radius) is necessary to explain the strong mid-infrared emission, and to fit the scattered light and continuum observations simultaneously. △ Less

Submitted 8 May, 2023; v1 submitted 3 May, 2023; originally announced May 2023.

Comments: 16 pages, 8 figures, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 677, A17 (2023)

Technical barriers for deep closed-loop geothermal

Authors: Mark McClure

Abstract: This is the most exciting time in my lifetime for geothermal. There are many, many innovative things happening. To name a few, promising new approaches to Enhanced Geothermal Systems, geothermal projects in sedimentary and lower enthalpy formations, new approaches for geothermal exploration, lithium extraction from produced brines, geothermal energy storage, integrations with CO2 storage and captu… ▽ More This is the most exciting time in my lifetime for geothermal. There are many, many innovative things happening. To name a few, promising new approaches to Enhanced Geothermal Systems, geothermal projects in sedimentary and lower enthalpy formations, new approaches for geothermal exploration, lithium extraction from produced brines, geothermal energy storage, integrations with CO2 storage and capture, and new technologies for producing energy from hot water that is coproduced with oil and gas. However, this post is about a concept about which I remain skeptical, deep closed-loop heat exchangers. These designs are sometimes called Advanced Geothermal Systems, AGS. To clarify, this blog post discusses the use of closed-loop heat exchangers in deep geothermal wells, 1000s of ft deep. This is a different application than the use of closed-loop heat exchangers in very shallow wells, 10s of ft deep, for ground source heat pumps. The fundamental challenge is that closed-loop heat exchanger designs rely on conduction, and sometimes free convection, to bring energy into the well. These processes are inherently much slower than forced convection, which is what drives energy transport into a conventional geothermal well. Slower energy transport means less revenue, and therefore, worse economics. In the past few years, there have been many creative new ideas for using closed-loop heat exchangers for geothermal energy production. Usually, these proposals recognize the challenge facing closed-loop, but introduce a wrinkle intended to overcome the problem. Despite these efforts, there are still daunting technical hurdles to widespread deployment of deep closed-loop geothermal. In this post, I go through some of the proposed designs and discuss key issues. △ Less

Submitted 22 March, 2023; originally announced March 2023.

Comments: This article was originally posted as a ResFrac Blog post

A new method for interpreting well-to-well interference tests and quantifying the magnitude of production impact: Theory and applications in a multi-basin case study

Authors: Mouin Almasoodi, Thad Andrews, Curtis Johnston, Ankush Singh, Mark McClure

Abstract: Interference tests are used in shale reservoirs to evaluate the strength of connectivity between wells. The results inform engineering decisions about well spacing. In this paper, we propose a new procedure for interpreting interference tests. We fit the initial interference response with the solution to the 1D diffusivity equation at an offset observation point. It is advantageous to use the init… ▽ More Interference tests are used in shale reservoirs to evaluate the strength of connectivity between wells. The results inform engineering decisions about well spacing. In this paper, we propose a new procedure for interpreting interference tests. We fit the initial interference response with the solution to the 1D diffusivity equation at an offset observation point. It is advantageous to use the initial interference response, rather than the subsequent trend, because it is less affected by nonlinearities, time-varying boundary conditions, and uncertainties about flow geometry and regime. From the curve fit, we estimate the hydraulic diffusivity and conductivity of the fractures connecting the wells. For engineering purposes, it would be useful to quantify the impact of interference on well production. Thus, we seek a relationship between the degree of production interference (DPI) and an appropriate dimensionless quantity that can be derived from the estimate of fracture conductivity. Using simulations run under a wide range of conditions, we find that the classical definition for dimensionless fracture conductivity does not achieve a consistent prediction of DPI. This occurs because the dimensionless fracture conductivity is derived assuming radial flow geometry, but the dominant flow geometry during shale production is linear. We also find that the CPG (Chow Pressure Group) metric does not yield consistently accurate predictions of DPI. As an alternative, we derive a dimensionless quantity similar to the classical dimensionless fracture conductivity, but derived for linear, flow geometry. Using this approach, we calculate a dimensionless interference length that collapses all cases onto a single curve that predicts DPI as a function of fracture conductivity, well spacing, and formation properties. We conclude by applying the new method to field cases from the Anadarko and Delaware Basins. △ Less

Submitted 3 February, 2023; originally announced February 2023.

An Ice Age JWST inventory of dense molecular cloud ices

Authors: M. K. McClure, W. R. M. Rocha, K. M. Pontoppidan, N. Crouzet, L. E. U. Chu, E. Dartois, T. Lamberts, J. A. Noble, Y. J. Pendleton, G. Perotti, D. Qasim, M. G. Rachid, Z. L. Smith, Fengwu Sun, Tracy L Beck, A. C. A. Boogert, W. A. Brown, P. Caselli, S. B. Charnley, Herma M. Cuppen, H. Dickinson, M. N. Drozdovskaya, E. Egami, J. Erkal, H. Fraser , et al. (17 additional authors not shown)

Abstract: Icy grain mantles are the main reservoir of the volatile elements that link chemical processes in dark, interstellar clouds with the formation of planets and composition of their atmospheres. The initial ice composition is set in the cold, dense parts of molecular clouds, prior to the onset of star formation. With the exquisite sensitivity of JWST, this critical stage of ice evolution is now acces… ▽ More Icy grain mantles are the main reservoir of the volatile elements that link chemical processes in dark, interstellar clouds with the formation of planets and composition of their atmospheres. The initial ice composition is set in the cold, dense parts of molecular clouds, prior to the onset of star formation. With the exquisite sensitivity of JWST, this critical stage of ice evolution is now accessible for detailed study. Here we show the first results of the Early Release Science program "Ice Age" that reveal the rich composition of these dense cloud ices. Weak ices, including, $^{13}$CO$_2$, OCN$^-$, $^{13}$CO, OCS, and COMs functional groups are now detected along two pre-stellar lines of sight. The $^{12}$CO$_2$ ice profile indicates modest growth of the icy grains. Column densities of the major and minor ice species indicate that ices contribute between 2 and 19% of the bulk budgets of the key C, O, N, and S elements. Our results suggest that the formation of simple and complex molecules could begin early in a water-ice rich environment. △ Less

Submitted 22 January, 2023; originally announced January 2023.

Comments: To appear in Nature Astronomy on January 23rd, 2023. 33 pages, 16 figures, 3 tables; includes extended and supplemental data sections. Part of the JWST Ice Age Early Release Science program's science enabling products. Enhanced spectra downloadable on Zenodo at the following DOI: 10.5281/zenodo.7501239

Disentangling protoplanetary disk gas mass and carbon depletion through combined atomic and molecular tracers

Authors: J. A. Sturm, A. S. Booth, M. K. McClure, M. Leemker, E. F. van Dishoeck

Abstract: The total disk gas mass and elemental C, N, O composition of protoplanetary disks are crucial ingredients for our understanding of planet formation. Measuring the gas mass is complicated, since H$_2$ cannot be detected in the cold bulk of the disk and the elemental abundances with respect to hydrogen are degenerate with gas mass in all disk models. We present new NOEMA observations of CO, $^{13}$C… ▽ More The total disk gas mass and elemental C, N, O composition of protoplanetary disks are crucial ingredients for our understanding of planet formation. Measuring the gas mass is complicated, since H$_2$ cannot be detected in the cold bulk of the disk and the elemental abundances with respect to hydrogen are degenerate with gas mass in all disk models. We present new NOEMA observations of CO, $^{13}$CO, C$^{18}$O and optically thin C$^{17}$O $J$=2-1 lines, and use additional high angular resolution Atacama Large Millimeter Array millimeter continuum and CO data to construct a representative model of LkCa 15. The transitions that constrain the gas mass and carbon abundance most are C$^{17}$O 2-1, N${_2}$H$^+$ 3-2 and HD 1-0. Using these three molecules we find that the gas mass in the LkCa 15 disk is $M_\mathrm{g}=0.01 ^{+0.01}_{-0.004} M_{\odot}$, a factor of six lower than estimated before. The carbon abundance is C/H = ($3 \pm 1.5) \times10^{-5}$, implying a moderate depletion of elemental carbon by a factor of 3-9. All other analyzed transitions also agree with these numbers, within a modeling uncertainty of a factor of two. Using the resolved \ce{C2H} image we find a C/O ratio of $\sim$1, which is consistent with literature values of H$_2$O depletion in this disk. The lack of severe carbon depletion in the LkCa 15 disk is consistent with the young age of the disk, but contrasts with the higher depletions seen in older cold transition disks. Combining optically thin CO isotopologue lines with N$_2$H$^+$ is promising to break the degeneracy between gas mass and CO abundance. The moderate level of depletion for this source with a cold, but young disk, suggests that long carbon transformation timescales contribute to the evolutionary trend seen in the level of carbon depletion among disk populations, rather than evolving temperature effects and presence of dust traps alone. △ Less

Submitted 19 September, 2022; originally announced September 2022.

Comments: 16 pages, 8 figures. Accepted for publication in A&A

LIDA - The Leiden Ice Database for Astrochemistry

Authors: W. R. M. Rocha, M. G. Rachid, B. Olsthoorn, E. F. van Dishoeck, M. K. McClure, H. Linnartz

Abstract: High quality vibrational spectra of solid-phase molecules in ice mixtures and for temperatures of astrophysical relevance are needed to interpret infrared observations toward protostars and background stars. Over the last 25 years, the Laboratory for Astrophysics at Leiden Observatory has provided more than 1100 spectra of diverse ice samples. Timely with the recent launch of the James Webb Space… ▽ More High quality vibrational spectra of solid-phase molecules in ice mixtures and for temperatures of astrophysical relevance are needed to interpret infrared observations toward protostars and background stars. Over the last 25 years, the Laboratory for Astrophysics at Leiden Observatory has provided more than 1100 spectra of diverse ice samples. Timely with the recent launch of the James Webb Space Telescope, we have fully upgraded the Leiden Ice Database for Astrochemistry (LIDA) adding recently measured spectra. The goal of this manuscript is to describe what options exist to get access to and work with a large collection of IR spectra, and the UV/vis to mid-infrared refractive index of H2O ice and astronomy-oriented online tools to support the interpretation of IR ice observations. LIDA uses Flask and Bokeh for generating the web pages and graph visualization, respectively, SQL for searching ice analogues within the database and Jmol for 3D molecule visualization. The infrared data in the database are recorded via transmission spectroscopy of ice films condensed on cryogenic substrates. The real UV/vis refractive indices of H2O ice are derived from interference fringes created from the simultaneous use of a monochromatic HeNe laser beam and a broadband Xe-arc lamp, whereas the real and imaginary mid-IR values are theoretically calculated. LIDA also offers online tools. The first tool, SPECFY, used to create a synthetic spectrum of ices towards protostars. The second tool aims at the calculation of mid-infrared refractive index values. LIDA allows to search, download and visualize experimental data of astrophysically relevant molecules in the solid phase, as well as to provide the means to support astronomical observations. As an example, we analyse the spectrum of the protostar AFGL 989 using the resources available in LIDA and derive the column densities of H2O, CO and CO2 ices. △ Less

Submitted 25 August, 2022; originally announced August 2022.

Comments: Accepted for publication in A&A. 29 pages, 9 figures. Online database at https://icedb.strw.leidenuniv.nl/

Journal ref: A&A 668, A63 (2022)

Ice Age : Chemo-dynamical modeling of Cha-MMS1 to predict new solid-phase species for detection with JWST

Authors: Mihwa Jin, Ka Ho Lam, Melissa K. McClure, Jeroen Terwisscha van Scheltinga, Zhi-Yun Li, Adwin Boogert, Eric Herbst, Shane W. Davis, Robin T. Garrod

Abstract: Chemical models and experiments indicate that interstellar dust grains and their ice mantles play an important role in the production of complex organic molecules (COMs). To date, the most complex solid-phase molecule detected with certainty in the ISM is methanol, but the James Webb Space Telescope (JWST) may be able to identify still larger organic species. In this study, we use a coupled chemo-… ▽ More Chemical models and experiments indicate that interstellar dust grains and their ice mantles play an important role in the production of complex organic molecules (COMs). To date, the most complex solid-phase molecule detected with certainty in the ISM is methanol, but the James Webb Space Telescope (JWST) may be able to identify still larger organic species. In this study, we use a coupled chemo-dynamical model to predict new candidate species for JWST detection toward the young star-forming core Cha-MMS1, combining the gas-grain chemical kinetic code MAGICKAL with a 1-D radiative hydrodynamics simulation using Athena++. With this model, the relative abundances of the main ice constituents with respect to water toward the core center match well with typical observational values, providing a firm basis to explore the ice chemistry. Six oxygen-bearing COMs (ethanol, dimethyl ether, acetaldehyde, methyl formate, methoxy methanol, and acetic acid), as well as formic acid, show abundances as high as, or exceeding, 0.01% with respect to water ice. Based on the modeled ice composition, the infrared spectrum is synthesized to diagnose the detectability of the new ice species. The contribution of COMs to IR absorption bands is minor compared to the main ice constituents, and the identification of COM ice toward the core center of Cha-MMS1 with the JWST NIRCAM/Wide Field Slitless Spectroscopy (2.4-5.0 micron) may be unlikely. However, MIRI observations (5-28 micron) toward COM-rich environments where solid-phase COM abundances exceed 1% with respect to the water ice column density might reveal the distinctive ice features of COMs. △ Less

Submitted 9 July, 2022; originally announced July 2022.

Comments: Accepted to ApJ (28 Pages, 11 figures)

Gemini-LIGHTS: Herbig Ae/Be and massive T-Tauri protoplanetary disks imaged with Gemini Planet Imager

Authors: Evan A. Rich, John D. Monnier, Alicia Aarnio, Anna S. E. Laws, Benjamin R. Setterholm, David J. Wilner, Nuria Calvet, Tim Harries, Chris Miller, Claire L. Davies, Fred C. Adams, Sean M. Andrews, Jaehan Bae, Catherine Espaillat, Alexandra Z. Greenbaum, Sasha Hinkley, Stefan Kraus, Lee Hartmann, Andrea Isella, Melissa McClure, Rebecca Oppenheimer, Laura M. Pérez, Zhaohuan Zhu

Abstract: We present the complete sample of protoplanetary disks from the Gemini- Large Imaging with GPI Herbig/T-tauri Survey (Gemini-LIGHTS) which observed bright Herbig Ae/Be stars and T-Tauri stars in near-infrared polarized light to search for signatures of disk evolution and ongoing planet formation. The 44 targets were chosen based on their near- and mid-infrared colors, with roughly equal numbers of… ▽ More We present the complete sample of protoplanetary disks from the Gemini- Large Imaging with GPI Herbig/T-tauri Survey (Gemini-LIGHTS) which observed bright Herbig Ae/Be stars and T-Tauri stars in near-infrared polarized light to search for signatures of disk evolution and ongoing planet formation. The 44 targets were chosen based on their near- and mid-infrared colors, with roughly equal numbers of transitional, pre-transitional, and full disks. Our approach explicitly did not favor well-known, "famous" disks or those observed by ALMA, resulting in a less-biased sample suitable to probe the major stages of disk evolution during planet formation. Our optimized data reduction allowed polarized flux as low as 0.002% of the stellar light to be detected, and we report polarized scattered light around 80% of our targets. We detected point-like companions for 47% of the targets, including 3 brown dwarfs (2 confirmed, 1 new), and a new super-Jupiter mass candidate around V1295 Aql. We searched for correlations between the polarized flux and system parameters, finding a few clear trends: presence of a companion drastically reduces the polarized flux levels, far-IR excess correlates with polarized flux for non-binary systems, and systems hosting disks with ring structures have stellar masses $<$ 3 Msun. Our sample also included four hot, dusty "FS CMa" systems and we detected large-scale ($>100$ au) scattered light around each, signs of extreme youth for these enigmatic systems. Science-ready images are publicly available through multiple distribution channels using a new FITS file standard jointly developed with members of the VLT/SPHERE team. △ Less

Submitted 12 June, 2022; originally announced June 2022.

Comments: 51 pages, 31 figures, 7 tables, accepted to AJ

The A to Z Guide to Accelerating Continuous Improvement with ResFrac

Authors: Mark McClure, Garrett Fowler, Chris Hewson, Charles Kang

Abstract: This document is designed to give a new ResFrac user the tools that they need to succeed. It describes what ResFrac does and how to use it. We cover much more than just button pressing. Using ResFrac successfully requires high-level engineering skills: how to design and execute a modeling workflow, how to think critically about problems, and how to communicate the results. This document is designed to give a new ResFrac user the tools that they need to succeed. It describes what ResFrac does and how to use it. We cover much more than just button pressing. Using ResFrac successfully requires high-level engineering skills: how to design and execute a modeling workflow, how to think critically about problems, and how to communicate the results. △ Less

Submitted 8 May, 2024; v1 submitted 29 May, 2022; originally announced May 2022.

Chemical Habitability: Supply and Retention of Life's Essential Elements During Planet Formation

Authors: Sebastiaan Krijt, Mihkel Kama, Melissa McClure, Johanna Teske, Edwin A. Bergin, Oliver Shorttle, Kevin J. Walsh, Sean N. Raymond

Abstract: Carbon, Hydrogen, Nitrogen, Oxygen, Phosphorus and Sulfur (CHNOPS) play key roles in the origin and proliferation of life on Earth. Given the universality of physics and chemistry, not least the ubiquity of water as a solvent and carbon as a backbone of complex molecules, CHNOPS are likely crucial to most habitable worlds. To help guide and inform the search for potentially habitable and ultimatel… ▽ More Carbon, Hydrogen, Nitrogen, Oxygen, Phosphorus and Sulfur (CHNOPS) play key roles in the origin and proliferation of life on Earth. Given the universality of physics and chemistry, not least the ubiquity of water as a solvent and carbon as a backbone of complex molecules, CHNOPS are likely crucial to most habitable worlds. To help guide and inform the search for potentially habitable and ultimately inhabited environments, we begin by summarizing the CHNOPS budget of various reservoirs on Earth, their role in shaping our biosphere, and their origins in the Solar Nebula. We then synthesize our current understanding of how these elements behave and are distributed in diverse astrophysical settings, tracing their journeys from synthesis in dying stars to molecular clouds, protoplanetary settings, and ultimately temperate rocky planets around main sequence stars. We end by identifying key branching points during this journey, highlighting instances where a forming planets' distribution of CHNOPS can be altered dramatically, and speculating about the consequences for the chemical habitability of these worlds. △ Less

Submitted 18 March, 2022; originally announced March 2022.

Comments: 39 pages, 7 figures; under review for publication as a chapter in Protostars and Planets VII (Editors: Shu-ichiro Inutsuka, Yuri Aikawa, Takayuki Muto, Kengo Tomida, and Motohide Tamura)

Commentary on Four New DFIT Papers: (a) In-Situ Measurements from EGS Collab; (b) Comparison of Techniques at the Bedretto Project; (c) Statistical Review of 62 DFITs; and (d) An Article Advocating the Tangent Method

Authors: Mark W. McClure

Abstract: This post provides commentary on recent four papers on diagnostic fracture injection testing (DFIT). The first paper uses in-situ deformation measurements to directly observe fractures opening and closing during fracture injection-falloff tests (Guglielmi et al., 2022). The second compares various stress measurement techniques in a series of fracture/injection tests from the Bedretto project (Brok… ▽ More This post provides commentary on recent four papers on diagnostic fracture injection testing (DFIT). The first paper uses in-situ deformation measurements to directly observe fractures opening and closing during fracture injection-falloff tests (Guglielmi et al., 2022). The second compares various stress measurement techniques in a series of fracture/injection tests from the Bedretto project (Broker and Ma, 2022). The third statistically reviews results from applying the interpretation procedure from McClure et al. (2019) to 62 DFITs across nine different shale plays (McClure et al., 2022). The fourth is an op-ed written in JPT (Journal of Petroleum Technology) by an advocate of the tangent method for estimating DFIT closure stress (Buijs, 2021; 2022). △ Less

Submitted 7 February, 2022; originally announced February 2022.

Tracing pebble drift and trapping using radial carbon depletion profiles in protoplanetary disks

Authors: J. A. Sturm, M. K. McClure, D. Harsono, S. Facchini, F. Long, M. Kama, E. A. Bergin, E. F. van Dishoeck

Abstract: The composition of planets may be largely determined by the chemical processing and accretion of icy pebbles in protoplanetary disks. Recent observations of protoplanetary disks hint at wide-spread depletion of gaseous carbon. The missing volatile carbon is likely frozen in CO and/or CO$_2$ ice on grains and locked into the disk through pebble trapping in pressure bumps or planetesimals. We presen… ▽ More The composition of planets may be largely determined by the chemical processing and accretion of icy pebbles in protoplanetary disks. Recent observations of protoplanetary disks hint at wide-spread depletion of gaseous carbon. The missing volatile carbon is likely frozen in CO and/or CO$_2$ ice on grains and locked into the disk through pebble trapping in pressure bumps or planetesimals. We present the results of the first successful ACA (Atacama Compact Array) [C I] $J$ = 1-0 mini-survey of seven protoplanetary disks. Using tailored azimuthally symmetric DALI (Dust And LInes) thermo-chemical disk models, supported by the [C I] $J$ = 1-0 and resolved CO isotopologue data, we determine the system-averaged elemental volatile carbon abundance in the outer disk of three sources. Six out of seven sources are detected in [C I] $J$ = 1-0 with ACA, four of which show a distinct disk component. Based on the modeling we find severe cold gaseous carbon depletion in the outer disk of DL Tau and moderate depletion in the outer disks of DR Tau and DO Tau. Combining the outer and inner disk carbon abundances, we demonstrate definitive evidence for radial drift in the disk of DL Tau, where the existence of multiple dust rings points to either short lived or leaky dust traps. We find dust locking in the compact and smooth disks of DO Tau and DR Tau, hinting at unresolved dust substructure. Comparing our results with stars of different ages and luminosities, we identify an observational evolutionary trend in gaseous carbon depletion that is consistent with dynamical models of CO depletion processes. Transport efficiency of solids in protoplanetary disks can significantly differ from what we expect based on the current resolved substructure in the continuum observations. This has important implications for our understanding of the impact of radial drift and pebble accretion on planetary compositions. △ Less

Submitted 11 January, 2022; originally announced January 2022.

Comments: 23 pages, 19 figures, accepted to A&A

Journal ref: A&A 660, A126 (2022)

Early planet formation in embedded protostellar disks: Setting the stage for the first generation of planetesimals

Authors: A. J. Cridland, G. P. Rosotti, B. Tabone, L. Tychoniec, M. McClure, E. F. van Dishoeck

Abstract: (Abridged) Recent surveys of young star formation regions have shown that the average Class II object does not have enough dust mass to make the cores of giant planets. Younger Class 0/I objects have enough dust in their embedded disk, which begs the questions: can the first steps of planet formation occur in these younger systems? The first step is building the first planetesimals, generally beli… ▽ More (Abridged) Recent surveys of young star formation regions have shown that the average Class II object does not have enough dust mass to make the cores of giant planets. Younger Class 0/I objects have enough dust in their embedded disk, which begs the questions: can the first steps of planet formation occur in these younger systems? The first step is building the first planetesimals, generally believed to be the product of the streaming instability. Hence the question can be restated: are the physical conditions of embedded disks conducive to the growth of the streaming instability? Here we model the collapse of a `dusty' proto-stellar cloud to show that if there is sufficient drift between the falling gas and dust, regions of the embedded disk can become sufficiently enhanced in dust to drive the streaming instability. We include four models, three with different dust grain sizes and one with a different initial cloud angular momentum to test a variety of collapse trajectories. We find a `sweet spot' for planetesimal formation for grain sizes of a few 10s of micron since they fall sufficiently fast relative to the gas to build a high dust-to-gas ratio along the disk midplane, but have slow enough radial drift speeds in the embedded disk to maintain the high dust-to-gas ratio. Unlike the gas, which is held in hydrostatic equilibrium for a time due to gas pressure, the dust can begin collapsing from all radii at a much earlier time. The streaming instability can produce at least between 7-35 M$_\oplus$ of planetesimals in the Class 0/I phase of our smooth embedded disks, depending on the size of the falling dust grains. This first generation of planetesimals could represent the first step in planet formation, and occurs earlier in the lifetime of the young star than is traditionally thought. △ Less

Submitted 13 December, 2021; originally announced December 2021.

Comments: Accepted for publication in A&A, 22 pages, 19 figures

Journal ref: A&A 662, A90 (2022)

The Scorpion Planet Survey: Wide-Orbit Giant Planets Around Young A-type Stars

Authors: Kevin Wagner, Daniel Apai, Markus Kasper, Melissa McClure, Massimo Robberto

Abstract: The first directly imaged exoplanets indicated that wide-orbit giant planets could be more common around A-type stars. However, the relatively small number of nearby A-stars has limited the precision of exoplanet demographics studies to $\gtrsim$10%. We aim to constrain the frequency of wide-orbit giant planets around A-stars using the VLT/SPHERE extreme adaptive optics system, which enables targe… ▽ More The first directly imaged exoplanets indicated that wide-orbit giant planets could be more common around A-type stars. However, the relatively small number of nearby A-stars has limited the precision of exoplanet demographics studies to $\gtrsim$10%. We aim to constrain the frequency of wide-orbit giant planets around A-stars using the VLT/SPHERE extreme adaptive optics system, which enables targeting $\gtrsim$100 A-stars between 100$-$200 pc. We present the results of a survey of 84 A-stars within the nearby $\sim$5$-$17 Myr-old Sco OB2 association. The survey detected three companions$-$one of which is a new discovery (HIP75056Ab), whereas the other two (HD 95086b and HIP65426b) are now-known planets that were included without a priori knowledge of their existence. We assessed the image sensitivity and observational biases with injection and recovery tests combined with Monte Carlo simulations to place constraints on the underlying demographics. We measure a decreasing frequency of giant planets with increasing separation, with measured values falling between 10$-$2% for separations of 30$-$100 au, and 95% confidence-level (CL) upper limits of $\lesssim$45$-$8% for planets on 30$-$100 au orbits, and $\lesssim$5% between 200$-$500 au. These values are in excellent agreement with recent surveys of A-stars in the solar neighborhood$-$supporting findings that giant planets at $\lesssim$100 au are more frequent around A-stars than around solar-type hosts. Finally, the relatively low occurrence rate of super-Jupiters on wide orbits, the positive correlation with stellar mass, and the inverse correlation with orbital separation are consistent with core accretion being their dominant formation mechanism. △ Less

Submitted 16 January, 2022; v1 submitted 3 December, 2021; originally announced December 2021.

Comments: Accepted for publication in AJ

Description of ResFrac automated history matching and optimization workflow

Authors: Charles A. Kang, Mark W. McClure, Somasekhar Reddy

Abstract: This document describes the functioning of the automated history matching and optimization workflow implemented with the ResFrac hydraulic fracturing and reservoir simulator. The purpose of the automated history matching and optimization workflow is enable automated calibration of simulator models to observed data and application to forward optimization. The automated workflow works by solving a f… ▽ More This document describes the functioning of the automated history matching and optimization workflow implemented with the ResFrac hydraulic fracturing and reservoir simulator. The purpose of the automated history matching and optimization workflow is enable automated calibration of simulator models to observed data and application to forward optimization. The automated workflow works by solving a formal mathematical optimization problem to minimize misfit with observations from any point in the lifecycle of a hydraulically fractured reservoir, or to maximize a quantity of interest associated such a reservoir, such as net present value. Objective function evaluations in the optimization problem consist of runs of the ResFrac simulator, which is a physics-based model of hydraulic fracturing and reservoir phenomena. The workflow employs a proxy model to improve computational speed and applies experimental design and Bayesian sampling techniques to generate points with which to train the proxy model. This document also provides an overview of the software infrastructure developed to support the automated workflow. △ Less

Submitted 25 November, 2021; originally announced November 2021.

Comments: 10 pages, one figure

Transverse Hydraulic Fracture Initiation: Insights from Intermediate-Scale Hydraulic Fracture Tests

Authors: Jeffrey Burghardt, Hunter Knox, Chris Strickland, Vincent Vermeul, Pengcheng Fu, Craig Ulrich, Mark McClure, Mathew Ingraham, Paul Schwering, Timothy Kneafsey, Doug Blankenship

Abstract: In a uniform stress field, a tensile hydraulic fracture will favor propagation in the plane formed by the maximum (most compressive) and intermediate principal stresses. For the propagation of hydraulic fractures in non-uniform stress fields, especially those where the principal stress directions change spatially, things become much more complex, are less well understood, and therefore are more di… ▽ More In a uniform stress field, a tensile hydraulic fracture will favor propagation in the plane formed by the maximum (most compressive) and intermediate principal stresses. For the propagation of hydraulic fractures in non-uniform stress fields, especially those where the principal stress directions change spatially, things become much more complex, are less well understood, and therefore are more difficult to predict. Such is often the case when a hydraulic fracture is initiated from a borehole drilled in the direction of the minimum principal stress. This is a well-known problem that has been the subject of several numerical and experimental studies. The purpose of this paper is to offer new insights from a series of hydraulic fracture tests conducted in sub-horizontal boreholes drilled in the Sanford Underground Research Facility. This project, called the Enhanced Geothermal Systems (EGS) Collab project, has developed a test bed at a depth of approximately 1478 meters underground consisting of sub-horizontal injection, production, and monitoring wells. This paper presents stimulation test data, and supporting geologic/geophysical characterization data. The test data are compared to previous numerical and experimental studies to determine how well this new experimental data compares with previous modeling, laboratory, and field tests. △ Less

Submitted 1 October, 2021; originally announced October 2021.

Report number: PNNL-SA-161540

Molecules with ALMA at Planet-forming Scales. XX. The Massive Disk Around GM Aurigae

Authors: Kamber R. Schwarz, Jenny K. Calahan, Ke Zhang, Felipe Alarcón, Yuri Aikawa, Sean M. Andrews, Jaehan Bae, Edwin A. Bergin, Alice S. Booth, Arthur D. Bosman, Gianni Cataldi, L. Ilsedore Cleeves, Ian Czekala, Jane Huang, John D. Ilee, Charles J. Law, Romane Le Gal, Yao Liu, Feng Long, Ryan A. Loomis, Enrique Macías, Melissa McClure, François Ménard, Karin I. Öberg, Richard Teague , et al. (3 additional authors not shown)

Abstract: Gas mass remains one of the most difficult protoplanetary disk properties to constrain. With much of the protoplanetary disk too cold for the main gas constituent, H2, to emit, alternative tracers such as dust, CO, or the H2 isotopolog HD are used. However, relying on disk mass measurements from any single tracer requires assumptions about the tracer's abundance relative to \hh\ and the disk tempe… ▽ More Gas mass remains one of the most difficult protoplanetary disk properties to constrain. With much of the protoplanetary disk too cold for the main gas constituent, H2, to emit, alternative tracers such as dust, CO, or the H2 isotopolog HD are used. However, relying on disk mass measurements from any single tracer requires assumptions about the tracer's abundance relative to \hh\ and the disk temperature structure. Using new Atacama Large Millimeter/submillimeter Array (ALMA) observations from the Molecules with ALMA at Planet-forming Scales (MAPS) ALMA Large Program as well as archival ALMA observations, we construct a disk physical/chemical model of the protoplanetary disk GM Aur. Our model is in good agreement with the spatially resolved CO isotopolog emission from eleven rotational transitions with spatial resolution ranging from 0.15'' to 0.46'' (24-73 au at 159 pc) and the spatially unresolved HD J=1-0 detection from Herschel. Our best-fit model favors a cold protoplanetary disk with a total gas mass of approximately 0.2 solar masses, a factor of 10 reduction in CO gas inside roughly 100 au and a factor of 100 reduction outside of 100 au. Despite its large mass, the disk appears to be on the whole gravitationally stable based on the derived Toomre Q parameter. However, the region between 70 and 100 au, corresponding to one of the millimeter dust rings, is close to being unstable based on the calculated Toomre Q of <1.7. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement. △ Less

Submitted 13 September, 2021; originally announced September 2021.

Comments: 18 pages, 12 figures

The formation of planetary systems with SPICA

Authors: I. Kamp, M. Honda, H. Nomura, M. Audard, D. Fedele, L. B. F. M. Waters, Y. Aikawa, A. Banzatti, J. E. Bowey, M. Bradford, C. Dominik, K. Furuya, E. Habart, D. Ishihara, D. Johnstone, G. Kennedy, M. Kim, Q. Kral, S. P. Lai, B. Larsson, M. McClure, A. Miotello, M. Momose, T. Nakagawa, D. Naylor , et al. (16 additional authors not shown)

Abstract: In this era of spatially resolved observations of planet forming disks with ALMA and large ground-based telescopes such as the VLT, Keck and Subaru, we still lack statistically relevant information on the quantity and composition of the material that is building the planets, such as the total disk gas mass, the ice content of dust, and the state of water in planetesimals. SPICA is an infrared spac… ▽ More In this era of spatially resolved observations of planet forming disks with ALMA and large ground-based telescopes such as the VLT, Keck and Subaru, we still lack statistically relevant information on the quantity and composition of the material that is building the planets, such as the total disk gas mass, the ice content of dust, and the state of water in planetesimals. SPICA is an infrared space mission concept developed jointly by JAXA and ESA to address these questions. The key unique capabilities of SPICA that enable this research are (1) the wide spectral coverage 10-220 micron, (2) the high line detection sensitivity of (1-2) 10-19 W m-2 with R~2000-5000 in the far-IR (SAFARI) and 10-20 W m-2 with R~29000 in the mid-IR (SMI, spectrally resolving line profiles), (3) the high far-IR continuum sensitivity of 0.45 mJy (SAFARI), and (4) the observing efficiency for point source surveys. This paper details how mid- to far-IR infrared spectra will be unique in measuring the gas masses and water/ice content of disks and how these quantities evolve during the planet forming period. These observations will clarify the crucial transition when disks exhaust their primordial gas and further planet formation requires secondary gas produced from planetesimals. The high spectral resolution mid-IR is also unique for determining the location of the snowline dividing the rocky and icy mass reservoirs within the disk and how the divide evolves during the build-up of planetary systems. Infrared spectroscopy (mid- to far-IR) of key solid state bands is crucial for assessing whether extensive radial mixing, which is part of our Solar System history, is a general process occurring in most planetary systems and whether extrasolar planetesimals are similar to our Solar System comets/asteroids. ... (abbreviated) △ Less

Submitted 25 June, 2021; originally announced June 2021.

Comments: accepted for publication in PASA

Infrared spectra of complex organic molecules in astronomically relevant ice matrices. III. Methyl formate and its tentative solid-state detection

Authors: Jeroen Terwisscha van Scheltinga, Giulia Marcandalli, Melissa K. McClure, Michiel R. Hogerheijde, Harold Linnartz

Abstract: Context. Infrared spectroscopy of star and planet forming regions is at the dawn of a new age with the upcoming James Webb Space Telescope. In support of these observations, laboratory spectra are required to identify complex organic molecules in the ices that cover the dust grains in these regions. Aims. This study aims to provide reference spectra to firmly detect icy methyl formate in the dif… ▽ More Context. Infrared spectroscopy of star and planet forming regions is at the dawn of a new age with the upcoming James Webb Space Telescope. In support of these observations, laboratory spectra are required to identify complex organic molecules in the ices that cover the dust grains in these regions. Aims. This study aims to provide reference spectra to firmly detect icy methyl formate in the different stages of star and planet forming regions. Methyl formate is mixed in astronomically relevant matrices, and the peak positions, FWHMs, and relative band intensities are characterized for different temperatures to provide an analytical tool for astronomers. Methods. Methyl formate is deposited at 15 K under high-vacuum conditions. Specifically, methyl formate is deposited pure and mixed with CO, H$_2$CO, CH$_3$OH, H$_2$O, and CO:H$_2$CO:CH$_3$OH combined. Throughout the experiment infrared spectra are acquired with a FTIR spectrometer in the range from 4000-500 cm$^{-1}$ (2.5-20 $μみゅー$m) at a spectral resolution of 0.5 cm$^{-1}$. Results. We present the characterization of five solid-state methyl formate vibrational modes in pure and astronomically relevant ice matrices. The five selected vibrational modes, namely the C=O stretch, C$-$O stretch, CH$_3$ rocking, O$-$CH$_3$ stretching, and OCO deformation, are best suited for a JWST identification of methyl formate. For each of these vibrational modes, and each of the mixtures the TvS heatmaps, peak position versus FWHM, and relative band intensities are given. Additionally, the acquired reference spectra of methyl formate are compared with Spitzer observations of HH 46. A tentative detection of methyl formate provides an upper limit to the column density of $1.7\times10^{17}$ cm$^{-2}$, corresponding to an upper limit relative to water of $\leq 2.2\%$ and $\leq 40\%$ with respect to methanol. △ Less

Submitted 5 May, 2021; originally announced May 2021.

Comments: Accepted for publication in A&A

Journal ref: A&A 651, A95 (2021)

Disk Evolution Study Through Imaging of Nearby Young Stars (DESTINYS): Late infall causing disk misalignment and dynamic structures in SU Aur

Authors: C. Ginski, S. Facchini, J. Huang, M. Benisty, D. Vaendel, L. Stapper, C. Dominik, J. Bae, F. Menard, G. Muro-Arena, M. Hogerheijde, M. McClure, R. G. van Holstein, T. Birnstiel, Y. Boehler, A. Bohn, M. Flock, E. E. Mamajek, C. F. Manara, P. Pinilla, C. Pinte, A. Ribas

Abstract: Gas-rich circumstellar disks are the cradles of planet formation. As such, their evolution will strongly influence the resulting planet population. In the ESO DESTINYS large program, we study these disks within the first 10 Myr of their development with near-infrared scattered light imaging. Here we present VLT/SPHERE polarimetric observations of the nearby class II system SU Aur in which we resol… ▽ More Gas-rich circumstellar disks are the cradles of planet formation. As such, their evolution will strongly influence the resulting planet population. In the ESO DESTINYS large program, we study these disks within the first 10 Myr of their development with near-infrared scattered light imaging. Here we present VLT/SPHERE polarimetric observations of the nearby class II system SU Aur in which we resolve the disk down to scales of ~7 au. In addition to the new SPHERE observations, we utilize VLT/NACO, HST/STIS and ALMA archival data. The new SPHERE data show the disk around SU Aur and extended dust structures in unprecedented detail. We resolve several dust tails connected to the Keplerian disk. By comparison with ALMA data, we show that these dust tails represent material falling onto the disk. The disk itself shows an intricate spiral structure and a shadow lane, cast by an inner, misaligned disk component. Our observations suggest that SU Aur is undergoing late infall of material, which can explain the observed disk structures. SU Aur is the clearest observational example of this mechanism at work and demonstrates that late accretion events can still occur in the class II phase, thereby significantly affecting the evolution of circumstellar disks. Constraining the frequency of such events with additional observations will help determine whether this process is responsible for the spin-orbit misalignment in evolved exoplanet systems. △ Less

Submitted 17 February, 2021; originally announced February 2021.

Comments: 18 pages, 12 figures, published in ApJL on 18-02-2021

Journal ref: ApJL, 908, L25 (2021)

Direct Imaging Discovery of a Young Brown Dwarf Companion to an A2V Star

Authors: Kevin Wagner, Dániel Apai, Markus Kasper, Melissa McClure, Massimo Robberto, Thayne Currie

Abstract: We present the discovery and spectroscopy of HIP 75056Ab, a companion directly imaged at a very small separation of 0.125 arcsec to an A2V star in the Scorpius-Centaurus OB2 association. Our observations utilized VLT/SPHERE between 2015$-$2019, enabling low-resolution spectroscopy (0.95$-$1.65 $μみゅーm$), dual-band imaging (2.1$-$2.25 $μみゅーm$), and relative astrometry over a four-year baseline. HIP 75056A… ▽ More We present the discovery and spectroscopy of HIP 75056Ab, a companion directly imaged at a very small separation of 0.125 arcsec to an A2V star in the Scorpius-Centaurus OB2 association. Our observations utilized VLT/SPHERE between 2015$-$2019, enabling low-resolution spectroscopy (0.95$-$1.65 $μみゅーm$), dual-band imaging (2.1$-$2.25 $μみゅーm$), and relative astrometry over a four-year baseline. HIP 75056Ab is consistent with spectral types in the range of M6$-$L2 and $T_{\rm eff}\sim$ 2000$-$2600 K. A comparison of the companion's brightness to evolutionary tracks suggests a mass of $\sim$20$-$30 M$_{Jup}$. The astrometric measurements are consistent with an orbital semi-major axis of $\sim$15$-$45 au and an inclination close to face-on (i$\lesssim$35$^o$). In this range of mass and orbital separation, HIP 75056Ab is likely at the low-mass end of the distribution of companions formed via disk instability, although a formation of the companion via core accretion cannot be excluded. The orbital constraints are consistent with the modest eccentricity values predicted by disk instability, a scenario that can be confirmed by further astrometric monitoring. HIP 75056Ab may be utilized as a low-mass atmospheric comparison to older, higher-mass brown dwarfs, and also to young giant planets. Finally, the detection of HIP 75056Ab at 0.125 arcsec represents a milestone in detecting low-mass companions at separations corresponding to the habitable zones of nearby Sun-like stars. △ Less

Submitted 28 September, 2020; v1 submitted 17 September, 2020; originally announced September 2020.

Comments: Accepted for publication in ApJL

A multi-frequency ALMA characterization of substructures in the GM Aur protoplanetary disk

Authors: Jane Huang, Sean M. Andrews, Cornelis P. Dullemond, Karin I. Oberg, Chunhua Qi, Zhaohuan Zhu, Tilman Birnstiel, John M. Carpenter, Andrea Isella, Enrique Macias, Melissa K. McClure, Laura M. Perez, Richard Teague, David J. Wilner, Shangjia Zhang

Abstract: The protoplanetary disk around the T Tauri star GM Aur was one of the first hypothesized to be in the midst of being cleared out by a forming planet. As a result, GM Aur has had an outsized influence on our understanding of disk structure and evolution. We present 1.1 and 2.1 mm ALMA continuum observations of the GM Aur disk at a resolution of ~50 mas (~8 au), as well as HCO$^+$ $J=3-2$ observatio… ▽ More The protoplanetary disk around the T Tauri star GM Aur was one of the first hypothesized to be in the midst of being cleared out by a forming planet. As a result, GM Aur has had an outsized influence on our understanding of disk structure and evolution. We present 1.1 and 2.1 mm ALMA continuum observations of the GM Aur disk at a resolution of ~50 mas (~8 au), as well as HCO$^+$ $J=3-2$ observations at a resolution of ~100 mas. The dust continuum shows at least three rings atop faint, extended emission. Unresolved emission is detected at the center of the disk cavity at both wavelengths, likely due to a combination of dust and free-free emission. Compared to the 1.1 mm image, the 2.1 mm image shows a more pronounced "shoulder" near R~40 au, highlighting the utility of longer-wavelength observations for characterizing disk substructures. The spectral index $αあるふぁ$ features strong radial variations, with minima near the emission peaks and maxima near the gaps. While low spectral indices have often been ascribed to grain growth and dust trapping, the optical depth of GM Aur's inner two emission rings renders their dust properties ambiguous. The gaps and outer disk ($R>100$ au) are optically thin at both wavelengths. Meanwhile, the HCO$^+$ emission indicates that the gas cavity is more compact than the dust cavity traced by the millimeter continuum, similar to other disks traditionally classified as "transitional." △ Less

Submitted 29 January, 2020; originally announced January 2020.

Comments: 31 pages, 15 figures, accepted by ApJ. Self-calibration and imaging scripts, images, and self-calibrated visibilities available at https://zenodo.org/record/3628656#.XjEXyRNKiYV

Measuring the atomic composition of planetary building blocks

Authors: M. K. McClure, C. Dominik, M. Kama

Abstract: Volatile molecules are critical to terrestrial planetary habitability, yet difficult to observe directly where planets form at the midplanes of protoplanetary disks. It is unclear whether the inner 1 AUえーゆー of disks are volatile-poor or if this region is resupplied with ice-rich dust from colder disk regions. Dust traps at radial pressure maxima bounding disk gaps can cut off the inner disk from such… ▽ More Volatile molecules are critical to terrestrial planetary habitability, yet difficult to observe directly where planets form at the midplanes of protoplanetary disks. It is unclear whether the inner 1 AUえーゆー of disks are volatile-poor or if this region is resupplied with ice-rich dust from colder disk regions. Dust traps at radial pressure maxima bounding disk gaps can cut off the inner disk from such volatile reservoirs. However, the trap retention efficiency and atomic composition of trapped dust have not been measured. We present a new technique to measure the absolute atomic abundances in gas accreting onto T Tauri stars and infer the bulk atomic composition and distribution of midplane solids retained in the disk around the young star TW Hya. We identify line emission from gas-phase material inside the dust sublimation rim of TW Hya. Gaussian decomposition of the strongest H Paschen lines isolates the inner disk hydrogen emission. We measure several key elemental abundances, relative to hydrogen, using a chemical photoionization model and infer dust retention in the disk. With a 1D transport model, we determine approximate radial locations and retention efficiencies of dust traps for different elements. Volatile and refractory elements are depleted from TW Hya's hot gas by factors of ~10^2 and up to 10^5, respectively. Dust traps beyond the CO and N2 snowline cumulatively sequester 96% of the total dust, while the trap at 2 AUえーゆー retains 3% of the initial dust mass. The high depletions of Si, Mg, and Ca are explained by a third trap at 0.3 AUえーゆー. TW Hya has a significant volatile reservoir rich in C- and N-ices in its outer ring structure. However, the lack of C resupply may leave the terrestrial planet-forming region dry and carbon-poor. Any planets that form within the silicate dust trap at 0.3 AUえーゆー may resemble Earth in terms of the degree of their volatile depletion. △ Less

Submitted 15 September, 2020; v1 submitted 16 October, 2019; originally announced October 2019.

Comments: Accepted for publication in A&A Letters. 6 pages, 3 figures, plus appendices

Journal ref: A&A 642, L15 (2020)

Carbon depletion observed inside T Tauri inner rims: Formation of icy, kilometer size planetesimals by 1 Myr

Authors: M. K. McClure

Abstract: The carbon content of protoplanetary disks is an important parameter to characterize planets formed at different disk radii. There is some evidence from far-infrared and submillimeter observations that gas in the outer disk is depleted in carbon, with a corresponding enhancement of carbon-rich ices at the disk midplane. Observations of the carbon content inside of the inner sublimation rim could c… ▽ More The carbon content of protoplanetary disks is an important parameter to characterize planets formed at different disk radii. There is some evidence from far-infrared and submillimeter observations that gas in the outer disk is depleted in carbon, with a corresponding enhancement of carbon-rich ices at the disk midplane. Observations of the carbon content inside of the inner sublimation rim could confirm how much carbon remains locked in kilometer size bodies in the disk. Using near-infrared spectra, I self-consistently determine the stellar, extinction, veiling, and accretion properties of the 26 stars in my sample. From these values, I extract the inner disk excess of the target stars and identify a series of C^0 recombination lines in 18 of these disks and use the CHIANTI atomic line database with an optically thin slab model to constrain the average n_e, T_e, and n_C for these lines in the five disks with a complete set of lines. By comparing these values with other slab models of the inner disk using the Cloudy photoionization code, I also constrain n_H and the carbon abundance, X_C, and hence the amount of carbon `missing' from the slab. The inner disks are extremely dense (n_H~10^{16} cm^{-3}), warm (T_e~4500 K), and moderately ionized (logX_e~3.3). Three of the five modeled disks show robust carbon depletion <=42 relative to the solar value. I discuss multiple ways in which the `missing' carbon could be locked out of the accreting gas. Given the high-density inner disk gas, evidence for radial drift, and lack of obvious gaps in these three systems, their carbon depletion is most consistent with the `missing' carbon being sequestered in kilometer size bodies. For DR Tau, nitrogen and silicon are also depleted by factors of 45 and 4, respectively, suggesting that the kilometer size bodies into which the grains are locked were formed beyond the N_2 snowline. △ Less

Submitted 14 October, 2019; originally announced October 2019.

Comments: 14 pages, 7 figures, appendices, accepted to A&A July 16th, 2019

A dust and gas cavity in the disc around CQ Tau revealed by ALMA

Authors: M. Giulia Ubeira Gabellini, Anna Miotello, Stefano Facchini, Enrico Ragusa, Giuseppe Lodato, Leonardo Testi, Myriam Benisty, Simon Bruderer, Nicolàs T. Kurtovic, Sean Andrews, John Carpenter, Stuartt A. Corder, Giovanni Dipierro, Barbara Ercolano, Davide Fedele, Greta Guidi, Thomas Henning, Andrea Isella, Woojin Kwon, Hendrik Linz, Melissa McClure, Laura Perez, Luca Ricci, Giovanni Rosotti, Marco Tazzari , et al. (1 additional authors not shown)

Abstract: The combination of high resolution and sensitivity offered by ALMA is revolutionizing our understanding of protoplanetary discs, as their bulk gas and dust distributions can be studied independently. In this paper we present resolved ALMA observations of the continuum emission ($λらむだ=1.3$ mm) and CO isotopologues ($^{12}$CO, $^{13}$CO, C$^{18}$O $J=2-1$) integrated intensity from the disc around the… ▽ More The combination of high resolution and sensitivity offered by ALMA is revolutionizing our understanding of protoplanetary discs, as their bulk gas and dust distributions can be studied independently. In this paper we present resolved ALMA observations of the continuum emission ($λらむだ=1.3$ mm) and CO isotopologues ($^{12}$CO, $^{13}$CO, C$^{18}$O $J=2-1$) integrated intensity from the disc around the nearby ($d = 162$ pc), intermediate mass ($M_{\star}=1.67\,M_{\odot}$) pre-main-sequence star CQ Tau. The data show an inner depression in continuum, and in both $^{13}$CO and C$^{18}$O emission. We employ a thermo-chemical model of the disc reproducing both continuum and gas radial intensity profiles, together with the disc SED. The models show that a gas inner cavity with size between 15 and 25 au is needed to reproduce the data with a density depletion factor between $\sim 10^{-1}$ and $\sim 10^{-3}$. The radial profile of the distinct cavity in the dust continuum is described by a Gaussian ring centered at $R_{\rm dust}=53\,$au and with a width of $σしぐま=13\,$au. Three dimensional gas and dust numerical simulations of a disc with an embedded planet at a separation from the central star of $\sim20\,$au and with a mass of $\sim 6\textrm{-} 9\,M_{\rm Jup}$ reproduce qualitatively the gas and dust profiles of the CQ Tau disc. However, a one planet model appears not to be able to reproduce the dust Gaussian density profile predicted using the thermo-chemical modeling. △ Less

Submitted 2 May, 2019; originally announced May 2019.

Discussion of the Paper SPE-187038-MS: Fracture Closure Stress: Reexamining Field and Laboratory Experiments of Fracture Closure Using Modern Interpretation Methodologies

Authors: Mark W. McClure

Abstract: In recent years, there has been discussion in the literature regarding methods of estimating the magnitude of the minimum principal stress from subsurface fracture injection tests, commonly called Diagnostic Fracture Injection Tests (DFITs). McClure et al. (2016) used modeling and mathematical analysis to propose that changes should be made to common interpretation techniques. Subsequently, Craig… ▽ More In recent years, there has been discussion in the literature regarding methods of estimating the magnitude of the minimum principal stress from subsurface fracture injection tests, commonly called Diagnostic Fracture Injection Tests (DFITs). McClure et al. (2016) used modeling and mathematical analysis to propose that changes should be made to common interpretation techniques. Subsequently, Craig et al. (2017) attempted to refute those findings using field and laboratory observations. This discussion paper reviews the interpretations from Craig et al. (2017) and concludes that their conclusions are unsupported by the data presented. △ Less

Submitted 15 April, 2019; originally announced April 2019.

The Disk Gas Mass and the Far-IR Revolution

Authors: Edwin A. Bergin, Klaus M. Pontoppidan, Charles M. Bradford, L. Ilsedore Cleeves, Neal J. Evans, Maryvonne Gerin, Paul F. Goldsmith, Quentin Kral, Gary J. Melnick, Melissa McClure, Karin Oberg, Thomas L. Roellig, Edward Wright, Richard Teague, Jonathan P. Williams, Ke Zhang

Abstract: The gaseous mass of protoplanetary disks is a fundamental quantity in planet formation. The presence of gas is necessary to assemble planetesimals, it determines timescales of giant planet birth, and it is an unknown factor for a wide range of properties of planet formation, from chemical abundances (X/H) to the mass efficiency of planet formation. The gas mass obtained from traditional tracers, s… ▽ More The gaseous mass of protoplanetary disks is a fundamental quantity in planet formation. The presence of gas is necessary to assemble planetesimals, it determines timescales of giant planet birth, and it is an unknown factor for a wide range of properties of planet formation, from chemical abundances (X/H) to the mass efficiency of planet formation. The gas mass obtained from traditional tracers, such as dust thermal continuum and CO isotopologues, are now known to have significant (1 - 2 orders of magnitude) discrepancies. Emission from the isotopologue of H2, hydrogen deuteride (HD), offers an alternative measurement of the disk gas mass. Of all of the regions of the spectrum, the far-infrared stands out in that orders of magnitude gains in sensitivity can be gleaned by cooling a large aperture telescope to 8 K. Such a facility can open up a vast new area of the spectrum to exploration. One of the primary benefits of this far-infrared revolution would be the ability to survey hundreds of planet-forming disks in HD emission to derive their gaseous masses. For the first time, we will have statistics on the gas mass as a function of evolution, tracing birth to dispersal as a function of stellar spectral type. These measurements have broad implications for our understanding of the time scale during which gas is available to form giant planets, the dynamical evolution of the seeds of terrestrial worlds, and the resulting chemical composition of pre-planetary embryos carrying the elements needed for life. Measurements of the ground-state line of HD requires a space-based observatory operating in the far-infrared at 112 microns. △ Less

Submitted 20 March, 2019; originally announced March 2019.

Comments: Science white paper submitted to the Astro2020 Decadal Survey

The trail of water and the delivery of volatiles to habitable planets

Authors: Klaus M. Pontoppidan, Andrea Banzatti, Edwin Bergin, Geoffrey A. Blake, Sean Brittain, Maryvonne Gerin, Paul Goldsmith, Quentin Kral, David Leisawitz, Dariusz Lis, Melissa McClure, Stefanie Milam, Gary Melnick, Joan Najita, Karin Öberg, Matt Richter, Colette Salyk, Martina Wiedner, Ke Zhang

Abstract: Water is fundamental to our understanding of the evolution of planetary systems and the delivery of volatiles to the surfaces of potentially habitable planets. Yet, we currently have essentially no facilities capable of observing this key species comprehensively. With this white paper, we argue that we need a relatively large, cold space-based observatory equipped with a high-resolution spectromet… ▽ More Water is fundamental to our understanding of the evolution of planetary systems and the delivery of volatiles to the surfaces of potentially habitable planets. Yet, we currently have essentially no facilities capable of observing this key species comprehensively. With this white paper, we argue that we need a relatively large, cold space-based observatory equipped with a high-resolution spectrometer, in the mid- through far-infrared wavelength range (25-600~$μみゅー$m) in order to answer basic questions about planet formation, such as where the Earth got its water, how giant planets and planetesimals grow, and whether water is generally available to planets forming in the habitable zone of their host stars. △ Less

Submitted 15 March, 2019; originally announced March 2019.

Comments: Science white paper submitted to the Astro2020 Decadal Survey

Multiple spiral arms in the disk around intermediate-mass binary HD 34700A

Authors: John D. Monnier, Tim J. Harries, Jaehan Bae, Benjamin R. Setterholm, Anna Laws, Alicia Aarnio, Fred C. Adams, Sean Andrews, Nuria Calvet, Catherine Espaillat, Lee Hartmann, Stefan Kraus, Melissa McClure, Chris Miller, Rebecca Oppenheimer, David Wilner, Zhaohuan Zhu

Abstract: We present the first images of the transition disk around the close binary system HD 34700A in polarized scattered light using the Gemini Planet Imager instrument on Gemini South. The J and H band images reveal multiple spiral-arm structures outside a large (R = 0.49" = 175 au) cavity along with a bluish spiral structure inside the cavity. The cavity wall shows a strong discontinuity and we clearl… ▽ More We present the first images of the transition disk around the close binary system HD 34700A in polarized scattered light using the Gemini Planet Imager instrument on Gemini South. The J and H band images reveal multiple spiral-arm structures outside a large (R = 0.49" = 175 au) cavity along with a bluish spiral structure inside the cavity. The cavity wall shows a strong discontinuity and we clearly see significant non-azimuthal polarization Uphi consistent with multiple scattering within a disk at an inferred inclination ~42deg. Radiative transfer modeling along with a new Gaia distance suggest HD 37400A is a young (~5 Myr) system consisting of two intermediate-mass (~2Msun) stars surrounded by a transitional disk and not a solar-mass binary with a debris disk as previously classified. Conventional assumptions of the dust-to-gas ratio would rule out a gravitational instability origin to the spirals while hydrodynamical models using the known external companion or a hypothetical massive protoplanet in the cavity both have trouble reproducing the relatively large spiral arm pitch angles (~30deg) without fine tuning of gas temperature. We explore the possibility that material surrounding a massive protoplanet could explain the rim discontinuity after also considering effects of shadowing by an inner disk. Analysis of archival Hubble Space Telescope data suggests the disk is rotating counterclockwise as expected from the spiral arm structure and revealed a new low-mass companion at 6.45" separation. We include an appendix which sets out clear definitions of Q, U, Qphi, Uphi, correcting some confusion and errors in the literature. △ Less

Submitted 8 January, 2019; originally announced January 2019.

Comments: Accepted by the Astrophysical Journal on Dec 15, 2018

A transitional disk around an intermediate mass star in the sparse population of the Orion OB1 Association

Authors: Alice Pérez-Blanco, Karina Maucó, Jesús Hernández, Nuria Calvet, Catherine Espaillat, Melissa McClure, Cesar Briceño, Connor Robinson, Daniel Feldman, Luis Villarreal, Paola D'Alessio

Abstract: We present a detailed study of the disk around the intermediate mass star SO 411, aiming to explain the spectral energy distribution of this star. We show that this is a transitional disk truncated at $\sim$11 au, with $\sim$0.03 lunar masses of optically thin dust inside the cavity. Gas also flows through the cavity, since we find that the disk is still accreting mass onto the star, at a rate of… ▽ More We present a detailed study of the disk around the intermediate mass star SO 411, aiming to explain the spectral energy distribution of this star. We show that this is a transitional disk truncated at $\sim$11 au, with $\sim$0.03 lunar masses of optically thin dust inside the cavity. Gas also flows through the cavity, since we find that the disk is still accreting mass onto the star, at a rate of $\sim 5x10^{-9}$ Msun/yr. Until now, SO 411 has been thought to belong to the $\sim$3 Myr old {$σしぐま$} Orionis cluster. However, we analyzed the second Gaia Data Release in combination with kinematic data previously reported, and found that SO 411 can be associated with an sparse stellar population located in front of the {$σしぐま$} Orionis cluster. If this is the case, then SO 411 is older and even more peculiar, since primordial disks in this stellar mass range are scarce for ages $>$5 Myr. Analysis of the silicate 10$μみゅー$m feature of SO 411 indicates that the observed feature arises at the edge of the outer disk, and displays a very high crystallinity ratio of $\sim$0.5, with forsterite the most abundant silicate crystal. The high forsterite abundance points to crystal formation in non-equilibrium conditions. The PAH spectrum of SO 411 is consistent with this intermediate state between the hot and luminous Herbig Ae and the less massive and cooler T Tauri stars. Analysis of the 7.7$μみゅー$m PAH feature indicates that small PAHs still remain in the SO 411 disk. △ Less

Submitted 25 September, 2018; originally announced September 2018.

Comments: Accepted in the Astrophysical Journal (17 pages, 9 Figures)

Magellan Adaptive Optics Imaging of PDS 70: Measuring the Mass Accretion Rate of a Young Giant Planet within a Gapped Disk

Authors: Kevin Wagner, Katherine B. Follette, Laird M. Close, Dániel Apai, Aidan Gibbs, Miriam Keppler, André Müller, Thomas Henning, Markus Kasper, Ya-Lin Wu, Joseph Long, Jared Males, Katie Morzinski, Melissa McClure

Abstract: PDS 70b is a recently discovered and directly imaged exoplanet within the wide ($\gtrsim$40 au) cavity around PDS 70 (Keppler et al. 2018, Müller et al. 2018). Ongoing accretion onto the central star suggests that accretion onto PDS 70b may also be ongoing. We present the first high contrast images at H$αあるふぁ$ (656 nm) and nearby continuum (643 nm) of PDS 70 utilizing the MagAO system. The combination… ▽ More PDS 70b is a recently discovered and directly imaged exoplanet within the wide ($\gtrsim$40 au) cavity around PDS 70 (Keppler et al. 2018, Müller et al. 2018). Ongoing accretion onto the central star suggests that accretion onto PDS 70b may also be ongoing. We present the first high contrast images at H$αあるふぁ$ (656 nm) and nearby continuum (643 nm) of PDS 70 utilizing the MagAO system. The combination of these filters allows for the accretion rate of the young planet to be inferred, as hot infalling hydrogen gas will emit strongly at H$αあるふぁ$ over the optical continuum. We detected a source in H$αあるふぁ$ at the position of PDS 70b on two sequential nights in May 2018, for which we establish a false positive probability of $<$0.1%. We conclude that PDS 70b is a young, actively accreting planet. We utilize the H$αあるふぁ$ line luminosity to derive a mass accretion rate of $\dot M= 10^{-8\pm1}$ M$_{Jup}/yr$, where the large uncertainty is primarily due to the unknown amount of optical extinction from the circumstellar and circumplanetary disks. PDS 70b represents the second case of an accreting planet interior to a disk gap, and is among the early examples of a planet observed during its formation. △ Less

Submitted 15 December, 2020; v1 submitted 27 July, 2018; originally announced July 2018.

Comments: 6 pages, 5 figures, accepted to ApJL

Revised SED of the triple protostellar system VLA 1623-2417

Authors: Nadia M. Murillo, Daniel Harsono, Melissa McClure, Shih-Ping Lai, Michiel R. Hogerheijde

Abstract: VLA 1623$-$2417 is a triple protostellar system deeply embedded in Ophiuchus A. Sources A and B have a separation of 1.1", making their study difficult beyond the submillimeter regime. Lack of circumstellar gas emission suggested that VLA 1623$-$2417 B has a very cold envelope and is much younger than source A, generally considered the prototypical Class 0 source. We explore the consequences of ne… ▽ More VLA 1623$-$2417 is a triple protostellar system deeply embedded in Ophiuchus A. Sources A and B have a separation of 1.1", making their study difficult beyond the submillimeter regime. Lack of circumstellar gas emission suggested that VLA 1623$-$2417 B has a very cold envelope and is much younger than source A, generally considered the prototypical Class 0 source. We explore the consequences of new ALMA Band 9 data on the spectral energy distribution (SED) of VLA 1623$-$2417 and their inferred nature. Using dust continuum observations spanning from centimeter to near-infrared wavelengths, the SED of each component in VLA 1623$-$2417 is constructed and analysed. The ALMA Band 9 data presented here show that the SED of VLA 1623$-$2417 B does not peak at 850 $μみゅー$m as previously expected, but instead presents the same shape as VLA 1623$-$2417 A at wavelengths shorter than 450 $μみゅー$m. The results presented here indicate that the previous assumption that the flux in $Herschel$ and Spitzer observations is solely dominated by VLA 1623$-$2417 A is not valid, and instead, VLA 1623$-$2417 B most likely contributes a significant fraction of the flux at $λらむだ~<$ 450 $μみゅー$m. These results, however, do not explain the lack of circumstellar gas emission and puzzling nature of VLA 1623$-$2417 B. △ Less

Submitted 2 July, 2018; originally announced July 2018.

Comments: 6 pages, 2 figures, accepted to A&A letters to the editor

Journal ref: A&A 615, L14 (2018)

The Evolution of Protoplanetary Disks: Probing the Inner Disk of Very Low Accretors

Authors: Thanawuth Thanathibodee, Nuria Calvet, Gregory Herczeg, Cesar Briceno, Catherine Clark, Megan Reiter, Laura Ingleby, Melissa McClure, Karina Mauco, Jesus Hernandez

Abstract: We report FUV, optical, and NIR observations of three T Tauri stars in the Orion OB1b subassociation with H$αあるふぁ$ equivalent widths consistent with low or absent accretion and various degrees of excess flux in the mid-infrared. We aim to search for evidence of gas in the inner disk in HST ACS/SBC spectra, and to probe the accretion flows onto the star using H$αあるふぁ$ and He I $λらむだ$10830 in spectra obtained… ▽ More We report FUV, optical, and NIR observations of three T Tauri stars in the Orion OB1b subassociation with H$αあるふぁ$ equivalent widths consistent with low or absent accretion and various degrees of excess flux in the mid-infrared. We aim to search for evidence of gas in the inner disk in HST ACS/SBC spectra, and to probe the accretion flows onto the star using H$αあるふぁ$ and He I $λらむだ$10830 in spectra obtained at the Magellan and SOAR telescopes. At the critical age of 5 Myr, the targets are at different stages of disk evolution. One of our targets is clearly accreting, as shown by redshifted absorption at free-fall velocities in the He I line and wide wings in H$αあるふぁ$; however, a marginal detection of FUV H$_2$ suggests that little gas is present in the inner disk, although the spectral energy distribution indicates that small dust still remains close to the star. Another target is surrounded by a transitional disk, with an inner cavity in which little sub-micron dust remains. Still, the inner disk shows substantial amounts of gas, accreting onto the star at a probably low, but uncertain rate. The third target lacks both a He I line or FUV emission, consistent with no accretion or inner gas disk; its very weak IR excess is consistent with a debris disk. Different processes occurring in targets with ages close to the disk dispersal time suggest that the end of accretion phase is reached in diverse ways. △ Less

Submitted 14 May, 2018; originally announced May 2018.

Comments: 13 pages, 3 figures. Accepted for publication in The Astrophysical Journal

Journal ref: ApJ (2018) 861, 73

ResFrac Technical Writeup

Authors: Mark McClure, Charles Kang, Soma Medam, Chris Hewson, Egor Dontsov, Ankush Singh, Carlo Peruzzo, Elizaveta Gordeliy

Abstract: ResFrac is a combined hydraulic fracturing, reservoir, and hydraulic fracturing simulator. It describes multiphase fluid flow (black oil or compositional), proppant transport, transport of non-Newtonian fluid additives, and thermal transport. It also includes stress shadowing from fracture propagation and porothermoelastic responses from pressure change in the matrix. It uses constitutive equation… ▽ More ResFrac is a combined hydraulic fracturing, reservoir, and hydraulic fracturing simulator. It describes multiphase fluid flow (black oil or compositional), proppant transport, transport of non-Newtonian fluid additives, and thermal transport. It also includes stress shadowing from fracture propagation and porothermoelastic responses from pressure change in the matrix. It uses constitutive equations that smoothly transition between equations for flow through an open crack to flow through a closed crack (with or without proppant). This document provides a detailed technical description of the code, along with validation simulations to confirm numerical accuracy. △ Less

Submitted 22 April, 2024; v1 submitted 5 April, 2018; originally announced April 2018.

The Orbit of the Companion to HD 100453A: Binary-Driven Spiral Arms in a Protoplanetary Disk

Authors: Kevin Wagner, Ruobing Dong, Patrick Sheehan, Daniel Apai, Markus Kasper, Melissa McClure, Katie M. Morzinski, Laird Close, Jared Males, Phil Hinz, Sascha P. Quanz, Jeffrey Fung

Abstract: HD 100453AB is a 10+/-2 Myr old binary whose protoplanetary disk was recently revealed to host a global two-armed spiral structure. Given the relatively small projected separation of the binary (1.05", or ~108 au), gravitational perturbations by the binary seemed to be a likely driving force behind the formation of the spiral arms. However, the orbit of these stars remained poorly understood, whic… ▽ More HD 100453AB is a 10+/-2 Myr old binary whose protoplanetary disk was recently revealed to host a global two-armed spiral structure. Given the relatively small projected separation of the binary (1.05", or ~108 au), gravitational perturbations by the binary seemed to be a likely driving force behind the formation of the spiral arms. However, the orbit of these stars remained poorly understood, which prevented a proper treatment of the dynamical influence of the companion on the disk. We observed HD 100453AB between 2015-2017 utilizing extreme adaptive optics systems on the Very Large Telescope and Magellan Clay Telescope. We combined the astrometry from these observations with published data to constrain the parameters of the binary's orbit to a=1.06"+/-0.09", e=0.17+/-0.07, and i=32.5+/- 6.5 degrees. We utilized publicly available ALMA CO data to constrain the inclination of the disk to i~28 degrees, which is relatively co-planar with the orbit of the companion and consistent with previous estimates from scattered light images. Finally, we input these constraints into hydrodynamical and radiative transfer simulations to model the structural evolution of the disk. We find that the spiral structure and truncation of the circumprimary disk in HD 100453 are consistent with a companion-dirven origin. Furthermore, we find that the primary star's rotation, its outer disk, and the companion exhibit roughly the same direction of angular momentum, and thus the system likely formed from the same parent body of material. △ Less

Submitted 11 January, 2018; originally announced January 2018.

Comments: 28 pages, 11 figures, Accepted to ApJ

An accurate and efficient method for calculating fluid exchange between fractures and matrix with a non-conforming mesh

Authors: Mark McClure

Abstract: I propose the 1D Subgrid Method for calculating fluid exchange between fractures and matrix with a non-conforming mesh. The method is demonstrated on linear and radial flow problems. The method has negligible computational cost and is accurate at all levels of mesh refinement. I propose the 1D Subgrid Method for calculating fluid exchange between fractures and matrix with a non-conforming mesh. The method is demonstrated on linear and radial flow problems. The method has negligible computational cost and is accurate at all levels of mesh refinement. △ Less

Submitted 10 September, 2017; v1 submitted 7 September, 2017; originally announced September 2017.

Measuring Catastrophic Forgetting in Neural Networks

Authors: Ronald Kemker, Marc McClure, Angelina Abitino, Tyler Hayes, Christopher Kanan

Abstract: Deep neural networks are used in many state-of-the-art systems for machine perception. Once a network is trained to do a specific task, e.g., bird classification, it cannot easily be trained to do new tasks, e.g., incrementally learning to recognize additional bird species or learning an entirely different task such as flower recognition. When new tasks are added, typical deep neural networks are… ▽ More Deep neural networks are used in many state-of-the-art systems for machine perception. Once a network is trained to do a specific task, e.g., bird classification, it cannot easily be trained to do new tasks, e.g., incrementally learning to recognize additional bird species or learning an entirely different task such as flower recognition. When new tasks are added, typical deep neural networks are prone to catastrophically forgetting previous tasks. Networks that are capable of assimilating new information incrementally, much like how humans form new memories over time, will be more efficient than re-training the model from scratch each time a new task needs to be learned. There have been multiple attempts to develop schemes that mitigate catastrophic forgetting, but these methods have not been directly compared, the tests used to evaluate them vary considerably, and these methods have only been evaluated on small-scale problems (e.g., MNIST). In this paper, we introduce new metrics and benchmarks for directly comparing five different mechanisms designed to mitigate catastrophic forgetting in neural networks: regularization, ensembling, rehearsal, dual-memory, and sparse-coding. Our experiments on real-world images and sounds show that the mechanism(s) that are critical for optimal performance vary based on the incremental training paradigm and type of data being used, but they all demonstrate that the catastrophic forgetting problem has yet to be solved. △ Less

Submitted 9 November, 2017; v1 submitted 7 August, 2017; originally announced August 2017.

Comments: To appear in AAAI 2018

An Incipient Debris Disk in the Chamaeleon I Cloud

Authors: C. C. Espaillat, Á. Ribas, M. K. McClure, J. Hernández, J. E. Owen, N. Avish, N. Calvet, R. Franco-Hernández

Abstract: The point at which a protoplanetary disk becomes a debris disk is difficult to identify. To better understand this, here we study the $\sim$40~AUえーゆー separation binary T~54 in the Chamaeleon I cloud. We derive a K5 spectral type for T~54~A (which dominates the emission of the system) and an age of $\sim$2~Myr. However, the dust disk properties of T~54 are consistent with those of debris disks seen aro… ▽ More The point at which a protoplanetary disk becomes a debris disk is difficult to identify. To better understand this, here we study the $\sim$40~AUえーゆー separation binary T~54 in the Chamaeleon I cloud. We derive a K5 spectral type for T~54~A (which dominates the emission of the system) and an age of $\sim$2~Myr. However, the dust disk properties of T~54 are consistent with those of debris disks seen around older and earlier-type stars. At the same time, T~54 has evidence of gas remaining in the disk as indicated by [Ne II], [Ne III], and [O I] line detections. We model the spectral energy distribution of T~54 and estimate that $\sim$3$\times$10$^{-3}$ Earth-masses of small dust grains ($<$0.25~$μみゅー$m) are present in an optically thin circumbinary disk along with at least $\sim$3$\times$10$^{-7}$ Earth-masses of larger ($>$10~$μみゅー$m) grains within a circumprimary disk. Assuming a solar-like mixture, we use Ne line luminosities to place a minimum limit on the gas mass of the disk ($\sim$3$\times$10$^{-4}$ Earth-masses) and derive a gas-to-dust mass ratio of $\sim$0.1. We do not detect substantial accretion, but we do see H$αあるふぁ$ in emission in one epoch, suggestive that there may be intermittent dumping of small amounts of matter onto the star. Considering the low dust mass, the presence of gas, and young age of T~54, we conclude that this system is on the bridge between the protoplanetary and debris disk stages. △ Less

Submitted 15 June, 2017; originally announced June 2017.

Comments: 15 pages, 9 figures, accepted to ApJ

Polarized disk emission from Herbig Ae/Be stars observed using Gemini Planet Imager: HD 144432, HD 150193, HD 163296, and HD 169142

Authors: John D. Monnier, T. J. Harries, A. Aarnio, F. Adams, S. Andrews, N. Calvet, C. Espaillat, L. Hartmann, S. Hinkley, S. Kraus, M. McClure, R. Oppenheimer, M. Perrin, D. Wilner

Abstract: In order to look for signs of on-going planet formation in young disks, we carried out the first J-band polarized emission imaging of the Herbig Ae/Be stars HD 150193, HD 163296, and HD 169142 using the Gemini Planet Imager (GPI), along with new H band observations of HD 144432. We confirm the complex "double ring" structure for the nearly face-on system HD 169142 first seen in H-band, finding the… ▽ More In order to look for signs of on-going planet formation in young disks, we carried out the first J-band polarized emission imaging of the Herbig Ae/Be stars HD 150193, HD 163296, and HD 169142 using the Gemini Planet Imager (GPI), along with new H band observations of HD 144432. We confirm the complex "double ring" structure for the nearly face-on system HD 169142 first seen in H-band, finding the outer ring to be substantially redder than the inner one in polarized intensity. Using radiative transfer modeling, we developed a physical model that explains the full spectral energy distribution (SED) and J- and H-band surface brightness profiles, suggesting that the differential color of the two rings could come from reddened starlight traversing the inner wall and may not require differences in grain properties. In addition, we clearly detect an elongated, off-center ring in HD 163296 (MWC 275), locating the scattering surface to be 18 AUえーゆー above the midplane at a radial distance of 77 AUえーゆー, co-spatial with a ring seen at 1.3mm by ALMA linked to the CO snow line. Lastly, we report a weak tentative detection of scattered light for HD 150193 (MWC 863) and a non-detection for HD 144432; the stellar companion known for each of these targets has likely disrupted the material in the outer disk of the primary star. For HD 163296 and HD 169142, the prominent outer rings we detect could be evidence for giant planet formation in the outer disk or a manifestation of large-scale dust growth processes possibly related to snow-line chemistry. △ Less

Submitted 15 February, 2017; originally announced February 2017.

Comments: Accepted Astrophysical Journal

Photophoretic Levitation and Trapping of Dust in the Inner Regions of Protoplanetary Disks

Authors: Colin P. McNally, Melissa K. McClure

Abstract: In protoplanetary disks, the differential gravity-driven settling of dust grains with respect to gas and with respect to grains of varying sizes determines the observability of grains, and sets the conditions for grain growth and eventually planet formation. In this work we explore the effect of photophoresis on the settling of large dust grains in the inner regions of actively accreting protoplan… ▽ More In protoplanetary disks, the differential gravity-driven settling of dust grains with respect to gas and with respect to grains of varying sizes determines the observability of grains, and sets the conditions for grain growth and eventually planet formation. In this work we explore the effect of photophoresis on the settling of large dust grains in the inner regions of actively accreting protoplanetary disks. Photophoretic forces on dust grains result from the collision of gas molecules with differentially heated grains. We undertake one dimensional dust settling calculations to determine the equilibrium vertical distribution of dust grains in each column of the disk. In the process we introduce a new treatment of the photophoresis force which is consistent at all optical depths with the representation of the radiative intensity field in a two-stream radiative transfer approximation. The levitation of large dust grains creates a photophoretic dust trap several scale heights above the mid-plane in the inner regions of the disk where the dissipation of accretion energy is significant. We find that differential settling of dust grains is radically altered in these regions of the disk, with large dust grains trapped in a layer below the stellar irradiation surface in where the dust to gas mass ratio can be enhanced by a factor of a hundred for the relevant particles. The photophoretic trapping effect has a strong dependence on particle size and porosity. △ Less

Submitted 9 January, 2017; v1 submitted 14 November, 2016; originally announced November 2016.

Comments: 15 pages, 5 figures, updated to match published version

Journal ref: The Astrophysical Journal, Volume 834, Issue 1, article id. 48, 13 pp. (2017)

Bed load proppant transport during slickwater hydraulic fracturing: insights from comparisons between published laboratory data and correlations for sediment and pipeline slurry transport

Authors: Mark W. McClure

Abstract: Bed load transport is the movement of particles along the top of a bed through rolling, saltation, and suspension created by turbulent lift above the bed surface. In recent years, there has been a resurgence of interest in the idea that bed load transport is significant for proppant transport during hydraulic fracturing. However, scaling arguments suggest that bed load transport is only dominant i… ▽ More Bed load transport is the movement of particles along the top of a bed through rolling, saltation, and suspension created by turbulent lift above the bed surface. In recent years, there has been a resurgence of interest in the idea that bed load transport is significant for proppant transport during hydraulic fracturing. However, scaling arguments suggest that bed load transport is only dominant in the laboratory and is negligible at the field scale. I review laboratory experiments that have been used to develop concepts of bed load transport in hydraulic fracturing. I also review the scaling arguments and laboratory results that have been used to argue that viscous drag, not bed load transport, is dominant at the field scale. I compare literature correlations for fluvial sediment transport and for pipeline slurry transport with published laboratory data on proppant transport in slot flow. The comparisons indicate that fluvial transport correlations are suitable for predicting the rate of sediment erosion due to flow of proppant-free fluid over a bed. The pipeline slurry correlations are suitable for predicting the ability of proppant that is already in suspension to flow in bed transport without settling, but only if the aspect ratio of the flowing region in the slot is close to unity. This can occur at the laboratory scale, but not in the field. The comparison indicates that at low rates of proppant flow, the equilibrium bed height in a laboratory proppant transport experiment can be predicted from pipeline slurry correlations. However, as the volumetric flow rate and the aspect ratio of the flowing region in the slot increase, the applicability of the pipeline slurry correlations breaks down. Laboratory experiments using higher volumetric flow rate indicate that bed load transport rates are too low to be significant at the field scale, even if the flow velocity is very high. △ Less

Submitted 12 November, 2016; originally announced November 2016.

Identifying potentially induced seismicity and assessing statistical significance in Oklahoma and California

Authors: Mark McClure, Riley Gibson, Kitkwan Chiu, Rajesh Ranganath

Abstract: We develop a statistical method for identifying induced seismicity from large datasets and apply the method to decades of wastewater disposal and seismicity data in California and Oklahoma. The method is robust against a variety of potential pitfalls. The study regions are divided into gridblocks. We use a longitudinal study design, seeking associations between seismicity and wastewater injection… ▽ More We develop a statistical method for identifying induced seismicity from large datasets and apply the method to decades of wastewater disposal and seismicity data in California and Oklahoma. The method is robust against a variety of potential pitfalls. The study regions are divided into gridblocks. We use a longitudinal study design, seeking associations between seismicity and wastewater injection along time-series within each gridblock. The longitudinal design helps control for non-random application of wastewater injection. We define a statistical model that is flexible enough to describe the seismicity observations, which have temporal correlation and high kurtosis. In each gridblock, we find the maximum likelihood estimate for a model parameter that relates induced seismicity hazard to total volume of wastewater injected each year. To assess significance, we compute likelihood ratio test statistics in each gridblock and each state, California and Oklahoma. Resampling is used to empirically derive reference distributions used to estimate p-values from the likelihood ratio statistics. In Oklahoma, the analysis finds with extremely high confidence that seismicity associated with wastewater disposal (or other related activities, such as reservoir depletion) has occurred. In California, the analysis finds that seismicity associated with wastewater disposal has probably occurred, but the result is not strong enough to be conclusive. We identify areas where temporal association between wastewater disposal and seismicity is apparent. Our method could be applied to other datasets, extended to identify risk factors that increase induced seismic hazard, or modified to test alternative statistical models for natural and induced seismicity. △ Less

Submitted 15 December, 2016; v1 submitted 10 November, 2016; originally announced November 2016.

Comments: Contact Mark McClure (mark.w.mccl@gmail.com) for access to the raw data, the code used to perform the analysis, or the Supporting Material

Journal ref: JGR Solid Earth 2017