Search NASA⌕ Search

SEARCH · Search NASA

Results for “climate model outputs”

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 469 records · Page 26

A data-driven global soil heterotrophic respiration dataset and the drivers of its inter-annual variability

Soil heterotrophic respiration (SHR), one of the primary carbon fluxes from terrestrial ecosystems to the atmosphere, is important for carbon-climate feedbacks because of its sensitivity to available litter and soil carbon, climatic conditions, and nutrient availability. However, until recently limited SHR data were available, and most published global SHR estimates have either a short time span, coarse spatial resolution, or reply on overly-simple model formulations. To better understand and quantify the global distribution of SHR and its sensitivity to climate variability, we produced a new global SHR dataset using Random Forest algorithms, up-scaling 455 point data from the Global Soil Respiration Database (SRDB 4.0) with gridded fields of climatic, edaphic and productivity as explanatory variables. We estimated a global total SHR of 46.8 Pg C yr-1 over 1985-2013 (95% confidence interval: 38.6-56.3 Pg C yr-1), with a significant increasing trend of 0.03 Pg C yr-2 during this period. We found that the choice of soil moisture datasets contributes more to the difference among these data-driven SHR members rather than that of productivity, temperature and precipitation data sources. We also analyzed the influence of climatic variables on the inter-annual variability (IAV) of our SHR product. Water availability was the dominant driver of IAV at global scales, although the inferred sensitivity depends on the choice of the soil moisture gridded dataset. At the ecosystem scale, temperature strongly controls the IAV of SHR in tropical forests, while water availability dominates in extra-tropical forest and semi-arid regions. Our machine-learning gridded SHR dataset and outputs from process-based land surface models (TRENDYv6) show agreement for a strong association between water variability and SHR IAV at the global scale, but the two approaches lead to different temporal trend globally and different controlling variables for IAV at the ecosystem scale. Our study provides evidence for the pervasive and important role of water availability in driving SHR, indicating both a direct effect limiting decomposition rates and an indirect effect through the amount of fresh organic matter made available to SHR from productivity. In consideration of potential limitations and uncertainties remaining in our data-driven SHR datasets, we call for a more scientifically designed observation network for SHR, more observation data compilation, and increased use of deep learning methods making maximum use of observation data in hand. This will benefit process-based models, and improve our understanding of SHR response to future anomalous environmental conditions.

Yao, Yitong↗

Weather, Climate, and Society: New Demands on Science and Services

A new algorithm has been constructed to estimate the path length of lightning channels for the purpose of improving the model predictions of lightning NOx in both regional air quality and global chemistry/climate models. This algorithm was tested and applied to VHF signals detected by the North Alabama Lightning Mapping Array (NALMA). The accuracy of the algorithm was characterized by comparing algorithm output to the plots of individual discharges whose lengths were computed by hand. Several thousands of lightning flashes within 120 km of the NALMA network centroid were gathered from all four seasons, and were analyzed by the algorithm. The mean, standard deviation, and median statistics were obtained for all the flashes, the ground flashes, and the cloud flashes. Channel length distributions were also obtained for the different seasons.

Source record↗

SCOPE: Spectral cloud optical property estimation using real-time GOES-R longwave imagery

The output of ground-based, solar power generation systems is strongly dependent on cloud cover, which is the main contributor to solar power variability and uncertainty. Cloud optical properties are typically over-simplified in forecasting applications due to the lack of realtime, accurate estimates. In this work, we introduce a method, the Spectral Cloud Optical Property Estimation (SCOPE), for estimating cloud optical properties directly from high-resolution (5-min, 2 km) imagery from Geostationary Operational Environmental Satellite (GOES)-R, which is the newest generation of the GOES system. The SCOPE method couples a two-stream, spectrally resolved radiative model with the longwave GOES-R sensor output to simultaneously estimate the cloud optical depth, cloud top height, and cloud thickness during both day and night at 5-min intervals. The accuracy of SCOPE is evaluated using one year (2018) of downwelling longwave (DLW) radiation measurements from the Surface Radiation Budget Network, which consists of seven sites spread across climatically diverse regions of the contiguous United States. During daytime clear-sky periods, SCOPE predicts DLW within instrument uncertainty (10 W m -2 ) for four of the seven locations, with the remaining locations yielding errors of the order of 11.2, 17.7, and 20.2 W m -2 . For daytime cloudy-sky, daytime all-sky (clear or cloudy), and nighttime all-sky periods, SCOPE achieves root mean square error values of 23.0–34.5 W m -2 for all seven locations. These results, together with the low-latency of the method (1 s per sample), show that SCOPE provides a viable solution to real-time, accurate estimation of cloud optical properties for both day and night.

