coolingFromNuclearWar: Global cooling from a nuclear war

In Ecdat: Data Sets for Econometrics

Description

Average surface temperature changes world wide and in the Northern Hemisphere 3 and 10 years after the injections of 5, 50 and 15 Tg (teragrams = millions of metric tons) of smoke into the upper troposphere, per Robock, Oman, and Stenchikov (2007).

These numbers are relative to the average for 1925-1975, which explains why the numbers are positive with smoke = 0.

Usage

data(coolingFromNuclearWar)

Format

A dataframe containing :

smoke

teragrams = millions of metric tons

dC3g, dC10g, dC3n, dC10n

average change in surface temperatature 3 and 10 years after injection of smoke into the upper troposphere globally (g) or in the Northern Hemisphers (n) in degrees Celsius.

Source

Alan Robock, Luke Oman, and Georgiy L. Stenchikov (2007) Nuclear winter revisited with a modern climate model and current nuclear arsenals: Still catastrophic consequences, Journal of Geophysical Research, 112

Examples

data(coolingFromNuclearWar)
matplot(coolingFromNuclearWar[, 'smoke'], 
    coolingFromNuclearWar[, 2:5], type='l')
(linFit <- lm(cbind(dC3g, dC10g, dC3n, dC10n)~smoke, 
      coolingFromNuclearWar))
      
# total change 
dC <- as.matrix(coolingFromNuclearWar[, 2:5] - 
        rep(unlist(coolingFromNuclearWar[1, -1]), e=4))
(linFit0 <- lm(dC~smoke, coolingFromNuclearWar))
summary(linFit0)

Ecdat documentation built on May 3, 2019, 1:24 p.m.