Search

We develop regionally optimized Gravity Recovery and Climate Experiment (GRACE) solutions to evaluate the mass balance of the drainage basins of Amery Ice Shelf, East Antarctica, and Getz Ice Shelf, West Antarctica. We find that the Amery region is near balance, while the Getz region is rapidly losing mass. We compare the results with the mass budget method (MBM) combining ice discharge along the periphery with surface mass balance derived from three regional climate models: (1) Regional Atmospheric Climate Model (RACMO) 2.3p1 and (2) 2.3p2, and (3) Modèle Atmosphérique Régional 3.6.41. For Amery, MBM/RACMO2.3p1 agrees with GRACE, while MBM/RACMO2.3p2 and MBM/MAR3.6.41 suggest a positive mass balance. For Getz, all estimates agree with a mass loss and the GRACE results are robust to uncertainties in glacial isostatic adjustment derived from an ensemble 128,000 forward models. Over the period April 2002 to November 2015, the mass loss of the Getz drainage basin is 22.9 ± 10.9 Gt/yr with an acceleration of 1.6 ± 0.9 Gt/yr2.

The continuous and precise mapping of glacier calving fronts is essential for monitoring and understanding rapid glacier changes in Antarctica and Greenland, which have the potential for significant sea level rise within the current century. This effort has been mostly restricted to the slow and painstaking manual digitalization of the calving front positions in thousands of satellite imagery products. Here, we have developed a machine learning toolkit to automatically detect glacier calving front margins in satellite imagery. The toolkit is based on semantic image segmentation using Convolutional Neural Networks (CNN) with a modified U-Net architecture to isolate the calving fronts from satellite images after having been trained with a dataset of images and their corresponding manually-determined calving fronts. As a case study we train our neural network on a varied set of Landsat images with lowered resolutions from Jakobshavn, Sverdrup, and Kangerlussuaq glaciers, Greenland and test the results on images from Helheim glacier, Greenland to evaluate the performance of the approach. The neural network is able to identify the calving front in new images with a mean deviation of 96.3 m from the true fronts, equivalent to 1.97 pixels on average, while the corresponding error for manually-determined fronts on the same resolution images is 92.5 m (1.89 pixels). We find that the trained neural network significantly outperforms common edge detection techniques, and can be used to continuously map out calving-ice fronts with a variety of data products.

Totten and Moscow University glaciers, in the marine-based sector of East Antarctica, contain enough ice to raise sea level by 5 m. Obtaining precise measurements of their mass balance is challenging owing to large area of the basins and the small mass balance signal compared to West Antarctic glaciers. Here we employ a locally optimized processing of Gravity Recovery and Climate Experiment (GRACE) harmonics to evaluate their mass balance at the sub-basin scale and compare the results with mass budget method (MBM) estimates using regional atmospheric climate model version 2.3 (RACMO2.3) or Modèle Atmosphérique Régional version 3.6.4 (MAR3.6.4). The sub-basin mass loss estimate for April 2002 to November 2015 is 14.8 ± 4.3 Gt/yr, which is weakly affected by glacial isostatic adjustment uncertainties (±1.4 Gt/yr). This result agrees with MBM/RACMO2.3 (15.8 ± 2.0 Gt/yr), whereas MBM/MAR3.6.4 underestimates the loss (6.6 ± 1.6 Gt/yr). For the entire drainage, the mass loss for April 2002 to August 2016 is 18.5 ± 6.6 Gt/yr, or 15 ± 4% of its ice flux. These results provide unequivocal evidence for mass loss in this East Antarctic sector.

Marine-terminating glaciers play a critical role in controlling Greenland's ice sheet mass balance. Their frontal margins interact vigorously with the ocean, but our understanding of this interaction is limited, in part, by a lack of bathymetry data. Here we present a multibeam echo sounding survey of 14 glacial fjords in the Uummannaq and Vaigat fjords, west Greenland, which extends from the continental shelf to the glacier fronts. The data reveal valleys with shallow sills, overdeepenings (>1300 m) from glacial erosion, and seafloor depths 100-1000 m deeper than in existing charts. Where fjords are deep enough, we detect the pervasive presence of warm, salty Atlantic Water (AW) (>2.5°C) with high melt potential, but we also find numerous glaciers grounded on shallow (

