SEARCH · Engineering Papers

Results for “parameterization”

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.

SourceResearch connectionFrom yearThrough year

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

At least 865 records · Page 48

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Assessment of Precipitating Marine Stratocumulus Clouds in the E3SMv1 Atmosphere Model: A Case Study from the ARM MAGIC Field Campaign

This paper presents a process-oriented evaluation of precipitating stratocumulus and its transition to cumulus in version 1 of the Energy Exascale Earth System Model (E3SMv1) using comprehensive case-study observations from a field campaign of the Atmospheric Radiation Measurement program (ARM). The E3SMv1 single-column model (SCM) of the marine boundary layer and its low clouds and precipitation are compared to observations including subcloud drizzle retrievals from a combination of Doppler radar and lidar backscatter measurements. Additionally, the SCM is also compared to a large-eddy simulation (LES) of the same case. The combination of advanced remote sensing observations and LES is a powerful framework to evaluate the physical parameterizations of large-scale models. Given the observed large-scale environment, the E3SMv1 SCM realistically represents the evolution of clouds and boundary layer structure during the stratocumulus-to-cumulus transition. The model well simulates the liquid water path and its diurnal cycle in the stratocumulus period as well as the two-layer vertical thermodynamic structure and lower cloud fraction in the transition period. E3SMv1’s success in simulating the cloud in the stratocumulus period permitted examination of its precipitation processes. Here problems were identified with E3SMv1 producing an unrealistically small subcloud precipitation fraction, an unrealistic double peak in the vertical profiles of precipitation mass, and drizzle that evaporates too close to the surface. Further model diagnostics determined that these unrealistic characteristics resulted from an overly long microphysics time step and an unrealistic parameterization of the precipitation fraction. These results imply that careful consideration of these issues is needed in order to better simulate precipitation processes in marine stratocumulus.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Fast Computational Algorithms for Partial Differential Equations and Uncertainty Quantifications

This project concerned the construction, testing and analysis of computational algorithms for solving parameterized and stochastic partial differential equations. The study and understanding of equations of this type is of fundamental importance in numerous engineering and scientific applications. Examples include simulation of plasma dynamics in models of electric propulsion and nuclear fusion, simulation of multiphase flows, such as the flow of water, gas and oil in reservoirs, and structural analysis of the dependence of structures on materials. Parametrization is used in such settings when properties of the models such as viscosity of fluids or electric resistivity of materials are not precisely understood and instead are treated as random variables. The resulting solutions are themselves random, and having such solutions will enable engineers to use probabilistic methods to assess the likelihood of events, for example, whether a pollutant in a liquid will exceed a limit, and to use such analyses to develop ways to ensure positive outcomes. Construction of accurate (high resolution) computational solutions is expensive, requiring significant computer time and computational resources, and there is need to reduce computational cost to make simulation useful and effective. The aim of the project was to construct algorithms to efficiently compute surrogate solutions to parameterized problems to allow for efficient and accurate simulation. The technical approach used focused on two related strategies, based on rank-reduction methods and reduced-order models. These methods construct surrogate solutions of parameter-dependent models by projection or interpolation into low-dimensional approximation spaces. Cost savings are achieved if the low-dimensional spaces can be identified and constructed efficiently and if the resulting low-dimensional algebraic systems can be solved cheaply. Accomplishments include: Theoretical and empirical demonstration of the effectiveness of fast multigrid solution strategies for computing low-rank representations of parameter-dependent solutions to discrete partial differential equations, including the first proof establishing so-called textbook convergence properties for low-rank methods. Development of efficient solution algorithms for solving nonlinear parameter-dependent partial differential equations used in models of fluid dynamics. Developent of efficient algorithms for low-rank representation of solutions of time-dependent simulations of fluid dynamics using multi-dimensional tensor representations of solutions.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Understanding Processes Controlling the Temporal and Spatial Variations of PBL Structures Over the ARM SGP Site

