Search NASA⌕ Search

SEARCH · Search NASA

Results for “AGE ESTIMATION”

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 955 records · Page 53

M dwarfs from the Einstein extended medium sensitivity survey

A complete sample of X-ray selected M dwarfs are presented which were serendipitously detected in the extended Einstein Observatory Medium Sensitivity Survey (MSS). The MSS detected only early M dwarfs (type MO through M5), 84 percent of which were emission stars (i.e., type Me V). It was calculated that the X-ray luminosity function for early M dwarfs with a log Lx greater than 27.6 and, by comparison to optical luminosity functions, find that about 25 percent of all early M dwarfs emit X-rays at a level greater than log Lx = 27.6. This result is compared to luminosity function estimates from optically selected samples and discussed in light of the sample's kinematic and mean age. Also presented is H alpha and Ca II K line fluxes for each member of the sample and it is shown that the surface H alpha and Ca II K luminosities do indeed correlate with Lx. The significance of this effect to theories of chromospheric and coronal heating is explained.

Fleming, Thomas A.↗

Geologic Evolution of the Martian Dichotomy in the Ismenius Area of Mars and Implications for Plains Magnetization

The origin of the Martian dichotomy, which divides highlands from lowlands, is unknown. We examine a section of the dichotomy ( 50 - 90E) defined by steep scarps and normal faults. Stratigraphy and age relationships preclude the formation of the 2.5 km high boundary via erosion. The abrupt disappearance of topographic knobs similar to 300 - 500 km to the northeast is interpreted as a buried fault. Alignment of the buried fault with grabens, stratigraphy, and age determinations using crater counts indicate that the lowland bench is down faulted highlands crust. The estimated local strain (3.5%) and fault pattern are broadly consistent with gravitational relaxation of a plateau boundary. Magnetic and gravity anomalies occur on either side of the buried fault. Admittance analysis indicates isostatic compensation. Although nonunique, a model with a 10 km thick intracrustal block under the lowland bench, a 20 km thick block under the plains, and an excess density of 200 kg/m(3) provides a good fit to the isostatic anomaly. A good fit to a profile of the magnetic field perpendicular to the dichotomy is produced using uniformly polarized intracrustal blocks 10 - 20 km thick, an intensity of 6 Am/m, a field inclination of -30 degrees, and gaps aligned with the isostatic anomalies. One interpretation is that high-density intrusions demagnetized the crust after dynamo cessation and that low-lying magnetized areas could be down faulted highlands crust. Another model (inclination of 30 degrees) has magnetized crust beneath the isostatic anomalies, separated by gaps. The gaps could result from hydrothermal alteration of the crust along fault zones.

faulting↗

Deriving brown carbon from multiwavelength absorption measurements: Method and application to AERONET and Aethalometer observations

The radiative impact of organic aerosols (OA) is a large source of uncertainty in estimating the global direct radiative effect (DRE) of aerosols. This radiative impact includes not only light scattering but also light absorption from a subclass of OA referred to as brown carbon (BrC). However the absorption properties of BrC are poorly understood leading to large uncertainties in modelling studies. To obtain observational constraints from measurements, a simple Absorption Ångström Exponent (AAE) method is often used to separate the contribution of BrC absorption from that of black carbon (BC). However, this attribution method is based on assumptions regarding the spectral dependence of BC that are often violated in the ambient atmosphere. Here we develop a new method that decreases the uncertainties associated with estimating BrC absorption. By applying this method to multi-wavelength absorption aerosol optical depth (AAOD) measurements at AERONET sites worldwide and surface aerosol absorption measurements at multiple ambient sites, we estimate that BrC globally contributes 6-40% of the absorption at 440nm. We find that the mass absorption coefficient of OA (OA-MAC) is positively correlated with BC/OA mass ratio. Based on the variability of BC properties and BC/OA emission ratio, we estimate a range of 0.05-1.2 m 2 /g for OA-MAC at 440nm. Using the combination of AERONET and OMI UV absorption observations we estimate that the AAE 388/440nm for BrC is generally ~4 world-wide, with a smaller value in Europe (< 2). Our analyses of two surface sites (Cape Cod, to the southeast of Boston, and the GoAmazon2014/5 T3 site, to the west of Manaus, Brazil) reveal no significant relationship between BrC absorptivity and photochemical aging in typical urban influenced conditions. However, the absorption of BrC measured during the biomass burning season near Manaus is found to decrease with photochemical aging with a lifetime of ~1 day. This lifetime is comparable to previous observations within a biomass burning plume but much slower than estimated from laboratory studies.