We use numerical climate simulations, paleoclimate data, and modern observations to study the effect of growing ice melt from Antarctica and Greenland. Meltwater tends to stabilize the ocean column, inducing amplifying feedbacks that increase subsurface ocean warming and ice shelf melting. Cold meltwater and induced dynamical effects cause ocean surface cooling in the Southern Ocean and North Atlantic, thus increasing Earth's energy imbalance and heat flux into most of the global ocean's surface. Southern Ocean surface cooling, while lower latitudes are warming, increases precipitation on the Southern Ocean, increasing ocean stratification, slowing deepwater formation, and increasing ice sheet mass loss. These feedbacks make ice sheets in contact with the ocean vulnerable to accelerating disintegration. We hypothesize that ice mass loss from the most vulnerable ice, sufficient to raise sea level several meters, is better approximated as exponential than by a more linear response. Doubling times of 10, 20 or 40 years yield multi-meter sea level rise in about 50, 100 or 200 years. Recent ice melt doubling times are near the lower end of the 10-40-year range, but the record is too short to confirm the nature of the response. The feedbacks, including subsurface ocean warming, help explain paleoclimate data and point to a dominant Southern Ocean role in controlling atmospheric CO2, which in turn exercised tight control on global temperature and sea level. The millennial (500-2000-year) timescale of deep-ocean ventilation affects the timescale for natural CO2 change and thus the timescale for paleo-global climate, ice sheet, and sea level changes, but this paleo-millennial timescale should not be misinterpreted as the timescale for ice sheet response to a rapid, large, human-made climate forcing. These climate feedbacks aid interpretation of events late in the prior interglacial, when sea level rose to +6-9 m with evidence of extreme storms while Earth was less than 1 °C warmer than today. Ice melt cooling of the North Atlantic and Southern oceans increases atmospheric temperature gradients, eddy kinetic energy and baroclinicity, thus driving more powerful storms. The modeling, paleoclimate evidence, and ongoing observations together imply that 2 °C global warming above the preindustrial level could be dangerous. Continued high fossil fuel emissions this century are predicted to yield (1) cooling of the Southern Ocean, especially in the Western Hemisphere; (2) slowing of the Southern Ocean overturning circulation, warming of the ice shelves, and growing ice sheet mass loss; (3) slowdown and eventual shutdown of the Atlantic overturning circulation with cooling of the North Atlantic region; (4) increasingly powerful storms; and (5) nonlinearly growing sea level rise, reaching several meters over a timescale of 50-150 years. These predictions, especially the cooling in the Southern Ocean and North Atlantic with markedly reduced warming or even cooling in Europe, differ fundamentally from existing climate change assessments. We discuss observations and modeling studies needed to refute or clarify these assertions.

The relationship between leading atmospheric teleconnection patterns and Greenland Ice Sheet (GrIS) temperature, precipitation, and surface mass balance (SMB) are investigated for the last 36 summers (1979-2014) based on Modern-Era Retrospective analysis for Research and Applications version 2 reanalyses. The results indicate that the negative phase of both the North Atlantic Oscillation (NAO) and Arctic Oscillation, associated with warm and dry conditions for the GrIS, lead to SMB decreases within 0-1 months. Furthermore, the positive phase of the East Atlantic (EA) pattern often lags the negative NAO, reflecting a dynamical linkage between these modes that acts to further enhance the warm and dry conditions over the GrIS, leading to a favorable environment for enhanced surface mass loss. The development of a strong negative NAO in combination with a strong positive EA in recent years leads to significantly larger GrIS warming compared to when the negative NAO occurs in combination with a negative or weak positive EA (0.69 K versus 0.13 K anomaly). During 2009 and 2011, weakened (as compared to conditions during the severe surface melt cases of 2010 and 2012) local high pressure blocking produced colder northerly flow over the GrIS inhibiting warming despite the occurrence of a strong negative NAO, reflecting an important role for the EA during those years. In particular, the EA acts with the NAO to enhance warming in 2010 and 2012, and weaken high pressure blocking in 2009 and 2011. In general, high pressure blocking primarily impacts the western areas of the GrIS via advective temperature increases, while changes in net surface radiative fluxes account for both western and eastern GrIS temperature changes.

Changing rainfall patterns have significant effect on water resources, agriculture output in many countries, especially the country like India where the economy depends on rain-fed agriculture. Rainfall over India has large spatial as well as temporal variability. To understand the variability in rainfall, spatial–temporal analyses of rainfall have been studied by using 107 (1901–2007) years of daily gridded India Meteorological Department (IMD) rainfall datasets. Further, the validation of IMD precipitation data is carried out with different observational and different reanalysis datasets during the period from 1989 to 2007. The Global Precipitation Climatology Project data shows similar features as that of IMD with high degree of comparison, whereas Asian Precipitation-Highly-Resolved Observational Data Integration Towards Evaluation data show similar features but with large differences, especially over northwest, west coast and western Himalayas. Spatially, large deviation is observed in the interior peninsula during the monsoon season with National Aeronautics Space Administration-Modern Era Retrospective-analysis for Research and Applications (NASA-MERRA), pre-monsoon with Japanese 25 years Re Analysis (JRA-25), and post-monsoon with climate forecast system reanalysis (CFSR) reanalysis datasets. Among the reanalysis datasets, European Centre for Medium-Range Weather Forecasts Interim Re-Analysis (ERA-Interim) shows good comparison followed by CFSR, NASA-MERRA, and JRA-25. Further, for the first time, with high resolution and long-term IMD data, the spatial distribution of trends is estimated using robust regression analysis technique on the annual and seasonal rainfall data with respect to different regions of India. Significant positive and negative trends are noticed in the whole time series of data during the monsoon season. The northeast and west coast of the Indian region shows significant positive trends and negative trends over western Himalayas and north central Indian region.

We use GRACE-derived terrestrial water storage (TWS) and ERA-interim air temperature, as proxy for available water and temperature constraints on vegetation productivity, inferred from MODIS satellite normalized difference vegetation index (NDVI), in Northern Eurasia during 2002-2011. We investigate how changes in TWS affect the correlation between NDVI and temperature during the non-frozen season. We find that vegetation growth exhibits significant spatial and temporal variability associated with varying trend in TWS and temperature. The largest NDVI gains occur over boreal forests associated with warming and wetting. The largest NDVI losses occur over grasslands in the Southwestern Ob associated with regional drying and cooling, with dominant constraint from TWS. Over grasslands and temperate forests in the Southeast Ob and South Yenisei, wetting and cooling lead to a dominant temperature constraint due to the relaxation of TWS constraints. Overall, we find significant monthly correlation of NDVI with TWS and temperature over 35% and 50% of the domain, respectively. These results indicate that water availability (TWS) plays a major role in modulating Eurasia vegetation response to temperature changes.

