Engineering topics
Gnanadesikan, Anand
Publications and source records attributed to Gnanadesikan, Anand.
How Well do Earth System Models Capture Apparent Relationships Between Phytoplankton Biomass and Environmental Variables?
Abstract As phytoplankton form the base of the marine food web, understanding the controls on their abundance is fundamental to understanding marine ecology and its sensitivity to global climate change. While many Earth System Models (ESMs) predict phytoplankton biomass, it is unclear whether they properly capture the mechanistic relationships that control this quantity in the real ocean. We used Random Forest analysis to analyze the output of 13 ESMs as well as two observational data sets. The target variable was phytoplankton carbon and the predictors included environmental parameters known to influence phytoplankton, including nutrients, light, mixed layer depth, salinity, temperature, and upwelling. We examined the following: (a) What fractions of variability in ESMs and observations can be linked to the large‐scale environmental variables simulated by ESMs? (b) What are the dominant predictors and relationships affecting phytoplankton biomass? (c) How well do ESMs simulate phytoplankton carbon and do they simulate the relationships we see in observations? About 88%–96% of the variability in observational data sets and greater than 98% in the ESMs was accounted for by environmental variables known to influence phytoplankton biomass. The dominant predictors in the observational data sets were shortwave radiation and dissolved iron, with temperature and ammonium also relatively important. All the ESMs show that shortwave radiation is the most important variable and most of them predict the right sign of sensitivity to most variables. However, the models predict that biomass reaches maximum levels at unrealistically low levels of iron and unrealistically high levels of light.
Improving the representation of isopycnal mixing in E3SM (Final Report)
The DOE-supported Energy Exascale Earth System Model (E3SM) is a major attempt by the Department of Energy to develop a new Earth System Model using unstructured grids, which allow for the model to concentrate resolution where it is most needed and avoid numerical artifacts experienced by previous Earth System Models (ESMs) using regular grids in which the lateral resolution becomes extremely fine as the resolution shrinks at the poles. This effort required rewriting many of the algorithms previously developed for regular grids. One algorithm that was problematic in the first version of E3SM was the representation of mixing due to turbulent ocean eddies with scales smaller than the model grid. These mesoscale eddies are the primary way by which tracers are stirred along density surfaces in the ocean interior. In particular, the process of isopycnal mixing exchanges fresher, more oxygenated waters from polar regions with saltier, nutrient-rich and oxygen poor waters from tropical regions. Previous work in our group has shown that this mechanism is vitally important for bringing oxygen into poorly ventilated tropical regions (Gnanadesikan et al., 2012, 2013) and plays a significant role in determining the uptake of anthropogenic carbon dioxide (Gnanadesikan et al., 2015). However, in E3SMv1 this process was turned off, as turning it on caused the model to become unstable. Additionally, the rate of isopycnal mixing is directly proportional to a mixing coefficient $A_{Redi}$ whose value varies from less than 400 m 2 /s to 2000 m 2 /s across contemporary climate models. The main thrusts of this proposal therefore were to: 1.) Improve the numerics of isopycnal mixing in E3SM; 2.) Identify ways in which climate models are sensitive to the isopycnal mixing coefficient; and 3.) Develop better representations of the isopycnal mixing coefficient.
Characterizing the Roles of Biogeochemical Cycling and Ocean Circulation in Regulating Marine Copper Distributions
Copper (Cu) is a key micronutrient for marine phytoplankton. Its oceanic biogeochemical cycle has elicited considerable attention due to dissolved Cu exhibiting a unique linear profile with depth. Several processes have been proposed for explaining this behavior. Here, in this study, we characterize the relationships between the observed Cu, PO 4 3– , and Si on a global scale. We find that the depth profiles of Cu resemble those of Si more than of PO 4 3– in the global ocean. To understand their relationships, we couple the biogeochemical and internal circulation processes in a model, fitting optimal Cu:PO 4 3– uptake ratios and remineralization length-scales to replicate the marine Cu distributions. The modeling results suggest that Cu uptake needs in the Southern Ocean are substantially higher than those in other oceanic regions. In addition, our modeling results indicate a deep Cu remineralization in the global ocean. We offer an alternative mechanism that relies on biogeochemical cycling and internal circulation to produce the linear depth profiles of dissolved Cu. Our results suggest that diatoms are likely the major phytoplankton dominating oceanic Cu cycling.
Dataset and scripts for manuscript "Using Neural Network Ensembles to Separate Ocean Biogeochemical and Physical Drivers of Phytoplankton Biogeography in Earth System Models"
Please note: The title of this version contains an updated title for the manuscript compared to the previous version of this dataset. This is only due to title updates during the peer review process for the manuscript. The zip file contains the scripts, functions, and source files for the manuscript titled "Using Neural Network Ensembles to Separate Ocean Biogeochemical and Physical Drivers of Phytoplankton Biogeography in Earth System Models." The manuscript has been submitted for peer review. Please consult the README file for information on the specifications of the files. These files may occasionally be updated to add annotations to the scripts to make them more user friendly and to correct any errors.
Earth System Model Development - Energy Exascale Earth System Model (2020 Principal Investigator Meeting Report)
As part of the U.S. Department of Energy (DOE) Office of Biological and Environmental Research (BER), the Earth System Model Development (ESMD) program area supports significant research in Earth system model development, including major efforts in the Energy Exascale Earth System Model (E3SM) project. Convened virtually on October 26–29, the 2020 ESMD-E3SM Principal Investigator (PI) Meeting provided a forum to survey ESMD sponsored research, recent progress made by the ESMD science team and the broader Earth system modeling community within and outside of DOE, foster interactions across ESMD projects, and discuss research gaps and opportunities for future advancements.
Scaling Global Warming Impacts on Ocean Ecosystems: Lessons From a Suite of Earth System Models
An important technique used by climate modelers to isolate the impacts of increasing greenhouse gasses on Earth System processes is to simulate the impact of an abrupt increase in carbon dioxide. The spatial pattern of change provides a “fingerprint” that is generally much larger than natural variability. Insofar as the response to radiative forcing is linear (the impact of quadrupling CO 2 is twice the impact of doubling CO 2 ) this fingerprint can then be used to estimate the impact of historical greenhouse gas forcing. However, the degree to which biogeochemical cycles respond linearly to radiative forcing has rarely been tested. In this paper, we evaluate which ocean biogeochemical fields are likely to respond linearly to changing radiative forcing, which ones do not, and where linearity breaks down. We also demonstrate that the representation of lateral mixing by mesoscale eddies, which varies significantly across climate models, plays an important role in modulating the breakdown of linearity. Globally integrated surface rates of biogeochemical cycling (primary productivity, particulate export) respond in a relatively linear fashion and are only moderately sensitive to mixing. By contrast, the habitability of the interior ocean (as determined by hypoxia and calcite supersaturation) behaves non-linearly and is very sensitive to mixing. This is because the deep ocean, as well as certain regions in the surface ocean, are very sensitive to the magnitude of deep wintertime convection. The cessation of convection under global warming is strongly modulated by the representation of eddy mixing.
The Impact of Parameterized Lateral Mixing on the Antarctic Circumpolar Current in a Coupled Climate Model
This study examines the impact of changing the lateral diffusion coefficient ARedi on the transport of the Antarctic Circumpolar Current (ACC). The lateral diffusion coefficient ARedi is poorly constrained, with values ranging across an order of magnitude in climate models. The ACC is difficult to accurately simulate, and there is a large spread in eastward transport in the Southern Ocean (SO) in these models. This paper examines how much of that spread can be attributed to different eddy parameterization coefficients. A coarse-resolution, fully coupled model suite was run with A Redi = 400, 800, 1200, and 2400 m 2 s -1 . Additionally, two simulations were run with two-dimensional representations of the mixing coefficient based on satellite altimetry. Relative to the 400 m 2 s -1 case, the 2400 m 2 s -1 case exhibits 1) an 11% decrease in average wind stress from 50° to 65°S, 2) a 20% decrease in zonally averaged eastward transport in the SO, and 3) a 14% weaker transport through the Drake Passage. The decrease in transport is well explained by changes in the thermal current shear, largely due to increases in ocean density occurring on the northern side of the ACC. In intermediate waters these increases are associated with changes in the formation of intermediate waters in the North Pacific. We hypothesize that the deep increases are associated with changes in the wind stress curl allowing Antarctic Bottom Water to escape and flow northward.
Interdecadal Trichodesmium Variability in Cold North Atlantic Waters
Studies of the nitrogen cycle in the ocean generally assume that the distribution of the marine diazotroph, Trichodesmium, is restricted to warm, tropical, and subtropical oligotrophic waters. Here we show evidence that Trichodesmium are widely distributed in the North Atlantic. We report an approximately vefold increase during the 1980s and 1990s in Trichodesmium presence near the British Isles with respect to the average over the last 50 years. A potential explanation is an increase in the Saharan dust source starting in the 1980s, coupled with changes in North Atlantic winds that opened a pathway for dust transport. Results from a coarse-resolution model in which winds vary but iron deposition is climatologically fixed suggest frequent nitrogen limitation in the region and reversals of the Portugal current, but it does not simulate the observed changes in Trichodesmium. Our results suggest that Trichodesmium may be capable of growth at temperatures below 20C and challenge assumptions about their latitudinal distribution. Therefore, we need to reevaluate assumptions about the temperature limitations of Trichodesmium and the dinitrogen (N2) xation capabilities of extratropical strains, which may have important implications for the global nitrogen budget.