Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “sub scale”

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 415 records · Page 23

The Lower Thermospheric Winter-To-Summer Meridional Circulation: 1: Driving Mechanism

In this study, the mechanism driving the narrow lower-thermospheric winter-to-summer meridional circulation is thoroughly investigated for the first time using the Specified Dynamics configuration runs of the Whole Atmosphere Community Climate Model eXtended (SD-WACCMX) simulations and the TIMED Doppler Interferometer (TIDI) observations. The mean meridional circulation in the SD-WACCMX is qualitatively consistent with the TIDI measurements, though the magnitude in the SD-WACCMX is about 50% weaker. The lower-thermospheric winter-to-summer circulation is mainly driven by the resolved wave forcing, including the tides and internally generated inertia gravity waves (GWs). The momentum forcing from the parameterized sub-grid scale GWs is not as significant as the resolved wave forcing in driving the lower-thermospheric meridional circulation. The GW parameterization scheme in the SD-WACCMX only includes GWs with phase velocities in the range of ±45 m/s, which might result in most of the parameterized sub-grid GWs dissipating and breaking in the mesosphere and hardly impacting the lower thermosphere. Only including slow GWs in the SD-WACCMX parameterization could potentially lead to the underestimation of the meridional wind in the model. Analysis also indicates the lower-thermospheric meridional circulation is stronger in the summer hemisphere, which is attributed to the hemispheric asymmetry in the resolved wave momentum forcing. This study underlines the importance of the whole atmosphere coupling through wave propagation and dissipation. This understanding can guide the model development with an accurate representation of underlying physical processes in the mesosphere and lower thermosphere which drives the lower-thermospheric circulation as well as the overall dynamics of this region.

Jack C. Wang↗

CO2 Flux Estimation Errors Associated with Moist Atmospheric Processes

Vertical transport by moist sub-grid scale processes such as deep convection is a well-known source of uncertainty in CO2 source/sink inversion. However, a dynamical link between vertical transport, satellite based retrievals of column mole fractions of CO2, and source/sink inversion has not yet been established. By using the same offline transport model with meteorological fields from slightly different data assimilation systems, we examine sensitivity of frontal CO2 transport and retrieved fluxes to different parameterizations of sub-grid vertical transport. We find that frontal transport feeds off background vertical CO2 gradients, which are modulated by sub-grid vertical transport. The implication for source/sink estimation is two-fold. First, CO2 variations contained in moist poleward moving air masses are systematically different from variations in dry equatorward moving air. Moist poleward transport is hidden from orbital sensors on satellites, causing a sampling bias, which leads directly to small but systematic flux retrieval errors in northern mid-latitudes. Second, differences in the representation of moist sub-grid vertical transport in GEOS-4 and GEOS-5 meteorological fields cause differences in vertical gradients of CO2, which leads to systematic differences in moist poleward and dry equatorward CO2 transport and therefore the fraction of CO2 variations hidden in moist air from satellites. As a result, sampling biases are amplified and regional scale flux errors enhanced, most notably in Europe (0.43+/-0.35 PgC /yr). These results, cast from the perspective of moist frontal transport processes, support previous arguments that the vertical gradient of CO2 is a major source of uncertainty in source/sink inversion.

moist atmospheric processes↗

Examples of data assimilation in mesoscale models

