Refine
Year of publication
- 2011 (141) (remove)
Document Type
- Article (111)
- Doctoral Thesis (24)
- Review (2)
- Habilitation Thesis (1)
- Other (1)
- Postprint (1)
- Preprint (1)
Keywords
- Holocene (6)
- Tibetan Plateau (6)
- Site effects (4)
- Holozän (3)
- Pollen (3)
- Wave propagation (3)
- climate change (3)
- Benzene (2)
- Climate change (2)
- Crustal structure (2)
Institute
- Institut für Geowissenschaften (141) (remove)
The Karaburun Peninsula, which is considered part of the Anatolide-Tauride Block of Turkey, contains clastic and carbonate sequences deposited on the northern margin of Gondwana. The Palaeozoic clastic sequence, which is intruded by the Early Triassic granitoid and tectonically overlies a Mesozoic melange sequence, can be divided into three subunits: a lower clastic subunit consisting of a sandstone-shale alternation, an upper clastic subunit consisting of black chert-bearing shales, sandstone and conglomerate, and a Permo-Carboniferous carbonate subunit. The lower Triassic Karaburun I-type granitoid has a high initial Sr-87/Sr-86 ratio (0.709021-0.709168), and low Nd-143/Nd-144 ratio (0.512004-0.512023) and epsilon Nd (-5.34 to -5.70) isotopic values. Geochronological data indicate a crystallization (intrusion) age of 247.1 +/- 2.0 Ma (Scythian). Geochemically, the acidic magmatism reflects a subduction-related continental-arc basin tectonic setting, which can be linked to the opening of the northern branch of Neo-Tethys as a continental back-arc rifting basin on the northern margin of Gondwana. This can be related to the closure through southward subduction of the Palaeotethys Ocean beneath Gondwana.
A sediment core from a closed basin lake (Lake Kuhai) from the semi-arid northeastern Tibetan Plateau was analysed for its pollen record to infer Lateglacial and post glacial vegetation and climatic change. At Lake Kuhai five major vegetation and climate shifts could be identified: (1) a change from cold and dry to relatively warmer and more moist conditions at 14.8 cal ka BP: (2) a shift to conditions of higher effective moisture and a stepwise warmer climate at 13.6 cal ka BP; (3) a further shift with increased moisture but colder conditions at 7.0 cal ka BP; (4) a return to a significantly colder and drier phase at 6.3 cal ka BP; (5) and a change back to relatively moist conditions at 2.2 cal ka BP. To investigate the response of lake ecosystems to climatic changes, statistical comparisons were made between the lake Kuhai pollen record and a formerly published ostracod and sedimentary record from the same sediment core. Furthermore, the pollen and lacustrine proxies from lake Kuhai were compared to a previously published pollen and lacustrine record from the nearby Lake Koucha. Statistical comparisons were done using non-metric multidimensional scaling and Procrustes rotation. Differences between lacustrine and pollen responses within one site could be identified, suggesting that lacustrine proxies are partly influenced by in-lake or local catchment processes, whereas the terrestrial (pollen) proxy shows a regional climate signal. Furthermore, we found regional differences in proxy response between lake Kuhai and Lake Koucha. Both pollen records reacted in similar ways to major environmental changes, with minor differences in the timing and magnitude of these changes. The lacustrine records were very similar in their timing and magnitude of response to environmental changes; however, the nature of change was at times very distinct. To place the current study in the context of Holocene moisture evolution across the Tibetan Plateau, we applied a five-scale moisture index and average link clustering to all available continuous palaeo-climate records from the Tibetan Plateau to possibly find general patterns of moisture evolution on the Plateau. However, no common regional pattern of moisture evolution during the Holocene could be detected. We assign this to complex responses of different proxies to environmental and atmospheric changes in an already very heterogeneous mountain landscape where minor differences in elevation can cause strong variation in microenvironments.
Characterizing the local site response in large cities is an important step towards seismic hazard assessment. To this regard, single station seismic noise measurements were carried out at 146 sites in the northern part of Santiago de Chile. This extensive survey allowed the fundamental resonance frequency of the sedimentary cover, derived from horizontal-to-vertical (H/V) spectral ratios, to be mapped. By inverting the spectral ratios under the constraint of the thickness of the sedimentary cover, known from previous gravimetric measurements, local S-wave velocity profiles have been retrieved. After interpolation between the individual profiles, the resulting high resolution 3D S-wave velocity model allows the entire area, as well as deeper parts of the basin, to be represented in great detail. Since one lithology shows a great scatter in the velocity values only a very general correlation between S-wave velocity in the uppermost 30 m (v(s)(30)) and local geology is found. Local S-wave velocity profiles can serve as a key factor in seismic hazard assessment, since they allow an estimate of the amplification potential of the sedimentary cover. Mapping the intensity distribution of the 27 February 2010 Maule, Chile, event (Mw = 8.8) the results indicate that local amplification of the ground motion might partially explain the damage distribution and encourage the use of the low cost seismic noise techniques for the study of seismic site effects.
The analysis of palaeoclimate time series is usually affected by severe methodological problems, resulting primarily from non-equidistant sampling and uncertain age models. As an alternative to existing methods of time series analysis, in this paper we argue that the statistical properties of recurrence networks - a recently developed approach - are promising candidates for characterising the system's nonlinear dynamics and quantifying structural changes in its reconstructed phase space as time evolves. In a first order approximation, the results of recurrence network analysis are invariant to changes in the age model and are not directly affected by non-equidistant sampling of the data. Specifically, we investigate the behaviour of recurrence network measures for both paradigmatic model systems with non-stationary parameters and four marine records of long-term palaeoclimate variations. We show that the obtained results are qualitatively robust under changes of the relevant parameters of our method, including detrending, size of the running window used for analysis, and embedding delay. We demonstrate that recurrence network analysis is able to detect relevant regime shifts in synthetic data as well as in problematic geoscientific time series. This suggests its application as a general exploratory tool of time series analysis complementing existing methods.
The aim of this paper is to estimate the Hurst parameter of Fractional Gaussian Noise (FGN) using Bayesian inference. We propose an estimation technique that takes into account the full correlation structure of this process. Instead of using the integrated time series and then applying an estimator for its Hurst exponent, we propose to use the noise signal directly. As an application we analyze the time series of the Nile River, where we find a posterior distribution which is compatible with previous findings. In addition, our technique provides natural error bars for the Hurst exponent.
Despite the importance of rhizosphere properties for water flow from soil to roots, there is limited quantitative information on the distribution of water in the rhizosphere of plants.
Here, we used neutron tomography to quantify and visualize the water content in the rhizosphere of the plant species chickpea (Cicer arietinum), white lupin (Lupinus albus), and maize (Zea mays) 12 d after planting.
We clearly observed increasing soil water contents (h) towards the root surface for all three plant species, as opposed to the usual assumption of decreasing water content. This was true for tap roots and lateral roots of both upper and lower parts of the root system. Furthermore, water gradients around the lower part of the roots were smaller and extended further into bulk soil compared with the upper part, where the gradients in water content were steeper.
Incorporating the hydraulic conductivity and water retention parameters of the rhizosphere into our model, we could simulate the gradual changes of h towards the root surface, in agreement with the observations. The modelling result suggests that roots in their rhizosphere may modify the hydraulic properties of soil in a way that improves uptake under dry conditions.
Climatic forcing on channel profiles in the eastern cordillera of the Coroico Region, Bolivia
(2011)
Orographic precipitation has a large impact on channel morphology and rock uplift via a positive feedback to erosion. We show that in the Eastern Cordillera of Bolivia, channel concavities reach their highest values where annual precipitation increases in the downstream direction, exceeding 3000 mm. The steepest channels are upstream of this zone of high concavity, where precipitation rates are <1000 mm yr(-1). Channels exhibit graded forms both upstream and downstream of this transient reach. We conclude that the prolonged effect of orographic erosion and related tectonic uplift is the preservation of channels with extreme concavities in the Eastern Cordillera.
Northwest Europe's largest heather-dominated sandy habitats are located in the nature reserve Luneburger Heide, Germany. Yet, even these appear to be losing their ability to support some of their stenotopic species such as the ladybird spider, Eresus kollari Rossi 1846, and are thus becoming increasingly important for the preservation of these species. The habitat requirements of this endangered spider species were investigated in order to obtain data that will help stabilize the last remnants of the species' population in northwest Germany. Several heathland habitats were surveyed by pitfall trapping during the mate-search period of the males. Two statistical methods were applied: logistic regression and boosted regression trees (BRT). Both methods showed that three habitat variables are of prime relevance in predicting the occurrence of E. kollari: a) thickness of the organic layer (a negative effect), b) soil temperature at a depth of 10 cm, and c) Calluna cover in the herb layer (both have positive effect). Our results show that choppering (removing above-ground biomass and most of O-layer) and burning are likely appropriate heathland management measures for the conservation of E. kollari. Such measures improve the species' habitat quality by creating a heterogenic (small-scaled) heathland structure with suitable microhabitats. As Calluna heathlands show a clear senescence of the dominant heather, it is essential that those habitat patches be conserved. Further measures, such as transfer experiments, are recommended.
Soil moisture at the plot or hill-slope scale is an important link between local vadose zone hydrology and catchment hydrology. However, so far only a few methods are on the way to close this gap between point measurements and remote sensing. One new measurement methodology that could determine integral soil moisture at this scale is the aboveground sensing of cosmic-ray neutrons, more precisely of ground albedo neutrons. The present study performed ground albedo neutron sensing (GANS) at an agricultural field in northern Germany. To test the method it was accompanied by other soil moisture measurements for a summer period with corn crops growing on the field and a later autumn-winter period without crops and a longer period of snow cover. Additionally, meteorological data and aboveground crop biomass were included in the evaluation. Hourly values of ground albedo neutron sensing showed a high statistical variability. Six-hourly values corresponded well with classical soil moisture measurements, after calibration based on one reference dry period and three wet periods of a few days each. Crop biomass seemed to influence the measurements only to minor degree, opposed to snow cover which has a more substantial impact on the measurements. The latter could be quantitatively related to a newly introduced field neutron ratio estimated from neutron counting rates of two energy ranges. Overall, our study outlines a procedure to apply the ground albedo neutron sensing method based on devices now commercially available, without the need for accompanying numerical simulations and suited for longer monitoring periods after initial calibration.
Quantification of subsurface water fluxes based on the one dimensional solution to the heat transport equation depends on the accuracy of measured subsurface temperatures. The influence of temperature probe setup on the accuracy of vertical water flux calculation was systematically evaluated in this experimental study. Four temperature probe setups were installed into a sand box experiment to measure temporal highly resolved vertical temperature profiles under controlled water fluxes in the range of +/- 1.3 md(-1). Pass band filtering provided amplitude differences and phase shifts of the diurnal temperature signal varying with depth depending on water flux. Amplitude ratios of setups directly installed into the saturated sediment significantly varied with sand box hydraulic gradients. Amplitude ratios provided an accurate basis for the analytical calculation of water flow velocities, which matched measured flow velocities. Calculated flow velocities were sensitive to thermal properties of saturated sediment and to temperature sensor spacing, but insensitive to thermal dispersivity equal to solute dispersivity. Amplitude ratios of temperature probe setups indirectly installed into piezometer pipes were influenced by thermal exchange processes within the pipes and significantly varied with water flux direction only. Temperature time lags of small sensor distances of all setups were found to be insensitive to vertical water flux.
The need to increase food production for a growing world population makes an assessment of global agricultural water productivities and virtual water flows important. Using the hydrology and agro-biosphere model LPJmL, we quantify at 0.5 degrees resolution the amount of blue and green water (irrigation and precipitation water) needed to produce one unit of crop yield, for 11 of the world's major crop types. Based on these, we also quantify the agricultural water footprints (WFP) of all countries, for the period 1998-2002, distinguishing internal and external WFP (virtual water imported from other countries) and their blue and green components, respectively. Moreover, we calculate water savings and losses, and for the first time also land savings and losses, through international trade with these products. The consistent separation of blue and green water flows and footprints shows that green water globally dominates both the internal and external WFP (84% of the global WFP and 94% of the external WFP rely on green water). While no country ranks among the top ten with respect to all water footprints calculated here, Pakistan and Iran demonstrate high absolute and per capita blue WFP, and the US and India demonstrate high absolute green and blue WFPs. The external WFPs are relatively small (6% of the total global blue WFP, 16% of the total global green WFP). Nevertheless, current trade of the products considered here saves significant water volumes and land areas (similar to 263 km(3) and similar to 41 Mha, respectively, equivalent to 5% of the sowing area of the considered crops and 3.5% of the annual precipitation on this area). Relating the proportions of external to internal blue/green WFP to the per capita WFPs allows recognizing that only a few countries consume more water from abroad than from their own territory and have at the same time above-average WFPs. Thus, countries with high per capita water consumption affect mainly the water availability in their own country. Finally, this study finds that flows/savings of both virtual water and virtual land need to be analysed together, since they are intrinsically related.
A poorly understood lag time of 15-20 m.y. exists between the initial Arabia-Eurasia continental collision in late Eocene to early Oligocene time and the acceleration of tectonic and sedimentary processes across the collision zone in the early to late Miocene. The late Eocene to Miocene-Pliocene clastic and shallow-marine sedimentary rocks of the Kond, Eyvanekey, and Semnan Basins in the Alborz Mountains (northern Iran) offer the possibility to track the evolution of this orogen in the framework of collision processes. A transition from volcaniclastic submarine deposits to shallow-marine evaporites and terrestrial sediments occurred shortly after 36 Ma in association with reversals in sediment provenance, strata tilting, and erosional unroofing. These events followed the termination of subduction arc magmatism and marked a changeover from an extensional to a contractional regime in response to initiation of continental collision with the subduction of stretched Arabian lithosphere. This early stage of collision produced topographic relief associated with shallow foreland basins, suggesting that shortening and tectonic loading occurred at low rates. Starting from the early Miocene (17.5 Ma), flexural subsidence in response to foreland basin initiation occurred. Fast sediment accumulation rates and erosional unroofing trends point to acceleration of shortening by the early Miocene. We suggest that the lag time between the initiation of continental collision (36 Ma) and the acceleration of regional deformation (20-17.5 Ma) reflects a two-stage collision process, involving the "soft" collision of stretched lithosphere at first and "hard" collision following the arrival of unstretched Arabian continental litho sphere in the subduction zone.
Hydrocarbons can be found in many different habitats and represent an important carbon source for microbes. As fossil fuels, they are also an important economical resource and through natural seepage or accidental release they can be major pollutants. DNA-specific stains and molecular probes bind to hydrocarbons, causing massive background fluorescence, thereby hampering cell enumeration. The cell extraction procedure of Kallmeyer et al. (2008) separates the cells from the sediment matrix. In principle, this technique can also be used to separate cells from oily sediments, but it was not originally optimized for this application. Here we present a modified extraction method in which the hydrocarbons are removed prior to cell extraction. Due to the reduced background fluorescence the microscopic image becomes clearer, making cell identification, and enumeration much easier. Consequently, the resulting cell counts from oily samples treated according to our new protocol are significantly higher than those treated according to Kallmeyer et al. (2008). We tested different amounts of a variety of solvents for their ability to remove hydrocarbons and found that n-hexane and in samples containing more mature oils methanol, delivered the best results. However, as solvents also tend to lyse cells, it was important to find the optimum solvent to sample ratio, at which hydrocarbon extraction is maximized and cell lysis minimized. A volumetric ratio of 1:2-1:5 between a formalin-fixed sediment slurry and solvent delivered highest cell counts. Extraction efficiency was around 30-50% and was checked on both oily samples spiked with known amounts of E. coli cells and oil-free samples amended with fresh and biodegraded oil. The method provided reproducible results on samples containing very different kinds of oils with regard to their degree of biodegradation. For strongly biodegraded oil MeOH turned out to be the most appropriate solvent, whereas for less biodegraded samples n-hexane delivered best results.
Ostracod shells in surface sediments from Ulungur Lake (Xinjiang, China) belong mainly to Limnocythere inopinata as the dominant species, and Candona neglecta and Darwinula stevensoni as accompanying, less abundant taxa. Shells of an additional nine species were recorded only sporadically. The three most abundant ostracods have wide tolerance ranges in terms of salinity, substrate and water depth. The similarly recorded bivalve Pisidium subtruncatum, and the gastropods Gyraulus chinensis and Radix auricularia belong to the most tolerant representatives of the genera. The bivalve and gastropods, in addition to the ostracod assemblage, reflect the fact that Ulungur Lake has experienced strong lake level and salinity variations due to water withdrawal in the catchment and the counteracting diversion of river waters to the lake in recent decades. The substrate in Ulungur Lake is typically fine-grained, apart from the delta region of the Ulungur River channel, which is marked by relatively coarse-grained detrital sediments barren of ostracod shells. This channel was created 40 years ago to divert water to Ulungur Lake and support its local fisheries and recreational facilities. A reassessed Holocene ostracod record from the lake shows that a significantly higher salinity and lower lake level existed in the early Holocene before 6.0 ka in response to the regional climate. In contrast, a higher lake level and lowest salinity is inferred for the late Holocene period between ca. 3.6 and 1.3 ka before present. Afterwards, the lake level declined and salinity increased in response to regional moisture reduction, although conditions similar to the early Holocene lake status were not re-established. Our surface-sediment-derived data provide a baseline for analysis of future environmental variations due to global climate change and regional water management.
In order to provide probabilistic projections of the future evolution of the Atlantic Meridional Overturning Circulation (AMOC), we calibrated a simple Stommel-type box model to emulate the output of fully coupled three-dimensional atmosphere-ocean general circulation models (AOGCMs) of the Coupled Model Intercomparison Project (CMIP). Based on this calibration to idealised global warming scenarios with and without interactive atmosphere-ocean fluxes and freshwater perturbation simulations, we project the future evolution of the AMOC mean strength within the covered calibration range for the lower two Representative Concentration Pathways (RCPs) until 2100 obtained from the reduced complexity carbon cycle-climate model MAGICC 6. For RCP3-PD with a global mean temperature median below 1.0 degrees C warming relative to the year 2000, we project an ensemble median weakening of up to 11% compared to 22% under RCP4.5 with a warming median up to 1.9 degrees C over the 21st century. Additional Greenland meltwater of 10 and 20 cm of global sea-level rise equivalent further weakens the AMOC by about 4.5 and 10 %, respectively. By combining our outcome with a multi-model sea-level rise study we project a dynamic sea-level rise along the New York City coastline of 4 cm for the RCP3-PD and of 8 cm for the RCP4.5 scenario over the 21st century. We estimate the total steric and dynamic sea-level rise for New York City to be about 24 cm until 2100 for the RCP3-PD scenario, which can hold as a lower bound for sea-level rise projections in this region, as it does not include ice sheet and mountain glacier contributions.
A parameterization for the motion of ice-shelf fronts on a Cartesian grid in finite-difference land-ice models is presented. The scheme prevents artificial thinning of the ice shelf at its edge, which occurs due to the finite resolution of the model. The intuitive numerical implementation diminishes numerical dispersion at the ice front and enables the application of physical boundary conditions to improve the calculation of stress and velocity fields throughout the ice-sheet-shelf system. Numerical properties of this subgrid modification are assessed in the Potsdam Parallel Ice Sheet Model (PISM-PIK) for different geometries in one and two horizontal dimensions and are verified against an analytical solution in a flow-line setup.
Coastal uplift and tsunami effects associated to the 2010 M(w)8.8 Maule earthquake in Central Chile
(2011)
On February 27, 2010 at 03:34:08 AM an M(w)8.8 earthquake, with epicenter located off Cobquecura (73.24 degrees W; 36.29 degrees S), severely hit Central Chile. The tsunami waves that followed this event affected the coastal regions between the cities of Valparaiso and Valdivia, with minor effects as far as Coquimbo. The earthquake occurred along the subduction of the Nazca oceanic plate beneath the South American plate. Coseismic coastal uplift was estimated through observations of bleached lithothamnioids crustose coralline algae, which were exposed after the mainshock between 34.13 degrees S and 38.34 degrees S, suggesting the latitudinal distribution of the earthquake rupture. The measured coastal uplift values varied between 240 +/- 20 cm at sites closer to the trench along the western coast of the Arauco peninsula and 15 +/- 10 cm at sites located farther east. A maximum value of 260 +/- 50 cm was observed at the western coast of Santa Maria Island, which is similar to the reported uplift associated with the 1835 earthquake at Concepcion. Land subsidence values on the order of 0.5 m to 1 m evidenced a change in polarity and position of the coseismic hinge at 110-120 km from the trench. In four sites along the coast we observed a close match between coastal uplift values deduced from bleached lithothamnioids algae and GPS measurements. According to field observations tsunami heights reached ea. 14 m in the coastal area of the Maule Region immediately north of the epicenter, and diminished progressively northwards to 4-2 m near Valparaiso. Along the coast of Cobquecura, tsunami height values were inferior to 2-4 m. More variable tsunami heights of 6-8 m were measured at Dichato-Talcahuano and Tirua-Puerto Saavedra, in the Biobio and Arauco regions, respectively, to the south of the epicenter. According to eyewitnesses, the tsunami reached the coast between 12 to 20 and 30 to 45 minutes in areas located closer and faraway from the earthquake rupture zone, respectively. Destructive tsunami waves arrived also between 2.5 and 4.5 hours after the mainshock, especially along the coast of the Biobio and Arauco regions. The tsunami effects were highly variable along the coast, as a result of geomorphological and bathymetric local conditions, besides potential complexities induced by the main shock.
Quantification of total cell abundance is one of the most fundamental parameters in the exploration of subsurface life. Despite all recent advances in molecular techniques, this parameter is usually determined by fluorescence microscopy. In order to obtain reliable and reproducible results, it is important not just to focus on the actual cell enumeration but also to consider the entire chain of processing. Starting with the retrieval of the sample, over subsampling and sample processing to the final step of fluorescence microscopy, there are many potential sources of contamination that have to be assessed and, if possible, avoided. Because some degree of sample contamination will always occur, it is necessary to employ some form of contamination control. Different tracers are available, each one with its specific advantages and drawbacks. In many cases, the problems arise not after the sample has arrived in a well-equipped laboratory with highly trained personnel, but much earlier at the drill site or in a field camp. In this review, I discuss the different aspects of cell enumeration in subsurface sediment, evaluating every step in the long process chain.
The Dead Sea Transform (DST) is a major left-lateral strike-slip fault that accommodates the relative motion between the African and Arabian plates, connecting a region of extension in the Red Sea to the Taurus collision zone in Turkey over a length of about 1100 km. The Dead Sea Basin (DSB) is one of the largest basins along the DST. The DSB is a morphotectonic depression along the DST, divided into a northern and a southern sub-basin, separated by the Lisan salt diapir. We report on a receiver function study of the crust within the multidisciplinary geophysical project, DEad Sea Integrated REsearch (DESIRE), to study the crustal structure of the DSB. A temporary seismic network was operated on both sides of the DSB between 2006 October and 2008 April. The aperture of the network is approximately 60 km in the E-W direction crossing the DSB on the Lisan peninsula and about 100 km in the N-S direction. Analysis of receiver functions from the DESIRE temporary network indicates that Moho depths vary between 30 and 38 km beneath the area. These Moho depth estimates are consistent with results of near-vertical incidence and wide-angle controlled-source techniques. Receiver functions reveal an additional discontinuity in the lower crust, but only in the DSB and west of it. This leads to the conclusion that the internal crustal structure east and west of the DSB is different at the present-day. However, if the 107 km left-lateral movement along the DST is taken into account, then the region beneath the DESIRE array where no lower crustal discontinuity is observed would have lain about 18 Ma ago immediately adjacent to the region under the previous DESERT array west of the DST where no lower crustal discontinuity is recognized.
Wavelet modelling of the gravity field by domain decomposition methods: an example over Japan
(2011)
With the advent of satellite gravity, large gravity data sets of unprecedented quality at low and medium resolution become available. For local, high resolution field modelling, they need to be combined with the surface gravity data. Such models are then used for various applications, from the study of the Earth interior to the determination of oceanic currents. Here we show how to realize such a combination in a flexible way using spherical wavelets and applying a domain decomposition approach. This iterative method, based on the Schwarz algorithms, allows to split a large problem into smaller ones, and avoids the calculation of the entire normal system, which may be huge if high resolution is sought over wide areas. A subdomain is defined as the harmonic space spanned by a subset of the wavelet family. Based on the localization properties of the wavelets in space and frequency, we define hierarchical subdomains of wavelets at different scales. On each scale, blocks of subdomains are defined by using a tailored spatial splitting of the area. The data weighting and regularization are iteratively adjusted for the subdomains, which allows to handle heterogeneity in the data quality or the gravity variations. Different levels of approximations of the subdomains normals are also introduced, corresponding to building local averages of the data at different resolution levels.
We first provide the theoretical background on domain decomposition methods. Then, we validate the method with synthetic data, considering two kinds of noise: white noise and coloured noise. We then apply the method to data over Japan, where we combine a satellite-based geopotential model, EIGEN-GL04S, and a local gravity model from a combination of land and marine gravity data and an altimetry-derived marine gravity model. A hybrid spherical harmonics/wavelet model of the geoid is obtained at about 15 km resolution and a corrector grid for the surface model is derived.