54 ENVIRONMENTAL SCIENCES↗

Molecular transformations of phenolic SOA during photochemical aging in the aqueous phase: competition among oligomerization, functionalization, and fragmentation

Abstract. Organic aerosol is formed and transformed in atmospheric aqueous phases (e.g., cloud and fog droplets and deliquesced airborne particles containing small amounts of water) through a multitude of chemical reactions. Understanding these reactions is important for a predictive understanding of atmospheric aging of aerosols and their impacts on climate, air quality, and human health. In this study, we investigate the chemical evolution of aqueous secondary organic aerosol (aqSOA) formed during reactions of phenolic compounds with two oxidants – the triplet excited state of an aromatic carbonyl (3C∗) and hydroxyl radical (•OH). Changes in the molecular composition of aqSOA as a function of aging time are characterized using an offline nanospray desorption electrospray ionization mass spectrometer (nano-DESI MS) whereas the real-time evolution of SOA mass, elemental ratios, and average carbon oxidation state (OSC) are monitored using an online aerosol mass spectrometer (AMS). Our results indicate that oligomerization is an important aqueous reaction pathway for phenols, especially during the initial stage of photooxidation equivalent to ∼ 2 h irradiation under midday winter solstice sunlight in Northern California. At later reaction times functionalization (i.e., adding polar oxygenated functional groups to the molecule) and fragmentation (i.e., breaking of covalent bonds) become more important processes, forming a large variety of functionalized aromatic and open-ring products with higher OSC values. Fragmentation reactions eventually dominate the photochemical evolution of phenolic aqSOA, forming a large number of highly oxygenated ring-opening molecules with carbon numbers (nC) below 6. The average nC of phenolic aqSOA decreases while average OSC increases over the course of photochemical aging. In addition, the saturation vapor pressures (C∗) of dozens of the most abundant phenolic aqSOA molecules are estimated. A wide range of C∗ values is observed, varying from < 10−20 µg m−3 for functionalized phenolic oligomers to > 10 µg m−3 for small open-ring species. The detection of abundant extremely low-volatile organic compounds (ELVOC) indicates that aqueous reactions of phenolic compounds are likely an important source of ELVOC in the atmosphere.

54 ENVIRONMENTAL SCIENCES↗

The Young L Dwarf 2MASS J11193254−1137466 Is a Planetary-mass Binary

We have discovered that the extremely red, low-gravity L7 dwarf 2MASS J11193254−1137466 is a 0.″14 (3.6 au) binary using Keck laser guide star adaptive optics imaging. 2MASS J11193254−1137466 has previously been identified as a likely member of the TW Hydrae Association (TWA). Using our updated photometric distance and proper motion, a kinematic analysis based on the BANYAN II model gives an 82% probability of TWA membership. At TWA’s 10 ± 3 Myr age and using hot-start evolutionary models, 2MASS J11193254−1137466AB is a pair of 3.7{sub −0.9}{sup +1.2} M{sub Jup} brown dwarfs, making it the lowest-mass binary discovered to date. We estimate an orbital period of 90{sub −50}{sup +80} years. One component is marginally brighter in K band but fainter in J band, making this a probable flux-reversal binary, the first discovered with such a young age. We also imaged the spectrally similar TWA L7 dwarf WISEA J114724.10−204021.3 with Keck and found no sign of binarity. Our evolutionary model-derived T{sub eff} estimate for WISEA J114724.10−204021.3 is ≈230 K higher than for 2MASS J11193254−1137466AB, at odds with the spectral similarity of the two objects. This discrepancy suggests that WISEA J114724.10−204021.3 may actually be a tight binary with masses and temperatures very similar to 2MASS J11193254−1137466AB, or further supporting the idea that near-infrared spectra of young ultracool dwarfs are shaped by factors other than temperature and gravity. 2MASS J11193254−1137466AB will be an essential benchmark for testing evolutionary and atmospheric models in the young planetary-mass regime.

79 ASTRONOMY AND ASTROPHYSICS↗