Using GPS radio occultation refractivity data collected over the period of 2002-2013, we present a new method for identification of the intertropical convergence zone (ITCZ). The ITCZ is identified by estimating the maximum in the monthly meridional refractivity and specific humidity field by applying a Gaussian fit at each longitude. The interannual variability and climatology of the ITCZ is presented from 12. years of refractivity data. This new method captures all the general features of ITCZ extent and its variability. We also examine the effects of the ITCZ vertically during different seasons. The ITCZ is observed mostly at eastern Pacific in May month, and it is zonally distributed in the September and October months of each year. The zonal variability is large between lower and higher levels, particularly over the Indian monsoon and western Pacific. The latitudinal difference in the vertical extent between 850. hPa and higher levels is larger during the northern hemisphere (NH) summer than NH winter.

The characteristics of the quasi-2-day wave (QTDW) in the upper stratosphere and lower mesospheric (USLM) altitudes over the northern hemisphere (NH) and southern hemisphere (SH) have been studied by using Global Positioning Radio Occultation (GPSRO) Constellation Observing System for Meteorology, Ionosphere, and Climate (COSMIC) temperature data from November 2006 to December 2010. We studied the seasonal, latitudinal, and interannual variability of the westward-propagating 2-day wave coincident with zonal wave number 3 in both hemispheres in the altitude range of 20-60 km. The Lomb-Scargle periodogram (LSP) analysis indicates the dominance of the QTDW in the USLM in both the hemispheres. The observed amplitude of the wave is maximum during the winter season in middle and higher latitudes, with monthly mean amplitudes being as high as ∼8 K. These amplitudes are found frequently during the late fall and peak to a maximum in the NH winter season. In the SH, QTDW amplitudes are found in the early winter season and appear till the early fall months. The QTDW varies from 49 ± 3 to 48 ± 1 h in the NH and SH, respectively, with westward-propagating wave number 3. The amplitudes of the wave are large during winter in both the hemispheres, and, comparatively, the NH amplitudes are larger than those of the SH in higher latitudes.

We compare four independent estimates of the mass balance of the Amundsen Sea Embayment of West Antarctica, an area experiencing rapid retreat and mass loss to the sea. We use ICESat and Operation IceBridge laser altimetry, Envisat radar altimetry, GRACE time-variable gravity, RACMO2.3 surface mass balance, ice velocity from imaging radars, and ice thickness from radar sounders. The four methods agree in terms of mass loss and acceleration in loss at the regional scale. Over 1992-2013, the mass loss is 83 ± 5 Gt/yr with an acceleration of 6.1 ± 0.7 Gt/yr2. During the common period 2003-2009, the mass loss is 84 ± 10 Gt/yr with an acceleration of 16.3 ± 5.6 Gt/yr2, nearly 3 times the acceleration over 1992-2013. Over 2003-2011, the mass loss is 102 ± 10 Gt/yr with an acceleration of 15.7 ± 4.0 Gt/yr2. The results reconcile independent mass balance estimates in a setting dominated by change in ice dynamics with significant variability in surface mass balance.

We use Gravity Recovery and Climate Experiment (GRACE) monthly gravity fields to determine the regional acceleration in ice mass loss in Greenland and Antarctica for 2003-2013. We find that the total mass loss is controlled by only a few regions. In Greenland, the southeast and northwest generate 70% of the loss (280±58 Gt/yr) mostly from ice dynamics, the southwest accounts for 54% of the total acceleration in loss (25.4±1.2 Gt/yr2) from a decrease in surface mass balance (SMB), followed by the northwest (34%), and we find no significant acceleration in the northeast. In Antarctica, the Amundsen Sea (AS) sector and the Antarctic Peninsula account for 64% and 17%, respectively, of the total loss (180±10 Gt/yr) mainly from ice dynamics. The AS sector contributes most of the acceleration in loss (11±4 Gt/yr2), and Queen Maud Land, East Antarctica, is the only sector with a significant mass gain due to a local increase in SMB (63±5 Gt/yr).

© 2015 Elsevier Inc. Using GPS radio occultation refractivity data collected over the period of 2002-2013, we present a new method for identification of the intertropical convergence zone (ITCZ). The ITCZ is identified by estimating the maximum in the monthly meridional refractivity and specific humidity field by applying a Gaussian fit at each longitude. The interannual variability and climatology of the ITCZ is presented from 12. years of refractivity data. This new method captures all the general features of ITCZ extent and its variability. We also examine the effects of the ITCZ vertically during different seasons. The ITCZ is observed mostly at eastern Pacific in May month, and it is zonally distributed in the September and October months of each year. The zonal variability is large between lower and higher levels, particularly over the Indian monsoon and western Pacific. The latitudinal difference in the vertical extent between 850. hPa and higher levels is larger during the northern hemisphere (NH) summer than NH winter.