42 ENGINEERING↗

Atmospheric Teleconnection over Eurasia Induced by Aerosol Radiative Forcing During Boreal Spring

The direct effects of aerosols on global and regional climate during boreal spring are investigated based on simulations using the NASA Global Modeling and Assimilation Office (GMAO) finite-volume general circulation model (fvGCM) with Microphyics of clouds in Relaxed Arakawa Schubert Scheme (McRAS). The aerosol loading are prescribed from three-dimensional monthly distribution of tropospheric aerosols viz., sulfate, black carbon, organic carbon, soil dust, and sea salt from output of the Goddard Ozone Chemistry Aerosol Radiation and Transport model (GOCART). The aerosol extinction coefficient, single scattering albedo, and asymmetric factor are computed as wavelength-dependent radiative forcing in the radiative transfer scheme of the fvGCM, and as a function of the aerosol loading and ambient relative humidity. We find that anomalous atmospheric heat sources induced by absorbing aerosols (dust and black carbon) excites a planetary scale teleconnection pattern in sea level pressure, temperature and geopotential height spanning North Africa through Eurasia to the North Pacific. Surface cooling due to direct effects of aerosols is found in the vicinity and downstream of the aerosol source regions, i.e., South Asia, East Asia, and northern and western Africa. Additionally, atmospheric heating is found in regions with large loading of dust (over Northern Africa, and Middle East), and black carbon (over South-East Asia). Paradoxically, the most pronounced feature in aerosol-induced surface temperature is an east-west dipole anomaly with strong cooling over the Caspian Sea, and warming over central and northeastern Asia, where aerosol concentration are low. Analyses of circulation anomalies show that the dipole anomaly is a part of an atmospheric teleconnection driven by atmospheric heating anomalies induced by absorbing aerosols in the source regions, but the influence was conveyed globally through barotropic energy dispersion and sustained by feedback processes associated with the regional circulations.

Kim, Maeng-Ki↗

Reducing uncertainty of polar to midlatitude linkages using DOE’s E3SM in a coordinated model-experiment setting

This project brought DOE’s climate modeling effort with the Energy Exascale Earth System Model into the Polar Amplification Model Intercomparison Project (PAMIP), which is part of the sixth and latest Coupled Model Intercomparison Project, CMIP6. PAMIP examines the causes and consequences of polar amplification, when external forcing results in a larger temperature increase in high latitudes than the global average, in a coordinated set of model experiments with a common modeling protocol. Our teams from UC Irvine and the University of Toronto have designed, carried out, analyzed, and disseminated PAMIP output from the Energy Exascale Earth System Model (E3SM) and the Community Earth System Model (CESM). PAMIP’s ongoing stream of significant new results have advanced progress in the community’s understanding and led to new outstanding research questions that have motivated further work. PAMIP has led to improved consensus on the atmospheric response to sea ice loss. The important finding is that for a similar sea ice anomaly forcing, the simulated atmospheric response in the troposphere is remarkably consistent among the 16 models’ runs analyzed. The zonal-mean tropospheric response consists of a very robust equatorward shift of the westerly flow in mid-latitudes. However, while the multi-model mean response is robust, it has a weak amplitude relative to internal variability. We identified a weakness in the models (including E3SM) in terms of their response to sea-ice forcing that is related to eddy forcing (or nonlinear dynamical effects) at mid- to high latitudes. In fact, E3SM is an outlier in terms of the models participating, and in that sense, it turned out to be a vital participant model. We found that reductions in energy transport due sea-ice loss and involving dry air only are compensated by increases in moist energy transport from warmer sea surface temperature in midlatitudes. This leads to a large spread in energy transport into the Arctic and is a potential source of spread in Arctic amplification. We identified an important role that climate modes, including tropical modes of variability (El Nino and the Southern Oscillation (ENSO); the Quasi-biennial Oscillation (QBO)) play in the response to sea-ice anomalies, including in ocean coupled experiments. Similarly, we identified and quantified the contribution of sea-ice thickness to the atmospheric response compared to the response to sea-ice extent only. We found that it is important to run large ensembles and even with an ensemble size of 100 simulations the response is largely influenced by internal variability. We demonstrated convincingly that Ural blocking, not sea-ice loss, provides the weakening of the stratospheric polar vortex in fall/early winter and a negative phase of the North Atlantic Oscillation that can last for up to two months. However, sea-ice anomalies can influence the background flow so that the response to Ural blocking is more persistent under low sea-ice conditions in the Barents/Kara Sea than high sea ice Atmospheric model hierarchies that progressively add individual processes have a long history in providing dynamical insight for modeling the atmosphere. Similarly, coupled model hierarchies that progressively add individual ocean processes can provide insights into the workings of the coupled climate system, however such hierarchies have not been available except for a non-dynamical slab ocean model. Because of the missing processes, surface flux corrections must be added to produce a target climate. In this project, we managed to overcome this problem and develop a globally coupled ocean model hierarchy in CESM that can turn on and off the processes of mixed-layer entrainment and Ekman flow. We used the hierarchy to study the impact of Arctic sea-ice loss on the climate system. We find that the effect of mixed-layer entrainment on ocean heat uptake influences the atmospheric circulation by shifting the latitudinal positions of the mid-latitude westerly jet and the Intertropical Convergence Zone (ITCZ). In quadrupled CO 2 experiments, we studied how air-sea coupling affects the response of tropical rainfall under global warming. In order to identify the importance of individual ocean processes, we used the hierarchy of ocean models to separate the effects of seasonal mixed-layer entrainment, wind-driven Ekman flows, and frictional flows. We showed that including Ekman and frictional flows allows our simulation to produce the Pacific Ocean's enhanced equatorial warming pattern and equatorward ITCZ contraction noted in previous climate simulations. We also showed that the frictional flow, which has yet to receive much attention, is as important as the Ekman flow in generating equatorial heat convergence.

