Gravity potential spherical harmonic series

Time variable gravity fields, reflecting variations of mass distribution in the system Earth is one of the key parameters to understand the changing Earth. Mass variations are caused either by redistribution of mass in, on or above the Earth's surface or by geophysical processes in the Earth's interior. The first set of observations of monthly variations of the Earth gravity field was provided by the US/German GRACE satellite mission beginning in 2002. This mission is still providing valuable information to the science community. However, as GRACE has outlived its expected lifetime, the geoscience community is currently seeking successor missions in order to maintain the long time series of climate change that was begun by GRACE. Several studies on science requirements and technical feasibility have been conducted in the recent years. These studies required a realistic model of the time variable gravity field in order to perform simulation studies on sensitivity of satellites and their instrumentation. This was the primary reason for the European Space Agency (ESA) to initiate a study on ''Monitoring and Modelling individual Sources of Mass Distribution and Transport in the Earth System by Means of Satellites''. The goal of this interdisciplinary study was to create as realistic as possible simulated time variable gravity fields based on coupled geophysical models, which could be used in the simulation processes in a controlled environment. For this purpose global atmosphere, ocean, continental hydrology and ice models were used. The coupling was performed by using consistent forcing throughout the models and by including water flow between the different domains of the Earth system. In addition gravity field changes due to solid Earth processes like continuous glacial isostatic adjustment (GIA) and a sudden earthquake with co-seismic and post-seismic signals were modelled. All individual model results were combined and converted to gravity field spherical harmonic series, which is the quantity commonly used to describe the Earth's global gravity field. The result of this study is a twelve-year time-series of 6-hourly time variable gravity field spherical harmonics up to degree and order 180 corresponding to a global spatial resolution of 1 degree in latitude and longitude. In this paper, we outline the input data sets and the process of combining these data sets into a coherent model of temporal gravity field changes. The resulting time series was used in some follow-on studies and is available to anybody interested.

The Atlantic Meridional Overturning Circulation during the last glacial

We reconstruct the geometry and strength of the Atlantic Meridional Overturning Circulation during Heinrich Stadial 2 and three Greenland interstadials of the 20-50 ka period based on the comparison of new and published sedimentary 231Pa/230Th data with simulated sedimentary 231Pa/230Th. We show that the deep Atlantic circulation during these interstadials was very different from that of the Holocene. Northern-sourced waters likely circulated above 2500 m depth, with a flow rate lower than that of the present day North Atlantic Deep Water (NADW). Southern-sourced deep waters most probably flowed northwards below 4000 m depth into the North Atlantic basin, and then southwards as a return flow between 2500 and 4000 m depth. The flow rate of this southern-sourced deep water was likely larger than that of the modern Antarctic Bottom Water (AABW). Our results further show that during Heinrich Stadial 2, the deep Atlantic was probably directly affected by a southern-sourced water mass below 2500 m depth, while a slow southward flowing water mass originating from the North Atlantic likely influenced depths between 1500 and 2500 m down to the equator.

Ostracoda in surface sediments of the Arctic Ocean

Ostracodes were studied from deep Arctic Ocean cores obtained during the Arctic 91 expedition of the Polarstern to the Nansen, Amundsen and Makarov Basins, the Lomonosov Ridge, Morris Jesup Rise and Yermak Plateau, in order to investigate their distribution in Arctic Ocean deep water (AODW) and apply these data to paleoceanographic reconstruction of bottom water masses during the Quaternary. Analyses of coretop assemblages from Arctic 91 boxcores indicate the following: ostracodes are common at all depths between 1000 and 4500 m, and species distribution is strongly influenced by water mass characteristics and bathymetry; quantitative analyses comparing Eurasian and Canada Basin assemblages indicate that distinct assemblages inhabit regions east and west of the Lomonosov Ridge, a barrier especially important to species living in lower AODW; deep Eurasian Basin assemblages are more similar to those living in Greenland Sea deep water (GSDW) than those in Canada Basin deep water; two upper AODW assemblages were recognized throughout the Arctic Ocean, one living between 1000 and 1500 m, and the other, having high species diversity, at 1500-3000 m. Downcore quantitative analyses of species' abundances and the squared chord distance coefficient of similarity reveals a distinct series of abundance peaks in key indicator taxa interpreted to signify the following late Quaternary deep water history of the Eurasian Basin. During the Last Glacial Maximum (LGM), a GSDW/AODW assemblage, characteristic of cold, well oxygenated deep water > 3000 m today, inhabited the Lomonosov Ridge to depths as shallow as 1000 m, perhaps indicating the influence of GSDW at mid-depths in the central Arctic Ocean. During Termination 1, a period of high organic productivity associated with a strong inflowing warm North Atlantic layer occurred. During the mid-Holocene, several key faunal events indicate a period of warming and/or enhanced flow between the Canada and Eurasian Basins. A long-term record of ostracode assemblages from kastenlot core PS2200-5 (1073 m water depth) from the Morris Jesup Rise indicates a quasi-cyclic pattern of water mass changes during the last 300 kyr. Interglacial ostracode assemblages corresponding to oxygen isotope stages 1, 5, and 7 indicate rapid changes in dissolved oxygen and productivity during glacial-interglacial transitions.

Heat flow measurements in the Weddell Sea

From January to March 1987, heat flow measurements were tried at four sites (Sites 689, 690, 695, and 696) during ODP Leg 113, in the Weddell Sea, Antarctica. At Site 690 (Maud Rise), a convex upward shaped temperature vs. depth profile was observed. This profile cannot be explained by steady-state conduction through solid materials only. We conclude that the minimum heat flow value at Site 690 is 45 mW/m2. A prominent bottom simulating reflector (BSR) was observed at 600 mbsf at Site 695. However, the observed temperature is too high to explain the BSR as a gas hydrate. The origin of the BSR remains unknown, although it is probably of biogenic origin as observed in the Bering Sea during DSDP Leg 19. After correcting for the effects of sedimentation, heat flow values at Sites 695 and 696 are 69 and 63 mW/m2, respectively. Furthermore, we compiled heat flow data south of 50°S. In the Weddell Sea region, the eastern part shows relatively low heat flow in comparison with the western part, with the boundary between them at about 15°W longitude.