This study investigates how warming and changes in precipitation may affect the cycling of carbon (C) in tundra soils, and between high Arctic tundra and the atmosphere. We quantified ecosystem respiration (Reco) and soil pore space CO2 in a polar semi-desert in northwestern Greenland under current and future climate conditions simulated by long-term experimental warming (+2 °C, +4 °C), water addition (+50% summer precipitation), and a combination of both (+4 °C × +50% summer precipitation). We also measured the 14C content of Reco and soil CO2 to distinguish young C cycling rapidly between the atmosphere and the ecosystem from older C stored in the soil for centuries to millennia. We identified changes in the amount and timing of precipitation as a key control of the magnitude, seasonality and sources of Reco in a polar semi-desert. Throughout each summer, small (4 mm), more winter snow and experimental irrigation were associated with higher Reco fluxes and the release of recently fixed (young) C. Warmer summers and experimental warming also resulted in higher Reco fluxes (+2 °C > +4 °C), but coincided with losses of older C. We conclude that in high Arctic, dry tundra systems, future magnitudes and patterns of old C emissions will be controlled as much by the summer precipitation regime and winter snowpack as by warming. The release of older soil C is of concern, as it may lead to net C losses from the ecosystem. Therefore, reliable predictions of precipitation amounts, frequency, and timing are required to predict the changing C cycle in the high Arctic. © Author(s) 2014.

The non-frozen (NF) season duration strongly influences the northern carbon cycle where frozen (FR) temperatures are a major constraint to biological processes. The landscape freeze-thaw (FT) signal from satellite microwave remote sensing provides a surrogate measure of FR temperature constraints to ecosystem productivity, trace gas exchange, and surface water mobility. We analysed a new global satellite data record of daily landscape FT dynamics derived from temporal classification of overlapping SMMR and SSM/I 37 GHz frequency brightness temperatures (Tb). The FT record was used to quantify regional patterns, annual variability, and trends in the NF season over northern (≥45°N) vegetated land areas. The ecological significance of these changes was evaluated against satellite normalized difference vegetation index (NDVI) anomalies, estimated moisture and temperature constraints to productivity determined from meteorological reanalysis, and atmospheric CO2 records. The FT record shows a lengthening (2.4 days decade-1; p < 0.005) mean annual NF season trend (1979-2010) for the high northern latitudes that is 26% larger than the Northern Hemisphere trend. The NDVI summer growth response to these changes is spatially complex and coincides with local dominance of cold temperature or moisture constraints to productivity. Longer NF seasons are predominantly enhancing productivity in cold temperature-constrained areas, whereas these effects are reduced or reversed in more moisture-constrained areas. Longer NF seasons also increase the atmospheric CO2 seasonal amplitude by enhancing both regional carbon uptake and emissions. We find that cold temperature constraints to northern growing seasons are relaxing, whereas potential benefits for productivity and carbon sink activity are becoming more dependent on the terrestrial water balance and supply of plant-available moisture needed to meet additional water use demands under a warming climate. © 2014 Taylor & Francis.

Long-term data available from Lidar systems located at three different locations namely São José dos Campos, Brazil (23.2° S, 45.8° W), Gadanki (13.5° N, 79.2° E) and Reunion (20.8° S, 55.5° E) have been used to investigate the long-term variations like Annual, Semi-annual, Quasi-biennial, El Nino Southern Oscillation and solar cycle. These oscillations are also extracted from simultaneous satellite borne measurements of HALogen Occultation Experiment (HALOE) instrument onboard UARS and SABER onboard TIMED over these stations making largest time series covering the entire middle atmosphere. A good agreement is found between the LIDAR and satellite-derived amplitudes and phases between 30 and 65 km altitude, which suggests that satellite measurements can be used to investigate the long-term trends globally. Latter measurements are extended to 80 km in order to further investigate these oscillations. Large difference in the amplitudes between the eastern pacific and western pacific is noticed in these oscillations. Changing from cooling trends in the stratosphere to warming trends in the mesosphere occurs more or less at altitude around 70 km altitude and this result agrees well with that observed by satellite measurements reported in the literature. The peak in the cooling trend does not occur at a fixed altitude in the stratosphere however maximum warming trend is observed around 75 km at all the stations. The observed long-term trends including various oscillations are compared with that reported with various techniques.

Glacial isostatic adjustment (GIA) represents a source of uncertainty for ice sheet mass balance estimates from the Gravity Recovery and Climate Experiment (GRACE) time-variable gravity measurements. We evaluate Greenland GIA corrections from Simpson et al (2009 Quat. Sci. Rev. 28 1631–57), A et al (2013 Geophys. J. Int. 192 557–72) and Wu et al (2010 Nature Geosci. 3 642–6) by comparing the spatial patterns of GRACE-derived ice mass trends calculated using the three corrections with volume changes from ICESat (Ice, Cloud, and land Elevation Satellite) and OIB (Operation IceBridge) altimetry missions, and surface mass balance products from the Regional Atmospheric Climate Model (RACMO). During the period September 2003–August 2011, GRACE ice mass changes obtained using the Simpson et al (2009 Quat. Sci. Rev. 28 1631–57) and Aet al (2013 Geophys. J. Int. 192 557–72) GIA corrections yield similar spatial patterns and amplitudes, and are consistent with altimetry observations and surface mass balance data. The two GRACE estimates agree within 2% on average over the entire ice sheet, and better than 15% in four subdivisions of Greenland. The third GRACE estimate corrected using the (Wuet al 2010 Nature Geosci. 3 642–6)) GIA shows similar spatial patterns, but produces an average ice mass loss for the entire ice sheet that is 64 − 67 Gt yr−1 smaller. In the Northeast the recovered ice mass change is 46–49 Gt yr−1 (245–270%) more positive than that deduced from the other two corrections. By comparing the spatial and temporal variability of the GRACE estimates with trends of volume changes from altimetry and surface mass balance from RACMO, we show that the Wuet al (2010 Nature Geosci. 3 642–6) correction leads to a large mass increase in the Northeast that is inconsistent with independent observations.