The keynote address was the problem of physical initialization of mesoscale models. The classic purpose of physical or diabatic initialization is to reduce or eliminate the spin-up error caused by the lack, at the initial time, of the fully developed vertical circulations required to support regions of large rainfall rates. However, even if a model has no spin-up problem, imposition of observed moisture and heating rate information during assimilation can improve quantitative precipitation forecasts, especially early in the forecast. The two key issues in physical initialization are the choice of assimilating technique and sources of hydrologic/hydrometeor data. Another example of data assimilation in mesoscale models was presented in a series of meso-beta scale model experiments with and 11 km version of the MASS model designed to investigate the sensitivity of convective initiation forced by thermally direct circulations resulting from differential surface heating to four dimensional assimilation of surface and radar data. The results of these simulations underscore the need to accurately initialize and simulate grid and sub-grid scale clouds in meso- beta scale models. The status of the application of the CSU-RAMS mesoscale model by the NOAA Forecast Systems Lab for producing real-time forecasts with 10-60 km mesh resolutions over (4000 km)(exp 2) domains for use by the aviation community was reported. Either MAPS or LAPS model data are used to initialize the RAMS model on a 12-h cycle. The use of MAPS (Mesoscale Analysis and Prediction System) model was discussed. Also discussed was the mesobeta-scale data assimilation using a triply-nested nonhydrostatic version of the MM5 model.

Carr, Fred↗

Using Open Innovation in Reducing Risk to Crews

In the exploration of destinations outside of Earth's neighborhood, specifically Mars, scientific and engineering inquiries have occurred by two means; observations from satellites and observations by landed spacecraft. Satellite observations (Mariner, MRO, Mars Odyssey, provide global-scale spatial and temporal data while landed spacecraft (Viking, Mars Pathfinder, Spirit, Opportunity, Phoenix Mars Lander) investigate highly localized areas of the surface of the planet. In preparation for human exploration, extensive knowledge of the surface and atmospheric environments should be known before the first human leaves Earth. The primary goal of performing reconnaissance on Mars on a sub-global scale is to know as much as possible about the environment to which crews will be subjected. At the current rate of launching and landing probes to Mars, it will take a very long time to understand the surface and atmospheric conditions associated with the regions where prospective crews may land. Meanwhile electronics and electrical systems are rapidly getting smaller. One can argue that to acquire the knowledge of the region, one must take hundreds, maybe thousands of measurements simultaneously. One means to perform such a task is to deploy a swarm of sensors. Such a swarm would perform an in-situ assessment of the region. Imagine a close flyby mission to Mars for example, where mini- to micro-sensors are deposited into the atmosphere over half an orbit or more. The sensors, captured by the atmospheric drag and Martian gravity slowly descend buffeted about by Martian winds and weather until they settle on the surface a great time later (think of how long dust takes to settle). As they descend they communicate a vast array of data; temperature, chemistry, pressure, radiation dose, electric or magnetic properties from a region of the planet and an individual sensor need not measure the same quantity as its neighbors. Initially, they could move at the whim of the environment but later versions could have locomotion or propulsion mechanisms. Humans wouldn't need to decide where the sensors go, the sensors do that for themselves. This is a key strength of a sensor swarm. The intelligence relies on the group not on a decision maker on earth. Real time sensor inputs direct what the swarm considers most interesting to investigate resulting in emergent behavior. We issued a $20,000 challenge to the global innovators to provide solutions as to how such a swarm could be initialized and by what protocols and methodologies by which they operate. Over 400 innovators from 49 countries took a look at the problem, with three receiving partial awards for solutions.

Mel Ferebee↗

Assessing within-Field Corn and Soybean Yield Variability from WorldView-3, Planet, Sentinel-2, and Landsat 8 Satellite Imagery

