Search NASA⌕ Search

SEARCH · Search NASA

Results for “data and proxy modeling”

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 433 records · Page 24

Dst Index in the 2008 GEM Modeling Challenge - Model Performance for Moderate and Strong Magnetic Storms

The GEM 2008 modeling challenge efforts are expanding beyond comparing in-situ measurements in the magnetosphere and ionosphere to include the computation of indices to be compared. The Dst index measures the largest deviations of the horizontal magnetic field at 4 equatorial magnetometers from the quiet-time background field and is commonly used to track the strength of the magnetic disturbance of the magnetosphere during storms. Models can calculate a proxy Dst index in various ways, including using the Dessler-Parker Sckopke relation and the energy of the ring current and Biot-Savart integration of electric currents in the magnetosphere. The GEM modeling challenge investigates 4 space weather events and we compare models available at CCMC against each other and the observed values of Ost. Models used include SWMF/BATSRUS, OpenGGCM, LFM, GUMICS (3D magnetosphere MHD models), Fok-RC, CRCM, RAM-SCB (kinetic drift models of the ring current), WINDMI (magnetosphere-ionosphere electric circuit model), and predictions based on an impulse response function (IRF) model and analytic coupling functions with inputs of solar wind data. In addition to the analysis of model-observation comparisons we look at the way Dst is computed in global magnetosphere models. The default value of Dst computed by the SWMF model is for Bz the Earth's center. In addition to this, we present results obtained at different locations on the Earth's surface. We choose equatorial locations at local noon, dusk (18:00 hours), midnight and dawn (6:00 hours). The different virtual observatory locations reveal the variation around the earth-centered Dst value resulting from the distribution of electric currents in the magnetosphere during different phases of a storm.

Rastaetter, Lutz↗

Quantitative Characterization of Hyper-Local Atmospheric Greenhouse Gas Sources

Atmospheric greenhouse gas (GHG) emissions are often characterized using stationary, tower-based sensors. Ground based sensors reside in the turbulent boundary layer and are subject to intense concentration impulses from hyper-local (<100m) point sources of emissions. These high frequency spikes are often filtered out in broader emission flux studies, losing valuable information about how hyper-local sources influence receptors. In this study, we investigated how empirical atmospheric data can be used to locate and quantify a concurrently measured hyper-local point source in a dense urban setting. An eddy covariance style tower and a low-cost sensor tower were deployed in various locations around an urban, hyper-local CO2/CH4 emissions source (a continuously measured restaurant exhaust vent). A model using different processing and statistical techniques was built to examine the most effective procedures for source isolation, directional location, and emission quantification. Using excess concentrations above a minimum baseline, we identify the source using bivariate polar plots and quantify the relationship between source size, receptor distance, and statistical proxies. Furthermore, we find that varying statistical thresholds allows for identification of less influential sources which are drowned out by larger or closer sources. Finally, we show that large sources can be effectively characterized using low-cost sensors, a valuable outcome informing how networks for monitoring larger areas could be implemented. This work may provide a basis for source identification and monitoring protocols for networks that feature sensors influenced by hyper-local point sources, subject to site-specific assumptions.

54 ENVIRONMENTAL SCIENCES↗

The CHUVA Lightning Mapping Campaign