The surface heat, moisture, and momentum fluxes are transferred to the atmosphere above through the planetary boundary layer (PBL), where vertical mixing due to turbulent eddies of different sizes plays critical roles. Therefore, reliably representing PBL processes in numerical models is critical for weather, climate, and air quality prediction. Currently, there are over ten PBL schemes that are selectable within the advanced research version of the Weather Research and Forecasting (WRF) model, indicative of the challenges in capturing the impacts of turbulence within the PBL in models. Further improvements in PBL parameterizations are needed for both weather and climate models, as emphasized in many recent national reports, but require an advanced understanding of the underlying boundary layer processes from observations. This project takes advantage of Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) investments in the atmospheric boundary layer observations and Large-Eddy Simulation (LES) ARM Symbiotic Simulation and Observation (LASSO) simulations to characterize PBL structures, understand key physical processes controlling the mixed layer development, and to evaluate PBL parameterization.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2019

Demeter calibration and uncertainty analysis

Demeter is a community spatial downscaling model that disaggregates land use and land cover changes projected by integrated human-Earth system models. Demeter has not been intensively calibrated, and we still lack a good knowledge about its sensitivity to key parameters and the parameter uncertainties. We used long-term global satellite-based land cover records to calibrate key Demeter parameters. The results identified the optimal parameter values and showed that the parameterization substantially improved the model’s performance. The parameters of intensification ratio and selection threshold were the most sensitive and needed to be carefully tuned, especially for regional applications. Further, small parameter uncertainties after calibration can be inflated when propagated into future scenarios, suggesting that users should consider the parameterization equifinality to better account for the uncertainties in the Demeter downscaled products. Our study provides a key reference for Demeter users, and ultimately contribute to reducing the uncertainties in Earth system model simulations.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Turbulent kinetic energy over large offshore wind farms observed and simulated by the mesoscale model WRF (3.8.1)

Abstract. Wind farms affect local weather and microclimates;hence, parameterizations of their effects have been developed for numerical weather prediction models.While most wind farm parameterizations (WFPs) include drag effects of wind farms, models differ on whether or not an additionalturbulent kinetic energy (TKE) source should be included in theseparameterizations to simulate the impact of wind farms on the boundary layer.Therefore, we use aircraft measurements above large offshore wind farms in stable conditions to evaluateWFP choices. Of the three case studies we examine, we find the simulated ambient backgroundflow to agree with observations of temperature stratification and winds. This agreement allows us toexplore the sensitivity of simulated wind farm effects with respect to modeling choicessuch as whether or not to include a TKE source, horizontal resolution, vertical resolutionand advection of TKE. For a stably stratified marine atmospheric boundary layer (MABL), a TKE source and a horizontal resolutionon the order of 5 km or finer are necessary to represent the impact of offshore wind farms on the MABL.Additionally, TKE advection results in excessively reduced TKE over the wind farms, which in turn causesan underestimation of the wind speed deficit above the wind farm. Furthermore, using fine verticalresolution increases the agreement of the simulated wind speed with satellite observations of surface wind speed.

17 WIND ENERGY↗

DOE OSTI2014

Perturbed Physics Ensemble Simulations of Cirrus on the Cloud System-resolving Scale

In this study, the effect of uncertainties in the parameterization of ice microphysical processes and initial conditions on the variability of cirrus microphysical and radiative properties are investigated in a series of cloud system-resolving perturbed physics ensemble (PPE) and initial condition ensemble (ICE) simulations. Three cirrus cases representative of mid-latitude, subtropical and tropical cirrus are examined. It is found that the variability in cirrus properties induced by perturbing uncertain parameters in ice microphysics parameterizations outweighs the variability induced by perturbing the initial conditions in midlatitude and subtropical cirrus. However, in tropical anvil cirrus the variability in the PPE and ICE simulations is about the same order of magnitude. The cirrus properties showing the largest sensitivity are ice water content (IWC) and cloud thickness whereas the averaged high cloud cover is only marginally affected. Changes in cirrus ice water path and outgoing longwave radiation are controlled primarily by changes in IWC and cloud thickness but not by changes is the averaged high cloud cover. The change in the vertical distribution of cloud fraction and cloud thickness is caused by changes in cirrus cloud base whereas cloud top is not sensitive to either perturbed physics or perturbed initial conditions. In all cirrus cases, the top three parameters controlling the microphysical variability and radiative impact of cirrus clouds are ice fall speeds, ice autoconversion size thresholds and heterogeneous ice nucleation. Changes in the ice deposition coefficient do not affect the ice water path and outgoing longwave radiation. Similarly, changes in the number concentration of aerosols available for homogeneous freezing have virtually no effect on the microphysical and radiative properties of midlatitude and subtropical cirrus but only little impact on tropical anvil cirrus. Overall, the sensitivity of cirrus microphysical and radiative properties to uncertainties in ice microphysics is largest for midlatitude cirrus and smallest for tropical anvil cirrus.