Crop yield monitoring is an important component in agricultural assessment. Multispectral remote sensing instruments onboard space-borne platforms such as Advanced Very High Resolution Radiometer (AVHRR), Moderate Resolution Imaging Spectroradiometer (MODIS), and Visible Infrared Imaging Radiometer Suite (VIIRS) have shown to be useful for efficiently generating timely and synoptic information on the yield status of crops across regional levels. However, the coarse spatial resolution data inherent to these sensors provides little utility at the management level. Recent satellite imagery collection advances toward finer spatial resolution (down to 1 m) alongside increased observational cadence (near daily) implies information on crops obtainable at field and within-field scales to support farming needs is now possible. To test this premise, we focus on assessing the efficiency of multiple satellite sensors, namely WorldView-3, Planet/Dove-Classic, Sentinel-2, and Landsat 8 (through Harmonized Landsat Sentinel-2 (HLS)), and investigate their spatial, spectral (surface reflectance (SR) and vegetation indices (VIs)), and temporal characteristics to estimate corn and soybean yields at sub-field scales within study sites in the US state of Iowa. Precision yield data as referenced to combine harvesters’ GPS systems were used for validation. We show that imagery spatial resolution of 3 m is critical to explaining 100% of the within-field yield variability for corn and soybean. Our simulation results show that moving to coarser resolution data of 10 m, 20 m, and 30 m reduced the explained variability to 86%, 72%, and 59%, respectively. We show that the most important spectral bands explaining yield variability were green (0.560 µm), red-edge (0.726 µm), and near-infrared (NIR - 0.865 µm). Furthermore, the high temporal frequency of Planet and a combination of Sentinel-2/Landsat 8 (HLS) data allowed for optimal date selection for yield map generation. Overall, we observed mixed performance of satellite-derived models with the coefficient of determination (R^2) varying from 0.21 to 0.88 (averaging 0.56) for the 30 m HLS and from 0.09 to 0.77 (averaging 0.30) for 3 m Planet. R^2 was lower for fields with higher yields, suggesting saturation of the satellite-collected reflectance features in those cases. Therefore, other biophysical variables, such as soil moisture and evapotranspiration, at similar fine spatial resolutions are likely needed alongside the optical imagery to fully explain the yields.

agriculture; yield; within-field; corn; soybean; r↗

Power spectrum constraints from spectral distortions in the cosmic microwave background

Using recent experimental limits on chemical potential distortions from Cosmic Background Explorer (COBE) Far Infrared Astronomy Satellite (FIRAS), and the large lever-arm spanning the damping of sub-Jeans scale fluctuations to the COBE DMR fluctuations, we set a constraint on the slope of the primordial power spectrum n. It is possible to analytically calculate the contribution over the full range of scales and redshifts, correctly taking into account fluctuation growth and damping as well as thermalization processes. Assuming conservatively that mu is less than 1.76 x 10(exp -4), we find that the 95% upper limit on n is only weakly dependent on other cosmological parameters, e.g., n is less than 1.60 (h=0.5) and n is less than 1.63 (h=1.0) for Omega(sub 0) = 1, with marginally weaker constraints for Omega(sub 0) is less than 1 in a flat model with a cosmological constant.

Hu, Wayne↗

An Evaluation of Remotely Sensed and In-Situ Data Sufficiency for SGMA-Scale Groundwater Studies in the Central Valley, California

California’s Central Valley aquifer is a critical freshwater resource for the state, providing drinking water to 6.5 million residents and irrigation water for more than half of the nation’s produce. Increasing demands have led to a steady decline of aquifer water levels, leading to irreversible compaction and land subsidence. The Sustainable Groundwater Management Act of California (SGMA), instituted in 2014, requires high-priority regions to establish sustainability plans by 2020-2022 that mitigate groundwater depletion and land subsidence. Many regions within the Central Valley lack groundwater monitoring wells or continuous well data, rendering in-situ monitoring challenging. When combined with other data, NASA’s Gravity Recovery and Climate Experiment (GRACE) and GRACE Follow On (GRACE-FO) missions provide valuable information about groundwater storage changes at a sub-basin scale. Interferometric Surface Aperture Radar (InSAR) measurements can map land subsidence, and Global Positioning System (GPS) can be used to estimate crustal uplift. Yet, all of these measurements have differences in resolution, coverage, discretization, and record length. Here we assess the various datasets that are potentially useful for Central Valley sub-basin groundwater monitoring and provide an analysis of gaps and other issues in regard to their suitability for SGMA-related analysis. Finally, we offer next steps and recommendations on data treatment and integration.

Kyra H Kim↗

