Simulated Greenland Surface Mass Balance in the GISS ModelE2 GCM: Role of the Ice Sheet Surface

The rate of growth or retreat of the Greenland and Antarctic ice sheets remains a highly uncertain component of future sea level change. Here we examine the simulation of Greenland ice sheet surface mass balance (GrIS SMB) in the NASA Goddard Institute for Space Studies (GISS) ModelE2 General Circulation Model (GCM). GCMs are often limited in their ability to represent SMB compared with polarregion Regional Climate Models (RCMs). We compare ModelE2 simulated GrIS SMB for presentday (19962005) simulations with fixed ocean conditions, at a spatial resolution of 2 latitude by 2.5 longitude (~200 km), with SMB simulated by the Modle Atmosphrique Rgionale (MAR) RCM (19962005 at a 25 km resolution). ModelE2 SMB agrees well with MAR SMB on the whole, but there are distinct spatial patterns of differences and large differences in some SMB components. The impact of changes to the ModelE2 surface are tested, including a subgridscale representation of SMB with surface elevation classes. This has a minimal effect on ice sheetwide SMB, but corrects local biases. Replacing fixed surface albedo with satellitederived values and an agedependent scheme has a larger impact, increasing simulated melt by 60100%. We also find that lower surface albedo can enhance the effects of elevation classes. Reducing ModelE2 surface roughness length to values closer to MAR reduces sublimation by ~50%. Further work is required to account for meltwater refreezing in ModelE2, and to understand how differences in atmospheric processes and model resolution influence simulated SMB

Potassium–carbonate co-substituted hydroxyapatite compositions : maximising the level of carbonate uptake for potential CO2 utilisation options

The authors would like to acknowledge the University of Aberdeen and the Royal Commission for the Exhibition of 1851 for providing financial support.CO2 utilisation is a rapidly growing area of interest aimed at reducing the magnitude of anthropogenic greenhouse gas emissions. We report the synthesis of potassium–carbonate (K–CO3) co-substituted hydroxyapatites with potassium and carbonate contents ranging from approximately 0.4–0.9 wt% and 3.4–13.0 wt% respectively via an aqueous precipitation reaction between calcium hydroxide, phosphoric acid and either potassium carbonate or potassium hydrogen–carbonate. The incorporated carbonate is situated on both hydroxyl and phosphate sites. A subsequent heat treatment in dry CO2 at 600 °C allowed for a K–CO3 co-substituted apatite containing approximately 16.9 wt% CO32− to be prepared, amongst the largest carbonate contents that have been reported for such a material to date. Although this work shows that K–CO3 co-substituted apatites with high levels of carbonate incorporation can be prepared using simple, room temperature, aqueous precipitation reactions with starting reagents unlikely to pose significant environmental risks, testing of these materials in prospective applications (such as solid fertilisers) is required before they can be considered a viable CO2 utilisation option. A preliminary assessment of the effect of potassium/carbonate substitution on the solubility of the as-prepared compositions showed that increasing carbonate substitution increased the solubility.Publisher PDFPeer reviewe

A System of Conservative Regridding for Ice-Atmosphere Coupling in a General Circulation Model (GCM)

The method of elevation classes, in which the ice surface model is run at multiple elevations within each grid cell, has proven to be a useful way for a low-resolution atmosphere inside a general circulation model (GCM) to produce high-resolution downscaled surface mass balance fields for use in one-way studies coupling atmospheres and ice flow models. Past uses of elevation classes have failed to conserve mass and energy because the transformation used to regrid to the atmosphere was inconsistent with the transformation used to downscale to the ice model. This would cause problems for two-way coupling. A strategy that resolves this conservation issue has been designed and is presented here. The approach identifies three grids between which data must be regridded and five transformations between those grids required by a typical coupled atmosphere-ice flow model. This paper develops a theoretical framework for the problem and shows how each of these transformations may be achieved in a consistent, conservative manner. These transformations are implemented in Glint2, a library used to couple atmosphere models with ice models. Source code and documentation are available for download. Confounding real-world issues are discussed, including the use of projections for ice modeling, how to handle dynamically changing ice geometry, and modifications required for finite element ice models

Future Antarctic bed topography and its implications for ice sheet dynamics

The Antarctic bedrock is evolving as the solid Earth responds to the past and ongoing evolution of the ice sheet. A recently improved ice loading history suggests that the Antarctic Ice Sheet (AIS) has generally been losing its mass since the Last Glacial Maximum. In a sustained warming climate, the AIS is predicted to retreat at a greater pace, primarily via melting beneath the ice shelves. We employ the glacial isostatic adjustment (GIA) capability of the Ice Sheet System Model (ISSM) to combine these past and future ice loadings and provide the new solid Earth computations for the AIS. We find that past loading is relatively less important than future loading for the evolution of the future bed topography. Our computations predict that the West Antarctic Ice Sheet (WAIS) may uplift by a few meters and a few tens of meters at years AD 2100 and 2500, respectively, and that the East Antarctic Ice Sheet is likely to remain unchanged or subside minimally except around the Amery Ice Shelf. The Amundsen Sea Sector in particular is predicted to rise at the greatest rate; one hundred years of ice evolution in this region, for example, predicts that the coastline of Pine Island Bay will approach roughly 45 mm yr−1 in viscoelastic vertical motion. Of particular importance, we systematically demonstrate that the effect of a pervasive and large GIA uplift in the WAIS is generally associated with the flattening of reverse bed slope, reduction of local sea depth, and thus the extension of grounding line (GL) towards the continental shelf. Using the 3-D higher-order ice flow capability of ISSM, such a migration of GL is shown to inhibit the ice flow. This negative feedback between the ice sheet and the solid Earth may promote stability in marine portions of the ice sheet in the future