cirrus↗

DOE OSTI2014

Intercomparison of the Cloud Water Phase among Global Climate Models

Mixed-phase clouds (clouds that consist of both cloud droplets and ice crystals) are frequently present in the Earth’s atmosphere and influence the Earth’s energy budget through their radiative properties, which are highly dependent on the cloud water phase. In this study, the phase partitioning of cloud water is compared among six global climate models (GCMs) and with Cloud and Aerosol Lidar with Orthogonal Polarization retrievals. It is found that the GCMs predict vastly different distributions of cloud phase for a given temperature, and none of them are capable of reproducing the spatial distribution or magnitude of the observed phase partitioning. While some GCMs produced liquid water paths comparable to satellite observations, they all failed to preserve sufficient liquid water at mixed-phase cloud temperatures. Our results suggest that validating GCMs using only the vertically integrated water contents could lead to amplified differences in cloud radiative feedback. The sensitivity of the simulated cloud phase in GCMs to the choice of heterogeneous ice nucleation parameterization is also investigated. The response to a change in ice nucleation is quite different for each GCM, and the implementation of the same ice nucleation parameterization in all models does not reduce the spread in simulated phase among GCMs. The results suggest that processes subsequent to ice nucleation are at least as important in determining phase and should be the focus of future studies aimed at understanding and reducing differences among the models.

cloud water phase↗

DOE OSTI2015

Impacts of a mushy‐layer thermodynamic approach in global sea‐ice simulations using the CICE sea‐ice model

Abstract We perform global simulations of the Los Alamos sea‐ice model, CICE, with a new thermodynamics component that has a fully prognostic, variable bulk salinity vertical profile based on mushy layer physics. The processes of gravity drainage, melt‐water flushing and snow‐ice formation are parameterized to allow the bulk salinity to evolve with time. We analyze the seasonal and spatial variation of sea‐ice bulk salinity, area, volume and thickness and compare these quantities to simulations using the previous thermodynamic component. Adjusting one of the gravity drainage parameters, we find good agreement between simulation results and fieldwork ice‐core observations of sea‐ice bulk salinity. As expected, bulk salinity is highest during periods of ice growth and lowest after periods of ice melt. In the northern hemisphere the new thermodynamics component produces thicker ice than the previous thermodynamics component. Of the nine major differences between the two models, differences in how salinities are calculated and how melt‐pond flushing is parameterized are the principal causes of this thickness difference. Thickness differences are smaller in the southern hemisphere than in the northern hemisphere since a greater fraction of ice melts, and differences cannot accumulate year‐on‐year. Model differences in how ice thickness changes and snow‐ice formation are calculated are the most important causes of the different thickness between the two thermodynamic components in the southern hemisphere. The melt‐pond area and volume are found to be highly sensitive to a parameter choice controlling drainage through macroscopic holes in the ice, in both hemispheres.

Turner, Adrian K.↗

DOE OSTI2014

Evaluation of convection‐permitting model simulations of cloud populations associated with the Madden‐Julian Oscillation using data collected during the AMIE/DYNAMO field campaign

Abstract Regional convection‐permitting model simulations of cloud populations observed during the 2011 Atmospheric Radiation Measurement (ARM) Madden‐Julian Oscillation Investigation Experiment/Dynamics of the Madden‐Julian Oscillation Experiment (AMIE/DYNAMO) field campaign are evaluated against ground‐based radar and ship‐based observations. Sensitivity of model simulated reflectivity, surface rain rate, and cold pool statistics to variations of raindrop breakup/self‐collection parameters in four state‐of‐the‐art two‐moment bulk microphysics schemes in the Weather Research and Forecasting (WRF) model is examined. The model simulations generally overestimate reflectivity from large and deep convective cells, and underestimate stratiform rain and the frequency of cold pools. In the sensitivity experiments, introduction of more aggressive raindrop breakup or decreasing the self‐collection efficiency increases the cold pool occurrence frequency in all of the simulations, and slightly reduces the reflectivity and precipitation statistics bias in some schemes but has little effect on the overall mean surface precipitation. Both the radar observations and model simulations of cloud populations show an approximate power law relationship between convective echo‐top height and equivalent convective cell radius.