In this paper, GPS RO stratospheric temperatures are compared with different satellite and different model data sets. The data sets used for comparison include Japanese 25-year Reanalysis (JRA-25), UK Met office (MetO), ERA-Interim, GEOS5, and satellite temperatures are AIRS_Aqua, HIRDLS, Aura_MLS and SABER. Apart from the comparison, we also studied the seasonal variation of temperature during summer and winter in both the hemispheres. The seasonal and latitudinal variation of temperature by GPS (CHAMP+COSMIC) with other reanalysis (JRA-25, GEOS5 and METO) and satellite measurements (AURA_MLS and SABER) for both the hemispheres show reasonably good agreement. The difference of about ±0.5 to ±0.75 K is observed at 20 km, ±1K to ±1.75K at 35 km while the high values of ~3 to 4 K are observed at upper stratosphere and lower mesosphere height regions (50 km to 60 km). The GPS RO temperatures above 35 km are noticed to be warmer than reanalysis and satellite data sets in SH region. This tendency increases with increasing height and reaches its maximum at 60 km, with magnitude of 3 K to 4 K with reanalysis data sets and 2 K to 3.5 K with satellite measurements. The calculated SAO and AO amplitudes based on GPS data are found to be comparable with the earlier results.

The evolution of ocean temperature measurement systems is presented with a focus on the development and accuracy of two critical devices in use today (expendable bathythermographs and conductivity-temperature-depth instruments used on Argo floats). A detailed discussion of the accuracy of these devices and a projection of the future of ocean temperature measurements are provided. The accuracy of ocean temperature measurements is discussed in detail in the context of ocean heat content, Earth's energy imbalance, and thermosteric sea level rise. Up-to-date estimates are provided for these three important quantities. The total energy imbalance at the top of atmosphere is best assessed by taking an inventory of changes in energy storage. The main storage is in the ocean, the latest values of which are presented. Furthermore, despite differences in measurement methods and analysis techniques, multiple studies show that there has been a multidecadal increase in the heat content of both the upper and deep ocean regions, which reflects the impact of anthropogenic warming. With respect to sea level rise, mutually reinforcing information from tide gauges and radar altimetry shows that presently, sea level is rising at approximately 3 mm yr-1 with contributions from both thermal expansion and mass accumulation from ice melt. The latest data for thermal expansion sea level rise are included here and analyzed. Key Points Oceanographic techniques and analysis have improved over many decadesThese improvements allow more accurate Earth-energy balance estimatesUnderstanding of ocean heat content and sea-level rise has also increased ©2013. American Geophysical Union. All Rights Reserved.

Time-variable gravity data from the Gravity Recovery and Climate Experiment (GRACE) mission have been available since 2002 to estimate the mass balance of the Greenland and Antarctic Ice Sheets. We analyze current progress and uncertainties in GRACE estimates of ice sheet mass balance. We discuss the impacts of errors associated with spherical harmonic truncation, spatial averaging, temporal sampling, and leakage from other time-dependent signals (e.g., glacial isostatic adjustment (GIA)). The largest sources of error for Antarctica are the GIA correction, the omission of l=1 terms, nontidal changes in ocean mass, and measurement errors. For Greenland, the errors come mostly from the uncertainty in the scaling factor. Using Release 5.0 (RL05) GRACE fields for January 2003 through November 2012, we find a mass change of −258 ± 41 Gt/yr for Greenland, with an acceleration of −31 ± 6 Gt/yr2, and a loss that migrated clockwise around the ice sheet margin to progressively affect the entire periphery. For Antarctica, we report changes of −83 ± 49 and −147 ± 80 Gt/yr for two GIA models, with an acceleration of −12 ± 9 Gt/yr2 and a dominance from the southeast pacific sector of West Antarctica and the Antarctic Peninsula.

The AURA-MLS daily mean temperatures and zonal wind from NASA-MERRA reanalysis for latitudes between 60° N and 80° N are used to investigate the planetary wave (PW) characteristics in the stratosphere and lower mesosphere during sudden stratospheric warming (SSW) (November 2008 to March 2009). Here, we used a novel method called empirical mode decomposition (EMD) to extract the PWs from the temperature data. The EMD is an interesting approach to decompose signals into locally periodic components, the intrinsic mode functions (IMFs), and will easily identify the embedded structures, even those with small amplitudes. The spectral analysis reveals prevailing planetary wave periods of ~6-day, ~8-day, ~15-day, and ~21–23-day in IMFs 1, 2, 3, and 4, respectively. Clear upward propagation of these waves (20–30 days) is observed, suggesting that sources for these oscillations are in the troposphere.

Within the last decade, the Greenland ice sheet (GrIS) and its surroundings have experienced record high surface temperatures (Mote, 2007; Box et al., 2010), ice sheet melt extent (Fettweis et al., 2011) and record-low summer sea-ice extent (Nghiem et al., 2007). Using three independent data sets, we derive, for the first time, consistent ice-mass trends and temporal variations within seven major drainage basins from gravity fields from the Gravity Recovery and Climate Experiment (GRACE; Tapley et al., 2004), surface-ice velocities from Inteferometric Synthetic Aperture Radar (InSAR; Rignot and Kanagaratnam, 2006) together with output of the regional atmospheric climate modelling (RACMO2/GR; Ettema et al., 2009), and surface-elevation changes from the Ice, cloud and land elevation satellite (ICESat; Sørensen et al., 2011). We show that changing ice discharge (D), surface melting and subsequent run-off (M/R) and precipitation (P) all contribute, in a complex and regionally variable interplay, to the increasingly negative mass balance of the GrIS observed within the last decade. Interannual variability in P along the northwest and west coasts of the GrIS largely explains the apparent regional mass loss increase during 2002-2010, and obscures increasing M/R and D since the 1990s. In winter 2002/2003 and 2008/2009, accumulation anomalies in the east and southeast temporarily outweighed the losses by M/R and D that prevailed during 2003-2008, and after summer 2010. Overall, for all basins of the GrIS, the decadal variability of anomalies in P, M/R and D between 1958 and 2010 (w.r.t. 1961-1990) was significantly exceeded by the regional trends observed during the GRACE period (2002-2011). © 2012 Elsevier B.V.