Diatom events and Neogene nannofossils of ODP Hole 201-1225A

Shipboard investigation of magnetostratigraphy and shore-based investigation of diatoms and calcareous nannofossils were used to identify datum events in sedimentary successions collected at Ocean Drilling Program (ODP) Leg 201 Site 1225. The goal was to extend the magnetic record previously studied at the same site, ODP Leg 138 Site 851, and provide a comprehensive age model for Site 1225. Two high-magnetic intensity zones at 0-70 and 200-255 meters below seafloor (mbsf) were correlated with lithologic Subunits IA and IC in Hole 1225A. Subunit IA (0-70 mbsf) contains the magnetic reversal record until the Cochiti Subchronozone (3.8 Ma) and has a sedimentation rate of 1.7 cm/k.y. This agrees with previous work done at Site 851. Subunit IC (200-255 mbsf) was not sampled at Site 851. Diatom and nannofossil biostratigraphy constrained this subunit, and we found it to contain the magnetic reversal record between Subchrons C4n.2r and C5n.2n (8.6-9.7 Ma), yielding a sedimentation rate of 2.7 cm/k.y. Biostratigraphy was used to establish the sedimentation rates within Subunits IB and ID (70-200 mbsf and 255-300 mbsf, respectively). These subunits had higher sedimentation rates (~3.4 cm/k.y.) and coincide with the late Miocene-early Pliocene biogenic bloom event (4.5-7 Ma) and the Miocene global cooling trend (10-15 Ma). High biogenic productivity associated with these subunits resulted in the pyritization of the magnetic signal. In lithologic Subunit ID, basement flow is another factor that may be altering the magnetic signal; however, the good correlation between the biostratigraphy and magnetostratigraphy indicates that the magnetic record was locked-in near the seafloor and suggests the age model is robust.

(Table 1) Characteristics of South Pole personnel according to forced expiratory flow (FEF25-75) response to altitude

The impact of acute altitude exposure on pulmonary function is variable. A large inter-individual variability in the changes in forced expiratory flows (FEFs) is reported with acute exposure to altitude, which is suggested to represent an interaction between several factors influencing bronchial tone such as changes in gas density, catecholamine stimulation, and mild interstitial edema. This study examined the association between FEF variability, acute mountain sickness (AMS) and various blood markers affecting bronchial tone (endothelin-1, vascular endothelial growth factor (VEGF), catecholamines, angiotensin II) in 102 individuals rapidly transported to the South Pole (2835 m). The mean FEF between 25 and 75% (FEF25-75) and blood markers were recorded at sea level and after the second night at altitude. AMS was assessed using Lake Louise questionnaires. FEF25-75 increased by an average of 12% with changes ranging from -26 to +59% from sea level to altitude. On the second day, AMS incidence was 36% and was higher in individuals with increases in FEF25-75 (41 vs. 22%, P = 0.05). Ascent to altitude induced an increase in endothelin-1 levels, with greater levels observed in individuals with decreased FEF25-75. Epinephrine levels increased with ascent to altitude and the response was six times larger in individuals with decreased FEF25-75. Greater levels of endothelin-1 in individuals with decreased FEF25-75 suggest a response consistent with pulmonary hypertension and/or mild interstitial edema, while epinephrine may be upregulated in these individuals to clear lung fluid through stimulation of beta2-adrenergic receptors.

Flow rates, hydraulic conductivities and permeabilities of sediments from ODP Holes 190-1173A and 190-1174B

Understanding the role of fluids in active accretionary prisms requires quantitative knowledge of parameters such as permeability. We report here the results of permeability tests on four samples from Ocean Drilling Program Leg 190 at the Nankai Trough accretionary prism-two from Site 1173 and two from Site 1174. Volcanic ash is present in one of the samples; otherwise, the material is hemipelagic mud. A constant-rate-of-flow technique was used at various effective pressures and rates of flow. The permeability of the four samples ranges between 10**-15 and 10**-18 m**2, with the ash-bearing sample showing the highest values.

Two-way traveltimes of internal layers and ice thickness between Dome C - Vostok - Dome A

Chinese scientists will start to drill a deep ice core at Kunlun station near Dome A in the near future. Recent work has predicted that Dome A is a location where ice older than 1 million years can be found. We model flow, temperature and the age of the ice by applying a three-dimensional, thermomechanically coupled full-Stokes model to a 70 × 70 km**2 domain around Kunlun station, using isotropic non-linear rheology and different prescribed anisotropic ice fabrics that vary the evolution from isotropic to single maximum at 1/3 or 2/3 depths. The variation in fabric is about as important as the uncertainties in geothermal heat flux in determining the vertical advection which in consequence controls both the basal temperature and the age profile. We find strongly variable basal ages across the domain since the ice varies greatly in thickness, and any basal melting effectively removes very old ice in the deepest parts of the subglacial valleys. Comparison with dated radar isochrones in the upper one third of the ice sheet cannot sufficiently constrain the age of the deeper ice, with uncertainties as large as 500 000 years in the basal age. We also assess basal age and thermal state sensitivities to geothermal heat flux and surface conditions. Despite expectations of modest changes in surface height over a glacial cycle at Dome A, even small variations in the evolution of surface conditions cause large variation in basal conditions, which is consistent with basal accretion features seen in radar surveys.