Glacial isostatic adjustment in 3-D earth models: Implications for the analysis of tide gauge records along the U.S. east coast

Tide gauge records of recent sea-level change along the U.S. east coast have received significant attention within the literature of glacial isostatic adjustment (GIA). Geographic trends in these tide gauge rates are not reduced by a GIA correction based on a commonly adopted radial viscosity profile (characterized, in particular, by a lower mantle viscosity 1−2×10²¹ Pa s), and this has led to speculation that the residual trends reflect contributions from neotectonics or oceanographic processes. While the trends can be significantly reduced by adopting an Earth model with a stiffer lower mantle, such a model appears to be incompatible with independent constraints from post-glacial decay times in Hudson Bay. We use a finite-element model of the GIA process to investigate whether 3-D viscosity variations superimposed onto the “common” radial viscosity profile may provide a route to reconciling the east coast sea-level trends. We find that the specific 3-D structure we impose has little impact on the geographic trends in the GIA-corrected rates. However, we do find that the imposed lateral variations in lower mantle viscosity introduce a nearly uniform upward shift of 0.5 mm/yr in GIA-induced sea-level rates along the U.S. east coast. Thus, inferences of regional (U.S. east coast) sea-level rise due to modern melting of ice reservoirs, based on tide gauge rates corrected using 1-D GIA models, may be significantly biased by this simplifying assumption.     

Use of GRACE determined secular gravity rates for glacial isostatic adjustment studies in North-America

Monthly geopotential spherical harmonic coefficients from the GRACE satellite mission are used to determine their usefulness and limitations for studying glacial isostatic adjustment (GIA) in North-America. Secular gravity rates are estimated by unweighted least-squares estimation using release 4 coefficients from August 2002 to August 2007 provided by the Center for Space Research (CSR), University of Texas. Smoothing is required to suppress short wavelength noise, in addition to filtering to diminish geographically correlated errors, as shown in previous studies. Optimal cut-off degrees and orders are determined for the destriping filter to maximize the signal to noise ratio. The halfwidth of the Gaussian filter is shown to significantly affect the sensitivity of the GRACE data (with respect to upper mantle viscosity and ice loading history). Therefore, the halfwidth should be selected based on the desired sensitivity.It is shown that increase in water storage in an area south west of Hudson Bay, from the summer of 2003 to the summer of 2006, contributes up to half of the maximum estimated gravity rate. Hydrology models differ in the predictions of the secular change in water storage, therefore even 4-year trend estimates are influenced by the uncertainty in water storage changes. Land ice melting in Greenland and Alaska has a non-negligible contribution, up to one-fourth of the maximum gravity rate.The estimated secular gravity rate shows two distinct peaks that can possibly be due to two domes in the former Pleistocene ice cover: west and south east of Hudson Bay. With a limited number of models, a better fit is obtained with models that use the ICE-3G model compared to the ICE-5G model. However, the uncertainty in interannual variations in hydrology models is too large to constrain the ice loading history with the current data span. For future work in which GRACE will be used to constrain ice loading history and the Earth's radial viscosity profile, it is important to include realistic uncertainty estimates for hydrology models and land ice melting in addition to the effects of lateral heterogeneity.     

Reconstruction of paleocoastlines for the northwestern South China Sea since the Last Glacial Maximum