Violation of Field Line Conservation and Associated Spatial Scales in Particle-In-Cell Simulations and MMS Data

For the first time, space flight technology exists to detect, in situ, violation of magnetic field line conservation. The violation of magnetic line conservation on scales smaller than the system size is a necessary and sufficient condition for finite magnetic field reconnection. We demonstrate that violation of line conservation produces a detectable, structured signature both in particle-in-cell simulations of reconnection and in data from the Magnetospheric Multi-Scale mission. In particle-in-cell simulations of asymmetric reconnection, the quantity-which we call M-that identifies this violation achieves significant values in electron skin depth-scale layers that extend away from the separator, with higher values emerging on the low-density, high-magnetic-field side of the current sheet. At the separator, M owes largely to perpendicular gradients in the parallel electric field, while it attains its highest values away from the separator in dispersed, layered structures associated with gradients in the perpendicular nonideal electric field and electron transport. Sub-ion scale bipolar forms of the quantity also appear further from the separator. In two MMS burst data intervals detecting the electron diffusion region, we find that M exceeds measurement uncertainties both at the separator and near the separatrices. One interval has highly sheared reconnecting fields and the other a stronger guide field. For one event, we determine the location and scale of M and the inner electron diffusion region relative to electron outflows and the magnetic separatrices. The measure can therefore serve as a potent diagnostic for magnetic reconnection in space measurements.

Wendel, Deirdre E.↗

Modeling Mars Cyclogenesis and Frontal Waves: Seasonal Variations and Implications on Dust Activity

Between late autumn through early spring,middle and high latitudes onMars exhibit strong equator-to-polemean temperature contrasts (i.e., "baroclinicity"). Data collected during the Viking era and observations from both the Mars Global Surveyor (MGS) and Mars Reconnaissance Orbiter (MRO) indicate that such strong baroclinicity supports vigorous, large-scale eastward traveling weather systems (i.e., transient synoptic period waves) [1, 2]. For a rapidly rotating, differentially heated, shallow atmosphere such as on Earth and Mars, these large-scale, extratropical weather disturbances are critical components of the global circulation. The wave-like disturbances serve as agents in the transport of heat and momentum between low and high latitudes of the planet. Through cyclonic/anticyclonic winds, intense shear deformations, contractions-dilatations in temperature and density, and sharp perturbations amongst atmospheric tracers (i.e., dust, volatiles (e.g., water vapor) and condensates (e.g., water-ice cloud particles)), Mars' extratropical weather systems have significant sub-synoptic scale ramifications by supporting atmospheric frontal waves (Fig. 1).

Mars Global Surveyor (MGS)↗

Dynamics and Stability of Capillary Surfaces: Liquid Switches at Small Scales

The dynamics and stability of systems of interfaces is central to a range of technologies related to the Human Exploration and Development of Space (HEDS). Our premise is that dramatic shape changes can be manipulated to advantage with minimal input, if the system is near instability. The primary objective is to develop the science base to allow novel approaches to liquid management in low-gravity based on this premise. HEDS requires efficient, reliable and lightweight technologies. Our poster will highlight our progress toward this goal using the capillary switch as an example. A capillary surface is a liquid/liquid or liquid/gas interface whose shape is determined by surface tension. For typical liquids (e.g., water) against gas on earth, capillary surfaces occur on the millimeterscale and smaller where shape deformation due to gravity is unimportant. In low gravity, they can occur on the centimeter scale. Capillary surfaces can be combined to make a switch a system with multiple stable states. A capillary switch can generate motion or effect force. To be practical, the energy barriers of such a switch must be tunable, its switching time (kinetics) short and its triggering mechanism reliable. We illustrate these features with a capillary switch that consists of two droplets, coupled by common pressure. As long as contact lines remained pinned, motions are inviscid, even at sub-millimeter scales, with consequent promise of low-power consumption at the device level. Predictions of theory are compared to experiment on i) a soap-film prototype at centimeter scale and ii) a liquid droplet switch at millimeter-scale.