The trends high-contrast imaging survey. VI. Discovery of a mass, age, and metallicity benchmark brown dwarf

The mass and age of substellar objects are degenerate parameters leaving the evolutionary state of brown dwarfs ambiguous without additional information. Theoretical models are normally used to help distinguish between old, massive brown dwarfs and young, low-mass brown dwarfs but these models have yet to be properly calibrated. We have carried out an infrared high-contrast imaging program with the goal of detecting substellar objects as companions to nearby stars to help break degeneracies in inferred physical properties such as mass, age, and composition. Rather than using imaging observations alone, our targets are pre-selected based on the existence of dynamical accelerations informed from years of stellar radial velocity (RV) measurements. In this paper, we present the discovery of a rare benchmark brown dwarf orbiting the nearby (d = 18.69 ± 0.19 pc), solar-type (G9V) star HD 4747 ([Fe/H] = −0.22 ± 0.04) with a projected separation of only ρ = 11.3 ± 0.2 au (θ = 0.″6). Precise Doppler measurements taken over 18 years reveal the companion’s orbit and allow us to place strong constraints on its mass using dynamics (msin i=55.3±1.9M{sub Jup}). Relative photometry (ΔK {sub s} = 9.05 ± 0.14, M{sub K{sub s}}=13.00±0.14, K{sub s}−L{sup ′}=1.34±0.46) indicates that HD 4747 B is most likely a late-type L-dwarf and, if near the L/T transition, an intriguing source for studying cloud physics, variability, and polarization. We estimate a model-dependent mass of m=72{sub −13}{sup +3} M{sub Jup} for an age of 3.3{sub −1.9}{sup +2.3} Gyr based on gyrochronology. Combining astrometric measurements with RV data, we calculate the companion dynamical mass (m=60.2±3.3M{sub Jup}) and orbit (e = 0.740 ± 0.002) directly. As a new mass, age, and metallicity benchmark, HD 4747 B will serve as a laboratory for precision astrophysics to test theoretical models that describe the emergent radiation of brown dwarfs.

79 ASTRONOMY AND ASTROPHYSICS↗

Characterization and Valuation of the Uncertainty of Calibrated Parameters in Microsimulation Decision Models

We evaluated the implications of different approaches to characterize the uncertainty of calibrated parameters of microsimulation decision models (DMs) and quantified the value of such uncertainty in decision making. We calibrated the natural history model of CRC to simulated epidemiological data with different degrees of uncertainty and obtained the joint posterior distribution of the parameters using a Bayesian approach. We conducted a probabilistic sensitivity analysis (PSA) on all the model parameters with different characterizations of the uncertainty of the calibrated parameters. We estimated the value of uncertainty of the various characterizations with a value of information analysis. We conducted all analyses using high-performance computing resources running the Extreme-scale Model Exploration with Swift (EMEWS) framework. The posterior distribution had a high correlation among some parameters. The parameters of the Weibull hazard function for the age of onset of adenomas had the highest posterior correlation of -0.958. When comparing full posterior distributions and the maximum-a-posteriori estimate of the calibrated parameters, there is little difference in the spread of the distribution of the CEA outcomes with a similar expected value of perfect information (EVPI) of $\$$653 and $\$$685, respectively, at a willingness-to-pay (WTP) threshold of $\$$66,000 per quality-adjusted life year (QALY). Ignoring correlation on the calibrated parameters’ posterior distribution produced the broadest distribution of CEA outcomes and the highest EVPI of $\$$809 at the same WTP threshold. Different characterizations of the uncertainty of calibrated parameters affect the expected value of eliminating parametric uncertainty on the CEA. Ignoring inherent correlation among calibrated parameters on a PSA overestimates the value of uncertainty.

97 MATHEMATICS AND COMPUTING↗

An Expanded Chemo-dynamical Sample of Red Giants in the Bar of the Large Magellanic Cloud