Hagos, Samson↗

DOE OSTI2015

Evaluation of cloud-resolving model simulations of midlatitude cirrus with ARM and A-train observations

Cirrus clouds are ubiquitous in the upper troposphere and still constitute one of the largest uncertainties in climate predictions. Our paper evaluates cloud-resolving model (CRM) and cloud system-resolving model (CSRM) simulations of a midlatitude cirrus case with comprehensive observations collected under the auspices of the Atmospheric Radiation Measurements (ARM) program and with spaceborne observations from the National Aeronautics and Space Administration A-train satellites. The CRM simulations are driven with periodic boundary conditions and ARM forcing data, whereas the CSRM simulations are driven by the ERA-Interim product. Vertical profiles of temperature, relative humidity, and wind speeds are reasonably well simulated by the CSRM and CRM, but there are remaining biases in the temperature, wind speeds, and relative humidity, which can be mitigated through nudging the model simulations toward the observed radiosonde profiles. Simulated vertical velocities are underestimated in all simulations except in the CRM simulations with grid spacings of 500 m or finer, which suggests that turbulent vertical air motions in cirrus clouds need to be parameterized in general circulation models and in CSRM simulations with horizontal grid spacings on the order of 1 km. The simulated ice water content and ice number concentrations agree with the observations in the CSRM but are underestimated in the CRM simulations. The underestimation of ice number concentrations is consistent with the overestimation of radar reflectivity in the CRM simulations and suggests that the model produces too many large ice particles especially toward the cloud base. Simulated cloud profiles are rather insensitive to perturbations in the initial conditions or the dimensionality of the model domain, but the treatment of the forcing data has a considerable effect on the outcome of the model simulations. Despite considerable progress in observations and microphysical parameterizations, simulating the microphysical, macrophysical, and radiative properties of cirrus remains challenging. Comparing model simulations with observations from multiple instruments and observational platforms is important for revealing model deficiencies and for providing rigorous benchmarks. But, there still is considerable need for reducing observational uncertainties and providing better observations especially for relative humidity and for the size distribution and chemical composition of aerosols in the upper troposphere.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2014

Precipitation characteristics of CAM5 physics at mesoscale resolution during MC3E and the impact of convective timescale choice

Abstract The physics suite of the Community Atmosphere Model version 5 (CAM5) has recently been implemented in the Weather Research and Forecasting (WRF) model to explore the behavior of the parameterization suite at high resolution and within the controlled setting of a limited area model. The initial paper documenting this capability characterized the behavior for a northern high‐latitude region. This paper characterizes the precipitation characteristics for continental, midlatitude, springtime conditions during the Midlatitude Continental Convective Clouds Experiment (MC3E) over the central United States. This period exhibited a range of convective conditions from those driven strongly by large‐scale synoptic regimes to more locally driven convection. The study focuses on the precipitation behavior at 32 km grid spacing to better anticipate how the physics will behave in a global model when used at similar grid spacing in the coming years. Importantly, one change to the Zhang‐McFarlane deep convective parameterization when implemented in WRF was to make the convective timescale parameter an explicit function of grid spacing. This study examines the sensitivity of the precipitation to the default value of the convective timescale in WRF, which is 600 s for 32 km grid spacing, to the value of 3600 s used for 2° grid spacing in CAM5. For comparison, a 1200 s and an infinite convective timescale are also used. The results show that the 600 s timescale gives the most accurate precipitation amount over the central United States. However, this setting has the worst precipitation diurnal cycle, with the convection too tightly linked to the daytime surface heating. Longer timescales greatly improve the diurnal cycle but result in less precipitation and produce a low bias. An analysis of rain rates shows the accurate precipitation amount with the shorter timescale is assembled from an over abundance of drizzle combined with too few heavy rain events. With longer timescales, one can improve the frequency distribution, particularly for the extreme rain rates. Ultimately, without changing other aspects of the physics, one must decide between accurate diurnal timing and rain amount when choosing an appropriate convective timescale.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2015