Steen, Paul H.↗

Vorticity imbalance and stability in relation to convection

A complete synoptic-scale vorticity budget was related to convection storm development in the eastern two-thirds of the United States. The 3-h sounding interval permitted a study of time changes of the vorticity budget in areas of convective storms. Results of analyses revealed significant changes in values of terms in the vorticity equation at different stages of squall line development. Average budgets for all areas of convection indicate systematic imbalance in the terms in the vorticity equation. This imbalance resulted primarily from sub-grid scale processes. Potential instability in the lower troposphere was analyzed in relation to the development of convective activity. Instability was related to areas of convection; however, instability alone was inadequate for forecast purposes. Combinations of stability and terms in the vorticity equation in the form of indices succeeded in depicting areas of convection better than any one item separately.

Read, W. L.↗

Part 2 of a Computational Study of a Drop-Laden Mixing Layer

This second of three reports on a computational study of a mixing layer laden with evaporating liquid drops presents the evaluation of Large Eddy Simulation (LES) models. The LES models were evaluated on an existing database that had been generated using Direct Numerical Simulation (DNS). The DNS method and the database are described in the first report of this series, Part 1 of a Computational Study of a Drop-Laden Mixing Layer (NPO-30719), NASA Tech Briefs, Vol. 28, No.7 (July 2004), page 59. The LES equations, which are derived by applying a spatial filter to the DNS set, govern the evolution of the larger scales of the flow and can therefore be solved on a coarser grid. Consistent with the reduction in grid points, the DNS drops would be represented by fewer drops, called computational drops in the LES context. The LES equations contain terms that cannot be directly computed on the coarser grid and that must instead be modeled. Two types of models are necessary: (1) those for the filtered source terms representing the effects of drops on the filtered flow field and (2) those for the sub-grid scale (SGS) fluxes arising from filtering the convective terms in the DNS equations. All of the filtered-sourceterm models that were developed were found to overestimate the filtered source terms. For modeling the SGS fluxes, constant-coefficient Smagorinsky, gradient, and scale-similarity models were assessed and calibrated on the DNS database. The Smagorinsky model correlated poorly with the SGS fluxes, whereas the gradient and scale-similarity models were well correlated with the SGS quantities that they represented.

Okongo, Nora↗

Mapping of the Ronda peridotite massif (Spain) from AVIRIS spectro-imaging survey: A first attempt

In both AVIRIS and ISM data, through the use of mixing models, geological boundaries of the Ronda massif are identified with respect to the surrounding rocks. We can also yield first-order vegetation maps. ISM and AVIRIS instruments give consistent results. On the basis of endmember fraction images, it is then possible to discard areas highly vegetated or not belonging to the peridotite massif. Within the remaining part of the mosaic, spectro-mixing analysis reveals spectral variations in the peridotite massif between the well-exposed areas. Spatially organized units are depicted, related to differences in the relative depth of the absorption band at 1 micron, and it may be due to a different pyroxene content. At this stage, it is worth noting that, although mineralogical variations observed in the rocks are at a sub-pixel scale for the airborne analysis, we see an emerging spatial pattern in the distribution of spectral variations across the massif which might be prevailingly related to mineralogy. Although it is known from fieldwork that the Ronda peridotite massif exhibits mineralogical variations at local scale in the content of pyroxene, and at regional scale in different mineral facies, ranging from garnet-, to spinel- to plagioclase-lherzolites, no attempt has been done yet to produce a synoptic map relating the two scales of analysis. The present work is a first attempt to reach this objective, though a lot more work is still required. In particular, for the purpose of mineralogical interpretation, it is critical to relate the airborne observation to field work and laboratory spectra of Ronda rocks already obtained, with the use of image endmembers and associated reference endmembers. Also, the pretty rough linear mixing model used here is taken as a 'black-box' process which does not necessarily apply correctly to the physical situation at the sub-pixel level. One may think of using the ground-truth observations bearing on the sub-pixel statistical characteristics (texture, structural pattern, surface distribution and vegetation contribution (grass,..)) to produce a more advanced mixing model, physically appropriate to the geologic and environmental contexts.