54 ENVIRONMENTAL SCIENCES↗

The Brewer–Dobson Circulation in CMIP6

The Brewer–Dobson circulation (BDC) is a key feature of the stratosphere that models need to accurately represent in order to simulate surface climate variability and change adequately. For the first time, the Climate Model Intercomparison Project includes in its phase 6 (CMIP6) a set of diagnostics that allow for careful evaluation of the BDC. Here, the BDC is evaluated against observations and reanalyses using historical simulations. CMIP6 results confirm the well-known inconsistency in the sign of BDC trends between observations and models in the middle and upper stratosphere. Nevertheless, the large uncertainty in the observational trend estimates opens the door to compatibility. In particular, when accounting for the limited sampling of the observations, model and observational trend error bars overlap in 40 % of the simulations with available output. The increasing CO2 simulations feature an acceleration of the BDC but reveal a large spread in the middle-to-upper stratospheric trends, possibly related to the parameterized gravity wave forcing. The very close connection between the shallow branch of the residual circulation and surface temperature is highlighted, which is absent in the deep branch. The trends in mean age of air are shown to be more robust throughout the stratosphere than those in the residual circulation.

Brewer–Dobson circulation↗

How Earth System Models Can Inform Key Dimensions of Marine Food Security in the Alaskan Arctic

The Arctic is home to several groups of Indigenous Peoples, each with distinct ways of interacting with their environment and ways of life. Arctic, Indigenous Peoples’ food sovereignty is tightly linked with food security. Subsistence harvesting activities provide nutritious and culturally vital foods for Alaska Native households and communities. Climate change is causing rapid and more unpredictable shifts in environmental conditions that impact three of the key aspects of food security, availability, stability, and accessibility. While communities monitor the abundance and health of food webs and environments as part of subsistence harvest practices, anticipating major transformations and changes in these systems is challenging. We explored the potential of Earth System Model output in helping anticipate or project physical or ecosystem changes relevant to Alaska Indigenous peoples’ food security needs. Through examples of model products, that provide measures of accessibility and availability of marine resources, we show that modern models, such as the Energy Exascale Earth System Model presented here, can provide estimates of a broad suite of variables relevant to food security. We investigate how Earth System Model output could contribute to exploring questions related to aspects of Arctic food security such as accessibility and availability and highlight present model shortcomings that, if addressed, would move Earth System Models closer to being a useful tool for understanding environmentally driven changes to the availability and accessibility of harvestable food resources. Our example model-derived food security indicators illustrate how Earth System Model output could be combined with relevant, non-model, information sources; These model products are meant only as a starting point and a tool for engaging community members and to present, in an accessible way, the model’s potential utility, or current lack thereof, to rights holders and stakeholders concerned about food security. We are hopeful that with example products in hand, additional model development efforts will have a higher likelihood of success in achieving an iterative discussion with stakeholders regarding feasible and desired products.