The range of relative sea level rise in the northwestern South China Sea since the Last Glacial Maximum was over 100 m. As a result, lowland regions including the Northeast Vietnam coast, Beibu Gulf, and South China coast experienced an evolution from land to sea. Based on the principle of reconstructing paleogeography and using recent digital elevation model, relative sea level curves, and sediment accumulation data, this paper presents a series of paleogeographic scenarios back to 20 cal. ka BP for the northwestern South China Sea. The scenarios demonstrate the entire process of coastline changes for the area of interest. During the late glacial period from 20 to 15 cal. ka BP, coastline slowly retreated, causing a land loss of only 1×10⁴ km², and thus the land-sea distribution remained nearly unchanged. Later in 15–10 cal. ka BP coastline rapidly retreated and area of land loss was up to 24×10⁴km², causing lowlands around Northeast Vietnam and South China soon to be underwater. Coastline retreat continued quite rapidly during the early Holocene. From 10 to 6 cal. ka BP land area had decreased by 9×10⁴km², and during that process the Qiongzhou Strait completely opened up. Since the mid Holocene, main controls on coastline change are from vertical crustal movements and sedimentation. Transgression was surpassed by regression, resulting in a land accretion of about 10×10⁴km². Supported by Key Laboratory of Marginal Sea Geology, Chinese Academy of Sciences (Grant No. MSGL0711), the Guangdong Natural Science Foundation (Grant No. 04001309) and Open Fund of the Key Laboratory of Marine Geology and Environment, Chinese Academy of Sciences (Grant No. MGE2007KG04)     

A soil chronosequence in Late Glacial and Neoglacial moraines, Humboldt Glacier, northwestern Venezuelan Andes

Late Glacial and Neoglacial (Little Ice Age) deposits on the Humboldt Massif were analyzed for relative-age dating parameters, including geomorphic and weathering characteristics, geochemical and soil properties. The soil chronosequence, formed in chemically uniform parent materials, provides an important database to study soil evolution in a tropical alpine environment. Extractable and total Fe and Al concentrations, examined to assess their use in relative-age determination, and as paleoenvironmental indicators, provide an important measure of the accumulation and downward profile movement over time of organically-bound Al, ferrihydrite and other crystalline forms (hematite and goethite) of extractable Fe. Ferrihydrite is particularly useful in determining former perched water levels in soils with relation to paleoclimate. The ratios of most Fe extracts are time dependent. The Fe_d/Fe_t ratio, within statistical limits, shows a slow increase from LIA (Little Ice Age) to Late Glacial soils, which closely correlates with other alpine soil studies in the middle latitudes and other tropical alpine locales. Values of Al_d (dithionite) and Al_o (oxalate extractable) generally do not correlate with time; however, Al_p (pyrophosphate extractable) measured against Al_t (total) provides insight on the downward translocation over time of organically-bound Al. Low leaching rates in this chronosequence are further supported by clay mineralogy trends and the geochemical data.     

The Stanley unconformity in Lake Huron basin: evidence for a climate-driven closed lowstand about 7900 <Superscript>14</Superscript>C BP,with similar implications for the Chippewa lowstand in Lake Michigan basin

J.L. Hough in 1962 recognized an erosional unconformity in the upper section of early postglacial lake sediments in northwestern Lake Huron. Low-level Lake Stanley was defined at 70 m below present water surface on the basis of this observation, and was inferred to follow the Main Algonquin highstand and Post-Algonquin lake phases about 10 ¹⁴C ka, a seminal contribution to the understanding of Great Lakes history. Lake Stanley was thought to have overflowed from the Huron basin through the Georgian Bay basin and the glacio-isostatically depressed North Bay outlet to Ottawa and St. Lawrence rivers. For this overflow to have occurred, Hough assumed that post-Algonquin glacial rebound was delayed until after the Lake Stanley phase. A re-examination of sediment stratigraphy in northwestern Lake Huron using seismic reflection and new core data corroborates the sedimentological evidence of Hough’s Stanley unconformity, but not its inferred chronology or the level of the associated lowstand. Erosion of previously deposited sediment, causing the gap in the sediment sequence down to 70 m present depth, is attributed to wave erosion in the shoreface of the Lake Stanley lowstand. Allowing for non-deposition of muddy sediment in the upper 20 m approximately of water depth as occurs in the present Great Lakes, the inferred water level of the Stanley lowstand is repositioned at 50 m below present in northwestern Lake Huron. The age of this lowstand is about 7.9 ± 0.3¹⁴C ka, determined from the inferred ¹⁴C age of the unconformity by radiocarbon-dated geomagnetic secular variation in six new cores. This relatively young age shows that the lowstand defined by Hough’s Stanley unconformity is the late Lake Stanley phase of the northern Huron basin, youngest of three lowstands following the Algonquin lake phases. Reconstruction of uplift histories for lake level and outlets shows that late Lake Stanley was about 25–30 m below the North Bay outlet, and about 10 m below the sill of the Huron basin. The late Stanley lowstand was hydrologically closed, consistent with independent evidence for dry regional climate at this time. A similar analysis of the Chippewa unconformity shows that the Lake Michigan basin also hosted a hydrologically closed lowstand, late Lake Chippewa. This phase of closed lowstands is new to the geological history of the Great Lakes. This is the ninth in a series of ten papers published in this special issue of Journal of Paleolimnology. These papers were presented at the 47th Annual Meeting of the International Association for Great Lakes Research (2004), held at the University of Waterloo, Waterloo, Ontario, Canada. P.F. Karrow and C.F.M Lewis were guest editors of this special issue.     