We report new spectroscopic observations obtained with the Michigan/Magellan Fiber System of 308 red giants (RGs) located in two fields near the photometric center of the bar of the Large Magellanic Cloud. This sample consists of 131 stars observed in previous studies (in one field) and 177 newly observed stars (in the second field) selected specifically to more reliably establish the metallicity and age distributions of the bar. For each star, we measure its heliocentric line-of-sight velocity, surface gravity, and metallicity from its high-resolution spectrum (effective temperatures come from photometric colors). The spectroscopic Hertzsprung–Russell diagrams—modulo small offsets in surface gravities—reveal good agreement with model isochrones. The mean metallicity of the 177-RG sample is [Fe/H] = −0.76 ± 0.02 with a metallicity dispersion σ = 0.28 ± 0.03. The corresponding metallicity distribution—corrected for selection effects—is well fitted by two Gaussian components: one metal-rich with a mean −0.66 ± 0.02 and a standard deviation 0.17 ± 0.01, and the other metal-poor with −1.20 ± 0.24 and 0.41 ± 0.06. The metal-rich and metal-poor populations contain approximately 85% and 15% of stars, respectively. We also confirm that the velocity dispersion in the bar center decreases significantly from 31.2 ± 4.3 to 18.7 ± 1.9 km s{sup −1} with increasing metallicity over the range −2.09 to −0.38. Individual stellar masses are estimated using the spectroscopic surface gravities and the known luminosities. We find that lower mass, hence older, RGs have larger metallicity dispersion and lower mean metallicity than the higher-mass, younger RGs. The estimated masses, however, extend to implausibly low values (∼0.1 M{sub ⊙}), making it impossible to obtain an absolute age–metallicity or age distribution of the bar.

79 ASTRONOMY AND ASTROPHYSICS↗

Low dose radiation risks for women surviving the a-bombs in Japan: generalized additive model

Background: Analyses of cancer mortality and incidence in Japanese A-bomb survivors have been used to estimate radiation risks, which are generally higher for women. Relative Risk (RR) is usually modelled as a linear function of dose. Extrapolation from data including high doses predicts small risks at low doses. Generalized Additive Models (GAMs) are flexible methods for modelling non-linear behaviour. Methods: GAMs are applied to cancer incidence in female low dose subcohorts, using anonymous public data for the 1958 – 1998 Life Span Study, to test for linearity, explore interactions, adjust for the skewed dose distribution, examine significance below 100 mGy, and estimate risks at 10 mGy. Results: For all solid cancer incidence, RR estimated from 0 – 100 mGy and 0 – 20 mGy subcohorts is significantly raised. The response tapers above 150 mGy. At low doses, RR increases with age-at-exposure and decreases with time-since-exposure, the preferred covariate. Using the empirical cumulative distribution of dose improves model fit, and capacity to detect non-linear responses. RR is elevated over wide ranges of covariate values. Results are stable under simulation, or when removing exceptional data cells, or adjusting neutron RBE. Estimates of Excess RR at 10 mGy using the cumulative dose distribution are 10 – 45 times higher than extrapolations from a linear model fitted to the full cohort. Below 100 mGy, quasipoisson models find significant effects for all solid, squamous, uterus, corpus, and thyroid cancers, and for respiratory cancers when age-at-exposure > 35 yrs. Results for the thyroid are compatible with studies of children treated for tinea capitis, and Chernobyl survivors. Results for the uterus are compatible with studies of UK nuclear workers and the Techa River cohort. Conclusion: Non-linear models find large, significant cancer risks for Japanese women exposed to low dose radiation from the atomic bombings. The risks should be reflected in protection standards.

60 APPLIED LIFE SCIENCES↗

Carbon fluxes and interannual drivers in a temperate forest ecosystem assessed through comparison of top-down and bottom-up approaches

