Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “global models”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,135 records · Page 63

Beta-Delayed Neutron Emission of 𝑁=84 132Cd

Using the time-of-flight technique, we measured the beta-delayed neutron emission of 132 Cd . From our large-scale shell model (LSSM) calculation using the N3⁢LO interaction [, Phys. Rev. Lett. 131, 022501 (2023)], we suggest the decay is dominated by the transformation of a neutron in the 𝑔7/2 orbital, deep below the Fermi surface, into a proton in the 𝑔9/2 orbital. We compare the beta-decay half-lives and neutron branching ratios of nuclei with 𝑍<50 and 𝑁 ≥82 obtained with our LSSM with those of leading “global” models such as finite-range droplet model (FRDM). Our calculations match known half-lives and neutron branching ratios well and suggest that current leading models overestimate the yet-to-be-measured half-lives. Our model, backed by the 132 Cd decay data presented here, offers robust predictive power for nuclei of astrophysical interest such as 𝑟-process waiting points.

Madurga, M [Department of Physics and Astronomy, U↗

Improving Predictability of Mixed-Phase Clouds and Aerosol Interactions in the Community Earth System Model (CESM) with ARM Measurements

The objective of this project is to improve the simulation and predictability of mixed-phase clouds and aerosol interactions in the Community Earth System Model (CESM) through comparisons with the DOE ARM observations. There are four major goals of the proposed study: (1) Improve the representation of ice microphysical processes in mixed-phase clouds; (2) Develop a long-term multi-sensor mixed-phase cloud dataset; (3) Test the performance of improved ice microphysics in CESM-CAM5 with the ARM data; and (4) Examine mixed-phase cloud microphysics-aerosol-dynamics-radiation interactions in CESM-CAM5. In this project, we have (1) Improved the representation of ice microphysical processes in mixed-phase clouds that include ice nucleation, ice depositional growth through the Wegener–Bergeron–Findeisen (WBF) process, feedbacks between cloud microphysics, dynamics and radiation, and provides links to aerosol properties in the Community Earth System Model (CESM)-Community Atmosphere Model version 5 (CAM5); (2) Developed a long-term multi-sensor mixed-phase cloud dataset. The mixed-phase cloud retrieval algorithms are improved by refining the phase determinations and ice number concentration retrievals. A multi-year mixed-phase cloud observation dataset (1999-2004; October 2013-February 2017) is regenerated at the ARM Barrow site with the improved algorithms. The dataset has been used in the global model evaluation to improve the model parameterizations; (3) Tested the performance of improved ice microphysics in CESM-CAM5 with the ARM data. CESM-CAM5 is run in the DOE Cloud-Associated Parameterizations Testbed (CAPT) and in the single-column model (SCM) mode to facilitate comparison with ARM observations. Multi-year CAPT simulations of seasonal variations of cloud microphysical properties, cloud liquid water path (LWP) and ice water path (IWP), cloud longwave and shortwave forcings, cloud occurrence frequency are evaluated against the ARM multi-sensor mixed-phase cloud retrievals; and (4) Examined mixed-phase cloud microphysics-aerosol-dynamics-radiation interactions in CESM-CAM5 that include the impacts of different ice nucleation mechanisms, sensitivity to dust IN concentration, treatment of the WBF process on mixed-phase cloud properties. Aerosol effect on mixed-phase clouds through the liquid phase (droplet activation) and ice phase processes (e.g., the glaciation indirect effect) is investigated.

54 ENVIRONMENTAL SCIENCES↗

A Numerical Model for the Dynamics of Pickup Ions Outside the Heliopause and IBEX “Ribbon” Observation

The leading mechanism for the origin of the Interstellar Boundary Explorer (IBEX) “ribbon” of enhanced energetic neutral atoms (ENAs) from the outer heliosphere is the so-called secondary ENA process. In this mechanism, neutralized ions escape the heliosphere, and after several times of charge exchange collision,s some propagate back toward Earth in directions nearly perpendicular to the local interstellar magnetic field. However, the physical processes governing the distribution of the pickup ions (PUIs) outside the heliopause are still unclear. In this study, we build a new global model of the IBEX ribbon, where a key component is to calculate the dynamics of the PUIs outside the heliopause by solving the gyrophase-averaged focused transport equations on top of the background magnetohydrodynamic-kinetic model for the heliosphere-interstellar medium interaction. We discuss how the properties of the simulated ribbon change with different scattering parameters and show simulation results of the ENA sky maps, ribbon centers, ribbon profiles at several azimuthal slices, the ribbon ENA source region, and the velocity distribution of PUIs in the ENA source region for different scattering frequencies. Our results show that the model can provide reasonable ribbon sky maps comparable to the observed ribbon when the PUI mean free path is on the order of several thousand astronomical units.

79 ASTRONOMY AND ASTROPHYSICS↗

Climate Scenarios for the NASA / USAID SERVIR Project: Challenges for Multiple Planning Horizons

SERVIR, an acronym meaning "to serve" in Spanish, is a joint venture between NASA and the U.S. Agency for International Development (USAID) which provides satellite-based Earth observation data, modeling, and science applications to help developing nations in Central America, East Africa and the Himalayas improve environmental decision making. Anticipating climate variability / climate change impacts has now become an important component of the SERVIR efforts to build capacity in these regions. Uncertainty in hydrometeorological components of climate variations and exposure to extreme events across scales from weather to climate are of particular concern. We report here on work to construct scenarios or outlooks that are being developed as input drivers for decision support systems (DSSs) in a variety of settings. These DSSs are being developed jointly by a broad array NASA Applied Science Team (AST) Investigations and user communities in the three SERVIR Hub Regions, Central America, East Africa and the Himalayas. Issues span hydrologic / water resources modeling, agricultural productivity, and forest carbon reserves. The scenarios needed for these efforts encompass seasonal forecasts, interannual outlooks, and likely decadal / multi-decadal trends. Providing these scenarios across the different AST efforts enables some level of integration in considering regional responses to climate events. We will discuss a number of challenges in developing this continuum of scenarios including the identification and "mining" of predictability, addressing multiple continental regions, issues of downscaling global model integrations to regional / local applications (i.e. hydrologic and crop modeling). We compare / contrast the role of the U.S. National Multi- Model Experiment initiative in seasonal forecasts and the CMIP-5 climate model experiments in supporting these efforts. Examples of these scenarios, their use, and an assessment of their utility as well as limitations will be presented.

Robertson, Franklin R.↗

Helioseismic Measurements of Convective Power in Solar Cycle 24

Constraining the parameters under which convection in the solar interior operates has important implications for describing how energy is transported by plasma motions, and various physical models have been employed to provide some theoretical estimates on the expected power spectrum. Past attempts to measure the convective power distributed among large spatial scales have, however, found differing and incompatible values. Here, we present measurements of the convective power spectrum in the upper convection zone for Carrington rotations in Solar Cycle 24 obtained from the helioseismic signal corresponding to East-West flows. We also perform calibration on synthetic data using the global acoustic GALE code to make assessments of the flow velocities at various length scales without the need for performing inversions. This allows us to derive the convective power spectrum from the flow maps and to compare with the results from the raw travel times. These results are compared against predictions of convective power in global models of convection produced by the EULAG code. Finally, we show how the steps in our analysis procedure (for example data segmentation, filtering, etc.) affect our estimates of the convective power by comparing with the synthetic data from the GALE code.

Heliophysics↗

Mapping Modeled Aerosol Species to Measured Lidar Ratios for the MIRA Project

Overarching goal: Create Climatological Lidar Ratio Maps for CALIPSO - Since CALIPSO is an elastic backscatter lidar, it is necessary to select a priori aerosol lidar ratios prior to computing extinction profiles. - Presently, the CALIPSO lidar ratio selection process uses a single lidar ratio for each of the 7 CALIPSO aerosol types in the troposphere. - Consequently, regional and seasonal variability of the lidar ratio is very limited in the present CALIPSO algorithm. - Here, we describe an approach for creating new climatological lidar ratio maps for each of the CALIPSO aerosol types using MODIS-CALIPSO constrained retrievals and global model simulations. - The purpose is to eventually replace the seven single-valued lidar ratios for the CALIPSO types with latitude- and longitude-dependent maps for each of the seven types.

G. L. Schuster↗

A New Paradigm for Assessing the Role of Agriculture in the Climate System and in Climate Change

This paper discusses the diverse climate forcings that impact agricultural systems, and contrasts the current paradigm of using global models downscaled to agricultural areas (a top-down approach) with a new paradigm that first assesses the vulnerability of agricultural activities to the spectrum of environmental risk including climate (a bottom-up approach). To illustrate the wide spectrum of climate forcings, regional climate forcings are presented including land-use/land-cover change and the influence of aerosols on radiative and biogeochemical fluxes and cloud/precipitation processes, as well as how these effects can be teleconnected globally. Examples are presented of the vulnerability perspective, along with a small survey of the perceived drought impacts in a local area, in which a wide range of impacts for the same precipitation deficits are found. This example illustrates why agricultural assessments of risk to climate change and variability and of other environmental risks should start with a bottom-up perspective.

Pielke, Roger A., Sr.↗

Recent Antarctic Sea Ice Extent Variability and Its Relation to Atmospheric Circulation in Reanalyses and Seasonal Forecasts

Sea ice is considered an important indicator of climate change. Interest has focused on the trend towards reduced sea ice cover in the Arctic Ocean, particularly during summer months. In the Southern Hemisphere, sea ice has also received attention due to a trend towards increasing cover, in opposition to the Northern Hemisphere and to climate model projections. This all changed in the austral spring 2016, when the Southern Ocean sea ice cover dropped dramatically, and has remained consistently at below average values up to the present. Here I review some of the issues and perceived causes of the 2016 sea ice reduction. Several studies have suggested a change in surface winds, which would allow warmer waters to move upwards in the water column and towards the surface. Results from the GMAO (Global Modeling and Assimilation Office) ocean data assimilation system are consistent with this hypothesis, and indicate a significant warming of subsurface waters occurred in the winter of 2016, and has remained in place.

Polar↗

Modeling the Frozen-In Anticyclone in the 2005 Arctic Summer Stratosphere

Immediately following the breakup of the 2005 Arctic spring stratospheric vortex, a tropical air mass, characterized by low potential vorticity (PV) and high nitrous oxide (N2O), was advected poleward and became trapped in the easterly summer polar vortex. This feature, known as a "Frozen-In Anticyclone (FrIAC)", was observed in Earth Observing System (EOS) Aura Microwave Limb Sounder (MLS) data to span the potential temperature range from approximately 580 to 1100 K (approximately 25 to 40 km altitude) and to persist from late March to late August 2005. This study compares MLS N2O observations with simulations from the Global Modeling Initiative (GMI) chemistry and transport model, the GEOS-5/MERRA Replay model, and the VanLeer Icosahedral Triangular Advection isentropic transport model to elucidate the processes involved in the lifecycle of the FrIAC which is here divided into three distinct phases. During the "spin-up phase" (March to early April), strong poleward flow resulted in a tight isolated anticyclonic vortex at approximately 70-90 deg N, marked with elevated N2O. GMI, Replay, and VITA all reliably simulted the spin-up of the FrIAC, although the GMI and Replay peak N2O values were too low. The FrIAC became trapped in the developing summer easterly flow and circulated around the polar region during the "anticyclonic phase" (early April to the end of May). During this phase, the FrIAC crossed directly over the pole between the 7th and 14th of April. The VITA and Replay simulations transported the N2O anomaly intact during this crossing, in agreement with MLS, but unrealistic dispersion of the anomaly occurred in the GMI simulation due to excessive numerical mixing of the polar cap. The vortex associated with the FrIAC was apparently resistant to the weak vertical hear during the anticyclonic phase, and it thereby protected the embedded N20 anomaly from stretching. The vortex decayed in late May due to diabatic processes, leaving the N2O anomaly exposed to horizontal and vertical wind shears during the "shearing phase" (June to August). The observed lifetime of the FrIAC during this phase is consistent with time-scales calculated from the ambient horizontal and vertical wind shear. Replay maintained the horizontal structure of the N2O anomaly similar to NILS well into August. The VITA simulation also captured the horizontal structure of the FrIAC during this phase, but VITA eventually developed fine-scale N2O structure not observed in MLS data.

Allen, D. R.↗

Methanediol from cloud-processed formaldehyde is only a minor source of atmospheric formic acid

Atmospheric formic acid is severely underpredicted by models. A recent study proposed that this discrepancy can be resolved by abundant formic acid production from the reaction ( 1 ) between hydroxyl radical and methanediol derived from in-cloud formaldehyde processing and provided a chamber-experiment-derived rate constant, k 1 = 7.5 × 10 −12 cm 3 s −1 . High-level accuracy coupled cluster calculations in combination with E,J -resolved two-dimensional master equation analyses yield k 1 = (2.4 ± 0.5) × 10 −12 cm 3 s −1 for relevant atmospheric conditions ( T = 260–310 K and P = 0–1 atm). We attribute this significant discrepancy to HCOOH formation from other molecules in the chamber experiments. More importantly, we show that reversible aqueous processes result indirectly in the equilibration on a 10 min. time scale of the gas-phase reaction HCHO + H 2 O ⇌ HOCH 2 OH (2) with a HOCH 2 OH to HCHO ratio of only ca . 2%. Although HOCH 2 OH outgassing upon cloud evaporation typically increases this ratio by a factor of 1.5–5, as determined by numerical simulations, its in-cloud reprocessing is shown using a global model to strongly limit the gas-phase sink and the resulting production of formic acid. Based on the combined findings in this work, we derive a range of 1.2–8.5 Tg/y for the global HCOOH production from cloud-derived HOCH 2 OH reacting with OH. The best estimate, 3.3 Tg/y, is about 30 times less than recently reported. The theoretical equilibrium constant K eq (2) determined in this work also allows us to estimate the Henry’s law constant of methanediol (8.1 × 10 5 M atm −1 at 280 K).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Electromagnetic Analysis of ITER Electron Cyclotron Emission Model

The ITER electron cyclotron emission (ECE) diagnostic system is located at diagnostics shielding module 2 (DSM2), equatorial port 9 (EP9), to measure electron temperature profile and electron temperature fluctuations, and also assess nonthermal electron distributions via the oblique view. Therefore, ECE has both radial and oblique views with two hot sources for calibration and two couples of mirrors for two different optical views. It shall receive that the large electromagnetic (EM) loads up to 100-MN/m3 force density due to the eddy currents on the copper mirrors generated by the short transient plasma disruption. The simplified EM analysis model using magnetic field ( B ) data and flux variations (dB/dt) method based on the worst case of plasma disruption, MD_DW_EXP16MS_CATIII, has been developed to calculate the force and moments for each ECE component and support structure for EP9 DSM2 Bay2 and Bay3. The B and dB/dt methods use the constant B and dB/dt assumed for each local bay space, but the real B and dB/dt are basically varied along the radial direction of each bay. The complex global model with the local ECE components and support structure using Maxwell transient has been developed to do EM analysis to compare the results of the B and dB/dt models. The volumetric EM force density of the whole ECE components and support structure can be used for the subsequent structural analysis. In conclusion, combined with the thermal and nuclear loads and seismic and initial loads, the ECE integration analysis can be finalized for the operation case.

Fang, J. [Princeton Plasma Physics Laboratory (PPP↗

Data Management as a Cluster Middleware Centerpiece

Through earth and space modeling and the ongoing launches of satellites to gather data, NASA has become one of the largest producers of data in the world. These large data sets necessitated the creation of a Data Management System (DMS) to assist both the users and the administrators of the data. Halcyon Systems Inc. was contracted by the NASA Center for Computational Sciences (NCCS) to produce a Data Management System. The prototype of the DMS was produced by Halcyon Systems Inc. (Halcyon) for the Global Modeling and Assimilation Office (GMAO). The system, which was implemented and deployed within a relatively short period of time, has proven to be highly reliable and deployable. Following the prototype deployment, Halcyon was contacted by the NCCS to produce a production DMS version for their user community. The system is composed of several existing open source or government-sponsored components such as the San Diego Supercomputer Center s (SDSC) Storage Resource Broker (SRB), the Distributed Oceanographic Data System (DODS), and other components. Since Data Management is one of the foremost problems in cluster computing, the final package not only extends its capabilities as a Data Management System, but also to a cluster management system. This Cluster/Data Management System (CDMS) can be envisioned as the integration of existing packages.

Zero, Jose↗

Decoupling in the Vertical Shape of HCHO During a Sea Breeze Event: the Effect on Trace Gas Satellite Retrievals and Column-to-Surface Translation

The effect of sea breeze circulation on stratification of the vertical formaldehyde (HCHO) concentration vertical profiles is explored using a regional atmospheric chemical transport model (CTM) for three synoptically stagnant days focused on the east coast of the U.S in June 2018. During this event, a significant thermal contrast between the Atlantic Ocean and the terrestrial regions (12–17 °C), observed by moderate resolution imaging spectroradiometer (MODIS) and well-captured by the WRF-CMAQ model (15–18 °C), is conducive to monsoon-like flow, perpendicular to the shorelines, carrying clean marine air masses over the land within a few hundreds of meters above the surface. In contrast, the westerly continental polluted air masses prevail in higher altitudes. These two conflicting flows result in atypical vertical shapes of HCHO concentrations increasing with altitude. This decoupling pattern is so pronounced that we observe total column HCHO negatively correlate with surface concentrations. Comparisons of an accredited global model, GEOS-CF, to surface wind measurements and MODIS skin temperature indicate its poor representation of the sea breeze timing and strength, resulting in GEOS-CF HCHO vertical shapes being drastically different from the WRF-CMAQ. Based on radiative transfer calculations, the differences in the vertical distribution of HCHO between the WRF-CMAQ and that of GEOS-CF in the first 3 km are sufficient to induce a 20–30% error in air mass factors (thus total vertical HCHO column abundances). Through an experiment involving converting HCHO total columns to surface mixing ratios, we demonstrate that GEOS-CF allocates noticeably more HCHO molecules (40–150%) to the surface layer due to the misrepresentation of the vertical shape of HCHO during the sea breeze event. It is known that a significant fraction of the human population lives in coastal areas prone to detrimental effects caused by air pollution, and elevated pollutant concentrations usually occur in synoptically stagnant atmospheric conditions when local circulation patterns come into play; accordingly, our experiments emphasize the importance of the effect a priori profiles can have on satellite-derived applications under such conditions. To ensure that the quantitative representation of satellite-based trace gas retrievals on a daily basis is trustworthy and useable for air quality applications, atmospheric models providing a priori profiles for satellite retrievals should be well-tuned to reproduce complex local circulation such as sea-land breezes.

sea breeze↗

A Surface Radiation Balance Data Set from Siple Dome in West Antarctica for Atmospheric and Climate Model Evaluation

A field campaign at Siple Dome in West Antarctica during the austral summer 2019-2020 offers an opportunity to evaluate climate model performance, particularly cloud microphysical simulation. Over Antarctic ice sheets and ice shelves, clouds are a major regulator of the surface energy balance, and in the warm season their presence occasionally induces surface melt that can gradually weaken an ice shelf structure. This dataset from Siple Dome, obtained using transportable and solar-poweredequipment, includes surface energy balance measurements, meteorology and cloud remote sensing. To demonstrate how these data can be used to evaluate model performance, comparisons are made with meteorological reanalysis known to give generally good performance over Antarctica (ERA5). Surface albedo measurements show expected variability with observed cloud amount, and can be used to evaluate a model's snowpack parameterization. One case study discussed involves a squall with northerly winds, during which ERA5 fails to produce cloud cover throughout one of the days. A second case study illustrates how shortwave spectroradiometer measurements that encompass the 1.6-micron atmospheric window reveal cloud phase transitions associated with cloud lifecycle. Here, continuously precipitating mixed-phase clouds become mainly liquid water clouds from local morning through the afternoon, not reproduced by ERA5. We challenge researchers to run their various regional or global models in a manner that has the large-scale meteorology follow the conditions of this field campaign, compare cloud and radiation simulations with this Siple Dome dataset, and potentially investigate why cloud microphysical simulations or other model components might produce discrepancies with these observations.

cloud remote sensing↗

The NASA/MSFC global reference atmospheric model: 1990 version (GRAM-90). Part 1: Technical/users manual

A technical description of the NASA/MSFC Global Reference Atmospheric Model 1990 version (GRAM-90) is presented with emphasis on the additions and new user's manual descriptions of the program operation aspects of the revised model. Some sample results for the new middle atmosphere section and comparisons with results from a three dimensional circulation model are provided. A programmer's manual with more details for those wishing to make their own GRAM program adaptations is also presented.

Justus, C. G.↗

Scenario development for safety assessment in deep geologic disposal of high-level radioactive waste and spent nuclear fuel: A review

Radiation and radioactive substances result in the production of radioactive wastes which require safe management and disposal to avoid risks to human health and the environment. To ensure permanent safe disposal, the performance of a deep geological repository for radioactive waste is assessed against internationally agreed risk-based standards. Assessing postclosure safety of the future system's evolution includes screening of features, events, and processes (FEPs) relevant to the situation, their subsequent development into scenarios, and finally the development and execution of safety assessment (SA) models. Global FEP catalogs describe important natural and man-made repository system features and identify events and processes that may affect these features into the future. By combining FEPs, many of which are uncertain, different possible future system evolution scenarios are derived. Repository licensing should consider both the reference or “base” evolution as well as alternative futures that may lead to radiation release, pollution, or exposures. Scenarios are used to derive and consider both base and alternative evolutions, often through production of scenario-specific SA models and the recombination of their results into an assessment of the risk of harm. Furthermore, while the FEP-based scenario development process outlined here has evolved somewhat since its development in the 1980s, the fundamental ideas remain unchanged. A spectrum of common approaches is given here (e.g., bottom–up vs. top–down scenario development, probabilistic vs. bounding handling of uncertainty), related to how individual numerical models for possible futures are converted into a determination as to whether the system is safe (i.e., how aleatoric uncertainty and scenarios are integrated through bounding or Monte Carlo approaches).

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Analysis of the Atmospheric Distribution, Sources, and Sinks of Oxygenated Volatile Organic Chemicals Based on Measurements over the Pacific during TRACE-P

Airborne measurements of a large number of oxygenated volatile organic chemicals (OVOC) were carried out in the Pacific troposphere (0.1-12 km) in winter/spring of 2001 (24 February to 10 April). Specifically, these measurements included acetone (CH3COCH3), methylethyl ketone (CH3COC2H5, MEK), methanol (CH3OH), ethanol (C2H5OH), acetaldehyde (CH3CHO), propionaldehyde (C2H5CHO), peroxyacylnitrates (PANs) (C(sub n)H(sub 2n+1)COO2NO2), and organic nitrates (C(sub n)H(sub 2n+1)ONO2). Complementary measurements of formaldehyde (HCHO), methyl hydroperoxide (CH3OOH), and selected tracers were also available. OVOC were abundant in the clean troposphere and were greatly enhanced in the outflow regions from Asia. Background mixing ratios were typically highest in the lower troposphere and declined toward the upper troposphere and the lowermost stratosphere. Their total abundance (Summation of OVOC) was nearly twice that of nonmethane hydrocarbons (Summation of C2-C8 NMHC). Throughout the troposphere, the OH reactivity of OVOC is comparable to that of methane and far exceeds that of NMHC. A comparison of these data with western Pacific observations collected some 7 years earlier (February-March 1994) did not reveal significant differences. Mixing ratios of OVOC were strongly correlated with each other as well as with tracers of fossil and biomass/biofuel combustion. Analysis of the relative enhancement of selected OVOC with respect to CH3Cl and CO in 12 plumes originating from fires and sampled in the free troposphere (3-11 km) is used to assess their primary and secondary emissions from biomass combustion. The composition of these plumes also indicates a large shift of reactive nitrogen into the PAN reservoir thereby limiting ozone formation. A three-dimensional global model that uses state of the art chemistry and source information is used to compare measured and simulated mixing ratios of selected OVOC. While there is reasonable agreement in many cases, measured aldehyde concentrations are significantly larger than predicted. At their observed levels, acetaldehyde mixing ratios are shown to be an important source of HCHO (and HO x ) and PAN in the troposphere. On the basis of presently known chemistry, measured mixing ratios of aldehydes and PANs are mutually incompatible. We provide rough estimates of the global sources of several OVOC and conclude that collectively these are extremely large (150-500 Tg C / yr) but remain poorly quantified.

Singh, H. B.↗

Fomalhaut's Debris Disk and Planet: Constraining the Mass of Formalhaut B from Disk Morphology

Following the optical imaging of exoplanet candidate Fomalhaut b (Fom b), we present a numerical model of how Fomalhaut's debris disk is gravitationally shaped by a single interior planet. The model is simple, adaptable to other debris disks, and can be extended to accommodate multiple planets. If Fom b is the dominant perturber of the belt, then to produce the observed disk morphology it must have a mass M(sub pl) < 3M(sub J), an orbital semimajor axis a(sub pl) > 101.5AU, and an orbital eccentricity e(sub pl) = 0.11 - 0.13. These conclusions are independent of Fom b's photometry. To not disrupt the disk, a greater mass for Fom b demands a smaller orbit farther removed from the disk; thus, future astrometric measurement of Fom b's orbit, combined with our model of planet-disk interaction, can be used to determine the mass more precisely. The inner edge of the debris disk at a approximately equals 133AU lies at the periphery of Fom b's chaotic zone, and the mean disk eccentricity of e approximately equals 0.11 is secularly forced by the planet, supporting predictions made prior to the discovery of Fom b. However, previous mass constraints based on disk morphology rely on several oversimplifications. We explain why our constraint is more reliable. It is based on a global model of the disk that is not restricted to the planet's chaotic zone boundary. Moreover, we screen disk parent bodies for dynamical stability over the system age of approximately 100 Myr, and model them separately from their dust grain progeny; the latter's orbits are strongly affected by radiation pressure and their lifetimes are limited to approximately 0.1 Myr by destructive grain-grain collisions. The single planet model predicts that planet and disk orbits be apsidally aligned. Fomalhaut b's nominal space velocity does not bear this out, but the astrometric uncertainties are difficult to quantify. Even if the apsidal misalignment proves real, our calculated upper mass limit of 3 M(sub J) still holds. Parent bodies are evacuated from mean-motion resonances with Fom b; these empty resonances are akin to the Kirkwood gaps opened by Jupiter. The belt contains at least 3M(sub Earth) of solids that are grinding down to dust, their velocity dispersions stirred so strongly by Fom b that collisions are destructive. Such a large mass in solids is consistent with Fom b having formed in situ.

Chiang, E.↗