The primary science objective for the CHUVA lightning mapping campaign is to combine measurements of total lightning activity, lightning channel mapping, and detailed information on the locations of cloud charge regions of thunderstorms with the planned observations of the CHUVA (Cloud processes of tHe main precipitation systems in Brazil: A contribUtion to cloud resolVing modeling and to the GPM (GlobAl Precipitation Measurement) field campaign. The lightning campaign takes place during the CHUVA intensive observation period October-December 2011 in the vicinity of S o Luiz do Paraitinga with Brazilian, US, and European government, university and industry participants. Total lightning measurements that can be provided by ground-based regional 2-D and 3-D total lightning mapping networks coincident with overpasses of the Tropical Rainfall Measuring Mission Lightning Imaging Sensor (LIS) and the SEVIRI (Spinning Enhanced Visible and Infrared Imager) on the Meteosat Second Generation satellite in geostationary earth orbit will be used to generate proxy data sets for the next generation US and European geostationary satellites. Proxy data, which play an important role in the pre-launch mission development and in user readiness preparation, are used to develop and validate algorithms so that they will be ready for operational use quickly following the planned launch of the GOES-R Geostationary Lightning Mapper (GLM) in 2015 and the Meteosat Third Generation Lightning Imager (LI) in 2017. To date there is no well-characterized total lightning data set coincident with the imagers. Therefore, to take the greatest advantage of this opportunity to collect detailed and comprehensive total lightning data sets, test and validate multi-sensor nowcasting applications for the monitoring, tracking, warning, and prediction of severe and high impact weather, and to advance our knowledge of thunderstorm physics, extensive measurements from lightning mapping networks will be collected in conjunction with electric field mills, field change sensors, high speed cameras and other lightning sensors, dual-polarimetric radars, and aircraft in-situ microphysics which will allow for excellent cross-network inter-comparisons, assessments, and physical understanding.

Goodman, Steven J.↗

Bryce Canyon Water Resources: Monitoring Vegetation Health and Water Availability in Bryce Canyon National Park for Drought Stress Mitigation Planning

Bryce Canyon National Park is home to groundwater-dependent ecosystems (GDEs) that are threatened by a multidecadal drought and increased groundwater extraction due to a spike in tourism. These ecosystems contain unique species that are only found in areas where near-surface groundwater is present, such as aspen groves and fens. These species contribute to the high biodiversity found in Bryce Canyon, which boosts an ecosystem’s productivity and the services it provides to the park. Unfortunately, many of these GDEs are too small to identify with traditional Earth observation platforms and are difficult to physically reach for monitoring purposes. This project partnered with the National Park Service to identify springs and seeps as a proxy for GDEs within Bryce Canyon from 2013–2022. Furthermore, this project tested the feasibility of various methods to detect and monitor springs and seeps and therefore facilitate the partner’s efforts to conserve these ecologically valuable GDEs in Bryce Canyon. The team mapped groundwater discharge with high resolution National Agriculture Imagery Program (NAIP) and assessed park vegetation trends with Landsat 8 Operational Land Imager (OLI) and PlanetScope imagery. In-situ precipitation data and the Western Land Data Assimilation System (WLDAS) were used to produce time series of climatic variables. Seeps and spring locations were predicted using random forest classification and maximum entropy machine learning models.

Groundwater dependent ecosystems↗

Aperture size distribution, length, and preferential location of bed-parallel veins in shale

Bed-parallel, calcite-filled veins (BPVs) are common in shale formations, and although they have been widely described in other studies, little is known about their population aperture size distribution. To address this knowledge gap, we analyzed BPV sizes in outcrops and cores from the Vaca Muerta Formation, Neuquén Basin, Argentina; in two cores from the Marcellus Formation, Appalachian Basin, northeast Pennsylvania; and in one core from the Wolfcamp Shale, Delaware Basin, West Texas. Nine out of ten aperture size populations follow a negative exponential distribution, with one following a weak power law. Bed-parallel vein size distribution and intensity vary among formations and within the same shale. We define three groups of distributions: (1) Vaca Muerta outcrops, with the highest BPV intensity and the largest BPVs (cumulative frequency of 4.9 BPVs per meter [BPVs/m] for apertures 0.265 mm to 8.7 cm); (2) Vaca Muerta cores with a similar BPV intensity overall but with no apertures wider than 1.2 cm; and (3) Vaca Muerta, Wolfcamp, and Marcellus cores with the fewest BPVs (cumulative frequency up to 0.63 BPVs/m) and very few wider than 1 cm. Aperture and length in two outcrop data sets are weakly positively correlated and follow power laws with exponents of 0.44 and 0.49. Mechanical interfaces at boundaries between different lithologies exert a strong control on BPV location, with 65–75% of observed interfaces having BPVs along them. Only 25–30% of the BPVs occur at observed material interfaces, however, and unless subtle, unobserved mechanical layering is present, other factors must also control location. BPV intensity and organic richness (TOC) from Vaca Muerta well logs are correlated in some instances but not in others, indicating TOC is not always a good proxy for BPV location or intensity. Furthermore, these findings provide useful information for modeling of hydraulic fracture treatments where BPVs may influence development of the stimulated fracture network, for example by limiting height growth.

58 GEOSCIENCES↗

ECP Report: Update on Proxy Applications and Vendor Interactions

The ExaLearn miniGAN team (Ellis and Rajamanickam) have released miniGAN, a generative adversarial network(GAN) proxy application, through the ECP proxy application suite. miniGAN is the first machine learning proxy application in the suite (note: the ECP CANDLE project did previously release some benchmarks) and models the performance for training generator and discriminator networks. The GAN's generator and discriminator generate plausible 2D/3D maps and identify fake maps, respectively. miniGAN aims to be a proxy application for related applications in cosmology (CosmoFlow, ExaGAN) and wind energy (ExaWind). miniGAN has been developed so that optimized mathematical kernels (e.g., kernels provided by Kokkos Kernels) can be plugged into to the proxy application to explore potential performance improvements. miniGAN has been released as open source software and is available through the ECP proxy application website (https://proxyapps.exascaleproject.ordecp-proxy-appssuite/) and on GitHub (https://github.com/SandiaMLMiniApps/miniGAN). As part of this release, a generator is provided to generate a data set (series of images) that are inputs to the proxy application.

97 MATHEMATICS AND COMPUTING↗

High-Dimensional Intrinsic Interpolation Using Gaussian Process Regression and Diffusion Maps

This article considers the challenging task of estimating geologic properties of interest using a suite of proxy measurements. The current work recast this task as a manifold learning problem. In this process, this article introduces a novel regression procedure for intrinsic variables constrained onto a manifold embedded in an ambient space. The procedure is meant to sharpen high-dimensional interpolation by inferring non-linear correlations from the data being interpolated. The proposed approach augments manifold learning procedures with a Gaussian process regression. It first identifies, using diffusion maps, a low-dimensional manifold embedded in an ambient high-dimensional space associated with the data. It relies on the diffusion distance associated with this construction to define a distance function with which the data model is equipped. This distance metric function is then used to compute the correlation structure of a Gaussian process that describes the statistical dependence of quantities of interest in the high-dimensional ambient space. The proposed method is applicable to arbitrarily high-dimensional data sets. Here, it is applied to subsurface characterization using a suite of well log measurements. The predictions obtained in original, principal component, and diffusion space are compared using both qualitative and quantitative metrics. Considerable improvement in the prediction of the geological structural properties is observed with the proposed method.

58 GEOSCIENCES↗

Mesoscale modeling of hypervelocity impacts using the CTH shock physics code

Material fragmentation after a hypervelocity impact is vital to predictive electro-optical and infrared (EO/IR) modeling. Successful comparisons with data require that hot, submicron fragments are generated in such impacts; yet, experimental data has so far been unable to produce fragments of this scale. The purpose of this work was to investigate how modeling assumptions of macro-scale, bulk materials might influence the generation of debris in hypervelocity impacts and ultimately the predicted EO/IR signatures of these debris clouds. Sphere-on-plate impact simulations simplified the comparison of different modeling approaches. In one set of simulations, materials were modeled with the traditional, bulk approach. Those findings were compared to simulations run with the mesoscale material grain structure explicitly modeled. This study focused on the comparison of two parameters that are tied directly to the EO/IR signature: strain rate at failure (a proxy for debris fragment size) and material temperature. Interfaces between grains, here due to void insertion, resulted in the most notable change in both the strain rate at failure and material temperature. Shock reflections from grain-void interfaces induced higher strain rates and material temperatures, and it is expected that similar effects may be produced from inclusions or dislocations in real materials. Thus, interfaces within a material may play an important role in producing smaller hot debris frag

42 ENGINEERING↗

Real-World Driving Features for Identifying Intelligent Driver Model Parameters

Driver behavior models play a significant role in representing different driving styles and the associated relationships with traffic patterns and vehicle energy consumption in simulation studies. The models often serve as a proxy for baseline human driving when assessing energy-saving strategies that alter vehicle velocity. Such models are especially important in connectivity-enabled energy-saving strategy research because they can easily adapt to changing driving conditions like posted speed limits or change in traffic light state. While numerous driver models exist, parametric driver models provide the flexibility required to represent variability in real-world driving through different combinations of model parameters. These model parameters must be informed by a representative set of parameter values for the driver model to adequately represent a real-world driver. It stands to reason that determining the parameter values from real-world driving data would serve the purpose of representing a real-world driver. Although the published literature is replete with techniques and consequences of using real-world driving data to determine parameter values for parametric driver models, none have explored them in the context of using shorter driving features where the parameter values may change over the course of a single trip for the same driver. In this study we consider the “intelligent driver model” (IDM) as our driver behavior model and use real-world driving data from the Transportation Secure Data Center (TSDC) maintained at the National Renewable Energy Laboratory (NREL). The TSDC includes real-world travel data from across the United States, from which NREL has created a wide range of driving routes consisting of road features such as speed limits, stop locations, and turn locations. The real-world driving data are categorized into different driving regimes and extracted into driving features. The driving features are then used to calibrate the parameter values for the IDM. The distribution of the parameters and the relationships among them are reported. The insights obtained from this study enable judicious usage of IDM or similar parametric driver models to represent baseline human driver behavior in simulations.

27 ARPA - Advanced Research Projects Agency-Energy↗

Real-World Driving Features for Identifying Intelligent Driver Model Parameters: Preprint

Driver behavior models play a significant role in representing different driving styles and the associated relationships with traffic patterns and vehicle energy consumption in simulation studies. The models often serve as a proxy for baseline human driving when assessing energy-saving strategies that alter vehicle velocity. Such models are especially important in connectivity-enabled energy-saving strategy research because they can easily adapt to changing driving conditions like posted speed limits or change in traffic light state. While numerous driver models exist, parametric driver models provide the flexibility required to represent variability in real-world driving through different combinations of model parameters. These model parameters must be informed by a representative set of parameter values for the driver model to adequately represent a real-world driver. It stands to reason that determining the parameter values from real-world driving data would serve the purpose of representing a real-world driver. While the published literature is replete with techniques and consequences of using real-world driving data to determine parameter values for parametric driver models, none have explored them in the context of using shorter driving features where the parameter values may change over the course of a single trip for the same driver. In this study we consider the “Intelligent Driver Model” (IDM) as the driver behavior model to explore and real-world driving data from the Transportation Secure Data Center (TSDC) maintained at the National Renewable Energy Laboratory (NREL). The TSDC includes real-world travel data from across the United States, from which NREL has created a wide range of driving routes consisting of road features such as speed limits, stop locations and turn locations. The real-world driving data are categorized into different driving regimes and extracted into driving features. The driving features are then used to calibrate the parameter values for the IDM. The distribution of the parameters and the relationships among them are reported. The insights obtained from this study enable judicious usage of IDM or similar parametric driver models to represent baseline human driver behavior in simulations.

27 ARPA - Advanced Research Projects Agency-Energy↗

Cosmological constraints from CODEX galaxy clusters spectroscopically confirmed by SDSS-IV/SPIDERS DR16

ABSTRACT This paper presents a cosmological analysis based on the properties of X-ray selected clusters of galaxies from the CODEX survey which have been spectroscopically followed up within the SPIDERS programme as part of the sixteenth data release (DR16) of SDSS-IV. The cosmological sub-sample contains a total of 691 clusters over an area of 5350 deg2 with newly measured optical properties provided by a reanalysis of the CODEX source catalogue using redMaPPer and the DESI Legacy Imaging Surveys (DR8). Optical richness is used as a proxy for the cluster mass, and the combination of X-ray, optical, and spectroscopic information ensures that only confirmed virialized systems are considered. Clusters are binned in observed redshift, $\tilde{z} \in \left[0.1, 0.6 \right)$ and optical richness, $\tilde{\lambda } \in \left[25, 148 \right)$ and the number of clusters in each bin is modelled as a function of cosmological and richness–mass scaling relation parameters. A high-purity sub-sample of 691 clusters is used in the analysis and best-fitting cosmological parameters are found to be $\Omega _{m_{0}}=0.34^{+0.09}_{-0.05}$ and $\sigma _8=0.73^{+0.03}_{-0.03}$. The redshift evolution of the self-calibrated richness–mass relation is poorly constrained due to the systematic uncertainties associated with the X-ray component of the selection function (which assumes a fixed X-ray luminosity–mass relation with h = 0.7 and $\Omega _{m_{0}}=0.30$). Repeating the analysis with the assumption of no redshift evolution is found to improve the consistency between both cosmological and scaling relation parameters with respect to recent galaxy cluster analyses in the literature.

59 BASIC BIOLOGICAL SCIENCES↗

Scalability of Several Asynchronous Many-Task Models for In Situ Statistical Analysis.

This report is a sequel to [PB16], in which we provided a first progress report on research and development towards a scalable, asynchronous many-task, in situ statistical analysis engine using the Legion runtime system. This earlier work included a prototype implementation of a proposed solution, using a proxy mini-application as a surrogate for a full-scale scientific simulation code. The first scalability studies were conducted with the above on modestly-sized experimental clusters. In contrast, in the current work we have integrated our in situ analysis engines with a full-size scientific application (S3D, using the Legion-SPMD model), and have conducted nu- merical tests on the largest computational platform currently available for DOE science ap- plications. We also provide details regarding the design and development of a light-weight asynchronous collectives library. We describe how this library is utilized within our SPMD- Legion S3D workflow, and compare the data aggregation technique deployed herein to the approach taken within our previous work.

97 MATHEMATICS AND COMPUTING↗

The Development of New Solar Indices for use in Thermospheric Density Modeling

New solar indices have been developed to improve thermospheric density modeling for research and operational purposes. Out of 11 new and 4 legacy indices and proxies, we have selected three (F10.7, S10.7, and M10.7) for use in the new JB2006 empirical thermospheric density model. In this work, we report on the development of these solar irradiance indices. The rationale for their use, their definitions, and their characteristics, including the ISO 21348 spectral category and sub-category, wavelength range, solar source temperature region, solar source feature, altitude region of terrestrial atmosphere absorption at unit optical depth, and terrestrial atmosphere thermal processes in the region of maximum energy absorption, are described. We also summarize for each solar index, the facility and instrument(s) used to observe the solar emission, the time frame over which the data exist, the measurement cadence, the data latency, and the research as well as operational availability. The new solar indices are provided in forecast (http://SpaceWx.com) as well as real-time and historical (http://sol.spacenvironment.net/~jb2006/) time frames. We describe the forecast methodology, compare results with actual data for active and quiet solar conditions, and compare improvements in F10.7 forecasting with legacy High Accuracy Satellite Drag Model (HASDM) and NOAA SEC forecasts.

Tobiska, W. Kent↗

Measurement of correlations among net-charge, net-proton, and net-kaon multiplicity distributions in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV

Correlations among conserved quantum numbers, such as the net-electric charge, the net-baryon, and the net-strangeness in heavy-ion collisions, are crucial for exploring the QCD phase diagram. In this paper, these correlations are investigated using net-proton number (as a proxy for the net-baryon), net-kaon number (for the net-strangeness), and net-charged particle number in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV with the ALICE detector. The observed correlations deviate from the Poissonian baseline, with a more pronounced deviation at LHC energies than at RHIC. Theoretical calculations of the Thermal-FIST hadron resonance gas model, HIJING, and EPOS LHC event generators are compared with experimental results, where a significant impact of resonance decays is observed. Thermal-FIST calculations under the grand canonical and canonical ensembles highlight significant differences, underscoring the role of local charge conservation in explaining the data. Recent lattice QCD studies have demonstrated that the magnetic field generated by spectator protons in heavy-ion collisions affects susceptibility ratios, in particular those related to the net-electric charge and the net-baryon numbers. The experimental findings are in qualitative agreement with the expectations of lattice QCD.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

X-Ray Emission from the Sun in Its Youth and Old Age

We have obtained ROSAT PSPC (Roentgen Satellite Position Sensitive Proportional Counter) pointed observations of two nearby G stars of ages 70 Myr and 9.5 Gyr that are of unique importance as proxies for the Sun at the two extremes of its main-sequence evolutionary lifetime. The younger star, HD 129333 (EK Dra; G0 V), a rapid rotator with a 2.7 day period, is a strong source with an X-ray luminosity L(x)(0.2-2.4 keV) = (7.5-11.5) x 10(exp 29) erg/s. Modeling suggests a two-temperature corona with T(1) = (2.0 +/- 0.3) x 10(exp 6) K and T(2) = (9.7 +/- 0.3) x 10(exp 6) K (formal uncertainties). A continuous emission measure distribution, increasing to higher temperatures and with a cutoff at (20-30) x 10(exp 6) K, yields even better fits to the data. The old star, beta Hyi (HR 98; G2 IV), represents the Sun in the future, near the end of its hydrogen-core burning stage, when it should be rotating more slowly (present P(rot) = 25.4 day) and should have lower levels of activity. The ROSAT measurements yield L(x) = (0.9-3.0) x 10(exp 27) ergs/s and a rather cool, single coronal temperature of T = (1.7 +/- 0.4) x 10(exp 6) K. For comparison, the Sun has L(x) approx. equal to 2 x 10(exp 27) ergs/s and a coronal temperature of about T = 2 x 10(exp 6) K. These stars provide information on the decline of the stellar (and specifically solar) magnetic activity from extreme youth to old age. HD 129333 is also important in that it yields an estimate of the solar soft X-ray flux in the early solar system at the epoch of the terminal stages of planetary accretion.

Dorren, J. D.↗

Iterative sampling of expensive simulations for faster deep surrogate training

Deep neural network (DNN) surrogates of expensive physics simulations are enabling a rapid change in the way that common experimental design and analysis tasks are approached. Surrogate models allow simulations to be performed in parallel and separately from downstream tasks, thereby enabling analyses that would be impossible with the simulation in-the-loop; surrogates based on DNNs can effectively emulate diverse non-scalar data of the types collected in fusion and laboratory-astrophysics experiments. The challenge is in training the surrogate model, for which large ensembles of physics simulations must be run, preferably without wasting computational effort on uninteresting simulations. Here, in this paper, we present an iterative sampling scheme that can preferentially propose simulations in interesting regions of parameter space without neglecting unexplored regions, allowing high-quality and wide-ranging surrogate models to be trained using 2–3 times fewer simulations compare to space-filling designs. Our approach uses an explicit importance function defined on the simulation output space, balanced against a measure of simulation density which serves as a proxy for surrogate accuracy. It is easy to implement and can be tuned to find interesting simulations early in the study, allowing surrogates to be trained quickly and refined as new simulations become available; this represents an important step towards the routine generation of deep surrogate models quickly enough to be truly relevant to experimental work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

O(N) Quantum Molecular Dynamics (QMD) Algorithms, Code Requirements, Release Notes, and Library Deployment

The Exascale Computing Project (ECP) Co-design Center for Particle Applications (CoPA) is addressing the challenges for particle applications to run on upcoming exascale computing architectures. This includes short- and long-range molecular dynamics (MD), particle-in-cell (PIC), and quantum molecular dynamics (QMD). This document addresses $\mathcal{O}(N)$ QMD. The goal of this CoPA team is to develop $\mathcal{O}(N)$ linear scaling methods for QMD optimized for new and next generation architectures. This will provide the capability to simulate larger systems and run for longer timeframes than is possible today. Focus is initially on the density matrix construction. This requires efficient matrix formats and linear algebra operations for distributed architectures. A significant reduction in the wallclock time for an MD timestep is necessary to achieve practical large scale QMD. Equally challenging is to develop scalable distributed implementations. This requires the exploration of efficient data formats (dense, sparse, etc.), data decomposition and distribution across compute nodes using block- and graph-based methods. Traditional programming models such as MPI/OpenMP, hybrid MPI/OpenMP/Accelerator, as well as asynchronous task-based models will be explored to achieve maximum parallelism. Parallelism at the level of a matrix operation, solver, and/or application needs to be explored to gain the most benefit on new and emerging architectures. This report serves as a “living document” of the additions and enhancements throughout the ECP CoPA project for $\mathcal{O}(N)$ QMD. The following sections address the applications we are working with, the quantum chemistry libraries (BML and PROGRESS), and our co-design exploration through the use of a proxy application.

36 MATERIALS SCIENCE↗

Basin-scale estimates of oceanic primary production by remote sensing - The North Atlantic

The monthly averaged CZCS data for 1979 are used to estimate annual primary production at ocean basin scales in the North Atlantic. The principal supplementary data used were 873 vertical profiles of chlorophyll and 248 sets of parameters derived from photosynthesis-light experiments. Four different procedures were tested for calculation of primary production. The spectral model with nonuniform biomass was considered as the benchmark for comparison against the other three models. The less complete models gave results that differed by as much as 50 percent from the benchmark. Vertically uniform models tended to underestimate primary production by about 20 percent compared to the nonuniform models. At horizontal scale, the differences between spectral and nonspectral models were negligible. The linear correlation between biomass and estimated production was poor outside the tropics, suggesting caution against the indiscriminate use of biomass as a proxy variable for primary production.

Platt, Trevor↗