Impacts of Shifts in Phytoplankton Community on Clouds and Climate via the Sulfur Cycle
Not Available
SEARCH · Search NASA
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.
Not Available
The Massachusetts Institute of Technology General Circulation Model (MITgcm) is widely used by the climate science community to simulate planetary atmosphere and ocean circulations. A defining feature of the MITgcm is that it has been developed to be compatible with an algorithmic differentiation (AD) tool, TAF, enabling the generation of tangent-linear and adjoint models. These provide gradient information which enables dynamics-based sensitivity and attribution studies, state and parameter estimation, and rigorous uncertainty quantification. Importantly, gradient information is essential for computing comprehensive sensitivities and performing efficient large-scale data assimilation, ensuring that observations collected from satellites and in-situ measuring instruments can be effectively used to optimize a large uncertain control space. As a result, the MITgcm forms the dynamical core of a key data assimilation product employed by the physical oceanography research community: Estimating the Circulation and Climate of the Ocean (ECCO) state estimate. Although MITgcm and ECCO are used extensively within the research community, the AD tool TAF is proprietary and hence inaccessible to a large proportion of these users. The new version 2 (MITgcm-AD v2) framework introduced here is based on the source-to-source AD tool Tapenade, which has recently been open-sourced. Another feature of Tapenade is that it stores required variables by default (instead of recomputing them) which simplifies the implementation of efficient, AD-compatible code. The framework has been integrated with the MITgcm model’s main branch and is now freely available.
Despite its importance for the global cycling of carbon, there are still large gaps in our understanding of the processes driving annual and seasonal carbon fluxes in the high-latitude Southern Ocean. This is due in part to a historical paucity of observations in this remote, turbulent, and seasonally ice-covered region. Here, we use autonomous biogeochemical float data spanning 6 full seasonal cycles and with circumpolar coverage of the Southern Ocean, complemented by atmospheric reanalysis, to construct a monthly climatology of the mixed layer budget of dissolved inorganic carbon (DIC). We investigate the processes that determine the annual mean and seasonal cycle of DIC fluxes in two different zones of the Southern Ocean—the Sea Ice Zone (SIZ) and Antarctic Southern Zone (ASZ). We find that, annually, mixing with carbon-rich waters at the base of the mixed layer supplies DIC which is, in the ASZ, either used for net biological production or outgassed to the atmosphere. In contrast, in the SIZ, where carbon outgassing and the biological pump are weaker, the surplus of DIC is instead advected northward to the ASZ. In other words, carbon outgassing in the southern Antarctic Circumpolar Current (ACC), which has been attributed to remineralized carbon from deep water upwelled in the ACC, is also due to the wind-driven transport of DIC from the SIZ. These results stem from the first observation-based carbon budget of the circumpolar Southern Ocean and thus provide a useful benchmark to evaluate climate models, which have significant biases in this region.
Abstract. The width of the Inter-Tropical Convergence Zone (ITCZ) affects tropical rainfall, Earth's albedo, large-scale circulation, and climate sensitivity. To better understand the ITCZ width and its effects on global climate, we present a protocol to force changes in ITCZ width in climate models. Starting from an aquaplanet configuration with a slab ocean, adding surface heat fluxes in the deep tropics forces the ITCZ to narrow, and subtracting them causes it to widen. The protocol successfully generates changes in ITCZ width in all four climate models used in this study. ITCZ width in each model responds linearly to forcing magnitude and sign. Comparing across the four climate models, a response to varying ITCZ width that is remarkably consistent among models is the ITCZ strength, which is greater the narrower the ITCZ. On the other hand, the effect of varying ITCZ width on climate sensitivity is divergent among our four models, varying even in sign. Results from this pilot study highlight the connections between surface fluxes, ITCZ width, and the wider climate. A comprehensive model intercomparison project (MIP) has the potential to advance understanding of both the physical processes shaping ITCZ width and its influence on remote atmospheric circulations and global climate.
The surface energy budget governs Arctic sea-ice growth/melt, yet observations are sparse, and reanalysis data sets suffer from systematic biases. Here, we train a neural network with observational data to bias-correct hourly ERA5 fluxes over Arctic ice-covered regions (≥70°N; sea-ice concentration >80%) for 1994–2024. Training data cover two full seasonal cycles and different sea-ice regimes. The neural network reduces RMSE for net shortwave radiation by ∼40%, downward longwave radiation by ∼16% and the total surface energy budget by ∼55%, eliminating the wintertime warm bias of ∼4 K in ERA5. Wintertime surface cooling is reduced by ∼50%, yielding thermodynamic ice-growth estimates of ∼80–120 cm, consistent with SMOS–CryoSat satellite thickness increases and in contrast to the 150–200 cm growth implied by ERA5. Our bias-corrected data capture the observed clear/cloudy states of the winter boundary layer and can be used to study Arctic climatology, evaluate climate models and drive sea-ice-ocean models.
A primary advantage to using reduced complexity climate models (RCMs) has been their ability to quickly conduct probabilistic climate projections, a key component of uncertainty quantification in many impact studies and multisector systems. Providing frameworks for such analyses has been a target of several RCMs used in studies of the future co-evolution of the human and Earth systems. In this paper, we present Matilda, an open-science R software package that facilitates probabilistic climate projection analysis, implemented here using the Hector simple climate model in a seamless and easily applied framework. The primary goal of Matilda is to provide the user with a turn-key method to build parameter sets from literature-based prior distributions, run Hector iteratively to produce perturbed parameter ensembles (PPEs), weight ensembles for realism against observed historical climate data, and compute probabilistic projections for different climate variables. This workflow gives the user the ability to explore viable parameter space and propagate uncertainty to model ensembles with just a few lines of code. The package provides significant freedom to select different scoring criteria and algorithms to weight ensemble members, as well as the flexibility to implement custom criteria. Additionally, the architecture of the package simplifies the process of building and analyzing PPEs without requiring significant programming expertise, to accommodate diverse use cases. We present a case study that provides illustrative results of a probabilistic analysis of mean global surface temperature as an example of the software application.
Model resolution plays a large role in accurately simulating the Southern Hemisphere circulation in both the ocean and atmosphere. Resolving the mesoscale field is important as it has been shown to have a significant impact on the large-scale climate in eddy-rich regions, which are regions of large CO2 absorption. The presence of ocean and atmospheric mesoscale features can affect sea surface temperatures, the strength and location of storm tracks, and many other air-sea processes. Additionally, with an improvement in resolution, the eddy kinetic energy in the ocean can be expected to change considerably. The significance model resolution has on the Southern Hemisphere is examined using the Community Climate System Model, Version 4, eddy-parameterizing and eddy-resolving simulations. The CO2 concentrations and ozone levels are specified independently to better understand how the mesoscale field responds to extreme changes in external forcing and the resulting climate impacts. Overall, in the eddy-parameterizing simulations, the ozone forcing is found to be more important than the changes in CO2 concentrations. However, in the case of the eddy-resolving simulations, the CO2 concentrations are found to be more dominant, especially in eddy-rich regions. These results demonstrate the need for an increase in model resolution for climate prediction.
Abstract Cumulus entrainment substantially regulates the earth's climate but remains poorly constrained in global climate models. Recent studies have shown that cumulus bulk entrainment (or dilution) is particularly sensitive to continentality, with the entrainment rate in simulated maritime cumuli nearly double that of continental cumuli. This study examines the impacts of such land–ocean entrainment contrasts on the current climate using 21‐year simulations with the Geophysical Fluid Dynamics Laboratory's High‐Resolution Atmospheric Model (HIRAM). In response to a 25% reduction in the HIRAM entrainment parameter c 0 over land, precipitation over tropical land regions increases by up to 40%. Along with directly facilitating enhanced convective precipitation, this c 0 reduction induces an increase in soil moisture, which may contribute to a further enhancement of convective precipitation over land. A 25% c 0 reduction over the oceans leads to more widespread modifications of convection patterns, with the strongest signal in the tropical Pacific. Deep convection shifts upstream (eastward) there, inducing enhanced large‐scale ascent over the central Pacific with compensating subsidence and reduced humidity and precipitation over the western Pacific (WP). Land–ocean variations in c 0 project onto the Pacific Walker circulation, with the 25% land reduction strengthening it by 4% and the 25% ocean reduction weakening it by 14%. These changes are driven by variations in convective and large‐scale stratiform heating over the Pacific. While reduced c 0 over land enhances diabatic heating in the Maritime Continent to strengthen the Walker circulation, reduced c 0 over the oceans decreases diabatic heating in the WP to weaken the Walker circulation.
Coupled ocean-atmospheric models suffer from the common bias of a spurious rain belt south of the central equatorial Pacific throughout the year. Observational constraints on key processes responsible for this bias are scarce. The recently available reanalysis from a coupled model system for the National Centers for Environmental Prediction (NCEP) Climate Forecast System Reanalysis (CFSR) data is a potential benchmark for climate models in this region. Its suitability for model evaluation and validation, however, needs to be established. This paper examines the mixed layer heat budget and the ocean surface currents - key factors for the sea surface temperature control in the double Inter-Tropical Convergence Zone in the central Pacific - from 5{sup o}S to 10{sup o}S and 170{sup o}E to 150{sup o}W. Two independent approaches are used. The first approach is through comparison of CFSR data with collocated station observations from field experiments; the second is through the residual analysis of the heat budget of the mixed layer. We show that the CFSR overestimates the net surface flux in this region by 23 W m{sup -2}. The overestimated net surface flux is mainly due to an even larger overestimation of shortwave radiation by 44 W m{sup -2}, which is compensated by a surface latent heat flux overestimated by 14 W m{sup -2}. However, the quality of surface currents and the associated oceanic heat transport in CFSR are not compromised by the surface flux biases, and they agree with the best available estimates. The uncertainties of the observational data from field experiments are also briefly discussed in the present study.
Marine heatwaves (MHWs), episodic periods of abnormally high sea surface temperature, severely affect marine ecosystems. Large marine ecosystems (LMEs) cover ~22% of the global ocean but account for 95% of global fisheries catches. Yet how climate change affects MHWs over LMEs remains unknown because such LMEs are confined to the coast where low-resolution climate models are known to have biases. Here, using a high-resolution Earth system model and applying a ‘future threshold’ that considers MHWs as anomalous warming above the long-term mean warming of sea surface temperatures, we find that future intensity and annual days of MHWs over the majority of the LMEs remain higher than in the present-day climate. Better resolution of ocean mesoscale eddies enables simulation of more realistic MHWs than low-resolution models. These increases in MHWs under global warming pose a serious threat to LMEs, even if resident organisms could adapt fully to the long-term mean warming.
Abstract The Pacific Decadal Oscillation (PDO), the leading mode of Pacific decadal sea surface temperature variability, arises mainly from combinations of regional air-sea interaction within the North Pacific Ocean and remote forcing, such as from the tropical Pacific and the Atlantic. Because of such a combination of mechanisms, a question remains as to how much PDO variability originates from these regions. To better understand PDO variability, the equatorial Pacific and the Atlantic impacts on the PDO are examined using several 3-dimensional partial ocean data assimilation experiments conducted with two global climate models: the CESM1.0 and MIROC3.2m. In these partial assimilation experiments, the climate models are constrained by observed temperature and salinity anomalies, one solely in the Atlantic basin and the other solely in the equatorial Pacific basin, but are allowed to evolve freely in other regions. These experiments demonstrate that, in addition to the tropical Pacific’s role in driving PDO variability, the Atlantic can affect PDO variability by modulating the tropical Pacific climate through two proposed processes. One is the equatorial pathway, in which tropical Atlantic sea surface temperature (SST) variability causes an El Niño-like SST response in the equatorial Pacific through the reorganization of the global Walker circulation. The other is the north tropical pathway, where low-frequency SST variability associated with the Atlantic Multidecadal Oscillation induces a Matsuno-Gill type atmospheric response in the tropical Atlantic-Pacific sectors north of the equator. These results provide a quantitative assessment suggesting that 12–29% of PDO variance originates from the Atlantic Ocean and 40–44% from the tropical Pacific. The remaining 27–48% of the variance is inferred to arise from other processes such as regional ocean-atmosphere interactions in the North Pacific and possibly teleconnections from the Indian Ocean.
Abstract Spatially compounding droughts over multiple regions pose amplifying pressures on the global food system, the reinsurance industry, and the global economy. Using observations and climate model simulations, we analyze the influence of various natural Ocean variability modes on the likelihood, extent, and severity of compound droughts across ten regions that have similar precipitation seasonality and cover important breadbaskets and vulnerable populations. Although a majority of compound droughts are associated with El Niños, a positive Indian Ocean Dipole, and cold phases of the Atlantic Niño and Tropical North Atlantic (TNA) can substantially modulate their characteristics. Cold TNA conditions have the largest amplifying effect on El Niño-related compound droughts. While the probability of compound droughts is ~3 times higher during El Niño conditions relative to neutral conditions, it is ~7 times higher when cold TNA and El Niño conditions co-occur. The probability of widespread and severe compound droughts is also amplified by a factor of ~3 and ~2.5 during these co-occurring modes relative to El Niño conditions alone. Our analysis demonstrates that co-occurrences of these modes result in widespread precipitation deficits across the tropics by inducing anomalous subsidence, and reducing lower-level moisture convergence over the study regions. Our results emphasize the need for considering interactions within the larger climate system in characterizing compound drought risks rather than focusing on teleconnections from individual modes. Understanding the physical drivers and characteristics of compound droughts has important implications for predicting their occurrence and characterizing their impacts on interconnected societal systems.
This project aimed to improve long term global climate simulations by resolving and enhancing the representation of the processes involved in the cycling of freshwater through estuaries and coastal regions. This was a collaborative multi-institution project consisting of physical oceanographers, climate model developers, and computational scientists. It specifically targeted the DOE objectives of advancing simulation and predictive capability of climate models through improvements in resolution and physical process representation. The main computational objectives were: 1. To develop computationally efficient, but physically based, parameterizations of estuary and continental shelf mixing processes for use in an Earth System Model (CESM). 2. To develop a two-way nested regional modeling framework in order to dynamically downscale the climate response of particular coastal ocean regions and to upscale the impact of the regional coastal processes to the global climate in an Earth System Model (CESM). 3. To develop computational infrastructure to enhance the efficiency of data transfer between specific sources and destinations, i.e., a point-to-point communication capability, (used in objective 1) within POP, the ocean component of CESM.
This project aims to better understand causes of ongoing changes in the Arctic climate system, particularly as decreasing sea ice trends have been observed in recent decades and are expected to continue in the future. As part of the Sea Ice Prediction Network, a multi-agency effort to improve sea ice prediction products on seasonal-to-interannual time scales, our team is studying sensitivity of sea ice to a collection of physical process and feedback mechanism in the coupled climate system. During 2017 we completed a set of climate model simulations using the fully coupled ACME-HiLAT model. The simulations consisted of experiments in which cloud, sea ice, and air-ocean turbulent exchange parameters previously identified as important for driving output uncertainty in climate models were perturbed to account for parameter uncertainty in simulated climate variables. We conducted a sensitivity study to these parameters, which built upon a previous study we made for standalone simulations (Urrego-Blanco et al., 2016, 2017). Using the results from the ensemble of coupled simulations, we are examining robust relationships between climate variables that emerge across the experiments. We are also using causal discovery techniques to identify interaction pathways among climate variables which can help identify physical mechanisms and provide guidance in predictability studies. This work further builds on and leverages the large ensemble of standalone sea ice simulations produced in our previous w14_seaice project.
In future climate simulations there is a pronounced region of reduced warming in the subpolar gyre of the North Atlantic Ocean known as the North Atlantic warming hole (NAWH). This study investigates the impact of the North Atlantic warming hole on atmospheric circulation and midlatitude jets within the Community Earth System Model (CESM). Specifically, a series of large-ensemble atmospheric model experiments with prescribed sea surface temperature (SST) and sea ice are conducted, in which the warming hole is either filled or deepened. Two mechanisms through which the NAWH impacts the atmosphere are identified: a linear response characterized by a shallow atmospheric cooling and increase in sea level pressure shifted slightly downstream of the SST changes, and a transient eddy forced response whereby the enhanced SST gradient produced by the NAWH leads to increased transient eddy activity that propagates vertically and enhances the midlatitude jet. The relative contributions of these two mechanisms and the details of the response are strongly dependent on the season, time period, and warming hole strength. Our results indicate that the NAWH plays an important role in midlatitude atmospheric circulation changes in CESM future climate simulations.
From weather to climate timescales, weeks to decades, societal needs require detailed predictions of the atmosphere, ocean, and land-surface to facilitate informed planning. Central to this scientific endeavor is the ability to accurately initialize the earth system models (ESMs), since frequently the skill in forecasts, predictions, and projections depends on the quality of the initialization. However, initialization of high-resolution ESMs remains an outstanding challenge, particularly for the climate modeling community where initialization has received relatively less attention and where longer-term information is needed. This need, after a review of previous workshop outcomes, led the U.S. Department of Energy (DOE) and National Oceanic and Atmospheric Administration (NOAA) to host a workshop that focused exclusively on the opportunities that data assimilation (DA) and, specifically, coupled data assimilation (CDA) might provide to the ESM enterprise with particular attention toward initialization. The workshop was held April 9-10, 2018, at the Rockville Hilton in Rockville, Maryland, with 41 participants. It was composed of seven sessions, each with a distinct topical theme that related to the broad challenge of initializing high-resolution ESMs. The first session was used to frame the workshop challenge. Two of the sessions (2 and 7) focused the current state of the art of CDA within laboratory and operational centers at both the short-term, subseasonal and long-term, climate time scales. Three of the sessions (3, 4, and 5) addressed new algorithmic approaches to CDA at high resolution. Session 6 focused on the computational science solutions for the management of large input/output streams required to carry out the assimilation procedure.
The sensitivity of the Atlantic Meridional Overturning Circulation (AMOC) to high–latitude freshwater input is one of the key uncertainties in the climate system. Considering the importance of the AMOC for global heat transports, and the vulnerability of the Greenland Ice Sheet (GrIS) to global warming, assessing this sensitivity is critical for climate change projections. Here we present a unique set of computational experiments to investigate the adjustment of the AMOC to enhanced melt water from the GrIS under present–day conditions. For the first time, the response in a global, strongly–eddying ocean model is systematically compared to that of an ocean model typical of IPCC–class climate models. We find that the overall decline of the AMOC on decadal time scales is quantitatively similar (<10%) in the two configurations. Furthermore, the transient response is significantly different, as the AMOC decline and reduction in wintertime convection is markedly more gradual and persistent in the strongly–eddying configuration.
In response to quadrupled CO 2 , the Southern Ocean primarily uptakes excess heat around 60°S, which is then redistributed by the northward ocean heat transport (OHT) and mostly stored in the ocean or released back to the atmosphere around 45°S. However, the relative roles of mean ocean circulation and ocean circulation change in the uptake and redistribution of heat in the Southern Ocean remain controversial. Here, in this study, a set of climate model experiments embedded with a novel partial coupling technique are used to separate the roles of mean ocean circulation (passive component) and ocean circulation change (active component). For the ocean heat uptake (OHU) response, the mean ocean circulation and ocean circulation change are of equal importance. The OHT response south of 50°S is mainly determined by mean ocean circulation, while the ocean circulation change generates an anomalous southward OHT north of 50°S. A heat budget analysis finds that the divergence of passive OHT acts to balance the passive surface heat gain to the south of ~50°S, while the convergence of active OHT acts to balance the active surface heat loss to the north of ~50°S. Intriguingly, all the increase in ocean heat storage (OHS) is attributable to the passive component, with the ocean circulation change playing almost no role. In the Southern Ocean, both the active and the passive ocean heat transports are overcompensated by the reverse atmospheric heat transport via the Bjerknes compensation.