% $Header$ % Purpose: Proposal for NSF Arctic Natural Sciences (ANS) Dirty Snow project % URLs: % http://dust.ess.uci.edu/prp/prp_ans-1.0/prp_ans.pdf % http://dust.ess.uci.edu/prp/prp_ans-1.0/prp_ans_fll.pdf % http://www.nsf.gov/awardsearch/showAward.do?AwardNumber=0431203 % CVS: cvs -d :ext:esmf.ess.uci.edu:/u/zender/cvs co -kk prp_ans % Usage (see also end of file): % cd ~/prp_ans-1.0;make -W prp_ans.tex prp_ans.dvi prp_ans.ps prp_ans.pdf prp_ans.txt;cd - % 2005 NSF Arctic Research Opportunities (ARO): % NSF 05-618 % Office of Polar Programs (OPP) % Arctic Sciences Section (ASS) % Arctic Natural Sciences (ANS) Program % Program Manager, ANS: Jane V. Dionne ``die-on'', (703) 292-7427, jdionne@nsf.gov % Deadlines: Full Proposal 20051216 % NSF FastLane Temporary Proposal #635594 PIN czen % NSF FastLane Proposal #0612954 (refer to as ANS-0612954) % Total 3-year budget request: $403289 % Project duration: 20060701--20090630 % Annual progress report deadlines: % Procurement control number (PCN): % UCI account number: 9-number-fund-sub-object = 9-123456-12345-1-1234 = 9-445925-12345-1-1234 % Physical Sciences budget code: 1234 % 2006 NSF OPP ARO ANS Round2: % Due November 10, 2006 % Usage: See end of file % fxm: Upload new bio_nsf \documentclass[12pt]{article} % Standard packages \usepackage{ifpdf} % Define \ifpdf \ifpdf % We are running PDFLaTeX \usepackage[pdftex]{graphicx} % Defines \includegraphics* \pdfcompresslevel=9 \usepackage{thumbpdf} % Generate thumbnails \usepackage{epstopdf} % Convert .eps, if found, to .pdf when required \else % We are not running PDFLaTeX \usepackage{graphicx} % Defines \includegraphics* \fi % endif PDFLaTeX \usepackage{amsmath} % \subequations, \eqref, \align \usepackage{array} % Table and array extensions, e.g., column formatting \usepackage{datetime} % \xxivtime, \ordinal \usepackage{longtable} % Multi-page tables, e.g., acronyms and symbols \usepackage{makeidx} % Index keyword processor: \printindex and \see \usepackage{mdwlist} % Compact list formats \itemize*, \enumerate* \usepackage{natbib} % \cite commands from aguplus \usepackage{times} % Postscript Times-Roman font KoD99 p. 375 \usepackage{tocbibind} % Add Bibliography and Index to Table of Contents \usepackage{url} % Typeset URLs and e-mail addresses % hyperref is last package since it redefines other packages' commands % hyperref options, assumed true unless =false is specified: % backref List citing sections after bibliography entries % baseurl Make all URLs in document relative to this % bookmarksopen Unknown % breaklinks Wrap links onto newlines % colorlinks Use colored text for links, not boxes % hyperindex Link index to text % plainpages=false Suppress warnings caused by duplicate page numbers % pdftex Conform to pdftex conventions % Colors used when colorlinks=true: % linkcolor Color for normal internal links % anchorcolor Color for anchor text % citecolor Color for bibliographic citations in text % filecolor Color for URLs which open local files % menucolor Color for Acrobat menu items % pagecolor Color for links to other pages % urlcolor Color for linked URLs \ifpdf % We are running PDFLaTeX \usepackage[backref,breaklinks,colorlinks,citecolor=blue,linkcolor=blue,urlcolor=blue,hyperindex,plainpages=false,pdftex]{hyperref} % Hyper-references \pdfcompresslevel=9 \else % We are not running PDFLaTeX \usepackage[backref=false,breaklinks,colorlinks=false,hyperindex,plainpages=false]{hyperref} % Hyper-references \fi % endif PDFLaTeX % Personal packages \usepackage{csz} % Library of personal definitions \usepackage{abc} % Alphabet as three letter macros \usepackage{dmn} % Dimensional units \usepackage{chm} % Commands generic to chemistry \usepackage{dyn} % Commands generic to fluid dynamics \usepackage{aer} % Commands specific to aerosol physics \usepackage{psd} % Particle size distributions \usepackage{rt} % Commands specific to radiative transfer \usepackage{hyp} % Hyphenation exception list \usepackage{jrn_agu} % AGU-sanctioned journal title abbreviations % Commands which must be executed in preamble % Commands which must be executed in preamble \newcounter{idxpg} % [nbr] Page number of Index %\allowdisplaybreaks[1] % AMSTeX permits eqnarray to wrap pages [1--4] %\listfiles % Print required files, versions \makeglossary % Glossary described on KoD95 p. 221 \makeindex % Index described on KoD95 p. 220 % Commands specific to this file \makeatletter \renewcommand{\fnum@table}{\textbf{\tablename~\thetable}} % Boldface ``Table: #'' \renewcommand{\fnum@figure}{\textbf{\figurename~\thefigure}} % Boldface ``Figure: #'' \makeatother % 1. Primary commands % 2. Derived commands % 3. Doubly-derived commands % Margins \oddsidemargin = 0.0in \evensidemargin = 0.0in \topmargin = -0.4in \textheight = 9in \textwidth = 6.5in \columnsep = 0.25in \headheight = 12pt \headsep = 12pt \footskip = 24pt % Vertical space % 20000717: None of these commands seem to have any effect here \setlength{\textfloatsep}{12pt} % Space between last top float or first bottom float and the text. \setlength{\intextsep}{12pt} % Space left on top and bottom of an in-text float. % \dbltextfloatsep is \textfloatsep for 2 column output. % \dblfloatsep is \floatsep for 2 column output. \setlength{\abovecaptionskip}{12pt} % Space above caption \setlength{\belowcaptionskip}{12pt} % Space below caption % Float placement % NB: Placement of figures is very sensitive to \textfraction \renewcommand\textfraction{0.0} % Minimum fraction of page that is text \setcounter{totalnumber}{73} % Maximum number of floats per page \setcounter{topnumber}{73} % Maximum number of floats at top of page \setcounter{dbltopnumber}{73} % Maximum number of floats at top of two-column page \setcounter{bottomnumber}{73} % Maximum number of floats at bottom of page \renewcommand\topfraction{1.0} % Maximum fraction of top of page occupied by floats \renewcommand\dbltopfraction{1.0} % Maximum fraction of top of two-column page occupied by floats \renewcommand\bottomfraction{1.0} % Maximum fraction of bottom of page occupied by floats \renewcommand\floatpagefraction{1.0} % Fraction of float page filled with floats \renewcommand\dblfloatpagefraction{1.0} % Fraction of double column float page filled with floats \begin{document} % Arctic Natural Sciences (ANS) \def\prpttl{Dirty Snow on Land, Glaciers, and Sea-Ice: Understanding Arctic Absorbing Aerosol Forcing and Feedbacks\\} % Cover page \pagenumbering{roman} \setcounter{page}{1} \pagestyle{headings} \thispagestyle{empty} {\noindent% On the Web at \url{http://dust.ess.uci.edu/prp/prp_ans-1.0/prp_ans.pdf}\\ NSF Arctic Natural Sciences (ANS) Proposal \hfill Submitted: December~22, 2005\\ Last modified: \today, \xxivtime \hfill Next Round Due: November~10, 2006} \begin{center} \textbf{\Large\prpttl} \bigskip Dr. Charles S. Zender \hfill \\ Department of Earth System Science \hfill \\ University of California, Irvine \hfill \\ \end{center} \vskip 0.5 cm \noindent\textbf{News/Preface:} This is NSF proposal 0612954. \noindent\textbf{Information for potential collaborators:} This NSF proposal responds to the 2005 NSF Arctic Research Opportunities (ARO) announcement, NSF 05-618. The proposale was submitted to the Arctic Natural Sciences (ANS) Program of the Arctic Sciences Section (ASS) in the Office of Polar Programs (OPP). The cognizant Program Manager is Jane V. Dionne \url{jdionne@nsf.gov}, (703)~292-7427. Suggestions for next round: \begin{enumerate} \item Don't write it in an all-nighter three days before Christmas \item Make satellite section less cheesey for NASA---e.g., include AMSR-E comparison \item Contact Dorothy Koch for NASA-GISS snow physics collaboration? \end{enumerate} \tableofcontents \listoffigures %\listoftables \clearpage %\onecolumn %\tableofcontents %\listoftables %\clearpage \pagenumbering{arabic} \setcounter{page}{0} \thispagestyle{empty} %\markleft{Dirty Snow} %\markright{} %\begin{center} %\textbf{\large\prpttl} %\end{center} \begin{center} \textbf{\Large\prpttl} \bigskip Dr. Charles S. Zender \hfill \\ Department of Earth System Science, University of California, Irvine \hfill \\ \end{center} \noindent{\large{\textbf{Project Summary.}}}\label{sxn:smr} %\enlargethispage*{0.5in} Surface and atmospheric concentrations of black carbon (BC), i.e, absorbing carbonaceous aerosol, are highly variable and slowly increasing in the Arctic. Current understanding Arctic BC climate impacts derives mostly from studies which focus on BC direct atmospheric radiative forcing and which neglect or drastically simplify surface BC interactions. However, the prevalence of bright surfaces (snow, glaciers, sea-ice, and clouds) make the Arctic uniquely susceptible to radiatively induced effects of surface BC and dust such as ice-albedo feedback amplification. Dirty snow feedbacks change throughout the aerosol lifecycle in the complex Arctic environment of snowfall, snowpack aging, snow-melt, drainage, and analogous sea-ice processes. Our goal is to assess absorbing aerosol interactions in the coupled Arctic climate system using models which represent the complex surface lifecycles of Arctic snow, BC, and dust, and which have been evaluated against satellite, in-situ, and laboratory measurements. The project builds upon, extends, and applies our existing, SNow, ICe, and Aerosol Radiative model, SNICAR. SNICAR has already helped us show that snowpack aging processes can trigger ice-albedo feedbacks which amplify direct surface forcing by an order of magnitude in mid-latitude regions. We hypothesize that Arctic BC and dust cause similar, or possibly greater, ice-albedo feedback amplification over glaciers and sea-ice. However, current coupled model systems cannot test this hypothesis. Hence, our project integrates BC, dust and snow lifecycles within SNICAR, and places SNICAR in existing coupled snow, ice sheet, and sea ice models. NASA MODIS, MISR, and AMSR-E retrievals will help us constrain model parameters, evaluate predictions, and interpret the regional and seasonal behavior of snowpack processes. We will support the IPY POLARCAT project by modeling Arctic soot concentration and snow reflectance following boreal fires. Snowpack specific surface area, albedo, crystal density, BC concentration, BC melt scavenging, and aging processes will (continue to) be evaluated against Arctic \textit{in~situ} measurements. \textbf{Scientific Merit:} Our studies of Arctic snow-ice-aerosol processes will improve understanding of ice-albedo feedbacks and polar climate amplification. Key scientific questions we will address include: 1)~Potential for Arctic BC to disrupt summertime sea-ice formation; 2)~Role of Boreal fire variability in Arctic reflectance, particularly Greenland; 3)~Relative roles of Arctic soot and dust in inducing surface melt, heating, and darkening. These questions will be addressed in the context of pre-industrial, present day, and next century timescales. \textbf{Broader Impacts:} Snow pervades the Arctic, and our improved snowpack model can and probably will improve understanding of Arctic hydrology. SNICAR will be freely available to run in both off-line and Community Climate System Model modes. We anticipate it will contribute to Arctic research areas including glacier mass balance, basin and catchment hydrology, sea-ice lifecycle, paleoclimate sensitivity, and snow chemistry. This project trains one graduate student in Arctic aerosol-climate interactions. PI Zender will incorporate this Arctic climate change research into a K--12 teacher training program and into presentations for lifelong learning students. \clearpage \begin{center} \textbf{\Large \prpttl} \end{center} \thispagestyle{empty} \section{Introduction}\label{sxn:ntr} Surface and atmospheric concentrations of black carbon (BC), i.e, absorbing carbonaceous aerosol, are highly variable and slowly increasing in the Arctic \cite[][]{PAA01}. Current understanding of Arctic BC climate impacts derives mostly from studies which focus on BC direct atmospheric radiative forcing and which neglect or drastically simplify surface BC interactions. However, bright surfaces (snow, glaciers, sea-ice, and clouds) make the Arctic uniquely susceptible to radiatively induced effects of surface BC and dust such as ice-albedo feedback amplification. Such feedbacks make dirty snow more efficacious than greenhouse gases at atmospheric temperature change \cite[][]{HaN04}. Dirty snow feedbacks change throughout the aerosol lifecycle in the complex Arctic environment of snowfall, snowpack aging, snow-melt, drainage, and analogous sea-ice processes \cite[e.g.,][]{LEM98,AHH03,FlZ06}. Our goal is to assess absorbing aerosol interactions in the coupled Arctic climate system using models which represent the complex surface lifecycles of Arctic snow, BC, and dust, and which have been evaluated against satellite, in-situ, and laboratory measurements. We use the terms soot and BC interchangeably to denote the light absorbing component of carbonaceous aerosol \cite[][]{BoB05}. Interest in BC effects on climate, particularly Arctic climate, has increased as general circulation model (GCM) aerosol capabilities and emission inventories have improved. Recent noteworthy studies suggest that anthropogenic soot may have caused one quarter of last century's observed warming \cite[][]{HaN04}, and significant reductions in Northern hemisphere albedo and sea-ice extent \cite[][]{Jac04}. A lesson learned from our mid-latitude snow studies is that such estimates are extremely sensitive to accurate treatment of snowpack aging and radiative transfer \cite[][]{FlZ05}, two areas where this project will devote significant attention. Dust can also play an important role in the Arctic, far from its dominant sources in the sub-tropical deserts \cite[][]{PGT02,ZBN03}. Dust embedded in Arctic ice cores records global climate change \cite[e.g.,][]{AAG98,RaK97,FWJ99,KoH01} and marks periods of abrupt climate change \cite[e.g.,][]{All00}. The extent to which this embedded dust may change Arctic and global climate is still uncertain, and potentially significant for understanding abrupt climate change mechanisms \cite[][]{AWL00,HKR01,MML06}. Ice-albedo feedback is arguably the most important positive feedback in the polar climate system \cite[e.g.,][]{Har94,HoB03,QuH06}. However, we are unaware of any fully interactive atmosphere-land-ocean-glacier-sea-ice model which accounts for soot and dust radiative interactions in the surface Arctic. Fortunately the necessary component models for an integrated model to assess, predict, and improve understanding of Arctic absorbing aerosol exist \cite[e.g.,][]{CBB06}. Coupled component models are commonly used to study Arctic climate \cite[e.g.,][]{HoB03,HBH06}. This project uses high quality component models with realistic snow-ice-aerosol physics \cite[][]{FlZ05,FlZ06} to ask questions about the integrated impacts of absorbing aerosols on Arctic climate. \csznote{ \begin{enumerate*} \item Interpreting historical changes in sea-ice extent and fraction \item Potential for Arctic BC to disrupt summertime sea-ice formation \item Role of extreme Boreal fire years in Arctic reflectance \item Use coordinated suite of NASA measurements \begin{enumerate*} \item MODIS \item MISR \item AMSR-E \end{enumerate*} \end{enumerate*} } % end csznote % Dusty snowpacks provide additional meltwater on Earth and % probably other planets \cite[e.g.,]{Clo87} \subsection{Organization}\label{sxn:org} This proposal is organized as follows. Section~\ref{sxn:prr} describes the results of our relevant, prior NSF-funded research. The project's historical and scientific context is in Section~\ref{sxn:bgr}. Our scientific objectives and specific hypotheses are in Section~\ref{sxn:obj}. Section~\ref{sxn:mth} describes the models and observations we will use to reach these goals. Section~\ref{sxn:crd} summarizes the project plan, personnel responsibilities, time-line, milestones, and travel. Projects related to ours, potential broader scientific impacts, and our education plan are in Section~\ref{sxn:mpc}. Four letters of support/collaboration and a list of acronyms and abbreviations appear as supplementary documents. \section{Results from Prior NSF Funding on Related Projects}\label{sxn:prr} PI Zender is PI on \href{http://www.nsf.gov/awardsearch/showAward.do?AwardNumber=0321380}{ATM-0321380}: ``Acquisition of an Earth System Modeling Facility (ESMF) for Coupled Climate, Chemistry, and Biogeochemistry Studies''. The \href{http://www.ess.uci.edu/esmf}{ESMF} opened in February 2004 and supports about thirty users. The ESMF is the main computational resource for UCI's Earth System Science (ESS) Department. All ESS modeling groups use the ESMF. This includes eight professors, four full-time researchers, and about a dozen graduate students performing global modeling/analysis as part of their dissertations. In order to utilize spare cycles when ESS jobs are not running, the ESMF provides free (but lower priority) access to researchers/students from any UCI department. A summary of (about one dozen) published and in-press papers which used and acknowledged the collaborative ESMF scientific computing facility is maintained \href{http://dust.ess.uci.edu/prp/prp_mri/prp_mri_apr_02.txt}{here}. Zender is/was PI on \href{http://www.nsf.gov/awardsearch/showAward.do?AwardNumber=0503148}{ATM-0503148}: ``SGER: Improving CCSM Snow/Ice Radiative and Heating Processes and Assessing the Importance of the Soot Albedo Effect''. This SGER grant funded graduate student Mark Flanner for one year, during which he published one paper \cite[]{FlZ05}, placed one manuscript in review \cite[]{FlZ06}, and gave three national meeting presentations. \cite{FlZ05} shows that resolving the vertical distribution of solar radiant heating within snowpack remediates significant climatological biases in the Tibetan Plateau region. \cite{FlZ06} study snowpack aging and its effect on albedo. This the first theoretical study to quantify the relative roles of initial size distribution, vertical temperature gradient, and snow density in snow albedo evolution (discussed further in Section~\ref{sxn:snicar}). This proposal requests one year of postdoctoral funding for Flanner. \section{Background}\label{sxn:bgr} \subsection{Relevance and Historic Trends}\label{sxn:hst} Bright surfaces (snow, glaciers, sea-ice and clouds) make the Arctic uniquely susceptible to radiatively induced effects of surface BC and dust such as ice-albedo feedback amplification \cite[e.g.,][]{HoB03,HaN04,QuH06} (Figure~\ref{fgr:sot_snw_fdb}). \begin{figure} \centering % \centering uses less vertical space than center-environment \includegraphics[width=0.50\hsize,angle=0,clip=true,trim=1.0in 1.0in 1.0in 1.0in]{/data/zender/fgr/snicar/sot_snw_fdb}% \caption[Snow-albedo feedback schematic]{ Absorbing aerosols like soot and dust mediate snow-albedo feedback via multiple paths. \label{fgr:sot_snw_fdb}} \end{figure} Snow-albedo feedback is triggered by any forcing mechanism which changes the areal extent of snow cover. A weaker, positive feedback associated with changes in net surface radiation is the change in growth rate of snow grains. Soot in the snowpack directly lowers snow albedo and increases the growth rate of snow grains, lowering albedo of the ice grains themselves. Furthermore, the instantaneous perturbation of soot is greater in larger-grained snowpack, effectively increasing the gain ($\GGG$) on feedback involving grain growth. Finally, a fourth mechanism perturbation may result from accumulation of hydrophobic impurities at the surface during melt events, as supported by observations from \cite{ClN85} and \cite{CGR96}. Absorbing aerosol also alters cloud reflectivity and lifecycle \cite[][]{CLV96,ATS00}. However, cloud-aerosol effects are not a focus of this project and will not be further mentioned. Surface and atmospheric concentrations of black carbon (BC), i.e, absorbing carbonaceous aerosol, are highly variable and slowly increasing in the Arctic. Most emission scenarios project an increase in anthropogenic BC emissions of 30--250\% in the 21st century \cite[][]{NAD00,KoH05}. Changes in fire regime also affect total BC emissions \cite[][]{VRC04}. Hence BC impacts on the Arctic will likely increase through the 21st century. Dust plays an important role in accelerating snow melt in some mid-latitude regions (Dr.~Tom Painter, NSIDC, personal communication). Dust-snow interactions may play a very important role in glacial climates due to increased aridity, equator-to-pole temperature gradients, exposed continental shelves, and peri-glacial dust generation \cite[][]{RSK97,MML05,MML06}. Estimates of past and future dust emissions contain large uncertainties including those due to land use change and to \COd\ fertilization of vegetative cover \cite[e.g.,][]{MaL03,TWH04}. Arctic atmospheric and snowpack BC measurements span a wide range of concentrations \cite[][]{ClN85,NoC88,HaN04}. Greenland concentrations are typically 1--4\,\ugxkg, and as high as 30\,\ugxkg\ \cite[][]{SCD02}. Models driven by satellite-derived emissions estimates suggest that boreal fires explain the largest component of interannual Arctic BC deposition variability \cite[][]{KoH05}. Our model \cite[][]{FZR05} uses prescribed fossil and biofuel BC emissions \cite[][]{BSY04}, and converts fire emissions \cite[][]{VRC03,VRC04,RVC05} to BC with estimated emissions factors \cite[]{AnM01}. We estimate that biomass burning BC emissions north of 40\dgrn\ increased from 0.1\,Tg to 0.58\,Tg between 1997, a weak fire year, and 1998, a strong fire year. We study Arctic climate-aerosol interactions within a coupled general circulation model (GCM) framework (Section~\ref{sxn:mdl}). Using this tool, we estimate Boreal fire emissions changes from 1997 to 1998 increase surface snowpack radiative forcing in the Arctic by about~50\% \cite[][]{FZR05} (Figure~\ref{fgr:frc_sfc}). \begin{figure} \centering % \centering uses less vertical space than center-environment \includegraphics[width=0.4\hsize,angle=270,clip=true,trim=0.0in 0.0in 0.0in 1.5in]{/data/zender/fgr/snicar/sotfrc_sfc_1997f11_JJA}% \includegraphics[width=0.4\hsize,angle=270,clip=true,trim=0.0in 0.0in 0.0in 0.0in]{/data/zender/fgr/snicar/sotfrc_sfc_1998f11_JJA}% \caption[Snowpack soot direct radiative forcing]{ Summertime mean surface direct radiative forcing [\wxmS] by soot in snowpack during weak (1997) and strong (1998) boreal burn years. \label{fgr:frc_sfc}} \end{figure} These estimates contain many uncertainties and potential Arctic aerosol-related biases including transport and deposition, size distribution, optical properties, aging, and cloud interactions. These longstanding issues are targets of ongoing research and International Polar Year (IPY) activities such as POLARCAT (Section~\ref{sxn:polarcat}). More to the point, recent investigations by us and others suggest that aerosol-snow interactions play a significant role in Arctic climate sensitivity via snow-ice-albedo feedback mechanisms on which rapid progress can be made. \subsection{Absorbing Arctic Aerosol in Models}\label{sxn:mdl} The fundamental physics of snowpack-aerosol radiative interactions have been understood for over 25~years \cite[][]{WiW80,WaW80,Clo87}. However, few GCMs explicitly account for these effects so the potentially important polar climate amplification of dirty snow is largely unstudied. Models which do not explicitly account for absorbing aerosol may prescribe snowpack darkening as a function of time since last snowfall \cite[e.g.,][]{ODB04}. This approach is consistent with observations that snowpacks darken with time due to decreasing snow grain specific surface area (SSA) and to increasing impurity concentration (e.g., BC deposition) \cite[e.g.,][]{WaW80,AHH03} (Figure~\ref{fgr:niwot}). \begin{figure} \centering % \centering uses less vertical space than center-environment % Rotate Mark's image by 270 degrees so baselines agree % cd ~/ppr_FlZ06;fl_foo='lggnx';cat ${fl_foo}.eps | epsffit -r 0 0 282 1020 | epsffit -r 0 0 282 1020 | epsffit -r 0 0 282 1020 > ${fl_foo}_270.eps \includegraphics[width=0.50\hsize,angle=0,clip=true,trim=0.0in 0.0in 0.0in 0.0in]{/home/zender/ppr_FlZ06/fgr_niwot2}% \includegraphics[width=0.50\hsize,angle=0,clip=true,trim=1.30in 0.0in 1.30in 0.0in]{/home/zender/ppr_FlZ06/lggnx_270}% \caption[Niwot Ridge albedo decay, Isothermal specific surface area]{ Left Panel: Observed and modeled albedo decay at Niwot Ridge following the January~2, 2001 snowfall event. Error bars represent on standard deviation of all measurements comprising each day's albedo change. Right Panel: Observed and modeled (SNICAR) isothermal snowpack specific surface area (SSA) \cite[]{LTD04,FlZ06}. \label{fgr:niwot}} \end{figure} The observed 0.05 snowpack albedo change within four days at Niwot Ridge illustrates the important role of snowpack aging in controlling surface energy budgets. Prescribed snowpack aging can capture this aging effect on timescales of a few days. More sophisticated physics are required to represent longer timescale albedo evolution \cite[][]{FlZ06}. Single layer models are significantly limited in their representation of aerosol transport and removal mechanisms and in correctly representing sub-surface radiative heating and melt. \cite{FlZ05} studied the vertical distribution of solar radiant heating on snowpack. Using accurate radiative transfer methods in a multi-layer snowpack model (Section~\ref{sxn:snicar}), we showed that 20--40\% of solar heating occurs beneath the top 2\,\cm\ of typical snowpack. This sub-surface heating leads to internal snowpack melt. Internal snowpack absorption can melt snow more efficiently per unit absorption than surface snowpack absorption which, in the Arctic, is typically balanced by strong sensible heat fluxes due to cold overlying air \cite[]{MPB04,FlZ05}. Snow insulation keeps the internal snowpack warmer than the surface so that radiatively-induced internal warming can trigger the moderate snow-grain size-feedback loop (Figure~\ref{fgr:sot_snw_fdb}) until melting triggers the stronger snow-cover feedback loop. %However, all GCMs because they confine solar absorption to the top %snow layer. Today's most sophisticated snow and ice models include Arctic soot, and, less frequently, dust \cite[][]{Jac04,HaN04,FZR05}. Interestingly, our preliminary simulations of Greenland show that soot can reduce total snow melt. Significant soot concentrations in the surface snowpack strongly absorb visible radiation some of which would otherwise penetrate into the snowpack (Figure~\ref{fgr:rfl_spc_rds_ffc}). \begin{figure} \centering % \centering uses less vertical space than center-environment \includegraphics[width=0.50\hsize,angle=0,clip=true,trim=0.0in 0.0in 0.0in 0.0in]{/data/zender/fgr/snicar/snw_sot_spectral_rfl}% % Rotate Mark's image by 270 degrees so baselines agree % epsffit coordinate arguments are: llx lly urx ury in Postscript units (points) % (lower-left x, lower-left y, upper-right x, upper-right y) % cd ${DATA}/fgr/snicar;fl_foo='rdsffc_1997f11-1998b11_JJA';cat ${fl_foo}.eps | epsffit -r 33 141 543 725 | epsffit -r 33 141 543 725 | epsffit -r 33 141 543 725 > ${fl_foo}_270.eps \includegraphics[width=0.50\hsize,angle=0,clip=true,trim=0.0in 0.0in 0.0in 0.0in]{/data/zender/fgr/snicar/rdsffc_1997f11-1998b11_JJA_270}% \caption[Spectral reflectance of snow; Response of effective radius]{ Left Panel: Spectral reflectance of pure snow and snow externally mixed with 200\,\ugxkg\ BC for different snow size distributions. Vertical lines show positions of MODIS Bands 3, 4, and~5 (left-to-right). Right panel: Change in summertime-mean effective radius $\rdsffc$\,[\um] of surface snowpack layer due to 1998. Cross-hatching indicates statistically significant changes ($\cnfstt < 0.05$) relative to simulations without boreal soot \cite[][]{FZR05}. \label{fgr:rfl_spc_rds_ffc}} \end{figure} Hence surface soot concentrations can cool the lower snowpack much as atmospheric soot cools the surface by reducing insolation. The screening effect of surface soot competes with the temperature-grain-size feedback (Figure~\ref{fgr:sot_snw_fdb}). The ramifications of these competing feedbacks for Arctic climate remain unknown. Nevertheless, it is clear that a prognostic approach which de-convolves SSA and aerosol effects on snowpack aging and heating is required for Arctic aerosol-climate impact studies such as this project. Previous efforts to understand absorbing aerosol impacts on the Arctic have made many approximations which disallow some important feedbacks that occur in Nature. This project, in collaboration with others, will result in a more complex and interactive Arctic system where \begin{enumerate*} \item Prognostic glaciers allow determination of net hydrology \item Sea-ice is sensitive to surface aerosol radiative effects \item Snow-ice-aerosol interactions are consistent across the cryosphere \item Soot is treated as fractal aggregates not spheres \item Snowpack aging reproduces observed features \item Glacial dust and boreal soot sources are accounted for % thermodynamic fixed Qflux sea-ice model in % slab ocean (lateral heat transport is prescribed) %\item Externally mixed soot not internally mixed soot %\item Provide improvement to global uniform snow grain size assumption % for future satellite surface reflectance retrievals. %\item Goal to monitor Arctic soot/dust concentration from satellites %No one has retrieved this yet \end{enumerate*} %Absorbing aerosol species like soot and dust are typically not %carried as prognostic, radiatively interactive tracers in %Arctic sea-ice models. %The Arctic Ocean Model Intercomparison Project %(\href{http://fish.cims.nyu.edu/project_aomip/overview.html}{AOMIP}) % \cite[e.g.,][]{UHM05}. \section{Scientific Objectives and Hypotheses}\label{sxn:obj} Our studies of Arctic snow-ice-aerosol processes will improve understanding of ice-albedo feedbacks and polar climate amplification in Nature, and their representations in models. Key scientific questions we will address include: \setcounter{enmrfr}{0} % Reset reference counter for this list \begin{enumerate} \item \enmrfrstp \label{idx_obj_sea_ice} \textbf{Objective}: Arctic absorbing aerosol impacts on polar climate sensitivity mediated by sea-ice\\ \textbf{Hypothesis}: Arctic BC can significantly disrupt summer sea-ice formation and extent during strong burn years\\ Multiple lines of evidence support this hypothesis: First, representation of thin sea-ice amplifies polar climate sensitivity \cite[][]{HoB03,HBH06}. Second, internal snowpack heating amplifies mid-latitude climate sensitivity \cite[][]{FlZ05}. Third, aging and absorbing aerosol content increase polar climate sensitivity \cite[][]{Jac04,HaN04,FZR05}. Moreover our preliminary investigations with slab ocean models and simple sea-ice models suggest a summertime Arctic sea-ice response to boreal soot in strong fire years. Summertime Arctic sea-ice extent has been declining in recent decades, likely related to greenhouse gas-induced warming \cite[][]{SMS03,SSF05}. We will quantify the extent to which snow-aerosol interactions may contribute to this trend. Boreal emissions estimates for recent strong fire years \cite[][]{RVC05} will allow us to search for connections to recent accelerations in Arctic sea-ice reduction \cite[g.g.,][]{SSF05}. \item \enmrfrstp \label{idx_obj_sot_dst} \textbf{Objective}: Quantify relative roles of Arctic soot and dust as polar climate amplifiers\\ \textbf{Hypothesis}: Boreal soot amplifies Arctic climate sensitivity more/less than dust in the present/LGM climate.\\ Soot is approximately an order of magnitude more absorptive than dust at solar wavelengths (assuming single scattering albedos of 0.5 and 0.95 for soot and dust, respectively). Dust deposition to Greenland in glacial periods is 2--20 times greater than present day \cite[e.g.,][]{RaK97,MML06}, enough to compete with BC as a snow-albedo trigger. However, differences between BC and dust deposition seasonality and variability will modulate the net solar forcing of these aerosols on Arctic surfaces. Dust of Asian provenance \cite[][]{BGR97} will likely deposit more continually than North American glacial dust \cite[][]{MML06}. The springtime maximum Asian dust export \cite[][]{ZBN03} is less efficacious for Arctic forcing than summertime aerosol events. How these spatio-temporal deposition patterns affect Arctic climate sensitivity is nearly completely unexplored. \item \enmrfrstp \label{idx_obj_grn_rfl} \textbf{Objective}: Role of Boreal fire variability in Arctic reflectance, particularly Greenland\\ \textbf{Hypothesis}: BC warms Greenland in strong fire years and cools Greenland in weak fire years.\\ Atmospheric BC cools the surface by backscattering and absorbing incident sunlight. Snowpack~BC heating compounded by snow-albedo feedback can exceed atmospheric~BC surface cooling in strong fire years \cite[][]{FZR05} (Figure~\ref{fgr:mlt_Grn}). Ice core analyses (Dr.~Eric Saltzman, UCI, personal communication) and model simulations \cite[][]{KoH05,FZR05} agree that boreal fires are the primary source of BC deposition to Greenland in strong fire years. We will convolve ice core records of historic BC deposition to Greenland with present day spatially explicit BC emissions data \cite[][]{RVC05} to study maximum changes in Greenland reflectance and melt due to boreal BC over the past 1000~years. \end{enumerate} \section{Methods: Arctic Models and Observations}\label{sxn:mth} It is important to emphasize that this project will not develop any Arctic climate model components from scratch. Our intellectual efforts are primarily directed toward uncovering the influence of previously neglected snow-ice-aerosol interactions in the Arctic system. All model development tasks outlined below involve improving physics in our in-house snow-ice-aerosol model (SNICAR) and/or merging these physics into high-quality Arctic system component models developed and maintained at National centers. \subsection{Community Climate System Model}\label{sxn:ccsm} An integrated Earth System Model which fully couples aerosols, snow, atmosphere, ocean, and land/sea-ice is required to test our hypotheses (Section~\ref{sxn:obj}). The NCAR CCSM is that model. CCSM polar climate simulations have been continually evaluated against meteorological analyses and satellite observations \cite[e.g.,][]{BrB981,HoB03,HBH06}. Arctic absorbing aerosols, soot and dust, are primarily emitted from non-frozen land surfaces at lower latitudes \cite[e.g.,][]{ZBN03,KoH05}, and, as such, cut across multiple climate ``spheres'' (biosphere-atmosphere-cryosphere). The project relies on our continuing external collaborations for realistic aerosol distributions and simulation codes. The CCSM BC/OC aerosol transport, deposition, and optics we use come from long time collaborators Drs.~Phil Rasch and Bill Collins (NCAR) \cite[][]{CRE01,CRE02}. Long time collaborator Dr.~Natalie Mahowald (NCAR, see attached letter of support) studies atmospheric biogeochemistry and dust-carbon-climate interactions on multiple timescales. Mahowald and PI~Zender are primary developers of the Dust Entrainment and Deposition (DEAD) mineral dust model \cite[][]{ZBN03,ZNT03,MLD03}. Ice core measurements are consistent with peri-glacial dust production during glacial periods \cite[][]{MML06}. Our simulations will account for this glacial dust production using the method and model of \cite{MML06}. \subsection{SNICAR}\label{sxn:snicar} The project builds upon, extends, and applies our existing, state-of-the-art, SNow, ICe, and Aerosol Radiative model, SNICAR \cite[][]{FlZ05,FlZ06}. SNICAR treats snowpack radiative transfer, aging, and aerosol interactions in a unified manner that allows for realistic feedbacks between solar radiation, snowpack temperature gradients, and aerosol concentration. The vertically resolved multi-layer radiative transfer component \cite[][]{WiW80,TMA89} treats snow as a collection of hexagonal prisms based on the equivalent surface area-to-volume approximation \cite[][]{GrW99,NGW03}. A lookup table (computed off-line) contains Mie parameters (single scattering albedo, extinction coefficient, and asymmetry parameter) for any lognormal size distribution of snow. SNICAR accounts for solar zenith angle, direct and diffuse incident radiation, bare surface reflectance \cite[][]{DZD03}, and vertically-resolved effective radius (\rdsffc), snow depth, density, and concentrations of absorbing impurities \cite[][]{WaW80}. An off-line version of SNICAR runs at high spectral resolution, 10\,\nm\ resolution in the solar spectrum from 0.3--5.0\,\um\ (470 bands). This is useful for simulating narrow-band satellite channels such as MODIS/MISR channels (Figure~\ref{fgr:rfl_spc_rds_ffc}). For climate simulations, SNICAR runs in a host snowpack model, typically the Community Land Model (CLM) \cite[][]{DZD03}. The CLM uses five vertical snowpack layers \cite[]{ODB04} and itself runs off-line forced by meteorological analyses or on-line in a GCM. We nest CLM/SNICAR in the Community Atmosphere Model (CAM) \cite[][]{CRB06} with modifications for prognostic soot driven by recent fossil fuel and biofuel emissions estimates \cite[][]{BSY04}, and biomass burning emissions \cite[][]{VRC03,VRC04,RVC05}. Our simulations suggest boreal soot causes seasonal net surface solar radiation forcings of 0.5--0.75\,\wxmS\ (Figure~\ref{fgr:frc_sfc}) in strong fire years. These forcings induce feedbacks such as larger snow grain size (Figure~\ref{fgr:rfl_spc_rds_ffc}) which together increase seasonal surface absorption by more than 1.5\,\wxmS\ \cite[][]{FZR05}. Soot-snow feedbacks in strong fire years appear to cause significant increases in meltwater production in Greenland snowpack (Figure~\ref{fgr:mlt_Grn}). \begin{figure} \centering % \centering uses less vertical space than center-environment \includegraphics[width=0.33\hsize,angle=270,clip=true,trim=0.5in 1.0in 0.0in 1.75in]{/data/zender/fgr/snicar/qmelt_1997f11-1998b11_JJA_GR}% \includegraphics[width=0.33\hsize,angle=270,clip=true,trim=0.5in 1.0in 0.0in 1.75in]{/data/zender/fgr/snicar/qmelt_1998f11-1998b11_JJA_GR}% \includegraphics[width=0.33\hsize,angle=270,clip=true,trim=0.5in 1.0in 0.0in 0.0in]{/data/zender/fgr/snicar/qmelt_1998g11-1998d11_JJA_GR}% \caption[Snow melt response to snowpack soot]{ Summertime mean change in Greenland snow melt [\mmxday] due to boreal soot during low (1997, left) and high (1998, middle and right) boreal burn years. Middle panel includes all feedbacks (soot in atmosphere and snowpack), while right panel includes atmospheric soot only. Cross-hatching indicates statistically significant changes ($\cnfstt < 0.05$) relative to simulations without boreal soot \cite[][]{FZR05}. \label{fgr:mlt_Grn}} \end{figure} Note that neglecting soot-snowpack interactions (and accounting only for atmospheric soot effects) eliminates or reverses the sign of most of the increased snow melt over Greenland. Hence, \textit{significant Arctic change is attributable to aerosol-snowpack feedbacks not represented in most GCMs}. This makes us eagerly anticipate results in year~3 when SNICAR is embedded in fully interactive sea-ice and glacier models which can respond to soot and to glacial dust sources \cite[][]{MML06}. \subsubsection{Snow Aging}\label{sxn:aging} Existing representation of snow aging in Arctic climate models is crude and empirical. Some models consider the role of temperature in albedo decay \cite[][]{Jor91,ODB04} though none consider the dominant role of temperature gradient (TG), which requires a multi-layer snow model such as SNICAR. High-latitude snowpack can have temperature gradients well in excess of 100\,\kxm\ owing to strong radiative cooling at the surface during night, and good thermal insulation of deeper snow. Vapor-density gradients induced by these temperature gradients cause rapid snow metamorphism \cite[][]{StB97}. Using first principles, we account for the roles of initial size distribution, temperature, temperature gradient, snow density, and inter-particle spacing in the evolution of snow specific surface area (SSA) \cite[][]{FlZ06}. Vapor diffusion from highly-curved surfaces characteristic of fresh, dendritic snow plays a smaller, but non-negligible role in grain growth via the Kelvin Effect in low TG environments \cite[e.g.,][]{Col80}. Recent SEM observations of well-sintered snow support the hypothesis that grain-boundary diffusion is an important mechanism in sintering \cite[][]{RoS06}, although the importance of this mechanism has been discounted in earlier studies. Meltwater flushing is the most important BC removal mechanism, since preferential gravitational settling only operates on external mixtures, and is likely extremely slow. Qualitative observations suggest that BC may become more concentrated in surface snow during melt events \cite[][]{WaW80,ClN85}. \cite{CGR96} spread hydrophobic and hydrophilic BC on top of snow, and noticed that hydrophobic BC remains in surface snow longer, maintaining lowered albedo for a longer time. Because of imprecise knowledge of vertical distributions of BC and meltwater formation in this experiment, we can only deduce that hydrophilic BC is scavenged about five times more efficiently than hydrophobic BC. Using a simple $\me$-folding removal model with this data, we estimate lower bounds for scavenging ratios of 0.01 and 0.002 for hydrophilic and hydrophobic BC, respectively. Even greater uncertainty exists for snow processes on sea-ice. Planned field studies from Warren and Grenfell (Section~\ref{sxn:insitu}) will help constrain these scavenging factors, while current studies by Dr.~Tom Painter (NSIDC) will help estimate scavenging efficiency for dust. In summary, this project will improve representation of these snowpack aging processes in SNICAR: \begin{enumerate*} \item Sintering and melting/re-freezing in temporal SSA decay. Mark Flanner and Dr.~Tom Painter (NSIDC) will collaborate on representing this process which may be very important for surface albedo in low-TG environments \item Frost deposition that brightens snow surfaces \cite[][]{Pir04}. We will explore means of parameterizing diurnal SSA increases due to hoar frost. \item Meltwater scavenging of soluble and insoluble aerosols \end{enumerate*} \subsubsection{Optics}\label{sxn:opt} Snow and aerosol optical properties link the snowpack microphysical properties (aerosol concentration, particle size distributions) to macroscopic net absorption (Figure~\ref{fgr:frc_sfc}a), reflectances (Figures~\ref{fgr:niwot} and~\ref{fgr:rfl_spc_rds_ffc}a), and heating rates that drive the snow melt and temperature change which trigger snow-albedo feedback. These responses are sensitive to optical property assumptions which this project will explore and improve, including \begin{enumerate*} \item BC indices of refraction: \cite{BoB05} question the OPAC properties \cite[][]{HKS98} (which we use) and recommend other measurements including \cite{ChC90} \item BC shape: Treating BC as spheres likely underestimates its single scattering albedo relative to more realistic shapes such as fractal aggregates \cite[][]{Sor01,BoB05} \item Aerosol mixing: BC and dust in remote regions such as the Arctic are primary deposited via wet scavenging \cite[][]{CCR01,CSK04,ZBN03} and so will often be internally mixed within snow grains. We will adopt an effective medium approximation \cite[e.g.,][]{BoH83} to represent internally mixed aerosols. This will increase snowpack absorption relative to our current externally mixed assumption. %\item Darkening by desert and peri-glacial dust sources \cite[][]{MML06}. \end{enumerate*} \subsection{Sea-Ice and Ice Sheets}\label{sxn:cice} Sea-ice is the fulcrum of ice-albedo feedbacks in the Arctic Ocean. Drs.~Bill Lipscomb and Elizabeth Hunke of LANL (see attached letter of support) are the principle developers of the Los Alamos sea-ice model (\href{http://climate.lanl.gov/Models/CICE}{CICE}), a primary component of the National Center for Atmospheric Research (NCAR) Community Climate System Model (CCSM). The CICE model contains an Ice Thickness Distribution (ITD) which maintains a half dozen prognostic categories of ice thickness in each grid cell \cite[][]{HBH06}. The thinnest ice category is most susceptible to changes in net solar radiation due to snowpack aging and aerosol concentration. Accounting for impurities in sea-ice is important to accurate radiative transfer throughout the atmosphere/ice/ocean system \cite[]{Gre91}. Currently, CICE uses a single layer snowpack upper boundary condition and neither the sea-ice nor the snowpack tracks absorbing impurities such as soot and dust as prognostic tracers. A multi-layer snowpack and prognostic aerosol tracer capability will be in CICE in 2006 (W.~Lipscomb, personal communication). Postdoc Mark Flanner will merge SNICAR physics (snow aging, radiative transfer, snow-aerosol optics) into this CICE configuration. The resulting CICE simulations will retain soot and dust deposited directly on bare sea-ice from the atmosphere and from melting snow cover. We hope to use sea-ice radiative transfer physics developed by Bruce Briegleb (NCAR) to account for aerosols embedded in the complex sea-ice-brine-pond matrix. Sea-ice extent and thickness may then respond to the full lifecycle of Arctic~BC. This will be a significant improvement to current models which remove BC from snow and sea-ice with rather ad~hoc mechanisms to prevent excessive concentrations from accumulating in multi-year land and sea-ice \cite[][]{Jac04,FZR05}. Residence time estimates for optically active Arctic soot will improve. This may play an important role in studies of Objective~\ref{idx_obj_sea_ice} (Section~\ref{sxn:obj}). Lipscomb and Hunke are developing an interactive ice sheet component for the CCSM. This ice sheet model is based on \href{http://wiki.nesc.ac.uk/read/glimmer-project}{GLIMMER} \cite[]{Pay99}. The snow and energy balance model which sits atop GLIMMER will be based on the Community Land Model (CLM). Once postdoc Flanner merges the SNICAR snow-aerosol physics into LANL's GLIMMER, the CCSM will have a glacier model component sensitive to realistic snowpack physics and aerosol interactions. This will complete the integrated treatment of clean and dirty snow in the cryosphere. The ice sheet-aerosol component is scheduled for completion in project year~3 but may lag due to its complexity and resource issues beyond our control. In any case, our numerical studies of the coupled Arctic land-ocean-sea-ice system with fixed ice sheets will proceed apace toward Objectives~\ref{idx_obj_sea_ice} and~\ref{idx_obj_sot_dst} (Section~\ref{sxn:obj}). \subsection{Numerical Experiment Strategy}\label{sxn:xpt} Our questions (Section~\ref{sxn:obj}) will be addressed in the context of pre-industrial, present day, and next century timescales when appropriate. Natural (i.e., unforced) interannual variability is quite large in the Arctic climate system \cite[e.g.,][]{BrB982,FWJ99}. Boreal fire variability is also quite large \cite[e.g.,][]{RVC05} and seems to explain most of the variability in Arctic BC deposition \cite[][]{FZR05}. Detecting and assessing the relatively small (though important) signal of aerosol-induced Arctic change (Figures~\ref{fgr:rfl_spc_rds_ffc}b and~\ref{fgr:mlt_Grn}) against the noisy background of natural Arctic variability is difficult. We will continue to employ an ensemble-based approach to increase the signal/noise ratio. The ensemble comprises multiple identical numerical experiments with slightly perturbed initial conditions. To obtain the climate responses presented in this proposal we conducted perennial 1997- and 1998-emissions experiments in separate 15\,yr. simulations. We used a Student's $\ttt$-test to quantify statistical significance of Arctic changes between the two ensembles. Significant ($\cnfstt < 0.05$) changes appear as cross-hatched regions in Figures~\ref{fgr:rfl_spc_rds_ffc}b and~\ref{fgr:mlt_Grn}. \subsection{Satellite Observations}\label{sxn:stl} NASA MODIS, MISR, and AMSR-E retrievals can constrain free model parameters and help us interpret the regional and seasonal behavior of snowpack processes. Figure~\ref{fgr:rfl_spc_rds_ffc}a shows simulated snow spectral reflectance expected in visible MODIS bands for various grain sizes. Soot concentration is most apparent in visible channels and particle sizes information is most distinguishable in the near infrared (NIR) \cite[][]{PDR03}, e.g., near MODIS channel~5. Greenland is an ideal location for evaluation of SNICAR from remote sensing platforms. Much of the ice sheet enjoys year round sub-freezing temperatures which remove the potentially confounding influence of liquid effects. Since soot concentrations rarely exceed 5\,\ugxkg\ in Greenland, surface snowpack effective radius~$\rdsffc$ is the most promising model parameters to retrieve and constrain. Remote sensing of snow grain size has been demonstrated with hyperspectral airborne instruments such as AVIRIS \cite[][]{PDR03,DoP04}. The same principles apply to MODIS and/or MISR radiances. However, these satellite products currently have problems associated with large zenith angles and topography. Once the MODIS/MISR spectral snow reflectance ($\rfl$) products reach robust operational status, we will use them to evaluate SNICAR spectral reflectance, and to try to infer~$\rdsffc$ (Figures~\ref{fgr:rfl_spc_rds_ffc}a and~b, respectively). (We will happily provide our predicted~$\rdsffc$ to any retrieval experts attempting to improve MODIS/MISR surface reflectance). In combination with temperature from meteorological analyses, retrieved~$\rfl$ and/or $\rdsffc$ will be used to evaluate SNICAR's snow aging physics \cite[][]{FlZ06} which predict significant temperature dependence for~$\rdsffc$ (Figure~\ref{fgr:niwot}). AMSR-E retrieves Snow Water Equivalent (SWE) over non-ice surfaces. SWE retrievals may provide useful constraints on SNICAR simulations of continental snowpack. In particular, we are interested in assessing the influence of extreme aerosol events on the duration of snow cover in seasonally snow covered regions which are most susceptible to snow-albedo feedback \cite[][]{FlZ05}. \subsection{In Situ Observations}\label{sxn:insitu} Arctic snowpack BC concentration is the key diagnostic which integrates aerosol source, transport, deposition, and melt processes. Drs.~Steve Warren (U.~Washington, see attached letter of support), Tom Grenfell and Tony Clarke (U.~Hawaii) proposed an ARO project ``Black carbon in Arctic snow and ice, and its effect on surface albedo'' which will measure BC in snow and ice in tundra regions of Alaska, Canada, and Russia, as well as on the Greenland Ice Sheet and the Arctic Ocean to update and improve the BC survey that Clarke conducted in 1983--1984. They will (continue) to share their measurements with us, including, potentially, BC measurements collected at Dr.~Konrad Steffen's automated meteorological sites in Greenland. Warren et~al.\ will also measure/estimate scavenging coefficients for removal of atmospheric BC by snow and removal of surface BC during snow melt (Section~\ref{sxn:aging}). These scavenging coefficients will provide important constraints for BC scavenging in CCSM and SNICAR, respectively. In addition to Warren and Clarke's earlier measurements, \cite{FlZ06} evaluated SNICAR against \textit{in situ} and laboratory measurements of snowpack specific surface area, crystal density, albedo, and curvature- and temperature-gradient growth processes. This project will enable us to continue these comparisons as new data become available. \subsection{IPY POLARCAT Participation}\label{sxn:polarcat} We are contributing to the IPY ``POLar study using Aircraft, Remote sensing, surface measurements and modeling of Climate, chemistry, Aerosols and Transport'' (\href{http://zardoz.nilu.no/~andreas/POLARCAT}{POLARCAT}) (see attached letter from IPY program office to POLARCAT PI~Stohl). Our contribution is to one of POLARCAT's main themes---the influence of boreal fire aerosol on Arctic surface properties. Although many observational aspects of POLARCAT are still pending, support for regular aerosol observations at Summit, Greenland appear to be in place. Commitments for aircraft campaigns extending from boreal forests to Greenland are likely in summer 2008. Using modeled/assimilated BC deposition from NCAR collaborator and POLARCAT Steering Committee member P.~Rasch, our group will estimate surface reflectance changes at Summit from significant boreal events upwind, and compare them to in~situ observations. The most important measurements we need to help reconcile discrepancies between our SNICAR model and broadband reflectance measurements are snow accumulation and vertical profiles of aerosol concentration, and snowpack temperature. Spectral reflectivity would be very valuable but is less likely to be measured. We will consider investigating other targets of opportunity which may arise during POLARCAT as instrument availability and fire events allow. \section{Project Coordination}\label{sxn:crd} \subsection{Personnel}\label{sxn:prs} PI~Zender will coordinate all project activities. Zender has extensive experience in global-scale aerosol modeling and ice cloud physics and radiative transfer. Zender will use DEAD \cite[][]{ZBN03} embedded in the CCSM/SNICAR with peri-glacial sources \cite[][]{MML06} to provide Arctic dust deposition fields. In addition, Zender will develop, test, and implement optical property improvements in SNICAR, including soot fractal aggregates \cite[][]{Sor01}, improved refractive indices \cite[][]{ChC90}, and internally mixed dust/soot/snow properties. Postdoc Mark Flanner (currently UCI graduate student in ESS) will merge SNICAR physics into the CICE sea-ice model and into the CCSM-compatible version of the GLIMMER ice sheet model provided by LANL. Scientific specialist Dr.~Chao Luo has specialized in atmospheric dust transport \cite[][]{LMD03,MaL03,LMJ04} and, more recently, cryospheric hydrology based on AMSR-E retrievals. Luo will assist with coupled model simulations and satellite data analysis. Zender will advise an ESS graduate student in Arctic aerosol studies. We will initially focus on satellite data analysis of Arctic surface reflectance. In the second and third years, we will use satellite products to help evaluate and constrain SNICAR, as described in Section~\ref{sxn:stl}. \subsection{Schedule and Milestones}\label{sxn:tm_ln} \noindent\textbf{Year~1}. \textit{Milestones}: 1a.~SNICAR interactive with CICE; 1b.~Fractal aggregate soot; 1c.~Snowpack can brighten diurnally\\ \textit{Tasks}: \begin{enumerate*} \item Couple SNICAR physics to sea-ice, examine impact on summer melt and extent \item Update soot optical properties to fractal aggregates \item Represent hoar frost, diamond dust snowpack brightening \end{enumerate*} \textit{Travel}: \begin{enumerate*} \item Zender and Postdoc Flanner to Boulder (supported by NCAR affiliate scientist and visitor funds) to collaborate with Mahowald and Rasch (NCAR) and Painter (NSIDC) \item Zender and Postdoc Flanner present results at AGU \end{enumerate*} \textbf{Year~2}. \textit{Milestones}: 2a.~Diurnal snowpack $\rfl$ matches observations; 2b.~SNICAR coupled to LANL glacier-model; 2c.~IPY POLARCAT participation\\ \textit{Tasks}: \begin{enumerate*} \item Represent sintering, especially in non-TG snowpacks \item Refine scavenging with Warren et~al.'s Greenland BC measurements; \item Compare SNICAR predictions with MODIS/MISR-inferred $\rfl$, $\rdsffc$ \item Quantify Arctic sea-ice sensitivity to soot and dust separately and together \end{enumerate*} \textit{Travel}: \begin{enumerate*} \item Zender to Norway/NILU for IPY POLARCAT team meeting \item Graduate student to Los Alamos to merge SNICAR into GLIMMER \item Zender presents results at AGU \end{enumerate*} \textbf{Year~3}. \textit{Milestones}: 3a.~SNICAR fully coupled in land/sea/glaciers; 3b.~POLARCAT event simulation\\ \textit{Tasks}: \begin{enumerate*} \item Estimate soot/dust effect on Greenland accumulation/ablation budgets \item Integrated absorbing aerosol impact on Arctic climate sensitivity \item Scale fire emissions from Randerson by ice core data (from Saltzman and McConnell) to estimate 1000\,year BC impacts on Greenland \end{enumerate*} \textit{Travel}: \begin{enumerate*} \item PI Zender and graduate student present results at AGU \end{enumerate*} \section{Related Projects, Broader Impacts and Education}\label{sxn:mpc} \subsection{Related Projects}\label{sxn:clb} This project is synergistic with multiple existing and newly-proposed projects. In addition to collaborative exchanges with LANL, NCAR, and U.~Washington described within the proposal text and endorsed in attached letters of support, the following projects will benefit from ours: \begin{enumerate*} \item Drs.~Steve Warren (U.~Washington, see attached letter of support and Section~\ref{sxn:insitu}), Tom Grenfell and Tony Clarke (U.~Hawaii). We will continue to provide Warren's group with BC deposition simulations to aid them in choosing BC measurement site locations. \item Drs.~Natalie Mahowald (see attached letter of support and Section~\ref{sxn:ccsm}) and Phil Rasch (NCAR). We will continue to make SNICAR code improvements available to NCAR CCSM component models. Mahowald and Dr. Peter Thornton are investigating changes in boreal fire regime in a dynamic vegetation and carbon/nitrogen cycling framework. The more realistic lower boundary condition SNICAR provides will improve these studies' realism. \item Drs.~Bill Lipscomb and Elizabeth Hunke, LANL (see attached letter of support and Section~\ref{sxn:cice}). \item Dr.~Tom Painter (NSIDC) studies snow aging, dust-driven heating of mid-latitude snowpacks, and retrieval of snowpack reflectance, size distribution and impurities from remote sensing instruments. Painter will use SNICAR physics within the SNTHERM snowpack model to account for dust-snowpack interactions on catchment/basin scale hydrology. \item Drs.~Jim Randerson and Yufang Jin (UC~Irvine) are our primary in-house collaborators on boreal soot impacts. Their integrated studies of C-cycling along boreal fire chronosequences provide the fire emissions estimates which drive our simulations \cite[e.g.,][]{RVC05}. We will continue to collaborate with them on integrated forcing estimates from boreal fires and to use CCSM/SNICAR to quantify soot indirect effects such snowpack-mediated radiative forcing \cite[][]{FZR05,RLF05}. \item Dr.~Eric Saltzman (UC~Irvine) measures trace gas and aerosol concentrations in ice cores \cite[e.g.,][]{SAD04} and Co-PIs a proposed ARO project ``High-Resolution, Biomass-Burning-Specific Tracers in Greenland Ice Cores over the Past 1000~Years''. In conjunction with Randerson's fire emission database, our project provides a method to quantify the impact of Saltzman's proxy measurements of biomass burning aerosol variability in Greenland over the last 1000\,yr. \end{enumerate*} \subsection{Improved Community Modeling Capabilities}\label{sxn:mdl_mpr} Snow pervades the Arctic, and improved snowpack models can and probably will improve understanding of Arctic hydrology. Zender has provided many useful research models to various communities including the Column Radiation Model (\href{http://www.cgd.ucar.edu/cms/crm}{CRM}), the Dust Entrainment and Deposition Model (\href{http://dust.ess.uci.edu/dead}{DEAD}), and the netCDF Operators (\href{http://nco.sf.net}{NCO}). SNICAR will be freely available both in off-line and in Community Climate System Model modes. We anticipate SNICARS and its physics will contribute to Arctic research areas including glacier mass balance (through realistic upper boundary radiation and melt conditions), basin and catchment hydrology (improved representation of snowpack sublimation vs.\ internal melt and surface percolation), sea-ice lifecycle (improved upper boundary condition), paleoclimate sensitivity (through improved accuracy of Arctic responsiveness to orbital and aerosol forcing), and snow chemistry (through improved representation of snowpack specific surface area). \subsection{Education}\label{sxn:edc} This project trains one graduate student in Arctic aerosol-climate interactions. Moreover, PI~Zender helps train Orange County K--12 teachers in Earth Science curricula. An NSF Math Science Partnership project called \href{http://focus.mspnet.org}{FOCUS} (Faculty Outreach Collaborations Uniting Scientists, Students and Schools) brings the teachers to UCI. Zender will incorporate Arctic climate change materials into \href{http://dust.ess.uci.edu/smn/smn_sun_focus_200508.pdf}{his FOCUS seminars}. Zender will also incorporate this Arctic research into his \href{http://dust.ess.uci.edu/smn/smn_sun_olli_200511.pdf}{seminars} to students of the Osher Lifelong Learning Institute (\href{http://unex.uci.edu/community/olli/index.asp?}{OLLI}). \newpage \section*{References}\label{sxn:rfr} % \section* does not number section \setcounter{page}{1} \thispagestyle{empty} % Bibliography \renewcommand\refname{(blue numbers show proposal section where citation occurs)} \vspace{-24.0pt} \newlength{\oldbaselineskip} \setlength{\oldbaselineskip}{\baselineskip} %\setlength{\baselineskip}{13.574pt} % 1.234 X 12pt \setlength{\baselineskip}{12.0pt} % 1.234 X 11pt \setlength{\bibsep}{4pt} % Space between natbib bibliography items \bibliographystyle{jas} \bibliography{bib} \setlength{\baselineskip}{\oldbaselineskip} % 1.851 X 11pt \newpage \printindex % Requires makeidx KoD95 p. 221 \addcontentsline{toc}{section}{Index} \newpage \subsection{Budget Justification}\label{sxn:bdg_jst} \setcounter{page}{1} \thispagestyle{empty} \begin{verbatim} % NB: Do not use LaTeX formatting in Budget Justification since must % upload into Liz's Word document Salaries and Wages One month of summer salary support for three years is requested for Prof. Charles Zender, the PI, who has primary responsibility for the proposed research. Salary support for Mark Flanner, a postdoctoral scholar is requested in year 1. Dr. Flanner is responsible for fully coupling the SNICAR aerosol/snowpack model to the glacier dynamics model and to the sea-ice model, and for performing initial fully coupled studies. Funds are requested to to support Dr. Chao Luo, Associate Specialist Step II, at a rate of 0.2 FTE for the duration of the project. Dr. Luo is the principal scientific programmer associated with UCI's Earth System Modeling facility. Dr. Luo has run the complex models involved in this project (SNICAR, DEAD, CLM, CAM, CCSM) and will devote 20% of his time to performing and to analyzing coupled model studies. A 2% cost of living increase was applied each year of this proposal as well as a 5% merit, where applicable. To be Named---Graduate Student Researcher III. Funds are requested to support one non-resident graduate student each year of the project. The graduate student support is requested at 50% for 9 months during the academic year and 100% for 3 months during the summer. The graduate student will work on isolating and evaluating dirty snowpack signatures in satellite data, modeling soot transport events in support of IPY POLARCAT, and, in year 3, performing fully coupled model studies. All salaries and wages were estimated using UCI's academic and staff salary scales. Employee Benefits Fringe Benefits were estimated using the composite rates agreed upon by the University of California Office of the President and the DHHS Audit Agency, the Cognizant Audit Agency for the University of California. Benefit rates used in this proposal are: Faculty - summer - 12.7% Academic (Specialist) - 17% Student employees - summer - 3% Student employees - academic year - 1.3% Fees are requested for one nonresident student for the duration of the project. University of California policy requires award payment of fees for any student with more than 25% support from a grant. In the first year, $26,690 is requested for non-resident fees and tuition, $29,258 in the second year. It is anticipated that the graduate student will advance to candidacy at the beginning of the third year. University policy provides a 75% reduction in nonresident tuition post advancement. Therefore, fees and tuition are reduced to $16,748 in year 3. Fees and tuition are excluded from indirect cost assessment. Equipment Equipment funds are requested for the first year only for one dual Opteron, dual core workstation at $6,000. This workstation will include adequate RAID'ed disk space (1 TB) for the graduate student to store and analyze satellite MODIS, MISR, and AMSR-E datasets. Travel Domestic: Round-trip travel at $1500 per trip is requested for the PI and Postdoc (year 1) or graduate student (years 2 and 3) to travel to national meetings (primarily AGU) to present results. Each trip includes roundtrip travel from Irvine to San Francisco or the East Coast, one-week hotel and per diem. These trips include estimated conference registration, abstract submission fees, RT airfare, lodging, meals and ground transportation. Travel estimates are based on historical usage. International: One round-trip at $3000 is requested for the PI to travel to Norway in Year 2 to participate in the IPY POLARCAT team meeting activities and to present results. This trip includes roundtrip travel from Irvine to Oslo, one-week hotel and per diem. These trips include estimated conference registration, abstract submission fees, RT airfare, lodging, meals and ground transportation. Travel estimates are based on historical usage. Other Direct Costs Charges for journals, photocopying, long distance phone, fax and postage charges pursuant to this project are requested each year. Included in these expenses are long-distance charges for usage directly related to the project, such as communication with colleagues, journals, and vendors. Photocopying of research materials including publications and results of this sponsored research project are requested as well as mail and shipping for materials related to this project. Support is requested in years 2 and 3 for publication costs pursuant to this project, which include utilization of expensive color figures. Costs were estimated according to historical usage. Indirect Costs Facilities and Administrative costs were estimated in accordance with UCI's approved indirect cost rate agreement. The indirect cost rate of 52.5% of MTDC effective 7/1/05 was based upon the nature and location of the work proposed. Graduate student fees and tuition and equipment are excluded from indirect cost assessment. UCI's indirect cost rate agreement was approved by DHHS, the Federal Cognizant Audit Agency for UCI on 12/5/01. \end{verbatim} \newpage \section{Facilities, Equipment, and Other Resources} \subsection{Computational Resources} \setcounter{page}{1} \thispagestyle{empty} Our ANS project is well-situated to take advantage of UCI's fast computing capabilities previously funded by NSF and other agencies. PI~Zender directs the Earth System Modeling Facility (ESMF), an NSF-supported MRI facility dedicated to coupled global climate, chemistry, and biogeochemistry simulations. The ESMF flagship machines is an 88-CPU Power4+ IBM supercomputer with 192\,\GB\ RAM and 16\,\TB\ of RAID storage. In Spring 2006, ESMF anticipates acquiring a new Beowulf cluster comprising approximately twenty two-way dual core Opteron nodes (80~CPUs) and about 5\,\TB\ of RAID storage. Since this ANS proposal is squarely fits the ESMF mission, the ESMF will host the primary modeling development and shorter simulations. Once this project is funded, we will request supplementary time at NCAR for long production simulations of the fully coupled CCSM/SNICAR code. \newpage \section{Acronyms and Abbreviations}\label{sxn:abb} \setcounter{page}{1} \thispagestyle{empty} %\begin{longtable}{ >{\raggedright}p{7.0em}<{} >{\raggedright}p{8.0em}<{} } \begin{longtable}{ r >{\raggedright}p{25.0em}<{} l } & \kill % NB: longtable requires caption as table entry \caption[Acronyms and Abbreviations]{\textbf{Acronyms and Abbreviations}% \label{tbl:abb}} \\ \hline\hline \rule{0.0ex}{\hlntblhdrskp}% Abbreviation & Description & \\[0.0ex] \hline \rule{0.0ex}{\hlntblntrskp}% \endfirsthead % Lines between and \endfirsthead appear at top of table \caption[]{(continued)} \\ % Set label for following pages Abbreviation & Description & \\[0.0ex] \hline \rule{0.0ex}{\hlntblntrskp}% \endhead % Previous block appears at top of every page \endlastfoot % Previous block appears at end of table AAA & Arctic Absorbing Aerosol & \\ AMSR-E & Advanced Microwave Scanning Radiometer (satellite instrument) & \\ ANS & Arctic Natural Sciences & \\ AOMIP & Arctic Ocean Model Intercomparison Project & \\ AR4 & Fourth Assessment Report & \\ ARF & Aerosol Radiative Forcing & \\ ATSR & Along Track Scanning Radiometer and Microwave Sounder & \\ AVIRIS & Airborne Visible/Infrared Imaging Spectrometer & \\ BC & Black Carbon (light-absorbing component of carbonaceous aerosol) & \\ BRDF & Bi-directional Reflectance Distribution Function & \\ CAM & Community Atmosphere Model & \\ CCSM & Community Climate System Model & \\ CFEP & Center for Educational Partnerships & \\ CICE & Los Alamos sea-ice model & \\ CLM & Community Land Model & \\ CRM & Column Radiation Model & \\ DEAD & Dust Entrainment And Deposition Model & \\ EMA & Effective Medium Approximation & \\ ESM & Earth System Model & \\ ESMF & Earth System Modeling Facility & \\ ESS & Earth System Science (Department) & \\ FOCUS & Faculty Outreach Collaborations Uniting Scientists, Students and Schools & \\ GCM & General Circulation Model & \\ GFED & Global Fire Emissions Database & \\ GHG & Greenhouse Gas & \\ GLIDE & General Land Ice Dynamic Elements (core of GLIMMER) & \\ GLIMMER & Ice Sheet Model of Payne et~al. & \\ GSFC & Goddard Space Flight Center & \\ IPCC & Intergovernmental Panel on Climate Change & \\ IPY & International Polar Year & \\ ITD & Ice Thickness Distribution & \\ LANL & Los Alamos National Laboratory & \\ LGM & Last Glacial Maximum & \\ MISR & Multi-angle Imaging Spectro-Radiometer (satellite instrument) & \\ MODIS & Moderate Resolution Imaging Spectroradiometer (satellite instrument) & \\ NASA & National Aeronautic and Space Administration & \\ NCAR & National Center for Atmospheric Research & \\ NCO & netCDF Operators & \\ NILU & Norwegian Institute for Air Research & \\ NIR & Near InfraRed & \\ NSIDC & National Snow and Ice Data Center & \\ OC & Organic Carbon & \\ OLLI & Osher Lifelong Learning Institute & \\ OPAC & Optical Properties of Aerosols and Clouds & \\ PI & Principle Investigator & \\ POLARCAT & POLar study using Aircraft, Remote sensing, surface measurements and modeling of Climate, chemistry, Aerosols and Transport (IPY project) & \\ RT & Radiative Transfer & \\ SEI & Science and Engineering Informatics & \\ SEM & Scanning Electron Microscopy & \\ SGER & Small Grant for Exploratory Research & \\ SNICAR & SNow, ICe, and Aerosol Radiative model & \\ SNTHERM & Snow Melt model & \\ SOM & Slab Ocean Model & \\ SSA & Specific Surface Area & \\ SWE & Snow Water Equivalent & \\ TG & Temperature Gradient & \\ UCI & University of California, Irvine & \\ \end{longtable} % end tbl:abb \newpage \section{Project-Wide Combined Collaborator and Advisor List}\label{sxn:prs_lst} \setcounter{page}{1} \thispagestyle{empty} % Currently includes: Zender All Personnel Associated with Proposal, Collaborators and Co-Editors of Project Senior Personnel, their Post-docs, and their Thesis Advisors: \begin{enumerate*} \item[] Ammann, C.~A. (NCAR) \item[] Bian,~H. (NASA/UMBC) \item[] Bonan, G.~B. (NCAR) \item[] Busacca,~A. (WSU) \item[] Cakmur,~R. (NASA GISS) \item[] Colarco,~P. (GSFC) \item[] Collins, W.~D. (NCAR) \item[] Cooper, W.~A. (NCAR) \item[] Famiglietti,~J. (UCI) \item[] Gaylord,~D. (WSU) \item[] Grini,~A. (U.~Oslo) \item[] Jenks,~S. (UCI) \item[] Khalsa, S.~J.~S. (NSIDC) \item[] Kiehl, J.~T. (NCAR) \item[] Kuester,~F. (UCI) \item[] Levin,~Z. (TAU, Israel) \item[] Mahowald, N.~M. (NCAR) \item[] Miller,~R. (NASA GISS) \item[] Moore, J.~K. (UCI) \item[] Okin,~G. (U.~Virginia) \item[] Painter,~T. (NSIDC) \item[] Pajarola,~R. (UCI) \item[] Papadopoulos,~P. (UCI) \item[] Rasch, P.~J. (NCAR) \item[] Tegen,~I. (IfT, Germany) \item[] Thomas, G.~T. (CU) \item[] Torres,~O. (NASA GSFC) \item[] Valero, F.~P.~J. (Scripps) \item[] Yu,~S. (Duke) \end{enumerate*} \newpage List is Alphabetical by Surname. \begin{enumerate*} \item[] Collaborators of Zender: \begin{enumerate*} \item[] Ammann, C.~A. (NCAR) \item[] Bian,~H. (NASA/UMBC) \item[] Bonan, G.~B. (NCAR) \item[] Busacca,~A. (WSU) \item[] Cakmur,~R. (NASA GISS) \item[] Colarco,~P. (GSFC) \item[] Collins, W.~D. (NCAR) \item[] Cooper, W.~A. (NCAR) \item[] Famiglietti,~J. (UCI) \item[] Gaylord,~D. (WSU) \item[] Grini,~A. (U.~Oslo) \item[] Jenks,~S. (UCI) \item[] Khalsa, S.~J.~S. (NSIDC) \item[] Kiehl, J.~T. (NCAR) \item[] Kuester,~F. (UCI) \item[] Levin,~Z. (TAU, Israel) \item[] Mahowald, N.~M. (NCAR) \item[] Miller,~R. (NASA GISS) \item[] Moore, J.~K. (UCI) \item[] Okin,~G. (U.~Virginia) \item[] Painter,~T. (NSIDC) \item[] Pajarola,~R. (UCI) \item[] Papadopoulos,~P. (UCI) \item[] Rasch, P.~J. (NCAR) \item[] Tegen,~I. (IfT, Germany) \item[] Thomas, G.~T. (CU) \item[] Torres,~O. (NASA GSFC) \item[] Valero, F.~P.~J. (Scripps) \item[] Yu,~S. (Duke) \end{enumerate*} \item[] Collaborators of Co-PI Dr.~Who: \begin{enumerate*} \item Hi \end{enumerate*} \end{enumerate*} \newpage \subsection{Supplementary Documents}\label{sxn:spl_doc} \setcounter{page}{1} \thispagestyle{empty} \begin{verbatim} 1. ${DATA}/prp/prp_ans/prp_ans_ltr_warren.pdf 2. ${DATA}/prp/prp_ans/prp_ans_ltr_lanl.pdf 3. ${DATA}/prp/prp_ans/prp_ans_ltr_mahowald.pdf 4. ${DATA}/prp/prp_ans/prp_ans_ltr_polarcat.pdf 5. ${DATA}/prp/prp_ans/prp_ans_clb.pdf 6. ${DATA}/prp/prp_ans/prp_ans_abb.pdf \end{verbatim} \clearpage \csznote{ % Text snippets for future proposals... } % end csznote \csznote{ # Dozier's pubs http://www.snow.ucsb.edu/REASoN-Publications.htm # Email distribution list Siri Jodha Singh Khalsa , Tom Painter , Jim Randerson , Mark Flanner , Jorge Talamantes , Charlie Zender % Letter sent to Bill Lipscomb 20051215: Dear Bill, Thanks for discussing collaborations on the phone yesterday. As I said, I am following up with a one page summary of our proposed NSF project, and some sample text should you choose to endorse our collaboration more formally with NSF. Please feel free to ruthlessly alter/edit! The below text is meant to save you time, not to constrain you or to obligate you to do work you're not already interested in. If you choose to write such a letter, it would be most helpful if you could e-mail me a PDF on your letterhead by Wednesday morning 12/21, although a simple e-mail will suffice if you are too busy for the deluxe PDF option. Thanks! Charlie "Dear Charlie, As one of the principle developers of the Los Alamos sea-ice model (CICE), a primary component of the National Center for Atmospheric Research (NCAR) Community Climate System Model (CCSM), I have a great interest in understanding and correctly representing the coupled cryosphere in Earth System Models. Your proposed NSF project "Dirty Snow on Land, Glaciers, and Sea-Ice: Understanding Arctic Absorbing Aerosol Forcing and Feedbacks" would continue to advance the state-of-the-art in snowpack radiative transfer which provides an important upper boundary condition in the cryosphere. We are happy to collaborate with you and Mark Flanner in two related areas during the NSF project. First, we are implementing a multi-layer snow thermodynamic model to improve the current single layer snow model in CICE. We will provide source code access and advice on technical details so that your project can merge SNICAR physics for snowpack aging and radiative processes into CICE. Second, we are developing an interactive ice sheet component, based on GLIMMER, for the CCSM. We envisage using the a form of the Community Land Model (CLM) and its snow model with SNICAR as the upper boundary condition. The integrated treatment of clean and dirty snow in the cryosphere which will result from our complementary efforts will provide a state-of-the-art toolkit with which to study the scientific questions posed in your proposal." % Letter sent to Steve Warren 20051215: Dear Steve, As I said, I am following up with some sample text should you choose to endorse our collaboration more formally with NSF. Please feel free to ruthlessly alter/edit! The below text is meant to save you time, not to constrain you or to obligate you to do work you're not already interested in. If you choose to write such a letter, it would be most helpful if you could e-mail me a PDF on your letterhead by Wednesday morning 12/21, although a simple e-mail will suffice if you are too busy for the deluxe PDF option. P.S. As I think you already know, SNICAR uses optical properties based on your and Grenfell's work on the surface-area-to-volume Mie approximation. Thanks again for developing this useful approximation! Charlie "Dear Charlie, As you know I have a great interest in understanding snow, ice, and Arctic aerosol properties and I approach these topics from a combination of laboratory and in situ measurements combined with radiative transfer modeling. Your proposed NSF project "Dirty Snow on Land, Glaciers, and Sea-Ice: Understanding Arctic Absorbing Aerosol Forcing and Feedbacks" would continue to advance the state-of-the-art in snowpack radiative transfer which provides an important upper boundary condition in the cryosphere. Moreover, your coupled snowpack aging and radiation model, SNICAR, is a very suitable tool for hind-casting the snowpack BC concentrations in Greenland which our group continues to measure. Our newly proposed NSF project "title goes here" will fund us to measure BC from snow samples collected at all Koni Steffen's AWS sites, which will provide good coverage of the Greenland ice sheet. We will be pleased to provide your project with these measurements and to collaborate on their interpretation. The synergy between our measurements and your modeling project to understand and represent Arctic snow-ice-aerosol feedbacks is exciting." % Letter sent to Phil Rasch 20051215: Hi Phil, We are writing an NSF Arctic proposal called "Dirty Snow on Land, Glaciers, and Sea-Ice: Understanding Arctic Absorbing Aerosol Forcing and Feedbacks". The proposal would fund a grad student, and a one year postdoc which Mark Flanner can use as bridge funding. The goal is to continue studying dirty snow, both BC and dust, and to move to the next level of feedbacks by really including sea-ice and glacier effects. We think we have land pretty well covered, and will work with Bill Lipscomb (and probably Bruce Briegleb, although I'm a bit unsure where that stands) on the sea-ice and glacier coupling to SNICAR. Mark is beginning to write a boreal soot manuscript on which you will be offered co-authorship. Not making you an official co-author on the study Mark presented at AGU was an oversight! Although we fixed it in the title slide your name was not in the program. Apologies. Back to the present. This proposal is due Thursday 12/22 (I got an extension to improve the collaboratory aspects of the proposal). This proposal will rely heavily on continued use of your BC aerosol model in CAM. I am verifying that we can rely on you for this, and that I can mention this in the proposal's collaborative portion. Is there anything I should know or that you'd like to add? Attached is the current one-page summary FYI. Thanks, Charlie % Letter sent to Natalie Mahowald 20051215: Hi Natalie, We are writing an NSF Arctic proposal called "Dirty Snow on Land, Glaciers, and Sea-Ice: Understanding Arctic Absorbing Aerosol Forcing and Feedbacks". The proposal would fund a grad student, and a one year postdoc which Mark Flanner can use as bridge funding. The goal is to continue studying dirty snow, both BC and dust, and to move to the next level of feedbacks by really including sea-ice and glacier effects. We think we have dirty land-snow pretty well covered, and will work with Bill Lipscomb (and probably Bruce Briegleb, although I'm a bit unsure where that stands) on the sea-ice and glacier coupling to SNICAR. This proposal is due Thursday 12/22 (I got an extension to improve the collaborative aspects of the proposal). This proposal will rely heavily on continued use of Phil's BC aerosol model in CAM and on your dust code, if you agree. May we rely on you for code access and advice as has been our habit? It would be especially cool to include your glaciogenic dust sources (pun intended). Of course, we'll offer you co-authorship on any of the dust-related results. If you agree, may I mention this in the proposal's collaborative portion? Is there anything I should know or that you'd like to add? Attached is the current one-page summary FYI. Thanks, Charlie } % end csznote \csznote{ % Usage: Place usage here at end of file so comment character % not needed cd ~/prp_ans-1.0;make -W prp_ans.tex prp_ans.dvi prp_ans.ps prp_ans.pdf prp_ans.txt;cd - scp ${HOME}/prp_ans-1.0/prp_ans.dvi ${DATA}/ps/prp_ans.pdf ${DATA}/ps/prp_ans.ps ${HOME}/prp_ans-1.0/prp_ans.tex ${HOME}/prp_ans-1.0/prp_ans.txt dust.ess.uci.edu:/var/www/html/prp/prp_ans-1.0 scp ${DATA}/ps/bio_nsf.pdf ${DATA}/prp/prp_ans/prp_ans_cv_zender.pdf # NB: latex2html works well on prp_ans.tex latex2html -dir /var/www/html/prp/prp_ans prp_ans.tex # NB: tth chokes on prp_ans.tex cd ${HOME}/prp_ans-1.0;tth -a -Lprp_ans -p./:${TEXINPUTS}:${BIBINPUTS} < ${HOME}/prp_ans-1.0/prp_ans.tex > prp_ans.html scp prp_ans.html dust.ess.uci.edu:/var/www/html/prp/prp_ans-1.0 # NB: tex4ht works well on prp_ans.tex cd ${HOME}/prp_ans-1.0;htlatex prp_ans.tex scp prp_ans*.css prp_ans*.html dust.ess.uci.edu:/var/www/html/prp/prp_ans-1.0 # NB: tex4moz works well on prp_ans.tex cd ${HOME}/prp_ans-1.0;/usr/share/tex4ht/mzlatex prp_ans.tex scp prp_ans*.css prp_ans*.html prp_ans*.xml dust.ess.uci.edu:/var/www/html/prp/prp_ans-1.0 # Proposal preparation: # Divvy up master PDF into FastLane components pdftk A=${DATA}/ps/prp_ans.pdf cat A3 output ${DATA}/prp/prp_ans/prp_ans_smr.pdf pdftk A=${DATA}/ps/prp_ans.pdf cat A4-18 output ${DATA}/prp/prp_ans/prp_ans_dsc.pdf pdftk A=${DATA}/ps/prp_ans.pdf cat A19-23 output ${DATA}/prp/prp_ans/prp_ans_rfr.pdf pdftk A=${DATA}/ps/prp_ans.pdf cat A24-26 output ${DATA}/prp/prp_ans/prp_ans_bdg_jst.pdf pdftk A=${DATA}/ps/prp_ans.pdf cat A27 output ${DATA}/prp/prp_ans/prp_ans_fcl.pdf pdftk A=${DATA}/ps/prp_ans.pdf cat A28-29 output ${DATA}/prp/prp_ans/prp_ans_abb.pdf pdftk A=${DATA}/ps/prp_ans.pdf cat A30 output ${DATA}/prp/prp_ans/prp_ans_clb.pdf # Add letters of support in one command rather than loop # pdftk does not allow input file to be output file pdftk A=${DATA}/prp/prp_ans/prp_ans_ltr_warren.pdf \ B=${DATA}/prp/prp_ans/prp_ans_ltr_lanl.pdf \ C=${DATA}/prp/prp_ans/prp_ans_ltr_mahowald.pdf \ D=${DATA}/prp/prp_ans/prp_ans_ltr_polarcat.pdf \ cat A B C D output ${DATA}/prp/prp_ans/prp_ans_ltr.pdf # Add supplementary files in one command rather than loop # pdftk does not allow input file to be output file pdftk A=${DATA}/ps/prp_ans.pdf \ B=${DATA}/prp/prp_ans/prp_ans_ltr.pdf \ C=${DATA}/prp/prp_ans/prp_ans_cv_zender.pdf \ D=${DATA}/prp/prp_ans/prp_ans_cp_zender.pdf \ cat A B C D output ${DATA}/ps/prp_ans_fll.pdf } % end csznote on usage \end{document}