We use monthly measurements of time-variable gravity from the GRACE (Gravity Recovery and Climate Experiment) satellite mission to quantify changes in terrestrial water storage (TWS) in the Lena river basin, Eurasia, during the period April 2002 to September 2010. We estimate a TWS increase of 32 10 km 3/yr for the entire basin, equivalent to an increase in water thickness of 1.3 0.4cm/yr over a basin of 2.4 million km2. We compare TWS estimates from GRACE with time series of precipitation (P) minus evapotranspiration (ET) from ERA-Interim reanalysis minus observational river discharge (R). We find an excellent agreement in annual and inter-annual variability between the two time series. Furthermore, we find that a bias of -20 10% in P-ET is sufficient to effectively close the water budget with GRACE. When we account for this bias, the time series of cumulative TWS from GRACE and climatological data agree to within 3.8 cm of water thickness, or 9% of the mean annual P. The TWS increase is not uniform across the river basin and exhibits a peak, over an area of 502,400 km2, centered at 118.5E, 62.5N, and underlain by discontinuous permafrost. In this region, we attribute the observed TWS increase of 68 19 km3 to an increase in subsurface water storage. This large subsurface water signal will have a significant impact on the terrestrial hydrology of the region, including increased baseflow and alteration of seasonal runoff. © 2012 by the American Geophysical Union.

In this study, global (50°S-50°N) distribution of water vapor is investigated using COSMIC GPS RO measurements. Detailed comparisons have been made between COSMIC and high resolution GPS radiosonde measurements across 13 tropical stations and model outputs (ERA-Interim, NCEP, and JRA-25 reanalyses data sets). In comparison with independent techniques like radiosonde (Väisälä), it is found that COSMIC GPS RO wet profiles are accurate up to 7-8. km (assuming radiosonde as standard technique). In general, comparisons with corresponding seasonal means of model outputs are qualitatively in good agreement, although they differ quantitatively especially over convective regions of South America, Africa, and Indonesia. In tropical latitudes, the COSMIC specific humidity values are higher than the model outputs. Among various model outputs, ERA-Interim data set show near realistic features to that observed by COSMIC GPS RO measurements. Large asymmetry in the specific humidity distribution is observed between northern and southern hemispheres. © 2011 Elsevier Ltd.

Ice sheet mass balance estimates have improved substantially in recent years using a variety of techniques, over different time periods, and at various levels of spatial detail. Considerable disparity remains between these estimates due to the inherent uncertainties of each method, the lack of detailed comparison between independent estimates, and the effect of temporal modulations in ice sheet surface mass balance. Here, we present a consistent record of mass balance for the Greenland and Antarctic ice sheets over the past two decades, validated by the comparison of two independent techniques over the last 8 years: one differencing perimeter loss from net accumulation, and one using a dense time series of time-variable gravity. We find excellent agreement between the two techniques for absolute mass loss and acceleration of mass loss. In 2006, the Greenland and Antarctic ice sheets experienced a combined mass loss of 475 ± 158 Gt/yr, equivalent to 1.3 ± 0.4 mm/yr sea level rise. Notably, the acceleration in ice sheet loss over the last 18 years was 21.9 ± 1 Gt/yr2 for Greenland and 14.5 ± 2 Gt/yr2 for Antarctica, for a combined total of 36.3 ± 2 Gt/yr2. This acceleration is 3 times larger than for mountain glaciers and ice caps (12 ± 6 Gt/yr2). If this trend continues, ice sheets will be the dominant contributor to sea level rise in the 21st century.

Since 2002, global snow water equivalent (SWE) estimates have been generated using Advanced Microwave Scanning Radiometer (AMSR-E)/Aqua data. Accurate estimates of SWE are important to improve monitoring and managing of water resources in specific regions. SWE and snow map product accuracy are functions of topography and of land cover type because landscape characteristics have a strong influence on redistribution and physical properties of snow cover, and influence the microwave properties of the surface. Here we evaluate the AMSR-E SWE and derived snow map products in the Mackenzie River Basin (MRB), Canada, which is characterized by complex topography and varying land cover types from tundra to boreal forest. We compare in situ snow depth observations and Moderate Resolution Imaging Spectroradiometer (MODIS) snow cover maps from January 2003 to December 2007 with passive microwave remotely sensed SWE from AMSR-E and derived snow cover maps. In the MRB the mean absolute error ranges from 12 mm in the early winter season to 50 mm in the late winter season and overestimations of snow cover maps based on a 1 mm threshold of AMSR-E SWE varies from 4% to 8%. The optimal threshold for AMSR-E SWE to classify the pixels as snow ranges from 6 mm to 9 mm. The overall accuracy of new snow cover maps from AMSR-E varies from 91% to 94% in different sub-basins in the MRB. © 2010 by the authors.

