Person:

Long, Michael

We describe the Harvard–NASA Emission Component version 1.0 (HEMCO), a stand-alone software component for computing emissions in global atmospheric models. HEMCO determines emissions from different sources, regions, and species on a user-defined grid and can combine, overlay, and update a set of data inventories and scale factors, as specified by the user through the HEMCO configuration file. New emission inventories at any spatial and temporal resolution are readily added to HEMCO and can be accessed by the user without any preprocessing of the data files or modification of the source code. Emissions that depend on dynamic source types and local environmental variables such as wind speed or surface temperature are calculated in separate HEMCO extensions.

HEMCO is fully compliant with the Earth System Modeling Framework (ESMF) environment. It is highly portable and can be deployed in a new model environment with only few adjustments at the top-level interface. So far, we have implemented HEMCO in the NASA Goddard Earth Observing System (GEOS-5) Earth system model (ESM) and in the GEOS-Chem chemical transport model (CTM).

By providing a widely applicable framework for specifying constituent emissions, HEMCO is designed to ease sensitivity studies and model comparisons, as well as inverse modeling in which emissions are adjusted iteratively. The HEMCO code, extensions, and the full set of emissions data files used in GEOS-Chem are available at http://wiki.geos-chem.org/HEMCO.

The GEOS-Chem global chemical transport model (CTM), used by a large atmospheric chemistry research community, has been re-engineered to also serve as an atmospheric chemistry module for Earth System Models (ESMs). This was done using an Earth System Modelling Framework (ESMF) interface that operates independently of the GEOS-Chem scientific code, permitting the exact same GEOS-Chem code to be used as an ESM module or as a stand-alone CTM. In this manner, the continual stream of updates contributed by the CTM user community is automatically passed on to the ESM module, which remains state-of-science and referenced to the latest version of the standard GEOS-Chem CTM. A major step in this re-engineering was to make GEOS-Chem grid-independent, i.e., capable of using any geophysical grid specified at run time. GEOS-Chem data "sockets" were also created for communication between modules and with external ESM code via the ESMF. The grid-independent, ESMF-compatible GEOS-Chem is now the standard version of the GEOS-Chem CTM. It has been implemented as an atmospheric chemistry module into the NASA GEOS-5 ESM. The coupled GEOS-5/GEOS-Chem system was tested for scalability and performance with a tropospheric oxidant-aerosol simulation (120 coupled species, 66 transported tracers) using 48–240 cores and MPI parallelization. Numerical experiments demonstrate that the GEOS-Chem chemistry module scales efficiently for the number of processors tested. Although inclusion of atmospheric chemistry in ESMs is computationally expensive, the excellent scalability of the chemistry module means that the relative cost goes down with increasing number of MPI processes.

Rivers discharge 28 ± 13 Mmol yr1 of mercury (Hg) to ocean margins, an amount comparable to atmospheric deposition to the global oceans. Most of the Hg discharged by rivers is sequestered by burial of benthic sediment in estuaries or the coastal zone, but some is evaded to the atmosphere and some is exported to the open ocean. We investigate the fate of riverine Hg by developing a new global 3-D simulation for Hg in the Massachusetts Institute of Technology ocean general circulation model. The model includes plankton dynamics and carbon respiration (DARWIN project model) coupled to inorganic Hg chemistry. Results are consistent with observed spatial patterns and magnitudes of surface ocean Hg concentrations. We use observational constraints on seawater Hg concentrations and evasion to infer that most Hg from rivers is sorbed to refractory organic carbon and preferentially buried. Only 6% of Hg discharged by rivers (1.8 Mmol yr1 ) is transported to the open ocean on a global basis. This fraction varies from a low of 2.6% in East Asia due to the barrier imposed by the Korean Peninsula and Japanese archipelago, up to 25% in eastern North America facilitated by the Gulf Stream. In the Arctic Ocean, low tributary particle loads and efficient degradation of particulate organic carbon by deltaic microbial communities favor a more labile riverine Hg pool. Evasion of Hg to the Arctic atmosphere is indirectly enhanced by heat transport during spring freshet that accelerates sea ice melt and ice rafting. Discharges of 0.23 Mmol Hg yr1 from Arctic rivers can explain the observed summer maximum in the Arctic atmosphere, and this magnitude of releases is consistent with recent observations. Our work indicates that rivers are major contributors to Hg loads in the Arctic Ocean.