重庆金佛山石笋的同位素年龄和古气候信息

对采自重庆金佛山梁天湾洞穴的一根长达295mm的石笋进行了高精度的TIM S-U系定年和碳氧稳定同位素测试分析,初步建立了末次冰期晚期( 31. 90～15. 21 ka B. P. )以来长江中游地区较高分辨率的古气候变化序列。石笋大约从31.90kaB. P.开始沉积,至15. 21kaB. P.沉积结束,属于末次冰期晚期的产物,其平均沉积速率为1. 77mm /100a。从石笋的氧同位素曲线中明显地检出了Heinrich事件,并且氧同位素曲线与周边地区的石笋记录有着良好的对比关系。      

GCM-simulated hydrology in the Arctic duringn the past 21,000 years

The Paleoclimates from Arctic Lakes and Estuaries (PALE) project has chosen to conduct high resolution data-model comparisons for the Arctic region at 21 and 10 (calendar) ka BP. The model simulations for 21, 10, and 0 ka BP were conducted with the GENESIS 2.0 GCM. The 10 ka BP simulation was coupled to the EVE vegetation model. The primary boundary conditions differing from present at 21 ka BP were the northern hemisphere ice sheets and lower CO2, and at 10 ka BP were the orbital insolation and smaller northern hemisphere ice sheets. The purpose of this article is to discuss the hydrological consequences of these simulations.At the Last Glacial Maximum (21 ka BP) the large ice sheets over North America and Eurasia and the lower CO2 levels produced a colder climate than present, with less precipitation throughout the Arctic, except where circulation was altered by the ice sheets. At 10 ka BP greater summer insolation resulted in a warmer and wetter Beringia, but conditions remained cold and dry in the north Atlantic sector, in the vicinity of the remnant ice sheets. Less winter insolation at 10 ka BP resulted in colder and drier conditions throughout the Arctic. Precipitation - evaporation generally correlated with precipitation except where changes in the surface type (ice sheets, vegetation at 10 ka BP, or sea level at 21 ka BP) caused large changes in the evaporation rate. The primary hydrological differences (from present) at 21 and 10 ka BP correlated with the temperature differences, which were a direct result of the large-scale boundary condition changes.     

末次盛冰期以来印度尼西亚穿越流对东印度洋上部水体结构的影响