Despite decades of research, gaining a comprehensive understanding of carbon (C) cycling in forests remains a considerable challenge. Uncertainties stem from persistent methodological limitations and the difficulty of resolving top-down estimates of ecosystem C exchange with bottom-up measurements of individual pools and fluxes. To address this, we derived estimates and associated uncertainties of ecosystem C fluxes for a 100–125 year old mixed temperate forest stand at the Bartlett Experimental Forest, New Hampshire, USA, using three different approaches: (1) tower-based eddy covariance, (2) a biometric approach involving C flux measurements of individual ecosystem subcomponents, and (3) an inventory approach involving changes in major C stocks over time. Our analysis made use of 13 years of data, collected over the period from 2004 to 2016. Estimates of mean annual net ecosystem production (NEP) ranged from 120 to 133 g C m -2 , demonstrating strong agreement among methods and suggesting that this aging forest acts as a moderate C sink. The use of multiple approaches to measure C fluxes and their uncertainties helped place constraints on difficult-to-measure processes such as aboveground contributions to ecosystem respiration and belowground allocation to mycorrhizal fungal biomass (which was estimated at 20% of net primary production). Analysis of interannual variability in C fluxes revealed a decoupling between annual wood growth and either current year or lagged NEP or GPP, suggesting that source limitation (C supply) is likely not controlling rates of wood production, at least on an interannual scale. Results also demonstrated a strong association between the maximum rate of C uptake during the growing season (Amax) and the length of the vernal window, defined as the period of time between soil thaw and the onset of photosynthesis. This suggests an important, but poorly understood, influence of winter and spring climate on mid-summer canopy physiology. Efforts to resolve the mechanisms responsible should be prioritized in light of ongoing and predicted changes in climate for the northeastern U.S. region, particularly during the winter and winter-spring transition period.

54 ENVIRONMENTAL SCIENCES↗

Improved Survival With Radiation Therapy in Stage I-II Primary Mediastinal B Cell Lymphoma: A Surveillance, Epidemiology, and End Results Database Analysis

Background: Primary mediastinal B cell lymphoma (PMBCL) is an uncommon lymphoma for which trials are few with small patient numbers. The role of radiation therapy (RT) after standard immunochemotherapy for early-stage disease has never been studied prospectively. We used the Surveillance, Epidemiology, and End Results (SEER) database to evaluate PMBCL and the impact of RT on outcomes. Methods and Materials: We queried the SEER database for patients with stage I-II PMBCL diagnosed from 2001 to 2011. Retrievable data included age, gender, race (white/nonwhite), stage, extranodal disease, year of diagnosis, and use of RT as a component of definitive therapy. Kaplan-Meier overall survival (OS) estimates, univariate (UVA) log-rank and multivariate (MVA) Cox proportional hazards regression analyses were performed. Results: Two hundred fifty patients with stage I-II disease were identified, with a median follow-up time of 39 months (range, 3-125 months). The median age was 36 years (range, 18-89 years); 61% were female; 76% were white; 45% had stage I disease, 60% had extranodal disease, and 55% were given RT. The 5-year OS for the entire cohort was 86%. On UVA, OS was improved with RT (hazard ratio [HR] 0.446, P=.029) and decreased in association with nonwhite race (HR 2.70, P=.006). The 5-year OS was 79% (no RT) and 90% (RT). On MVA, white race and RT remained significantly associated with improved OS (P=.007 and .018, respectively). The use of RT decreased over time: 61% for the 67 patients whose disease was diagnosed from 2001 to 2005 and 53% in the 138 patients treated from 2006 to 2010. Conclusion: This retrospective population-based analysis is the largest PMBCL dataset to date and demonstrates a significant survival benefit associated with RT. Nearly half of patients treated in the United States do not receive RT, and its use appears to be declining. In the absence of phase 3 data, the use of RT should be strongly considered for its survival benefit in early-stage disease.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

CALPHAD-based ICME design of single-step aging to enhance mechanical strength of WAAM Haynes 282

To match the strength of wire-arc additive manufactured Haynes 282 to its wrought counterpart via a single-step aging heat treatment, the CALPHAD (Calculation of Phase Diagrams) method is integrated with physics-based process-structure-property models and experimental validation. The integrated computational materials engineering (ICME) framework simulates the effects of aging on γ′ and M 23 C 6 precipitation and the resulting yield strength. To improve simulation reliability, the interfacial energies between γ/γ′ and γ/M 23 C 6 carbides were estimated by comparison with precipitation kinetic modeling and measured precipitate sizes. γ′ and M23C6 were found to precipitate simultaneously between 640 and 860 °C, producing microstructures similar to those produced by two-step aging. The optimal γ′ size for peak yield stress was calculated to be 20–23 nm. WAAM Haynes 282 aged at 780 °C for 50 h exceeded the mechanical performance of its wrought counterpart subjected to two-step aging, though desired properties can also be achieved at 800 °C for 16 h or less. The error in yield strength is less than 20 MPa, demonstrating good agreement between the modeling framework and experiments. Creep studies showed that WAAM Haynes 282 exceeded the calculated rupture time, reaching 481 h. This proposed methodology can accelerate the design of aging heat treatments for any γ′-strengthened nickel-base alloy, minimizing the resources required for trial-and-error experiments.