Impact of the dynamical core on the direct simulation of tropical cyclones in a high-resolution global model

Our paper examines the impact of the dynamical core on the simulation of tropical cyclone (TC) frequency, distribution, and intensity. The dynamical core, the central fluid flow component of any general circulation model (GCM), is often overlooked in the analysis of a model's ability to simulate TCs compared to the impact of more commonly documented components (e.g., physical parameterizations). The Community Atmosphere Model version 5 is configured with multiple dynamics packages. This analysis demonstrates that the dynamical core has a significant impact on storm intensity and frequency, even in the presence of similar large-scale environments. In particular, the spectral element core produces stronger TCs and more hurricanes than the finite-volume core using very similar parameterization packages despite the latter having a slightly more favorable TC environment. Furthermore, these results suggest that more detailed investigations into the impact of the GCM dynamical core on TC climatology are needed to fully understand these uncertainties. Key Points The impact of the GCM dynamical core is often overlooked in TC assessments The CAM5 dynamical core has a significant impact on TC frequency and intensity A larger effort is needed to better understand this uncertainty

58 GEOSCIENCES↗

DOE OSTI2015

How does increasing horizontal resolution in a global climate model improve the simulation of aerosol-cloud interactions?

The Community Atmosphere Model Version 5 is run at horizontal grid spacing of 2, 1, 0.5, and 0.25 degrees, with the meteorology nudged towards the Year Of Tropical Convection analysis, and cloud simulators and the collocated A-Train satellite observations are used to explore the resolution dependence of aerosol-cloud interactions. The higher-resolution model produces results that agree better with observations, showing an increase of susceptibility of cloud droplet size, indicating a stronger first aerosol indirect forcing (AIF), and a decrease of susceptibility of precipitation probability, suggesting a weaker second AIF. The resolution sensitivities of AIF are attributed to those of droplet nucleation and precipitation parameterizations. The annual average AIF in the northern hemisphere mid-latitudes (where most anthropogenic emissions occur) in the 0.25° model is reduced by about 1 W m⁻² (-30%) compared to the 2° model, leading to a 0.26 W m⁻² reduction (-15%) in the global annual average AIF.

Ma, Po -Lun↗

DOE OSTI2015

Quantifying sources and sinks of reactive gases in the lower atmosphere using airborne flux observations

Atmospheric composition is governed by the interplay of emissions, chemistry, deposition, and transport. Substantial questions surround each of these processes, especially in forested environments with strong biogenic emissions. Utilizing aircraft observations acquired over a forest in the southeast U.S., we calculate eddy covariance fluxes for a suite of reactive gases and apply the synergistic information derived from this analysis to quantify emission and deposition fluxes, oxidant concentrations, aerosol uptake coefficients, and other key parameters. Evaluation of results against state-of-the-science models and parameterizations provides insight into our current understanding of this system and frames future observational priorities. As a near-direct measurement of fundamental process rates, airborne fluxes offer a new tool to improve biogenic and anthropogenic emissions inventories, photochemical mechanisms, and deposition parameterizations.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2016

An ensemble constrained variational analysis of atmospheric forcing data and its application to evaluate clouds in CAM5: Ensemble 3DCVA and Its Application

Large-scale atmospheric forcing data can greatly impact the simulations of atmospheric process models including Large Eddy Simulations (LES), Cloud Resolving Models (CRMs) and Single-Column Models (SCMs), and impact the development of physical parameterizations in global climate models. This study describes the development of an ensemble variationally constrained objective analysis of atmospheric large-scale forcing data and its application to evaluate the cloud biases in the Community Atmospheric Model (CAM5). Sensitivities of the variational objective analysis to background data, error covariance matrix and constraint variables are described and used to quantify the uncertainties in the large-scale forcing data. Application of the ensemble forcing in the CAM5 SCM during March 2000 intensive operational period (IOP) at the Southern Great Plains (SGP) of the Atmospheric Radiation Measurement (ARM) program shows systematic biases in the model simulations that cannot be explained by the uncertainty of large-scale forcing data, which points to the deficiencies of physical parameterizations. The SCM is shown to overestimate high clouds and underestimate low clouds. These biases are found to also exist in the global simulation of CAM5 when it is compared with satellite data.