The annual exchange of water between the continents and oceans is observed by GPS, gravimetry, and altimetry. However, the global average amplitude of this annual cycle (observed amplitude of ∼8 mm) is not representative of the effects that would be observed at individual tide gauges or at ocean bottom pressure recorders because of self-attraction and loading effects (SAL). In this paper, we examine the spatial variation of sea level change caused by the three main components that load the Earth and contribute to the water cycle: hydrology (including snow), the atmosphere, and the dynamic ocean. The SAL effects cause annual amplitudes at tide gauges (modeled here with a global average of ∼9 mm) to vary from less than 2 mm to more than 18 mm. We find a variance reduction (global average of 3 to 4%) after removing the modeled time series from a global set of tide gauges. We conclude that SAL effects are significant in places (e.g., the south central Pacific and coastal regions in Southeast Asia and west central Africa) and should be considered when interpreting these data sets and using them to constrain ocean circulation models.

Widespread glacier acceleration has been observed in Greenland in the past few years associated with the thinning of the lower reaches of the glaciers as they terminate in the ocean. These glaciers thin both at the surface, from warm air temperatures, and along their submerged faces in contact with warm ocean waters. Little is known about the rates of submarine melting and how they may affect glacier dynamics. Here we present measurements of ocean currents, temperature and salinity near the calving fronts of the Eqip Sermia, Kangilerngata Sermia, Sermeq Kujatdleq and Sermeq Avangnardleq glaciers in central West Greenland, as well as ice-front bathymetry and geographical positions. We calculate water-mass and heat budgets that reveal summer submarine melt rates ranging from 0.7±0.2 to 3.9±0.8 m d -1. These rates of submarine melting are two orders of magnitude larger than surface melt rates, but comparable to rates of iceberg discharge. We conclude that ocean waters melt a considerable, but highly variable, fraction of the calving fronts of glaciers before they disintegrate into icebergs, and suggest that submarine melting must have a profound influence on grounding-line stability and ice-flow dynamics. © 2010 Macmillan Publishers Limited. All rights reserved.

Greenland's main outlet glaciers have more than doubled their contribution to global sea level rise over the last decade. Recent work has shown that Greenland's mass loss is still increasing. Here we show that the ice loss, which has been well-documented over southern portions of Greenland, is now spreading up along the northwest coast, with this acceleration likely starting in late 2005. We support this with two lines of evidence. One is based on measurements from the Gravity Recovery and Climate Experiment (GRACE) satellite gravity mission, launched in March 2002. The other comes from continuous Global Positioning System (GPS) measurements from three long-term sites on bedrock adjacent to the ice sheet. The GRACE results provide a direct measure of mass loss averaged over scales of a few hundred km. The GPS data are used to monitor crustal uplift caused by ice mass loss close to the sites. The GRACE results can be used to predict crustal uplift, which can be compared with the GPS data. In addition to showing that the northwest ice sheet margin is now losing mass, the uplift results from both the GPS measurements and the GRACE predictions show rapid acceleration in southeast Greenland in late 2003, followed by a moderate deceleration in 2006. Because that latter deceleration is weak, southeast Greenland still appears to be losing ice mass at a much higher rate than it was prior to fall 2003. In a more general sense, the analysis described here demonstrates that GPS uplift measurements can be used in combination with GRACE mass estimates to provide a better understanding of ongoing Greenland mass loss; an analysis approach that will become increasingly useful as long time spans of data accumulate from the 51 permanent GPS stations recently deployed around the edge of the ice sheet as part of the Greenland GPS Network (GNET). Copyright © 2010 by the American Geophysical Union.

Widespread glacier acceleration has been observed inGreenland in the past few years associated with the thinningof the lower reaches of the glaciers as they terminate inthe ocean. These glaciers thin both at the surface, fromwarm air temperatures, and along their submerged facesin contact with warm ocean waters. Little is known aboutthe rates of submarine melting and how they may affectglacier dynamics. Here we present measurements of oceancurrents, temperature and salinity near the calving fronts ofthe Eqip Sermia, Kangilerngata Sermia, Sermeq Kujatdleq andSermeq Avangnardleq glaciers in central West Greenland, aswell as ice-front bathymetry and geographical positions. Wecalculate water-mass and heat budgets that reveal summersubmarine melt rates ranging from 0.7±0.2 to3.9±0.8md−1.These rates of submarine melting are two orders of magnitudelarger than surface melt rates, but comparable to rates oficeberg discharge. We conclude that ocean waters melt aconsiderable, but highly variable, fraction of the calving frontsof glaciers before they disintegrate into icebergs, and suggestthat submarine melting must have a profound influence ongrounding-line stability and ice-flow dynamics.

We use monthly measurements of time-variable gravity from the GRACE (Gravity Recovery and Climate Experiment) satellite gravity mission to determine the ice mass-loss for the Greenland and Antarctic Ice Sheets during the period between April 2002 and February 2009. We find that during this time period the mass loss of the ice sheets is not a constant, but accelerating with time, i.e., that the GRACE observations are better represented by a quadratic trend than by a linear one, implying that the ice sheets contribution to sea level becomes larger with time. In Greenland, the mass loss increased from 137 Gt/yr in 2002–2003 to 286 Gt/yr in 2007–2009, i.e., an acceleration of −30 ± 11 Gt/yr2 in 2002–2009. In Antarctica the mass loss increased from 104 Gt/yr in 2002–2006 to 246 Gt/yr in 2006–2009, i.e., an acceleration of −26 ± 14 Gt/yr2 in 2002–2009. The observed acceleration in ice sheet mass loss helps reconcile GRACE ice mass estimates obtained for different time periods.