Gibson, Georgina A.↗

Time-Variable Gravity and Ocean Connections

Any large mass transport in the Earth system produces changes in the gravity field. Low harmonic degree components of such gravity variations have been observed by the satellite-laser-ranging (SLR) technique for the past quarter century, particularly in 52, the Earth's dynamic oblateness. 52 undergoes a slight decrease due primarily to the post-glacial rebound of the mantle, but large interannual anomalies have been observed, notably during 1998-2002. Intriguing evidences for the cause of the latter have been found in the extratropical Pacific basins, especially related to the Pacific Decadal Oscillation, and perhaps in related land hydrology. We will examine the latest results based on ocean altimetry, sea-surface temperature, and ocean and hydrology model outputs. Without firm estimates for the steric effects (which have no gravity signal), we will point out possible underestimation of OGCMs with respect to temporal variabilities. Besides J2, SLR also derived time series for other low-degree gravity components. While the formal uncertainty of these terms is significantly higher, some of these series have significant signal that show correlation to various climatic signals. For example, there is a significant correlation of the sectoral S2,2 with the Southern Oscillation Index signifying the influence of El Nino/La Nina. Cases such as these demonstrate the utility of assessing the mass component of climate variations, and anticipate the utility of GRACE-type space gravity observations with much higher spatial resolution.

Chao, Benjamin F.↗

Optimized Umkehr Profile Algorithm for Ozone Trend Analyses

The long-term record of Umkehr measurements from four NOAA Dobson spectrophotometers was reprocessed after updates to the instrument calibration procedures. In addition, a new data quality-control tool was developed for the Dobson automation software (WinDobson). This paper presents a comparison of Dobson Umkehr ozone profiles from NOAA ozone network stations (Boulder, OHP, MLO, Lauder) against several satellite records, including Aura Microwave Limb Sounder (MLS; ver. 4.2), and combined SBUV and OMPS records (NASA AGG and NOAA COH). A subset of satellite data is selected to match Dobson Umkehr observations at each station spatially (distance less than 200 km) and temporally (within 24 hours). Umkehr Averaging Kernels (AKs) are applied to vertically smooth all overpass satellite profiles prior to comparisons. The station Umkehr record consists of several instrumental records, which have different optical characterizations, and thus instrument-specific stray light contributes to the data processing errors and creates step changes in the record. This work evaluates the overall quality of Umkehr long-term measurements at NOAA ground-based stations and assesses the impact of the instrumental changes on the stability of the Umkehr ozone profile record. This paper describes a method designed to correct biases and discontinuities in the retrieved Umkehr profile that originate from the Dobson calibration process, repair, or optical realignment of the instrument. The M2GMI and GMI CTM ozone profile model output matched to station location and date of observation is used to evaluate instrumental step changes in the Umkehr record. Homogenization of the Umkehr record and discussion of the apparent stray light error in retrieved ozone profiles are the focus of this paper. Homogenization of ground-based records is of great importance for studies of long-term ozone trends and climate change.

Umkehr↗

Update and Status of the Mars Climate Modeling Center at NASA Ames Research Center

NASA Ames Research Center has a long-established history in the study and numerical modeling of Mars’ atmosphere and climate. Such studies began with the late Jim Pollack and have continued under the leaderships of Bob Haberle, Jeff Hollingsworth, and now Melinda Kahre. Our group has grown to more than 10 members, including civil servants, research scientists, postdocs, and students. While we still prioritize science and the model development that supports it, we have recently enhanced efforts to make our codes and output publicly available and boost community engagement through the hosting of modeling tutorials, etc. Here our goal is to present the status of the Ames Mars GCM, our ongoing science projects, the tools we have recently made publicly available, and our plans for releases and continued community engagement.

M. A. Kahre↗

The stratopause semiannual oscillation in the NCAR Community Climate Model