58 GEOSCIENCES↗

DOE OSTI2016

Potential sea salt aerosol sources from frost flowers in the pan-Arctic region

In order to better represent observed wintertime aerosol concentrations at Barrow, Alaska, we implemented an observationally-based parameterization for estimating sea salt production from frost flowers in the Community Earth System Model (CESM). In this work, we evaluate the potential influence of this sea salt source on the pan-Arctic (60ºN-90ºN) climate. Results show that frost flower salt emissions substantially increase the modeled surface sea salt aerosol concentration in the winter months when new sea ice and frost flowers are present. The parameterization reproduces both the magnitude and seasonal variation of the observed submicron sea salt aerosol concentration at surface in Barrow during winter much better than the standard CESM simulation without a frost-flower salt particle source. Adding these frost flower salt particle emissions increases aerosol optical depth by 10% and results in a small cooling at surface. Here, the increase in salt particle mass concentrations of a factor of 8 provides nearly two times the cloud condensation nuclei concentration, as well as 10% increases in cloud droplet number and 40% increases in liquid water content near coastal regions adjacent to continents. These cloud changes reduce longwave cloud forcing by 3% and cause a small surface warming, increasing the downward longwave flux at the surface by 2 W m -2 in the pan-Arctic under the present-day climate.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2015

Energy considerations in the Community Atmosphere Model (CAM)

An error in the energy formulation in the Community Atmosphere Model (CAM) is identified and corrected. Ten year AMIP simulations are compared using the correct and incorrect energy formulations. Statistics of selected primary variables all indicate physically insignificant differences between the simulations, comparable to differences with simulations initialized with rounding sized perturbations. The two simulations are so similar mainly because of an inconsistency in the application of the incorrect energy formulation in the original CAM. CAM used the erroneous energy form to determine the states passed between the parameterizations, but used a form related to the correct formulation for the state passed from the parameterizations to the dynamical core. If the incorrect form is also used to determine the state passed to the dynamical core the simulations are significantly different. In addition, CAM uses the incorrect form for the global energy fixer, but that seems to be less important. The difference of the magnitude of the fixers using the correct and incorrect energy definitions is very small.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2015

Complex functionality with minimal computation: Promise and pitfalls of reduced‐tracer ocean biogeochemistry models

Abstract Earth System Models increasingly include ocean biogeochemistry models in order to predict changes in ocean carbon storage, hypoxia, and biological productivity under climate change. However, state‐of‐the‐art ocean biogeochemical models include many advected tracers, that significantly increase the computational resources required, forcing a trade‐off with spatial resolution. Here, we compare a state‐of‐the art model with 30 prognostic tracers (TOPAZ) with two reduced‐tracer models, one with 6 tracers (BLING), and the other with 3 tracers (miniBLING). The reduced‐tracer models employ parameterized, implicit biological functions, which nonetheless capture many of the most important processes resolved by TOPAZ. All three are embedded in the same coupled climate model. Despite the large difference in tracer number, the absence of tracers for living organic matter is shown to have a minimal impact on the transport of nutrient elements, and the three models produce similar mean annual preindustrial distributions of macronutrients, oxygen, and carbon. Significant differences do exist among the models, in particular the seasonal cycle of biomass and export production, but it does not appear that these are necessary consequences of the reduced tracer number. With increasing CO 2 , changes in dissolved oxygen and anthropogenic carbon uptake are very similar across the different models. Thus, while the reduced‐tracer models do not explicitly resolve the diversity and internal dynamics of marine ecosystems, we demonstrate that such models are applicable to a broad suite of major biogeochemical concerns, including anthropogenic change. These results are very promising for the further development and application of reduced‐tracer biogeochemical models that incorporate “sub‐ecosystem‐scale” parameterizations.

58 GEOSCIENCES↗