海洋上部水体垂向结构变化对于理解热带海区在全球气候变化中的作用有着重要意义。通过分析印度尼西亚穿越流(ITF)出口处东印度洋帝汶海区SO18480-3孔中的浮游有孔虫表层种Globigerinoides ruber和温跃层种Pulleniatina obliquiloculata壳体氧碳同位素,并借助12个AMS¹⁴C测年数据重建了末次盛冰期(LGM)以来该区温跃层深度和营养盐水平的演化序列。壳体氧同位素(δ¹⁸O)记录表明温跃层古海洋学特征的变化幅度要大于表层海水,其差值(Δδ¹⁸O_(P-G))有效地反映了温跃层深度的变化,即冰消期和晚全新世温跃层较浅,LGM和早中全新世温跃层较深;并揭示出与全新世相比,LGM期间ITF总流量未显著减小,ITF对该区上部水体结构的影响受到了东西太平洋之间不对称性的调节。碳同位素(δ¹³C)记录则表明该区的古海洋学变化在不同程度上受到了南大洋的影响,并受本区上部水体垂向结构的控制,其差值(Δδ¹³_(G-P))在一定程度上反映了该区上部水体营养盐水平的变化。     

末次盛冰期东亚气候的成因检测

在国际古气候模拟比较计划设置的标准试验方案下,首先利用中国科学院大气物理研究所的全球大气环流模式（IAP-AGCM）模拟了末次盛冰期东亚气候状况,然后通过4组数值敏感性试验逐一模拟了大气CO₂浓度、海洋表面温度（SST）和海冰、陆地冰盖和地形、东亚植被变化4项强迫因子的单独气候效应,进而对末次盛冰期东亚气候的成因进行了检测。结果表明,末次盛冰期除华南局部略有升温外,中国年均地表气温显著降低,降温幅度总体上向北增大,青藏高原处存在一个降温中心。其中,SST和海冰变化是华南局部略偏暖的主因,它同时导致了东亚其他区域地表气温的显著降低,特别是在东北亚地区;陆地冰盖和地形变化对于东亚地表气温的显著冷却作用主要体现在东亚的西北部;大气CO₂浓度降低会引起东亚地区0.2～0.9℃的普遍降温;相对而言,东亚植被的降温作用（0.5～1.0℃）主要显现在中国40°N以南的区域。与此同时,SST和海冰变化能引起中国东部年均降水一定程度的减少,而大气CO₂浓度、陆地冰盖和地形、东亚植被单独变化均不会显著影响东亚年均降水的分布状况,然而,上述四项因子的共同变化会通过协同作用引起中国东部年均降水的显著减少,西部地区降水则与现在差别不大。此外,末次盛冰期东亚夏季风的显著减弱源于SST和海冰变化,冬季风变化则可归因于SST和海冰、陆地冰盖和地形的变化。     

南海北部天然气水合物远景区末次冰期以来底栖有孔虫稳定同位素特征及其影响因素

底栖有孔虫碳同位素负偏移是地质记录中天然气水合物释放的重要证据之一.对南海北部西沙海槽和东沙陆坡等天然气水合物远景区XH-27PC和DS-4PC柱状样分别进行顶空气甲烷含量分析、有机碳含量分析、粒度分析和有孔虫氧碳同位素分析.结合碳酸盐含量及AMS 14C测年,揭示研究区末次冰期以来底栖有孔虫的稳定同位素特征.结果显示:西沙海槽BSR区沉积物中甲烷含量较低;底栖有孔虫碳同位素负偏不明显,与顶空气甲烷含量呈弱正相关(R=0.32),与有机碳含量有强负相关(R=-0.82),说明低通量甲烷不足以引起底栖有孔虫碳同位素显著偏移,在无甲烷或甲烷轻微渗漏的环境中有机碳的厌氧氧化是影响底栖有孔虫碳同位素组成的主要因素.东沙陆坡BSR区沉积物中含有大量的甲烷气体;底栖有孔虫氧同位素记录在末次冰期异常偏重,可能与天然气水合物的分解释放有关;同时可识别出多期碳同位素快速负偏事件,其成因很可能是末次冰期海平面下降导致海底沉积物的温度、压力条件发生变化,从而引发水合物甲烷失稳分解,底栖有孔虫吸收富12C的甲烷源碳致使壳体碳同位素负偏移.