Pinet, P. C.↗

A Priori Subgrid Scale Modeling for a Droplet Laden Temporal Mixing Layer

Subgrid analysis of a transitional temporal mixing layer with evaporating droplets has been performed using a direct numerical simulation (DNS) database. The DNS is for a Reynolds number (based on initial vorticity thickness) of 600, with droplet mass loading of 0.2. The gas phase is computed using a Eulerian formulation, with Lagrangian droplet tracking. Since Large Eddy Simulation (LES) of this flow requires the computation of unfiltered gas-phase variables at droplet locations from filtered gas-phase variables at the grid points, it is proposed to model these by assuming the gas-phase variables to be given by the filtered variables plus a correction based on the filtered standard deviation, which can be computed from the sub-grid scale (SGS) standard deviation. This model predicts unfiltered variables at droplet locations better than simply interpolating the filtered variables. Three methods are investigated for modeling the SGS standard deviation: Smagorinsky, gradient and scale-similarity. When properly calibrated, the gradient and scale-similarity methods give results in excellent agreement with the DNS.

Okongo, Nora↗

Linear theory of tearing in a high-beta plasma

The linear dispersion relation for a collisionless plasma in a sheared one-dimensional current sheet is calculated with reference to conditions in the daytime magnetopause. Calculations are extended to include plasmas with beta approximately equal to 1. It is found that the tearing mode eigenstructure and temporal growth rate are a sensitive function of the ratios l sub s/l sub G, l sub s/l sub T, and l sub s/l sub n, where l sub s is the shearing length of the magnetic field, l sub G is the gradient scale length, and l sub T is the temperature gradient scale length. In particular, if beta is approximately equal to 1, and l sub s is less than l sub G, l sub n, and l sub T, then the thickness of the layer over which particles are resonantly accelerated by the induction magnetic field is approximately rho, a thermal gyroradius. If the above conditions are not satisfied, plasma gradients may electrostatically stabilize the mode.

Quest, K. B.↗

Distributed Disdrometer and Rain Gauge Measurement Infrastructure Developed for GPM Ground Validation

Global Precipitation Mission (GPM)retrieval algorithm validation requires datasets characterizing the 4-D structure, variability, and correlation properties of hydrometeor particle size distributions (PSD) and accumulations over satellite fields of view (FOV;<10 km). Collection of this data provides a means to assess retrieval errors related to beam filling and algorithm PSD assumptions. Hence, GPM Ground Validation is developing a deployable network of precipitation gauges and disdrometers to provide fine-scale measurements of PSD and precipitation accumulation variability. These observations will be combined with dual-frequency, polarimetric, and profiling radar data in a bootstrapping fashion to extend validated PSD measurements to a large coverage domain. Accordingly, a total of 24 Parsivel disdrometers(PD), 5 3rd-generation 2D Video Disdrometers (2DVD), 70 tipping bucket rain gauges (TBRG),9 weighing gauges, 7 Hot-Plate precipitation sensors (HP), and 3 Micro Rain Radars (MRR) have been procured. In liquid precipitation the suite of TBRG, PD and 2DVD instruments will quantify a broad spectrum of rain rate and PSD variability at sub-kilometer scales. In the envisioned network configuration 5 2DVDs will act as reference points for 16 collocated PD and TBRG measurements. We find that PD measurements provide similar measures of the rain PSD as observed with collocated 2DVDs (e.g., D0, Nw) for rain rates less than 15 mm/hr. For heavier rain rates we will rely on 2DVDs for PSD information. For snowfall we will combine point-redundant observations of SWER distributed over three or more locations within a FOV. Each location will contain at least one fenced weighing gauge, one HP, two PDs, and a 2DVD. MRRs will also be located at each site to extend the measurement to the column. By collecting SWER measurements using different instrument types that employ different measurement techniques our objective is to separate measurement uncertainty from natural variability in SWER and PSD. As demonstrated using C3VP polarimetric radar, gauge, and 2DVD/PD datasets these measurements can be combined to bootstrap an area wide SWER estimate via constrained modification of density-diameter and radar reflectivity-snowfall relationships. These data will be combined with snowpack, airborne microphysics, radar, radiometer, and tropospheric sounding data to refine GPM snowfall retrievals. The gauge and disdrometer instruments are being developed to operate autonomously when necessary using solar power and wireless communications. These systems will be deployed in numerous field campaigns through 2016. Planned deployment of these systems include field campaigns in Finland (2010), Oklahoma (2011), Canada (2012) and North Carolina (2013). GPM will also deploy 20 pairs of TBRGs within a 25 km2 region along the Virginia coast under NASA NPOL radar coverage in order to quantify errors in point-area rainfall measurements.