Fluctuations of an Atomic Ledge Bordering a Crystalline Facet

When a high symmetry facet joins the rounded part of a crystal, the step line density vanishes as sqrt(r) with r denoting the distance from the facet edge. This means that the ledge bordering the facet has a lot of space to meander as caused by thermal activation. We investigate the statistical properties of the border ledge fluctuations. In the scaling regime they turn out to be non-Gaussian and related to the edge statistics of GUE multi-matrix models.Comment: Version with major revisions -- RevTeX, 4 pages, 2 figure

Pattern Formation of Glioma Cells: Effects of Adhesion

We investigate clustering of malignant glioma cells. \emph{In vitro} experiments in collagen gels identified a cell line that formed clusters in a region of low cell density, whereas a very similar cell line (which lacks an important mutation) did not cluster significantly. We hypothesize that the mutation affects the strength of cell-cell adhesion. We investigate this effect in a new experiment, which follows the clustering dynamics of glioma cells on a surface. We interpret our results in terms of a stochastic model and identify two mechanisms of clustering. First, there is a critical value of the strength of adhesion; above the threshold, large clusters grow from a homogeneous suspension of cells; below it, the system remains homogeneous, similarly to the ordinary phase separation. Second, when cells form a cluster, we have evidence that they increase their proliferation rate. We have successfully reproduced the experimental findings and found that both mechanisms are crucial for cluster formation and growth.Comment: 6 pages, 6 figure

Twisted Classical Poincar\'{e} Algebras

We consider the twisting of Hopf structure for classical enveloping algebra $U(\hat{g})$, where $\hat{g}$ is the inhomogenous rotations algebra, with explicite formulae given for $D=4$ Poincar\'{e} algebra $(\hat{g}={\cal P}_4).$ The comultiplications of twisted $U^F({\cal P}_4)$ are obtained by conjugating primitive classical coproducts by $F\in U(\hat{c})\otimes U(\hat{c}),$ where $\hat{c}$ denotes any Abelian subalgebra of ${\cal P}_4$, and the universal $R-$matrices for $U^F({\cal P}_4)$ are triangular. As an example we show that the quantum deformation of Poincar\'{e} algebra recently proposed by Chaichian and Demiczev is a twisted classical Poincar\'{e} algebra. The interpretation of twisted Poincar\'{e} algebra as describing relativistic symmetries with clustered 2-particle states is proposed.Comment: \Large \bf 19 pages, Bonn University preprint, November 199

Development of a telescope for medium-energy gamma-ray astronomy

The Advanced Energetic Pair Telescope (AdEPT) is being developed at GSFC as a future NASA MIDEX mission to explore the medium-energy (5–200 MeV) gamma-ray range. The enabling technology for AdEPT is the Three- Dimensional Track Imager (3-DTI), a gaseous time projection chamber. The high spatial resolution 3-D electron tracking of 3-DTI enables AdEPT to achieve high angular resolution gamma-ray imaging via pair production and triplet production (pair production on electrons) in the medium-energy range. The low density and high spatial resolution of 3-DTI allows the electron positron track directions to be measured before they are dominated by Coulomb scattering. Further, the significant reduction of Coulomb scattering allows AdEPT to be the first medium-energy gamma-ray telescope to have high gamma-ray polarization sensitivity. We review the science goals that can be addressed with a medium-energy pair telescope, how these goals drive the telescope design, and the realization of this design with AdEPT. The AdEPT telescope for a future MIDEX mission is envisioned as a 8 m3 active volume filled with argon at 2 atm. The design and performance of the 3-DTI detectors for the AdEPT telescope are described as well as the outstanding instrument challenges that need to be met for the AdEPT mission

Antarctic firn compaction rates from repeat-track airborne radar data: I. Methods

While measurements of ice-sheet surface elevation change are increasingly used to assess mass change, the processes that control the elevation fluctuations not related to ice-flow dynamics (e.g. firn compaction and accumulation) remain difficult to measure. Here we use radar data from the Thwaites Glacier (West Antarctica) catchment to measure the rate of thickness change between horizons of constant age over different time intervals: 2009–10, 2010–11 and 2009–11. The average compaction rate to ∼25 m depth is 0.33 m a–1, with largest compaction rates near the surface. Our measurements indicate that the accumulation rate controls much of the spatio-temporal variations in the compaction rate while the role of temperature is unclear due to a lack of measurements. Based on a semi-empirical, steady-state densification model, we find that surveying older firn horizons minimizes the potential bias resulting from the variable depth of the constant age horizon. Our results suggest that the spatiotemporal variations in the firn compaction rate are an important consideration when converting surface elevation change to ice mass change. Compaction rates varied by up to 0.12 m a–1 over distances 20% larger during the 2010–11 interval than during 2009–10