The middle atmospheric version of the NCAR Community Climate Model (CCM2) has been used to study the development of the equatorial semiannual oscillation (SAO) in the stratosphere. The model domain extends from the ground to about 80 km, with a vertical resolution of 1 km. Transport of nitrous oxide (N2O) with simplified photochemistry is included in the calculation to illustrate the influence of tropical circulations on the distribution of trace species. Diagnosis of model output reveals two distinct phases in the evolution of the zonal mean state on the equator. In early December, a strong and broad easterly jet appears near the stratopause in connection with a midlatitude wave event (sudden stratospheric warming) that reverses the winter westerlies of the Northern Hemisphere throughout the upper stratosphere. When the wave forcing dies out, the radiative drive allows the westerlies to recover at midlatitudes, while easterlies persist in the tropics. The resulting strong meridional gradient of the zonal mean wind provides favorable conditions for the development of inertial instability at lower latitudes. The meridional circulation associated with the instability shapes the 'nose' of the easterly jet, reducing the extension of the unstable region. In equinoctial conditions, a jet of westerlies appears in the lower equatorial mesosphere and descends to lower altitudes; positive accelerations associated with the descending westerlies are due primarily to Kelvin waves. The descent of the westerly jet does not reproduce well the observed behavior of the SAO westerly phase, either in amplitude or in the extent of downward propagation. As a consequence, the model does not simulate the 'double peak' observed in the tropical distribution of N2O. Comparison of wave amplitudes in the model with those derived from satellite observations shows that the calculated amplitudes are larger than observed in the upper stratosphere. It follows that inadequate Kelvin wave forcing is not the cause of the weak westerly phase in the model, and that some other mechanism must be responsible for the generation of the strong westerly phase observed.

Sassi, Fabrizio↗

Groundwater and Terrestrial Water Storage

Groundwater, soil moisture, surface water, snow, and ice are the components of terrestrial water storage (TWS). On multi-annual timescales groundwater typically controls TWS variability, except in permanently frozen regions (Li et al. 2015). Even on an annual basis TWS changes are a reasonable proxy for groundwater storage changes, the latter being insufficiently observed in most of the world. From 2002 to 2017 the Gravity Recovery and Climate Experiment (GRACE; Tapley et al. 2004) and since 2018 the GRACE Follow On (GRACE-FO) satellite missions have enabled estimation of TWS anomalies (departures from the long-term mean) based on precise observations of variations in Earth’s gravity field. To bridge the 11-month gap between GRACE and GRACE-FO we make use of output from a land surface model that assimilates data from both missions (Li et al. 2019)

Matthew Rodell↗

Canopy spectral reflectance along the NGEE-Arctic tram transect, Utqiagvik (Barrow), Alaska, 2015-2017

Full-spectrum (350-2500 nm) canopy spectral reflectance measured in July 2015, 2016 and 2017 along the NGEE-Arctic tram located on the Barrow Environmental Observatory (BEO), Barrow, Alaska. Measurements were made using a SVC HR1024i spectroradiometer mounted on a manually operated instrument cart. Spectra were collected at 137 stop points at 0.5 m intervals along the tram transect, with each measurement calibrated to variable illumination conditions using a Spectralon white plate. Data were collected in July, 2015-2016, and in June, July and September in 2017. Processed spectral data are presented in .csv and .xlsx files for each year and the complete SVC instrument output (.sig) are also included in .zip files. The data package includes some photographs showing the tram transect and spectroradiometer setup and a written protocol with details of equipment used and methods.The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research.The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska.Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗

NASA Downscaling Project: Final Report

A team of researchers from NASA Ames Research Center, Goddard Space Flight Center, the Jet Propulsion Laboratory, and Marshall Space Flight Center, along with university partners at UCLA, conducted an investigation to explore whether downscaling coarse resolution global climate model (GCM) predictions might provide valid insights into the regional impacts sought by decision makers. Since the computational cost of running global models at high spatial resolution for any useful climate scale period is prohibitive, the hope for downscaling is that a coarse resolution GCM provides sufficiently accurate synoptic scale information for a regional climate model (RCM) to accurately develop fine scale features that represent the regional impacts of a changing climate. As a proxy for a prognostic climate forecast model, and so that ground truth in the form of satellite and in-situ observations could be used for evaluation, the MERRA and MERRA - 2 reanalyses were used to drive the NU - WRF regional climate model and a GEOS - 5 replay. This was performed at various resolutions that were at factors of 2 to 10 higher than the reanalysis forcing. A number of experiments were conducted that varied resolution, model parameterizations, and intermediate scale nudging, for simulations over the continental US during the period from 2000 - 2010. The results of these experiments were compared to observational datasets to evaluate the output.

dynamical downsizing↗

NASA Downscaling Project