We use 22 monthly GRACE (Gravity Recovery and Climate Experiment) gravity fields to estimate the linear trend in Greenland ice mass during 2002-2004. We recover a decrease in total ice mass of 82 ± 28 km3 of ice per year, consistent with estimates from other techniques. Our uncertainty estimate is dominated by the effects of GRACE measurement errors and errors in our post glacial rebound (PG) correction. The main advantages of GRACE are that it is sensitive to the entire ice sheet, and that it provides mass estimates with only minimal use of supporting physical assumptions or ancillary data. Copyright 2005 by the American Geophysical Union.

There is considerable controversy regarding the long-term strength of continents (Te). While some authors obtain both low and high Te estimates from the Bouguer coherence and suggest that both crust and mantle contribute to lithospheric strength, others obtain estimates of only

Eleven monthly GRACE gravity field solutions are now available for analyses. We show those fields can be used to recover monthly changes in water storage, both on land and in the ocean, to accuracies of 1.5 cm of water thickness when smoothed over 1000 km. The amplitude of the annually varying signal can be determined to 1.0 cm. Results are 30% better for a 1500 km smoothing radius, and 40% worse for a 750 km radius. We estimate the annually varying component of water storage for three large drainage basins (the Mississippi, the Amazon, and a region draining into the Bay of Bengal), to accuracies of 1.0–1.5 cm.

The NASA/DLR satellite gravity mission GRACE, launched in March, 2002, will map the Earth's gravity field at scales of a few hundred km and greater, every 30 days for five years. These data can be used to solve for time-variations in the gravity field with unprecedented accuracy and resolution. One of the many scientific problems that can be addressed with these time-variable gravity estimates, is post glacial rebound (PGR): the viscous adjustment of the solid Earth in response to the deglaciation of the Earth's surface following the last ice age. In this paper we examine the expected sensitivity of the GRACE measurements to the PGR signal, and explore the accuracy with which the PGR signal can be separated from other secular gravity signals. We do this by constructing synthetic GRACE data that include contributions from a PGR model as well as from a number of other geophysical processes, and then looking to see how well the PGR model can be recovered from those synthetic data. We conclude that the availability of GRACE data should result in improved estimates of the Earth's viscosity profile.

Measurements of ice elevation from the Geoscience Laser Altimeter System (GLAS) aboard the Ice, Cloud, and Land Elevation Satellite can be combined with time-variable geoid measurements from the Gravity Recovery and Climate Experiment (GRACE) satellite mission to learn about ongoing changes in polar ice mass and viscoelastic rebound of the lithosphere under the ice sheet. We estimate the accuracy in recovering the spatially varying ice mass trend and postglacial rebound signals for Antarctica, from combining 5 years of simulated GRACE and GLAS data. We obtain root-mean square accuracies of 5.3 and 19.9 mm yr-1 for postglacial rebound and ice mass trend, respectively, when smoothed over 250 km scales. The largest source of error in the combined signals is the effect of the unknown time-variable accumulation on the density of the ice column. To estimate this contribution and so obtain better estimates of ice mass trend and postglacial rebound, we add Global Positioning System (GPS) measurements of vertical velocities as additional constraints. Using an empirical relation between the errors in postglacial rebound and ice mass trend that result from the unknown density variation within the ice column, we are able to solve for all three unknowns in the problem: ice mass trend, postglacial rebound, and the snow compaction trend. The addition of a plausible distribution of GPS measurements reduces the errors in estimates of postglacial rebound and ice mass trend to 3.4 and 15.9 mm yr-1, respectively.

The Gravity Recovery and Climate Experiment (GRACE) satellite mission will resolve temporal variations in gravity orders of magnitude more accurately and with considerably higher resolution than any existing satellite. Effects of atmospheric mass over land will be removed prior to estimating the gravitational field, using surface pressure fields generated by global weather forecast centers. To recover the continental hydrological signal with an accuracy of 1 cm of equivalent water thickness down to scales of a few hundred kilometers, atmospheric pressure must be known to an accuracy of 1 mbar or better. We estimate errors in analyzed pressure fields and the impact of those errors on GRACE surface mass estimates by comparing analyzed fields with barometric surface pressure measurements in the United States and North Africa/Arabian peninsula. We consider (1) the error in 30-day averages of the pressure field, significant because the final GRACE product will average measurements collected over 30-day intervals, and (2) the short-period error in the pressure fields which would be aliased by GRACE orbital passes. Because the GRACE results will average surface mass over scales of several hundred kilometers, we assess the pressure field accuracy averaged over those same spatial scales. The atmospheric error over the 30-day averaging period, which will map directly into GRACE data, is generally < 0.5 mbar. Consequently, analyzed pressure fields will be adequate to remove the atmospheric contribution from GRACE hydrological estimates to subcentimeter levels. However, the short-period error in the pressure field, which would alias into GRACE data, could potentially contribute errors equivalent to 1 cm of water thickness. We also show that given sufficiently dense barometric coverage, an adequate surface pressure field can be constructed from surface pressure measurements alone. Copyright 2001 by the American Geophysical Union.

The Earth's structure and the geod) nam ic processes that originate it are strictly interconnected: an accurate knowledge of the Earth·s structure and its physical properties is fundamental for the understanding of the dynamics of our planet and critical for realistic simulations of seismic ground motion. On the other side the separation and the definition of the different dynamic processes acting in a certain area help to better understand the structure and the physical properties of the investigated region. In recent years much emphasis has been given to the inversion or observed phase and waveform data for the best fitting set or model parameters at a variety of scale-lengths.