Petersen, Walter A.↗

Modeling fluvial erosion on regional to continental scales

The fluvial system is a major concern in modeling landform evolution in response to tectonic deformation. Three stream bed types (bedrock, coarse-bed alluvial, and fine-bed alluvial) differ in factors controlling their occurrence and evolution and in appropriate modeling approaches. Spatial and temporal transitions among bed types occur in response to changes in sediment characteristics and tectonic deformation. Erosion in bedrock channels depends upon the ability to scour or pluck bed material; this detachment capacity is often a power function of drainage area and gradient. Exposure of bedrock in channel beds, due to rapid downcutting or resistant rock, slows the response of headwater catchments to downstream baselevel changes. Sediment routing through alluvial channels must account for supply from slope erosion, transport rates, abrasion, and sorting. In regional landform modeling, implicit rate laws must be developed for sediment production from erosion of sub-grid-scale slopes and small channels.

Howard, Alan D.↗

Formation of Large Regular Satellites of Giant Planets in an Extended Gaseous Nebula: Subnebula Model and Accretion of Satellites

We model the subnebulae of Jupiter and Saturn wherein satellite accretion took place. We expect a giant planet subnebula to be composed of an optically thick (given gaseous opacity) inner region inside of the planet's centrifugal radius (located at r(sub c, sup J) = l5R(sub J) for Jupiter and r(sub c, sup S) = 22R(sub S) for Saturn), and an optically thin, extended outer disk out to a fraction of the planet's Roche lobe, which we choose to be R(sub roche)/5 (located at approximately 150R(sub J) near the inner irregular satellites for Jupiter, and approximately 200R(sub S) near Phoebe for Saturn). This places Titan and Ganymede in the inner disk, Callisto and Iapetus in the outer disk, and Hyperion in the transition region. The inner disk is the leftover of the gas accreted by the protoplanet. The outer disk results from the solar torque on nebula gas flowing into the protoplanet during the time of giant planet gap opening. For the sake of specificity, we use a cosmic mixture 'minimum mass' model to constrain the gas densities of the inner disks of Jupiter and Saturn (and also Uranus). For the total mass of the outer disk we use the simple scaling M(sub disk) = M(sub P)tau(sub gap)/tau(sub acc), where M(sub P) is the mass of the giant planet, tau(sub gap) is the gap opening timescale, and tau(sub acc) is the giant planet accretion time. This gives a total outer disk mass of approximately 100M(sub Callisto) for Jupiter and possibly approximately 200M(sub Iapetus) for Saturn (which contain enough condensables to form Callisto and Iapetus respectively). Our model has Ganymede at a subnebula temperature of approximately 250 K and Titan at approximately 100 K. The outer disks of Jupiter and Saturn have constant temperatures of 130 K and 90 K respectively.

Mosqueira, I.↗