A team of researchers from NASA Ames Research Center, Goddard Space Flight Center, the Jet Propulsion Laboratory, and Marshall Space Flight Center, along with university partners at UCLA, conducted an investigation to explore whether downscaling coarse resolution global climate model (GCM) predictions might provide valid insights into the regional impacts sought by decision makers. Since the computational cost of running global models at high spatial resolution for any useful climate scale period is prohibitive, the hope for downscaling is that a coarse resolution GCM provides sufficiently accurate synoptic scale information for a regional climate model (RCM) to accurately develop fine scale features that represent the regional impacts of a changing climate. As a proxy for a prognostic climate forecast model, and so that ground truth in the form of satellite and in-situ observations could be used for evaluation, the MERRA and MERRA-2 reanalyses were used to drive the NU-WRF regional climate model and a GEOS-5 replay. This was performed at various resolutions that were at factors of 2 to 10 higher than the reanalysis forcing. A number of experiments were conducted that varied resolution, model parameterizations, and intermediate scale nudging, for simulations over the continental US during the period from 2000-2010. The results of these experiments were compared to observational datasets to evaluate the output.

Ferraro, Robert↗

Changing Characteristics of Tropical Extreme Precipitation–Cloud Regimes in Warmer Climates

In this study, we investigated the changing characteristics of climatic scale (monthly) tropical extreme precipitation in warming climates using the Energy Exascale Earth System Model (E3SM). The results are from Atmospheric Model Intercomparison Project (AMIP)-type simulations driven by (a) a control experiment with the present-day sea surface temperature (SST) and CO 2 concentration, (b) P4K, the same as in (a) but with a uniform increase of 4K in the SST globally, and (c) the same as in (a), but with an imposed SST and CO 2 concentration from the outputs of the coupled E3SM forced by a 4xCO 2 concentration. We found that as the surface warmed under P4K and 4xCO 2 , both convective and stratiform rain increased. Importantly, there was an increasing fractional contribution of stratiform rain as a function of the precipitation intensity, with the most extreme but rare events occurring preferentially over land more than the ocean, and more so under 4xCO 2 than P4K. Extreme precipitation was facilitated by increased precipitation efficiency, reflecting accelerated rates of recycling of precipitation cloud water (both liquid and ice phases) in regions with colder anvil cloud tops. Changes in the vertical profiles of clouds, condensation heating, and vertical motions indicate increasing precipitation–cloud–circulation organization from the control and P4K to 4xCO 2 . The results suggest that large-scale ocean warming, that is, P4K, was the primary cause contributing to an organization structure resembling the well-known mesoscale convective system (MCS), with increased extreme precipitation on shorter (hourly to daily) time scales. Additional 4xCO 2 atmospheric radiative heating and dynamically consistent anomalous SST further amplified the MCS organization under P4K. Analyses of the surface moist static energy distribution show that increases in the surface moisture (temperature) under P4K and 4xCO 2 was the key driver leading to enhanced convective instability over tropical ocean (land). However, a fast and large increase in the land surface temperature and lack of available local moisture resulted in a strong reduction in the land surface relative humidity, reflecting severe drying and enhanced convective inhibition (CIN). It is argued that very extreme and rare “record-breaking” precipitation events found over land under P4K, and more so under 4xCO 2 , are likely due to the delayed onset of deep convection, that is, the longer the suppression of deep convection by CIN, the more severe the extreme precipitation when it eventually occurs, due to the release of a large amount of stored surplus convective available potential energy in the lower troposphere during prolonged CIN.

54 ENVIRONMENTAL SCIENCES↗

A bespoke model of Arctic river basins based on hillslope delineation: Model Archive