CALPHAD↗

Status Report - Cane Fiberboard Properties And Degradation Rates For Storage Of The 9975 Shipping Package In KAC

Thermal, mechanical and physical properties have been measured on cane fiberboard samples following accelerated aging for up to approximately 10 years. The aging environments have included elevated temperature < 250 ºF (the maximum allowed service temperature for fiberboard in 9975 packages) and elevated humidity. The results from this testing have been analyzed, and aging models fit to the data. Correlations relating several properties (thermal conductivity, energy absorption, weight, dimensions and density) to their rate of change in potential storage environments have been developed. Combined with an estimate of the actual conditions the fiberboard experiences in KAC, these models allow development of service life predictions.

42 ENGINEERING↗

Accelerated Thermal Aging of Harvested Hypalon Jacket for Remaining Useful Life Determination and Diagnosis

For nuclear power plants (NPPs) considering second license renewal for their operation beyond 60 years, knowledge of long-term operation, condition monitoring, and viability for the reactor components including reactor pressure vessel, concrete structures, and cable systems is essential. Such knowledge will provide NPPs with a game plan to predict performance and to estimate costs that are associated with monitoring or replacement programs for the affected systems. For cable systems that encompass a wide variety of materials, manufacturers, and in-plant location, accelerated aging of harvested cable jacket and insulation can provide insight as to remaining useful life and methods for monitoring. Accelerated thermal aging in air at temperatures between 80°C and 120°C was carried out on a multi-conductor control rod cable that had been in-service for over 30 years and was made by Boston Insulated Wire with Hypalon™ cable jacket and ethylene-propylene rubber insulation. From elongation at break (EAB) measurements and supporting Arrhenius analysis of the jacket material, an activation energy of 97.84 kJ/mol was estimated and the time to degradation, as represented by 50% EAB at the expected maximum operating temperature of 45°C, was estimated to be 80 years. These values were slightly below previous measurements on similar BIW Hypalon cable jacket and could either be attributed to in-service degradation or variations in material properties from process variations. In addition, results from indenter modulus measurements and Fourier transform infrared spectroscopy suggest possible markers that could be beneficial to monitor cable conditions.

cable aging, cable insulation, Hypalon, second lic↗

Estimates of total organic and inorganic chlorine in the lower stratosphere from in situ and flask measurements during AASE 2

Aircraft sampling has provided extensive in situ and flask measurements of organic chlorine species in the lower stratosphere. The recent Airborne Arctic Stratospheric Expedition 2 (AASE 2) included two independent measurements of organic chlorine species using whole air sample and real-time techniques. From the whole air sample measurements we derive directly the burden of total organic chlorine (CCl(y)) in the lower stratosphere. From the more limited real-time measurements we estimate the CCl(y) burden using mixing ratios and growth rates of the principal CCl(y) species in the troposphere in conjunction with results from a two-dimensional photochemical model. Since stratospheric chlorine is tropospheric in origin and tropospheric mixing ratios are increasing, it is necessary to establish the average age of a stratospheric air parcel to assess its total chlorine (Cl(sub Total)) abundance. Total inorganic chlorine (Cl(y)) in the parcel is then estimated by the simple difference, Cl(y) = Cl(sub Total) - CCl(y). The consistency of the results from these two quite different techniques suggests that we can determine the CCl(y) and Cl(y) in the lower stratosphere with confidence. Such estimates of organic and inorganic chlorine are crucial in evaluating the photochemistry controlling chlorine partitioning and hence ozone loss processes in the lower stratosphere.

Woodbridge, E. L.↗

Volcanic contribution to crustal growth in the Central Andes: A new estimate and a discussion of uncertainties

Volcanism above subduction zones is a major mechanism for crustal growth and compared to some other proposed processes calculation of growth rates is relatively easy given accurate volumes and ages of volcanic material. Francis and Rundie (1976) first used this approach in a small region of the Central Andes, and extrapolated their result to the entire Central Andean arc. The derived rate of 3 to 4.2 x 10 to the minus 6th cubic km yr/km of arc length is compared with an independent estimate based upon a reconnaissance census of all major volcanoes in the Central Andes.

Wood, C. A.↗