300 likes | 505 Views
Attributing direct radiative forcing to specific emissions using adjoint sensitivities. Daven K Henze, Drew T. Shindell , Robert J. D. Spurr. g -con. Radiative Forcing Transfer Functions. How to calculate the radiative forcing change for a given change in emissions?. IPCC, 2007.
E N D
Attributing direct radiative forcing to specific emissions using adjoint sensitivities Daven K Henze, Drew T. Shindell, Robert J. D. Spurr g-con
Radiative Forcing Transfer Functions How to calculate the radiativeforcing change for a given change in emissions? IPCC, 2007 Global contributions to aerosol direct RF
Radiative Forcing Transfer Functions How to calculate the radiative forcing change for a given change in emissions? Using transfer function T:
Radiative Forcing Transfer Functions How to calculate the radiative forcing change for a given change in emissions? Using transfer function T: Approximate Tusing adjoint: Calculated using GC adjoint(Henze et al., 2007 and LIDORT (Spurr, 2002)
Radiative Forcing Transfer Functions The % change in radiative forcing per change in BC emission: note: per change in any BC emission. This shows variation in efficiency of BC emissions forcing.
Radiative Forcing Transfer Functions The % change in radiative forcing per change in SO2 emission:
Applying to MFR 2030 – 2000 inventories BC, Total SO2, Total
Applying to MFR 2030 – 2000 inventories BC, Total SO2, Total
Applying toCLE 2030 – 2000 inventories BC, Total SO2, Total
Applying toCLE 2030 – 2000 inventories BC, RESALL SO2, RESALL
Applying toCLE 2030 – 2000 inventories BC, POWER SO2, POWER
Validation: BC CLE 2030 – 2000 perturbation Looks good
Validation: SO2 CLE 2030 – 2000 perturbation Looks OK, but adjoint-approach biased?
Validation: SO2 Check: does reducing perturbation reduce nonlinearity? CLE 2030 – 2000 perturbation 10% perturbations Yes. The adjoint code is accurate.
Validation: SO2 ESO2 E’’2030 E2000 E’2030 RF |ADJ| <|FD| |ADJ| >|FD| Can we anticipate bias?
Validation: SO2 Yes, bias can be anticipated. Also, overall ordering remains the same. E2030 < E2000 (Europe ) E2030 > E2000 (China, India) CLE 2030 – 2000 perturbation Conclusion: adjoint sensitivities provide a rapid means of exploring the effect of specific emissions changes on aerosol DRF.
The end Thanks to: Columbia Univ. Earth Institute Fellowship Drew Shindell, Rob Spurr, Nadine Unger, John Seinfeld NASA GSFC: NCCS NASA JPL: SCC
Radiative Forcing with GEOS-Chem Forward model Sensitivity calculation: GEOS-Chem GC Adj Henze et al., 2007 [SIA], [BC], RH, Ddry Mie Code derivate mode Mie Code Grainger et al., 2004 Following Martin et al., 2004, Drury et al., 2008 Radiative Transfer Code Radiative Transfer Code LIDORT (Spurr, 2002) Weighting functions TOA upward SW flux
Radiative Forcing with GEOS-Chem GEOS-Chem [SIA], RH Literature Ddry N, Dwet Mie Code Grainger et al., 2004 (tabulate as )
Validating Radiative Forcing Sensitivity Phase function coefficients for SIA(Dwet) Dmax ignore Dmin
Validating Radiative Forcing Sensitivity Mie results for extinction at discrete mode diameters:
Radiative Forcing (forward calculation) Chemical Transport Model Aerosol concentrations Mie Code optical properties
Radiative Forcing Sensitivity GEOS-Chem Adj Henze et al., 2007 LIDORT Spurr, 2002 Jacobian calculation [SIA]*, [BC]* Mie Derivative Grainger et al., 2004
Radiative Forcing (forward calculation) Chemical Transport Model Radiative Transfer Code Aerosol concentrations TOA upward SW fluxes Mie Code optical properties
Next Steps • Validate the transfer functions • Apply to various emissions perturbations of interest
Radiative Forcing Sensitivity 1 week 24 hr Radiative Transfer Code LIDORT (Spurr 2002) Jacobian calculation GEOS-Chem Adj [SIA]*, [BC]* Mie Derivative
Radiative Forcing with GEOS-Chem GEOS-Chem Radiative Transfer Code [BC] LIDORT (Spurr, 2002) [SIA], RH Ddry TOA upward SW flux N, Dwet Mie Code Grainger et al., 2004 Following Martin et al., 2004; Drury et al,. 2008 (external mixture)
Radiative Forcing with GEOS-Chem Forward model Sensitivity calculation GEOS-Chem GEOS-Chem Adj 1 wk Henze et al., 2007 [SIA], [BC], RH, Ddry [SIA]*, [BC]* Mie Code Grainger et al., 2004 Mie Code derivate mode Radiative Transfer Code Radiative Transfer Code LIDORT (Spurr, 2002) Following Martin et al., 2004, Drury et al., 2008 Weighting functions TOA upward SW flux 24 hr
Applying to CLE 2030 inventories note: this takes about 10 seconds