NSF Major Research Instrumentation Proposal February 28, 2004
Dr. Charles S. Zender Dr. Susan E. Trumbore
Department of Earth System Science IGPP Director
University of California at Irvine University of California at Irvine
1. Recommended proposal structure: (a) Research Activities: 10 pages (All) (b) Research Instrumentation and Needs: 2 pages (Wessel/Zender) (c) Impact of Infrastructure Projects: 2 pages (Prather) (d) Project and management Plans: 1 page (Zender) 2. Strategic questions: (a) Email Sandy Shore and Russ Kelz seeking support for ATM proposal (b) According to Schmoltner, “big project” threshold is $800-900k, we are at $876,887
Project Summary.
This proposal asks for funds to purchase a high-performance computer and storage system as the center of the UCI Earth System Modeling Facility (ESMF). The UCI ESMF will be devoted to the integration, synthesis, and analysis of large models and datasets required to advance fundamental understanding of the coupled physical climate, chemistry, and biogeochemical cycles of the Earth system. The ESMF represents a major enhancement in computational capability over the workstation-based and older shared-memory resources currently in use at UCI. It is tailored for the merging of ESM components (e.g., atmospheric chemistry, ocean biology, land hydrology) that normally consume the available computing resources of individual research projects. Although primarily a development facility for faculty and graduate researchers, the ESMF will produce the decade-long simulations of the coupled system that are needed for basic scientific studies. In terms of computing power, the ESMF fits between our individual PI and school-based resources and national modeling facilities (e.g., NCAR, SDSC), but in terms of development it is unique. The UCI ESMF fills a niche, allowing faculty, post-graduates, and students to work with high-performance computing in an environment where they control the computer and the code development.
Intellectual Merits. The process of coupling existing ESM components will identify some of the more important feedbacks and help identify aspects of the models that need attention or further development. The ESMF will also be used operationally (i.e., 100-year simulations) with component subsets of the full ESM, to address some outstanding questions about feedbacks anticipated with global change in the 21st century. Three subset-ESMs are planned for the early development: aerosol/gas chemistry, physical climate, and land+ocean biogeochemistry; catchment-scale riverine nutrient transport, Aeolian erosion, and marine ecosystems; and ocean and land surface memory effects on the predictability of the hydrologic cycle. These subset-ESMs will provide new insights (and problems) as we explore scientific questions that highlight the new capability of coupling these components. The initial subset-ESMs will produce new scientific results which will be pursued in operational studies with the full ESM. The ESMF will foster interdisciplinary collaborations between research groups which bring together an extremely wide breadth of scientific expertise in developing, analyzing and improving a coupled Earth System Model. These improvements will feed back into national modeling efforts. The ESMF will be maintained as a multi-user resource as part of the UCI Department of Earth System Science (ESS). The primary users of the ESMF are UCI ESS Professors Famiglietti, Magnusdottir, Moore, Prather, Primeau, Yu, and Zender, and their associated researchers, post-doctoral fellows, students and collaborators.
Broader Impacts. Development of ESMs is essential if the public is to understand the full implications of climate change and make the best informed political decisions regarding adaptation and mitigation. Integration of the human dimensions side of global change needs the details provided by an ESM rather than the reduced dimensions of integrated assessment models. Work at UCI’s ESMF will contribute to ESMs around the globe. As an academic center for ESM development, the ESMF will train the next generation of scientists and research staff in applying high-performance computing to large, coupled models and datasets in the Earth sciences. The ESS Department’s undergraduate curriculum encourages research projects, and the ESMF will sustain computing projects for several undergraduates. The ESMF will establish projects with UCI’s Center for Educational Partnerships (CFEP), which links Southern California agencies including school districts, parents, and community colleges. ESMF projects will also bring in students from a nearby Minority Serving Institution, California State University at Bakersfield (CSUB), giving these students access to environmental science models, data, and mentors unavailable at CSUB.
Faculty of the UCI Department of Earth System Science (ESS) propose to acquire a high performance computer and storage system to become the UCI Earth System Modeling Facility (ESMF). The UCI ESMF will be devoted to the integration, synthesis, and analysis of large geophysical models and datasets required to advance fundamental understanding of the coupled climate, chemical, and biogeochemical cycles of the Earth System. The ESMF represents a major enhancement in computational capability over the workstation-based and older shared-memory resources currently in use at UCI. It is tailored for the merging of ESM components (e.g., atmospheric chemistry, ocean biology, land hydrology) that—within each PI’s research—normally consume the available computing resources of individual research projects. The ESMF is optimized for providing highly cached, shared-memory computing for each component, with more relaxed requirements for parallel computation and communication between the different ESM components.
In our recent experience with national and international assessments of climate, key questions were unanswerable because of the lack of capability in integrating coupled climate system models, specifically the capability that the ESMF will provide. The scientific expertise in these individual areas of Earth system science exists within the ESS department, and its members are primary contributors to coupled Earth system models being assembled in many places. Yet there is no computational facility in the US where a group of researchers with the breadth and depth of the PIs can easily collaborate on science arising from the coupling of Earth system climate, chemical, and biogeochemical cycles. In terms of computing power, the proposed ESMF fits between our individual PI and school-based resources and national modeling facilities (e.g., NCAR, SDSC), but in terms of development it is unique. The UCI ESMF fills a niche, allowing faculty, post-graduates, and students to work with high-performance computing in an environment where they control the computer and the code development.
Although primarily a development facility for faculty and graduate researchers, the ESMF will produce the decade-long simulations of the coupled system that are needed for basic scientific studies. An important component of the ESMF is the dedicated disk space for data synthesis and analysis. In the next three years ESS researchers will require uninterrupted access to 30-100 terabytes (TB) of data to keep pace with the increases in temporal, spatial, and spectral resolutions of remote sensing measurements and coupled climate system models. This ESMF will allow ESS researchers to pursue data-intensive research utilizing the large geophysical datasets from current and next generation numerical models and satellite observations.
The principle researchers making use of the ESMF are the research groups of UCI Earth System Science Professors Magnusdottir, Moore, Prather, Primeau, Yu, Zender, and Famiglietti (joint with Civil Engineering): 1. Graduate students: J. Abatzoglou, A. Berg, S. Bortz, L.-M. Chen, M. Flanner, E. Kwon, Y. Hu, J.-W. Lee, D. Ryu, H. Syed, C.-C. Wang 2. Post-Docs: H. Bian, J. Hsu 3. Researchers: S. Holl, X. Ma, D. Newman, S. Tyler, X. Zhu 4. Undergraduates: S. Bamattre and others unnamed
The PIs have committed to coordinate modeling and experimental research around hypotheses and experiments testable using a unified Earth System Model (ESM). These faculty (see attached bio-sketches) have identified specific research projects (Section 2.2, below) that take advantage of the ESMF.
Other faculty anticipate using the ESMF in support their research activities in the near future. These include Professors Blake (Chemistry, Using ESMF simulations to evaluate aircraft-measured trace gases in support of field programs), Dabdub (Mechanical & Aerospace Engineering, Using ESMF to provide boundary conditions for nested regional air quality modeling), Druffel (ESS, Comparing ESMF simulations with long term coral 14C measurements), Friehe (MAE, Comparing measured and ESMF-modeled air sea exchange), Smecker-Hane (Physics and Astronomy), and Smyth (Information and Computer Science, Fitting hidden Markov models to observations and ESMF simulations of precipitation and storm tracks), and Trumbore (ESS, Understanding terrestrial C cycling using ESMF-simulations and soil 14C measurements). The ESMF will be maintained as a multi-user ESS departmental resource.
Section 2 describes the research opportunities the PIs will take advantage of with the ESMF. Section 3 describes why our current computational resources are inadequate to perform this research. Then we describe the computational requirements of an ESMF suitable for solving the coupled climate-biogeochemistry problems that are our highest priority. Section 4 describes the broader impacts of the ESMF on society, research training, and minority opportunities at UCI. Section 5 concludes with the plans for ESMF administration and future maintenance. Appendices contain the Budget Justification, two vendor quotes, a parts list, a Research Facilities statement, letters of support from CSUB, the UCI IGPP and VCR, and a list of Acronyms and Abbreviations.
None of the PIs has received previous NSF funding for instrumentation.
To enhance coupled Earth System research at UCI, ESS researchers have agreed to coordinate modeling and experimental research around hypotheses and experiments testable using a unified Earth System Model (ESM). Many ESM components will be based on the Community Climate System Model (CCSM) coordinated by the National Center for Atmospheric Research (NCAR) Blackmon et al. (2001). The ESMF will be used to study and extend climate related tracers and biogeochemical sub-models in the existing CCSM. This allows us to focus on the cutting edge science while building on a state-of-the-art climate system model. The depth of Earth System modeling at UCI has reached a critical mass where we can now leverage eachother’s expertise to make fundamental advances the coupled Earth system.
The PIs will use the ESMF to coordinate and direct their autonomous modeling research into a fully coupled Earth System Model effort at a faster pace than can be done through collaborating through external groups such as Community Climate System Model (CCSM) working groups (WGs) (Blackmon et al., 2001). These WGs decide where improvements and extensions are needed in the CCSM, coordinate scientific strategies to meet these needs, allocate computational resources to test competing physics packages, and evaluate the results. The ESMF is dedicated to the design and initial scientific application of next generation couplings (cf. Table 2) that are led by UCI. ESMF researchers are active members of several CCSM WGs, including the Atmospheric Model (AMWG), Biogeochemistry (BGCWG), Climate Variability (CVWG), and Land Model (LMWG) groups. The ESMF will facilitate progress and ease student involvement and training in UCI focus areas across these WGs, so the ESMF is complementary to CCSM WGs. Our improvements to ESM components will feed back into national modeling efforts such as CCSM and NASA’s Global Modeling Initiative (GMI) through our strong associations with these efforts. Thus any negative impact of the ESMF on existing national collaborations of the PIs is expected to be minimal and outweighed by the positive impact of ESMF contributions to these efforts.
Why are we not pursuing the research outlined below at a remote supercomputer location (such as NCAR or SDSC) open to proposals for large scale Earth system computational projects? We will address this question directly and state our reasons for establishing the ESMF: 1. UCI’s research niche: The ESMF is for focused research on coupled next generation biogeochemistry-climate problems. The PIs have unique strength in linking chemistry, biogeochemistry, and climate cycles together. 2. Scale of research: ESMF experiments will leverage control simulations from external centers. 3. Unfettered collaboration: The ESMF is devoted to quick turnaround (order one week) projects to enable university researchers to dynamically pursue interesting hypotheses without excessive administrative overhead. 4. Different resolution, more free parameters, and quicker turnaround 5. Responsiveness: Being on-site, the ESMF will be responsive to our specific science and teaching needs. Students and researchers will gain valuable experience implementing and analyzing their own experiments.
In contrast to external centers such as NCAR and SDSC, the ESMF is designed as a dedicated facility for a group of close-knit collaborators. Our mission is to develop the physical and chemical couplings required to address unanswered questions.
One way we will use the ESMF is studying climate and biogeochemical feedbacks of a particular component of the Earth System where the remainder of the system is represented by the standard NCAR CCSM. The CCSM provides an efficient, well-documented base upon which to pursue innovative modeling studies. CCSM components run in both stand-alone and coupled mode, and already contain with many biogeochemical, hydrologic, and radiative improvements developed and contributed by UCI researchers. Table 1 shows the relations between the ESMF models, research areas, and researchers.
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Most ESMF researchers will use CCSM components. Named in the table are researchers with interests in specific physical processes described by the model in question. Often these researchers are significant contributors to the associated models. Clearly UCI researchers have developed (and continue to maintain and develop) many of the core physical process models in the CCSM, including CRM, DEAD, OEM, and RTM.
Other UCI models, such as the UCI CTM and the Catchment RTM represent processes (prognostic tropospheric chemistry and catchment hydrology, respectively) that are currently not represented in most Earth System Models (including CCSM). Studying the effects of these new interactions will be a prime focus of the ESMF. Studies have already begun on the next generation of processes we intend to link to our Earth system model. These includes riverine erosion, and multi-component aerosols.
The ESMF is not designed for the fully-coupled multi-century simulations that many are pursuing at NCAR and SDSC. The ESMF will be a staging area for new collaborations among UCI faculty. As these projects reach maturity, we will wish to perform multi-decadal ensemble “production runs”, e.g., to analyze decadal variability or feedback behavior assuming a variety of initial conditions (e.g., IPCC scenarios). These type of ensemble runs are best suited to external locations or separately dedicated facilities such as NCAR, NERSC, or SDSC.
We believe ESS researchers would make important scientific advances on understanding the coupled Earth system if we could join our current models and expertise in a single effort. The ESMF will provide this opportunity. We are poised for leaps in understanding in at least three distinctive aspects of the coupled Earth system. Thus three subset-ESMs are planned for early development: (1) Aerosol/gas chemistry, physical climate, and land+ocean biogeochemistry; (2) Catchment-scale riverine nutrient transport, Aeolian erosion, and terrestrial plus marine ecosystems; (3) Ocean and land surface memory effects on the predictability of the hydrologic cycle. These subset-ESMs will provide new insights (and problems) as we explore scientific questions that highlight the new capability of coupling these components.
Sections 2.2.1- 2.2.3 describe these cross-cutting research projects which are of interest to many ESS researchers and constitute the primary collaborative foci of ESMF research. By the end of this proposal, these initial subset-ESMs will have produced new scientific results, and the development phase of the more complete ESM will have culminated with initial operational simulations. Sections 2.2.4- 2.2.9 describe related research projects that will benefit from the ESMF.
Aerosol chemistry plays important roles in the chemical composition and oxidative capacity of the atmosphere (e.g., Tie et al., 2001; Bian et al., 2003). Recent progress in fast inorganic aerosol thermodynamic equilibrium models (Zhang et al., 2000) makes on-line aerosol nucleation, coagulation, and growth feasible in coupled models. Aerosol indirect effects on climate (e.g., Kaufman and Fraser, 1997; Penner et al., 2001) are strongly sensitive to cloud distribution and formation processes (Zender and Kiehl, 1997; Ackerman et al., 2000) which are in turn sensitive to climate (Magnusdottir et al., 2003). Our work has focused on competition between photochemical and heterogeneous chemistry forcing due to aerosols that can lead to significant non-linear chemical feedbacks in and downwind of aerosol source regions (Bian and Zender, 2003; Bian et al., 2003). The UCI CTM is now a full tropospheric aerosol-chemistry model which accounts for aerosol chemistry and uptake on internally mixed dust, sea salt, sulfate, nitrate, and carbonaceous particles. We are poised to examine links between climate and chemistry modes (Prather, 1996; Yu and Mechoso, 1999), air quality, and the distribution and direct and indirect forcing of aerosols and greenhouse gases (Prather and Ehhalt, 2001; Yu et al., 2001; Collins et al., 2002).
ESS researchers are poised to address key questions concerning the roles of biogeochemically significant elements in ecosystem and climate change. Deposition of aeolian mineral dust is thought to provide limiting nutrients such as soluble iron to remote ecosystems such as the High Nutrient Low Chlorophyll (HNLC) regions of the Southern Oceans (Martin and Fitzwater, 1988; Martin, 1990). Our work has addressed preliminary questions concerning mineral dust transport and distributions (Zender et al., 2003a; Mahowald et al., 2002; Zender et al., 2003b), the impact of atmospheric iron deposition on ocean biogeochemical cycling (Moore et al., 2002), and aerosol effects on atmospheric chemistry (Bian and Zender, 2003). Future work will examine important land-ocean links by coupling models of atmospheric aerosol transport and riverine transport with a global ocean circulation/biogeochemistry model. We seek to quantify the roles of atmospheric and riverine inputs of Fe, Si, N, and P in driving ocean biogeochemistry and air-sea carbon dioxide flux at regional to global spatial scales. UCI is uniquely suited to the development and coupling of the models that simulate these important biogeochemical cycles.
The regional hydrological cycle in the North American is influenced locally by land surface conditions and remotely by ocean conditions (sea surface temperatures—SSTs) associated with El Niño-Southern Oscillation (ENSO) (Higgins et al., 1998; Farrara and Yu, 2003). Both the land surface and ocean possess long-term memory and are capable of prolonging the useful lead time for predictions of North American precipitation beyond the limit imposed by atmospheric internal variability. A better understanding of the role of ocean and land-surface memory in determining interannual variations in North American precipitation is essential for improving extended-range forecasts and water resources planning. The land surface memory is not independent of the ocean memory. Global-scale atmosphere-ocean interactions in the Tropical Pacific can alter the paths of the Pacific storm track and lead to land surface anomalies in snow cover, soil moisture and surface vegetation in the western United States. With the proposed ESMF, we will integrate the existing coupled atmosphere-ocean GCM (CGCM), land surface model, and regional-scale hydrological model at UCI to 1) assess the individual contributions of ocean and land-surface memory and their inter-dependence to the hydrological cycle in North America; and 2) downscale global ENSO predictions for regional-scale hydrology forecasts in the North America. For task (1), series of long-term (on the order of 100 years) experiments will be performed with the CGCM interactively coupled with/without the land surface model. For task (2), the regional hydrological model will be nested within the CGCM to downscale its extended-range ENSO predictions for ensemble seasonal forecasts of regional hydrological cycle in North America.
Our group addresses the hydrologic cycle at regional and global scales and its interaction within the climate system. An ESMF would provide the computational resources to explore links between the water cycle and other Earth system components that are at the forefront of climate and hydrologic research. Given the diverse group of faculty participants in this proposal, and in the ESS department as well, the ESMF would enable frontier analyses and simulations that are simply not possible at other institutions. Active and planned projects that fall into this category include 1) regional, high-resolution modeling of coupled ecologic-hydrologic change in the Sierra Nevadas in response to fire suppression and climate change; 2) coupled model simulations of the role of terrestrial water in land-ocean-atmosphere interaction, including coupled biogeochemical process; and 3) development of a high resolution global river network model that can assimilate satellite and ground-based data for improved climate and hydrologic prediction. These types of studies typically involve processing large amounts of satellite data, running higher resolution and longer time-scale climate simulations, and analysis of large volumes of model output. The computational resources of the ESMF would significantly improve our capabilities to conduct such cutting-edge studies.
Extra-tropical low-frequency variability in the atmosphere (beyond the two weeks that are considered predictable) exhibits strongly preferred spatial scales with almost no preferred time scales. Superposed on this “natural variability” are strong trends in certain climatic fields over the last couple of decades. Our research group tries to understand the forcing of the unprecedented trends. To accomplish this, the natural variability must be understood. In atmospheric general circulation models we examine forcing from below, in terms of SST and sea-ice anomalies, from above or the stratospheric polar vortex, and from within or internal forcing, which is substantial in the nonlinear chaotic system (Deser et al., 2003). Even while restricting ourselves to the atmosphere, we have been limited in computational resources. The ESMF would allow us to examine various forcing mechanisms in unison, and to do coupled atmosphere/ocean climate runs. We would process both high temporal resolution reanalysis data derived from observations, and model output to produce eddy statistics, which is essential in any regional study such as one focusing on winter rain in Southern California under different climate scenarios.
The research of our group focuses on biogeochemical cycling in the oceans and the linkages between ocean, land, and atmosphere. We developed a state of the art marine ecosystem model that includes explicit iron cycling and several key phytoplankton functional groups (Moore et al., 2002). Ongoing efforts seek to quantify the roles of these functional groups in elemental cycling in the oceans and to determine the influence of atmospheric and riverine inputs on ocean biogeochemistry. Coupled simulations of ocean circulation, ecosystem dynamics, and biogeochemical cycling have recently been developed in the context of the NCAR CCSM POP ocean model (Doney et al., 2003; Moore et al., 2003). Future research will involve coupling this global ocean biogeochemistry model with models of atmospheric and riverine transport from the land surface to the oceans for important biogeochemical elements (i.e. Iron, Silicon, Nitrogen. . . ). The model will also be used in coupled simulations of the global carbon cycle and climate change over decadal to century timescales. A second research track involves the analysis of oceanic remote sensing data from multiple sensors to study physical-biological interactions at regional to global spatial scales (Moore and Abbott, 2000). Both of these research foci would benefit greatly from the proposed ESMF. Currently ocean simulations must be done off-site and only small portions of the model output can be analyzed locally.
The demands for greatly expanded computational capability in atmospheric chemistry come from several research foci; two are specific to developing and improving the chemistry models; and the other to coupling across the Earth system. For one, the basic development of chemistry-transport models (CTMs) now depends on very high resolution models for comparison with detailed campaign measurements. UCI participated in the NASA TRACE-P campaign that collected a massive set of trace gas and aerosol measurements to study the export of ozone and other pollutants from east Asia (Wild et al., 2003). Use of this data set for model validation and testing requires the CTM to be run globally at a resolution of 1.8o latitude or better, for more than forty chemical species, and for a wide range of sensitivity tests to evaluate the uncertainty in emissions. A second focus is chemical-mode studies: it is becoming clear that perturbations to global atmospheric chemistry must be studied in terms of the long-lived perturbations generated by short-lived species (Prather, 1996; Wild and Prather, 2000; Olsen et al., 2000). This requires many sensitivity and perturbation calculations made at modest resolution (128 × 64 × 24) but for about a decade to identify the modes. A major computational need within these mode studies is to run the CTM in a linearized mode to generate the necessary sensitivity functions in order to derive a reduced dimension model. Such calculations require about 100 effective years to be calculated per single model year. The third focus involves the coupling of atmospheric chemistry with the physical climate and biogeochemical models. In this case, we anticipate that a modest resolution CTM needs to be run in parallel (with hourly data exchange) with the atmospheric GCM and the terrestrial biosphere model.
Our research focuses on understanding and predicting climate variations in the coupled atmosphere-ocean system, ranging from global-scale ENSO (Yu and Mechoso, 2001; Yu et al., 2002); to regional-scale monsoon rainfall cycle (Farrara and Yu, 2003). Our investigations rely greatly on long-term simulations performed with state-of-the-art, coupled atmosphere-ocean general circulation models (CGCMs) (Yu and Malone, 2001) that have global coverage and high horizontal resolutions. High-performance computational facility is essential to those model simulations. With the ESMF, we will address the following research issues: 1. Mechanisms of decadal ENSO modulation, with a focus on the inter-basin interactions between the tropical Pacific and Indian Oceans. 2. Contributions of ocean and land-surface memory and their inter-dependence to the year-to-year variations of summer precipitation in North America. 3. Air-sea exchanges in the coastal climate system of California. Atmospheric GCMs coupled with global oceanic models, coastal oceanic models, and land surface models with be used for this research.
Our research centers on predicting the radiative, dynamic, chemical, and biogeochemical interactions of natural and anthropogenic aerosols and trace gases on global scales. We create, maintain, and supply the research community an open source global mineral dust aerosol model used by UCI, NCAR, UCSB, U. Oslo, GFDL, and GSFC to study the radiative and biogeochemical interactions of mineral dust and climate. Since this work demands simulations of complex radiative-dynamical interactions over global domains, nearly all of our research and publications depends on supercomputer simulations. The research problems that our group will attack with the ESMF in the next three years are 1. Past, Present, and Future changes in climate and chemical composition due to direct and indirect radiative, chemical, and biogeochemical forcing of mineral dust aerosol. 2. Assessment of the relative magnitudes of global wind and riverine erosion to long range transport of terrigenic biogeochemically significant minerals. 3. Quantification of the roles of geomorphology, surface hydrology, mesoscale dynamics in making semi-arid landscapes vulnerable to wind erosion. Ensemble decadal offline transport and coupled global climate model numerical simulations to are required to advance research in all of these areas. The computational resources of the ESMF would allow us to discretize important aerosol radiative processes at finer temporal resolution and thus substantially reduce temporal aliasing errors in our simulations.
A major component of improving and testing models that simulate the Earth System is (1) the inclusion of the best understanding of physical, chemical and biological controls on surface fluxes of energy, water, and nutrients; and (2) testing of model predictions with appropriately scaled field measurements. Many UCI researchers make laboratory and field measurements of processes simulated by an ESM, and all of these participate in the UCI Institute of Geophysics and Planetary Physics (IGPP) branch. The UCI IGPP strongly endorses the founding of the ESMF (see attached letter of support from S. Trumbore) since the ESMF effort will complement IGPP process studies and production of comparison data sets.
Isotopic tracers are particularly useful in studying the global hydrological and carbon cycles because they integrate processes across a number of scales. They are a major emphasis of measurement programs at the UCI IGPP. The distribution of trace gases in the atmosphere as predicted by tracer transport models depends critically on the sources and sinks of those gases or their precursors at Earth’s surface. Transfers of trace gases, energy and water at the land-air and ocean-air interface are controlled by complex interactions of atmospheric turbulence, boundary-layer development, and biological and chemical modification of air in contact with the surface. Major advances in modeling earth surface fluxes will continue to come from insight into these processes gleaned from direct eddy covariance observations of exchange (IGPP researchers Friehe, Goulden, Saltzman), in studying the factors controlling the rate of production and consumption of specific trace gases at the land and ocean surfaces (Cicerone, Druffel, Famiglietti, Goulden, Reeburgh, Saltzman, Tyler, Trumbore), and reactive chemistry of those gases in the atmosphere (Dabdub, Finlayson-Pitts, Prather, Saltzman).
ESM evaluation will be facilitated by comparison of predicted spatial and temporal distributions with observations for a suite of trace gases in the atmosphere and ocean (Blake, Saltzman) and their isotopic composition (Druffel, Saltzman, Tyler, Trumbore). Specific examples where the confluence of process investigations, and large scale observations within the IGPP at UCI can make major contributions when combined with ESMF modeling are in the areas of (1) high latitude climate change feedbacks to terrestrial ecosystems (observations by Goulden, Trumbore and postdocs); (2) methane distribution in the atmosphere, especially why growth rates have changed in recent decades (Reeburgh, Tyler, Blake, and Rowland); (3) the carbon cycle (Druffel, Goulden, Prather, Southon, Trumbore); (4) links between hydrology, trace gas fluxes, and ocean circulation (Famiglietti, Goulden, Trumbore, Primeau, Yu); and (5) links between ocean circulation, dust inputs, and chemistry of surface waters (Druffel, Primeau, Yu, Zender). This synergy between modeling and observation programs is another unique aspect of the proposed ESMF at UCI. While the data produced by IGPP facilities (the stable isotope trace gas facility and UCI Keck Carbon Cycle AMS facility) are widely available to all researchers, the presence of the ESMF in the same institution will provide impetus for increased collaboration and planning.
The computational requirements of the coupled Earth system research described above are significant. At least three of the ten most powerful national and international supercomputing centers are primarily devoted to Earth system simulation (Meuer, 2002). The proposed UCI ESMF will run many of the same models heavily employed by these national and international centers.
Experience has shown that the best metric of computational requirements for our research is coupled model throughput in terms of simulated years per real (wall-clock) day. This “throughput” is always substantially less than the raw computational speed (i.e., floating point performance) of high performance computing systems because coupled Earth system models have global communication requirements (e.g., spectral transforms) which perform better with shared memory than with massively parallel systems (e.g., Beowulf clusters).
The prototypical research project using the ESMF requires decadal simulations of the coupled climate system in experimental configurations (e.g., various aerosol interactions turned “on” and “off”) and a control configuration (e.g., no aerosols). The computational requirements of the subset ESMs described in Sections 2.2.1- 2.2.3 are summarized in Table 2.
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
We expect our subset ESMs to run 2-4 times slower than the uncoupled CAM, and we need to simulate 60-120 years per experiment. We would like the ESMF turnaround time to be about one week per experiment. Thus we require an ESMF throughput of 150-300 years of CAM per week.
Storage requirements of ESMF researchers are driven by three main tasks: (1) Storage of model input (forcing) datasets; (2) Archival and analysis of model simulations; (3) Archival and analysis of observations, notably satellite-derived datasets. Accounting for these three needs, we estimate that ESMF investigators will require 30-100 TB storage within 3 years.
The past decade has seen a tremendous advancement in data assimilation techniques for atmospheric and oceanic data which are used to drive models. The met fields that the UCI CTM uses to run with in off-line mode typically are more than 30 GB per year, and we anticipate that 40 years need to be stored on-line at the ESMF. The NCEP/NCAR reanalysis data represent the most reliable and complete atmospheric dataset available (Kalnay, 1996). The entire NCEP reanalysis from 1948-2002 (T62L28, 6 hour resolution) is about 500 GB. Due to its large size, many ESMF researchers can access only specific portions of this dataset. This entire datasest will be stored on-line at the ESMF, improving our potential research while avoiding duplication.
A typical atmospheric circulation simulation might run for 40 years and generate a minimum of 22 gigabytes (GB) of data. A decadal coupled climate-chemistry-biogeochemistry scenario might generate 500 GB of data. When we compare the results of just two different such climate scenarios our raw data is already in the terabyte range.
ESMF researchers currently analyze, or will soon analyze, data from the following satellite instruments: Nimbus 7 and Earth Probe TOMS, NOAA 7-14 AVHRR, SeaStar SeaWiFS, Terra and Aqua CERES and MODIS, Terra MOPITT, and Aqua AMSR-E and GRACE. These instruments collect gigabytes of processed data every year. Analysis of climate variability on interannual to decadal scales therefore can require processing 100 GB or more of data. Cross-comparisons of measurements from different instruments (such as correlating and merging optical depth estimates from several radiometers) can demand more storage than personal workstations are currently able to provide. The multi-terabyte storage facility (a RAID disk farm) we are requesting will allow us to perform integrated research on these large datasets in ways that are unavailable to us at the present time.
Over the past few years there has been a substantial increase in the number of Earth system modeling groups at UCI, but no increase in centralized computational resources. In mid-1999 the UCI acquired a 32 CPU SGI O2000 computer (”krein”). The O2000 is shared by the whole School of Physical Sciences (SPS) and cannot be dedicated to ESM modeling. Even if krein were dedicated as an ESMF, its throughput and storage system are inadequate for our purposes by factors of about 10 and 100, respectively. Two ESS groups (Magnusdottir and Zender) use krein for uncoupled atmospheric simulations. Other ESS groups use their own workstations (Primeau, Prather) and/or external supercomputer facilities, chiefly GSFC (Famiglietti), NCAR (Magnusdottir, Moore, and Zender), and NPACI/SDSC (Yu). This Balkanization of our simulation environments impedes effective local collaboration. The local siting of the ESMF will allow ESS faculty, post-graduates, and students to control the computer administration and code development.
The scientific problems the ESMF will tackle determine the optimal topology of the computational components. Figure 1 shows the organization of computations, information flow, and storage in an ESM experiment.
In a typical experiment, five to seven computational nodes would be dedicated to specific component models (Atmosphere, Chemistry, Ocean, Land, and Sea Ice) that run independently of the other nodes. In practice, the load-balancing node would be used for additional throughput or higher resolution a particular model of interest, or for a new sub-model, e.g., catchment hydrology or aerosol thermodynamics. It is also possible to run two or more of the faster models (e.g., land and sea ice) on one node, freeing up nodes for models, such as CAM, that scale well with MPI. For experiments where the load balancing node would not improve throughput it would remain available for interactive work or use as a backup in case another node fails.
Figure 1 illustrates a coupling configuration where all models communicate indirectly via the flux coupler. The flux coupler handles the details of regridding fluxes among the various model resolutions using a mass and energy conserving algorithm. In the CCSM, fluxes are accumulated and then exchanged between models (nodes) about once per simulation hour to day. The bandwidth requirements for this rate of coupling are not extreme but still must be taken into account.
The proposed ESMF must be modular and expandable, meaning that the simulation throughput could be increased by adding multi-CPU nodes. Figure 1 shows why the optimal configuration for running an Earth System Model has between five and seven computational nodes. The administrative node handles overhead associated with system utilities, queue management, and interactive work. We conclude that the optimal design for the minimal ESMF is an eight node system. Two simulations running concurrently (e.g., control and experiment) would require, on average, 14 total nodes; three simulations would require 20, etc. The cost difference between highly distributed (Beowulf) and highly shared memory realization of this configuration are significant. There we give more detail on the rationale for our decision in the attached budget justification.
IBM Corporation, an independent and reputable provider of high-speed scientific computing environments, can meet all these design requirements. IBM suggests a 64 processor p655 Power4 system similar to, but (much) smaller than, the fastest production systems at NCAR, SDSC, and LLNL. Their quote includes eight 8-way shared memory nodes based on the upcoming 1.5 GHz Power4 chips, and GB Ethernet interconnects for each node. The theoretical peak performance of this machine is 384 GFlops, and we expect throughput of about 50 GFlops. The inter-node bandwidth demands of the coupled ESM are met by Gigabit Ethernet switches. More expensive, proprietary IBM interconnects are available but not deemed necessary for ESMF applications. The attached quote and parts list contain complete hardware details.
The ESMF must be capable of simulating 60-120 yr wk-1 of the coupled system (Section 2.4). IBM has converted CAM benchmarks performed on a 32-processor IBM p690 at Oak Ridge National Lab (ORNL) (Worley, 2002) to estimate the CAM throughput of their proposed p655-based ESMF. The ORNL p690 system, using 32 processors, simulates about 70 yr wk-1 of CAM (T42×L26, Eulerian dynamical core). IBM estimates that the proposed ESMF, using 6 compute nodes (48 processors) will simulate at least 120 yr wk-1 of CAM (about 180 yr wk-1 for the SLD dynamical core). This estimate accounts for the differences in interconnect and chip speeds between the ORNL and ESMF systems. Based on our experience, 120 yr wk-1 of CAM is equivalent to at least 60 yr wk-1 of a typical ESMF simulation from Table 2.
The models in Table 1 perform well on these IBM systems. Most ESS modelers are already familiar with the IBM run-time environments NCAR and SDSC, and using the similarly configured ESMF will be easy. Our selection of IBM Corporation equipment for the ESMF quote is influenced by the success of other Universities and National laboratories running coupled climate system models on shared memory IBM systems. However, the optimal price/performance/maintenance system for high performance scientific computing (HPSC) solutions for climate problems changes rapidly. If this MRI is funded, we will solicit updated and competing bids for systems from other vendors (e.g., SGI, HP, Beowulf) which can demonstrate comparable simulation throughput as the quoted IBM p655 system.
Selecting an external storage system is relatively easy since high quality, fault-tolerant storage systems are rapidly commoditizing. Western Scientific will provide an array of ten Tornado F4 3.2 TB IDE-SCSI devices—a total initial capacity of 32 TB. The largest logical partition on this system is about 3 TB. This is enough archival space for any one of our experiments. Western Scientific has worked closely with IBM to ensure that the their systems systems will mesh.
The ESMF will build upon an ongoing collaboration with the Department of Physics and Geology at California State University, Bakersfield (CSUB) coordinated by Dr. Jorge Talamantes and PI Zender. Faculty in the Physics and Geology Department (and possibly other departments such as Chemistry and Computer Science) at CSUB will provide mentorship for CSUB undergraduates interested in environmental science and modeling, geophysical data analysis and statistics, and HPC. These students will get computer accounts on the ESMF and advice from relevant ESMF researchers. These Cal State students have no comparable opportunity to explore research projects in Earth system science and HPC. CSUB has been designated by the Department of Education as a Minority Serving Institution, with an undergraduate Latino population of 36.6%, and an overall minority undergraduate population of 53% which closely mirrors the regional population. Thus our existing UCI-CSUB collaboration will evolve into a major outreach component of this proposal.
The ESMF will also establish links with UCI’s Center for Educational Partnerships (CFEP), the campus outreach coordinator. We will provide new opportunities for undergraduate research, teacher education, mentoring, and K-12 education in the geosciences with a perspective on treating the Earth as a coupled system. CFEP links Southern California agencies whose goals are to increase academic success in students. CFEP supports a large number of programs specifically designed for students in elementary, intermediate, and high schools. These programs assist students to become eligible or competitively eligible for UC schools. Through established programs and community outreach CFEP fosters community links to improve the preparation of all students, particularly those from under-served groups, for success in higher education.
One program, the UCI Summer Science Institute (SSI), is an excellent outreach to high school science teachers. The SSI offers professional development programs for teachers and additional learning opportunities or all students to ensure equal access to high quality, effective and motivating science education. ESS faculty have participated in the institute, and the funding of the ESMF would provide a focus for summer course with some hands-on experience for teachers to understand the complexities and rewards of modeling the Earth system. Another program to which the ESMF will contribute is the California State Summer School in Mathematics and Science (COSMOS), established by the state legislature in 1999. COSMOS motivates the most creative minds of the new generation of prospective scientists, engineers, and mathematicians to actively participate in the business and higher educational sectors of California. COSMOS serves California high school students in grades 9-12 who are gifted and talented in mathematics and science. The program, which offers a curriculum not traditionally offered in high school, exposes them to an intensive learning experience that will enhance their academic development and shape their educational and career goals. The curriculum was designed by UCI faculty in eight subject areas: astronomy, biological sciences, isotope and atmospheric chemistry, cognitive science, computer science, engineering, mathematics and physics. We will expand the atmospheric chemistry curriculum of COSMOS in the broader area of environmental modeling and analysis.
The ESMF will benefit investigators who are not directly involved in a number of ways. The ESMF will make extensive data sets, pre-processing, and storage capacity readily available to students and postdocs without current access to large computational facilities. This will allow these researchers to carry out interdisciplinary studies using coupled models, or to compare in situ measurements with larger global satellite and numerical model records. Such innovative studies will likely in some case greatly expand the scope of their advisors’ current work. ESS faculty plan to enhance the graduate curricula by developing a set of ESM modeling exercises to which a small portion of the ESMF resources would be dedicated. Finally, the existence of a terabyte-scale data processing and archiving facility will help build the scientific infrastructure of the ESS department. This will benefit the ESS department in future recruitment of graduate students, post-docs and faculty, and improve the competitiveness of our research proposals.
Campus-wide the capability of the ESMF will draw a range of non-geoscientists—such as engineers, social and political scientists, legal scholars, and others—who will want to address a wide range of technological and societal options for the future. As the capability of the ESM develops, and we are able to evaluate future impacts of global change, the interest in merging Earth system science with the human dimensions side will become stronger and benefit many academic programs on campus.
National and international assessments are the chief mechanism for providing policy-makers with consensus scientific input on the on the past and future of Earth’s climate. ESS researchers are key contributors to the IPCC and other assessments (Prather et al., 1995; Olson et al., 1997; Prather and Sausen, 1999; Prather and Ehhalt, 2001). Many UCI-based models (DEAD, RTM, and OEM, see Table 1) are already integrated in CCSM and will be used in NCAR’s upcoming IPCC simulations for 2006. Given the strong links between ESS researchers and development and use of models contributing to national and international climate assessments, it is clear that ESMF research will influence these assessments broadly and for many years.
The ESMF will be maintained as a multi-user resource as part of the UCI ESS Department. The ESMF Director will coordinate ESMF scientific use and promote outreach linkages. PI Zender will serve as ESMF director. Zender has 10+ years experience with HPC Earth system modeling. Dr. Frank J. Wessel of the UCI Research Computing Support (RCS) group will serve as ESMF technical director. Wessel has 15 years experience in HPC procurement, use, system administration, and management. Dr. Wessel will oversee the system administrator, coordinate vendor relations, and determine appropriate facility upgrade and maintenance strategies to ensure maximum ESMF uptime and security. The Director, technical Director, and three to five ESMF investigators will form an ESMF advisory committee responsible for facility oversight and planning. We will accommodate all reasonable requests UCI researchers who desire to use the ESMF for purposes consistent with this proposal. This includes the individual purchase of additional nodes on the ESMF using startup or project funds.
Based on experience and discussions with UCI Network and Computing Services (NACS) and the Research Computing Support (RCS) Group, we have budgeted 1/2-FTE Programmer Analyst II for ESMF system administration and support. The system Administrator will have primary responsibility for day-to-day maintenance of the ESMF, accounts, and software upgrades. Our request includes vendors’ fees which fully cover the first three years of maintenance. These maintenance costs go to the vendor and are in addition to the time spent by ESMF personnel.
The initial purchase of the ESMF covers the first three years of software and hardware maintenance and system administration costs. After year 3 (late 2006), we will cover these costs entirely with UCI funds. We estimate the annual costs of operation as $60,000 for the 0.5 FTE system administrator, $45,000 for hardware maintenance, and $25,000 for software maintenance (see attached quote from IBM). The total annual cost of operation is therefore about $130,000.
The ESMF will be the nucleus of future ESS and UCI HPC activity. The University is committed to the maintenance and support of the ESMF beyond the third year through a combination of income derived from a recharge facility, support from IGPP, possible integration into the support structure of a larger scientific HPC facility, and financial support from SPS and NACS. After three years, significant users (those consuming more than 1% of the ESMF resources) will be charged to help defray operational costs. A similar system has been used successfully for the UCI SGI Origin 2000, purchased in 1999. The UCI IGPP has committed $20,000 per year to the ESMF facility after year 3 (see attached letter of support from Sue Trumbore). Growth in the number of ESS modelers, compounded by the increasing complexity of their models, will strongly motivate us to maintain and enhance the ESMF.
Ackerman, A. S., O. B. Toon, D. E. Stevens, A. J. Heymsfield, V. Ramanathan and E. J. Welton, 2000: Reduction of tropical cloudiness by soot. Science, 288, 1042-1047.
Bian, H., M. Prather and T. Takemura, 2003: Tropospheric aerosol impacts on trace-gas budgets through photolysis. J. Geophys. Res., 108(D8), 4242, doi:10.1029/2002JD002743.
Bian, H. and C. S. Zender, 2003: Mineral dust and global tropospheric chemistry: Relative roles of photolysis and heterogeneous uptake. J. Geophys. Res., 108(D21), 4672, doi:10.1029/2002JD003143.
Blackmon, M., B. Boville, F. Bryan, R. Dickinson, P. Gent, J. Kiehl, R. Moritz, D. Randall, J. Shukla, S. Solomon, G. Bonan, S. Doney, I. Fung, J. Hack, E. Hunke, J. Hurrell, J. Kutzbach, J. Meehl, B. Otto-Bliesner, R. Saravanan, E. K. Schneider, L. Sloan, M. Spall, K. Taylor, J. Tribbia and W. Washington, 2001: The community climate system model. Bull. Am. Meteorol. Soc., 82(11), 2357-2376.
Collins, W. D., P. J. Rasch, B. E. Eaton, D. W. Fillmore, J. T. Kiehl, C. T. Beck and C. S. Zender, 2002: Simulation of aerosol distributions and radiative forcing for INDOEX: Regional climate impacts. J. Geophys. Res., 107(D19), 8028, doi:10.1029/2000JD000032.
Deser, C., G. Magnusdottir, R. Saravanan and A. Phillips, 2003: The effect of North Atlantic SST and sea-ice anomalies on the winter circulation in CCM3. Part II: Direct and indirect components of the response. Submitted to J. Clim., Available from ftp://ess1.ess.uci.edu/pub/gudrun/NAO/Part2.
Doney, S. C., K. Lindsay and J. K. Moore, 2003: Global ocean carbon cycle modeling. in M. Fasham, J. Field, T. Platt and B. Zeitschel, editors, Ocean Biogeochemistry: A JGOFS synthesis. Springer-Verlag, New York, NY.
Farrara, J. D. and J.-Y. Yu, 2003: Interannual variations in the North American monsoon and SST anomalies: A general circulation model study. In Press in J. Clim.
Higgins, R. W., K. C. Mo and Y. Yao, 1998: Interannual variability of the United States summer precipitation regime with emphasis on the southwestern monsoon. J. Clim., 11, 2582-2606.
Kalnay, E., 1996: The NCEP/NCAR 40-year reanalysis project. Bull. Am. Meteorol. Soc., 77(3), 437-471.
Kaufman, Y. J. and R. S. Fraser, 1997: The effect of smoke particles on clouds and climate forcing. Science, 277, 1636-1639.
Magnusdottir, G., C. Deser and R. Saravanan, 2003: The effect of North Atlantic SST and sea-ice anomalies on the winter circulation in CCM3. Part I: Main features and storm-track characteristics of the response. Submitted to J. Clim., Available from ftp://ess1.ess.uci.edu/pub/gudrun/NAO/Part1.
Mahowald, N. M., C. S. Zender, C. Luo, D. Savoie, O. Torres and J. del Corral, 2002: Understanding the 30 year Barbados desert dust record. J. Geophys. Res., 107(D21), 4561, doi:10.1029/2002JD002097.
Martin, J. H., 1990: Glacial-interglacial CO2 change: The iron hypothesis. Paleoceanography, 5, 1-13.
Martin, J. H. and S. E. Fitzwater, 1988: Iron-deficiency limits phytoplankton growth in the northeast Pacific subarctic. Nature, 331(6154), 341-343.
Meuer, H., 2002: TOP500 supercomputer sites. “ http://www.top500.org/list/2002/11”.
Moore, J. K. and M. R. Abbott, 2000: Phytoplankton chlorophyll distributions and primary production in the Southern Ocean. J. Geophys. Res., 108(28), 709-728,722.
Moore, J. K., S. Doney and K. Lindsay, 2003: Upper ocean ecosystem dynamics and iron cycling in a global 3d model. Manuscript in preparation, Available from http://www.ess.uci.edu/ jkmoore.
Moore, J. K., S. C. Doney, D. M. Glover and I. Y. Fung, 2002: Iron cycling and nutrient-limitation patterns in surface waters of the World Ocean. Deep-Sea Res. II, 49, 463-507.
Olsen, S. C., B. J. Hannegan, X. Zhu and M. J. Prather, 2000: Evaluating ozone depletion from very short-lived halocarbons. Geophys. Res. Lett., 27, 1475-1478.
Olson, J., M. Prather, T. Berntsen, G. Carmichael, R. Chatfield, P. Connell, R. Derwent, L. Horowitz, S. Jin, M. Kanakidou, P. Kasibhatla, R. Kotamarthi, M. Kuhn, K. Law, J. Penner, L. Perliski, S. Sillman, F. Stordal, A. Thompson and O. Wild, 1997: Results from the Intergovernmental Panel on Climatic Change Photochemical Model Intercomparison (PhotoComp). J. Geophys. Res., 102(D5), 5979-5991.
Penner, J. E., M. Andreae, H. Annegarn, L. Barrie, J. Feichter, D. Hegg, A. Jayaraman, R. Leaitch, D. Murphy, J. Nganga and G. Pitari, 2001: Aerosols, their direct and indirect effects. in J. T. Houghton, Y. Ding, D. J. Griggs, M. Noguer, P. J. van der Linden, X. Dai, K. Maskell and C. A. Johnson, editors, Climate Change 2001: The Scientific Basis. Contribution of Working Group I to the Third Assessment Report of the Intergovernmental Panel on Climate Change, Chap. 5, pp. 291-336. Cambridge Univ. Press, Cambridge, UK, and New York, NY, USA.
Prather, M., R. Derwent, D. Ehhalt, P. Fraser, E. Sanhueza and X. Zhou, 1995: Other trace gases and atmospheric chemistry. in J. T. Houghton, L. G. M. Filho, J. P. Bruce, H. Lee, B. A. Callander and E. F. Haites, editors, Climate Change 1994: Radiative Forcing of Climate and an Evaluation of the IPCC IS92 Emission Scenarios, Chap. 2, pp. 73-126. Cambridge Univ. Press, Cambridge, UK, and New York, NY, USA.
Prather, M. and D. Ehhalt, 2001: Atmospheric chemistry and greenhouse gases. in J. T. Houghton, Y. Ding, D. J. Griggs, M. Noguer, P. J. van der Linden, X. Dai, K. Maskell and C. A. Johnson, editors, Climate Change 2001: The Scientific Basis. Contribution of Working Group I to the Third Assessment Report of the Intergovernmental Panel on Climate Change, Chap. 4, pp. 239-287. Cambridge Univ. Press, Cambridge, UK, and New York, NY, USA.
Prather, M. and R. Sausen, 1999: Potential climate change from aviation. in J. E. Penner, D. H. Lister, D. J. Griggs, D. J. Dokken and M. McFarland, editors, Aviation and the Global Atmosphere, Chap. 6, pp. 185-215. Cambridge Univ. Press, Cambridge, UK, and New York, NY, USA.
Prather, M. J., 1996: Natural modes and time scales in atmospheric chemistry: theory, GWPs for CH4 and CO, and runaway growth. Geophys. Res. Lett., 23, 2597-2600.
Tie, X., G. Brasseur, L. Emmons, L. Horowitz and D. Kinnison, 2001: Effects of aerosols on tropospheric oxidants: A global model study. J. Geophys. Res., 106(D19), 22931-22964.
Wild, O. and M. J. Prather, 2000: Excitation of the primary tropospheric chemical mode in a global three-dimensional model. J. Geophys. Res., 105, 24647-24660.
Wild, O., J. K. Sundet, M. J. Prather, I. S. A. Isaksen, H. Akimoto, E. V. Browell and S. J. Oltmans, 2003: CTM ozone simulations for spring 2001 over the Western Pacific: Comparisons with TRACE-P lidar, ozonesondes and TOMS columns. Submitted to J. Geophys. Res.
Worley, P., 2002: Community Atmospheric Model (CAM) performance benchmarks. “ http://www.csm.ornl.gov/~worley/evaluation/CAM/benchmark.html”.
Yu, J.-Y. and R. C. Malone, 2001: A global high-resolution UCLA-POP coupled atmosphere-ocean GCM for assessing climate and variability in the ocean. in Proceedings of the Symposium on Climate Variability, the Ocean, and Social Impacts, 14-18 Jan., Albuquerque, New Mexico, pp. 160-162.
Yu, J.-Y. and C. Mechoso, 1999: Links between annual variations of Peruvian stratocumulus clouds and of SSTs in the eastern equatorial Pacific. J. Clim., 12, 3305-3318.
Yu, J.-Y. and C. R. Mechoso, 2001: A coupled atmosphere-ocean GCM study of the ENSO cycle. J. Clim., 14, 2329-2350.
Yu, J.-Y., C. R. Mechoso, J. C. McWilliams and A. Arakawa, 2002: Impacts of the Indian Ocean on the ENSO cycle. Geophys. Res. Lett., 29, 461-464.
Yu, S., C. S. Zender and V. K. Saxena, 2001: Direct radiative forcing and atmospheric absorption by boundary layer aerosol in the southeastern US: model estimates on the basis of new observations. Atmos. Environ., 35(23), 3967-3977.
Zender, C. S., H. Bian and D. Newman, 2003a: Mineral Dust Entrainment And Deposition (DEAD) model: Description and 1990s dust climatology. J. Geophys. Res., 108(D14), 4416, doi:10.1029/2002JD002775.
Zender, C. S. and J. T. Kiehl, 1997: Sensitivity of climate simulations to radiative effects of tropical anvil structure. J. Geophys. Res., 102(D20), 23793-23803.
Zender, C. S., D. J. Newman and O. Torres, 2003b: Spatial heterogeneity in aeolian erodibility: Uniform, topographic, geomorphic, and hydrologic hypotheses. J. Geophys. Res., 108(D17), 4543, doi:10.1029/2002JD003039.
Zhang, Y., C. Seigneur, J. H. Seinfeld, M. Jacobson, S. L. Clegg and F. S. Binkowski, 2000: A comparative review of inorganic aerosol thermodynamic equilibrium modules: Similarities, differences, and their likely causes. Atmos. Environ., 34, 117-137.