This dataset is a model archive of the paper A bespoke model of Arctic river basins based on hillslope delineation (in prep), which introduces a watershed decomposition and parameterization method for large scale permafrost hydrology simulation. With this dataset, this study aims to address the research question: whether a computationally efficient hillslope-based modeling framework can reliably simulate discharge at Arctic river-basin scales. This dataset contains model input and output data for five modeling scenarios at a study site located in the Sagavanirktok River basin. The five modeling scenarios include three modeling cases under temperate conditions using full 3D, decomposed 3D, and decomposed 2D modeling strategies; and two modeling cases under actual Arctic conditions with permafrost using full 3D and decomposed 2D modeling strategies. Simulations were performed using the Advanced Terrestrial Simulator (ATS, v1.6 for three temperate scenarios and v1.5 for two Arctic scenarios), a physics-rich integrated surface–subsurface hydrologic model with cryo-hydrology features. For the three temperate models, simulations were conducted for the period of 10/01/1993 - 09/30/2002; and for the two Arctic models, simulations were conducted for the period of 01/01/1994 - 12/31/2002. To facilitate reproducibility of simulations, all datasets are organized hierarchically. The dataset contains: (1) Mesh files (.exo) for full 3D model, decomposed 3D models, and decomposed 2D models, located in huc/190604020802_gauge15906000/mesh/. Mesh files can be visualized through Paraview or read by Python. (2) Climate forcings (.h5) for full 3D model and decomposed 3D/2D models are located in huc/190604020802_gauge15906000/daymet_onePiece/, and huc/190604020802_gauge15906000/vp_pr_revised_daymet_1980_2006_with_wind/ separately. Accessible by Python. (3) Raw measured gage discharge (.csv) from USGS, located in huc/190604020802_gauge15906000/gaged_basin15906000_discharge_usgs/. Accessible by Python. (4) Delineated subdomain raster (.tif) and shape files (.shp), and the final parameterized results (.npy) for decomposed models, located in huc/190604020802_gauge15906000/data_preprocessed-meshing. Accessible by Python. (5) Temperate models are located in nonpermaf_huc190604020802_gauge15906000/, which includes three cases: decomposed 2D models (inside model_0*-hillslope_*), decomposed 3D models (inside model_1*-subcatchment_*), and full 3D model (inside model_2*-onepiece_*). Two step spin-up results (checkpoint_final.h5) are located in model_*1-*_spinup_steadystate and model_*2-*_spinup_cycle, separately, which are used to initialize real transient models. The input files (.xml) and output results (.dat) of the real transient models are located in model_*3-*_transient/. Especially, for two example hillslope models (ID=-11 and 11), additional h5py files are included in model_03-hillslope_transient/hillslope-11/, model_03-hillslope_transient/hillslope11, model_13-subcatchment_transient/subcatchment-11/, model_13-subcatchment_transient/subcatchment/11, respectively, which are used to plot the saturation figure (Figure 5) in the manuscript. Accessible by Python. (6) Arctic models are located in huc190604020802_gauge15906000/, which includes two cases: decomposed 2D models (inside model_04-hillslope_transient), and full 3D model (inside model_05-onepiece_transient_mannp1_ra). Three step spin-up results (checkpoint_final.h5) are located in model_01-column_freezeup/, model_02-column_spinup/, model_03-hillslope_spinup/, respectively, which are used to initialize real 2D transient hillslope models. The input files (.xml) and output results (.dat) of transient 2D hillslope models are located in model_04-hillslope_transient/. The input files (.xml) and output results (.dat) of the full 3D transient model is located in model_05-onepiece_transient_mannp1_ra/. The full 3D transient model is initialized by model_02-column_spinup/. Accessible by Python. (7) The MOSART routed discharge results (.csv) under Arctic conditions is located in huc190604020802_gauge15906000/MOSART/. Accessible by Python. (8) All Python codes (.py) used to parameterize full 3D model to decomposed 2D models are located in script/. These codes fit with watershed workflow (a watershed delineation tool) v1.4 under the branch gaob/v1.4 from https://github.com/gaobhub/watershed-workflow.git.

EARTH SCIENCE > CRYOSPHERE↗

Carbon dioxide emissions from Deccan volcanism and a K/T boundary greenhouse effect

A greenhouse warming caused by increased emissions of carbon dioxide from the Deccan Traps volcanism has been suggested as the cause of the terminal Cretaceous extinctions on land and in the sea. Total eruptive and noneruptive CO2 output by the Deccan eruptions (from 6 to 20 x 10 to the 16th moles) over a period of several hundred thousand years is estimated based on best estimates of the CO2 weight fraction of the original basalts and basaltic melts, the fraction of CO2 degassed, and the volume of the Deccan Traps eruptions. Results of a model designed to estimate the effects of increased CO2 on climate and ocean chemistry suggest that increases in atmospheric pCO2 due to Deccan Traps CO2 emissions would have been less than 75 ppm, leading to a predicted global warming of less than 1 C over several hundred thousand years. It is concluded that the direct climate effects of CO2 emissions from the Deccan eruptions would have been too weak to be an important factor in the end-Cretaceous mass extinctions.

Caldeira, Ken↗