Monitoring absolute vertical land motions and absolute sea-level changes from GPS and tide gauges data over French Polynesia

In this study,we estimate the absolute vertical land motions at three tidal stations with collocated Global Navigation Satellite System(GNSS)receivers over French Polynesia during the period 2007-2020,and obtain,as ancillary results,estimates of the absolute changes in sea level at the same locations.To verify our processing approach to determining vertical motion,we first modeled vertical motion at the International GNSS Service(IGS)THTI station located in the capital island of Tahiti and compared our estimate with previous independent determinations,with a good agreement.We obtained the following estimates for the vertical land motions at the tide gauges:Tubuai island,Austral Archipelago-0.92±0.17 mm/yr,Vairao village,Tahiti Iti:-0.49±0.39 mm/yr,Rikitea,Gambier Archipelago-0.43±0.17 mm/yr.The absolute variations of the sea level are:Tubuai island,Austral Archipelago 5.25±0.60 mm/yr,Vairao village,Tahiti Iti:3.62±0.52 mm/yr,Rikitea,Gambier Archipelago 1.52±0.23 mm/yr.We discuss these absolute values in light of the values obtained from altimetric measurements and other means in French Polynesia.

Exploring Variability in Sea Level at a Tide Gauge Station through Control Charts

Monitoring temporal changes in sea level is important in assessing coastal risk.Sea level anomalies at a tide gauge station,if kinematically conceived,include systematic variations such as trend,acceleration,periodic oscillations,and random disturbances.Among them,the non-stationary nature of the random sea level variations of known or unknown origin at coastal regions has been long recognized by the sea level community.This study proposes the analyses of subgroups of random residual statistics of a rigorously formulated kinematic model solution of tide gauge variations using X-bar and S control charts.The approach is demonstrated using Key West,Florida tide gauge records.The mean and standard errors of 5-year-long subgroups of the residuals revealed that sea level changes at this location have been progressively intensifying from 1913 to the present.Increasing oscillations in sea level at this locality may be attributed partly to the thermal expansion of seawater with increasing temperatures causing larger buoyancy-related sea level fluctuations as well as the intensification of atmospheric events including wind patterns and the impact of changes in inverted barometer effects that will alter coastal risk assessments for the future.

The Impact of Sea Level Rise on Roadway Design and Evacuation Routes in Delaware

As the global temperature continues to increase, the sea level continues to rise at a rapid rate that has never been seen before. This becomes an issue for many facets of life but one of the most impacted is the transportation infrastructure. Many people living in low elevation coastal areas can become trapped by flooding with no way in or out. With Delaware being a coastal state, this would affect a large portion of the population and will have detrimental effects over time if nothing is done to combat sea level rise. The issue with sea level rise in transportation is that once the roads become flooded, they become virtually unusable and detour routes would be needed. If all the roads in a coastal area were to be affected by sea level rise, the options for detours would become limited. This article looks at direct solutions to combat sea level rise and indirect solutions that would specifically help transportation infrastructure and evacuation routes in Delaware. There is not one solution that can fix every problem, so many solutions are laid out to see what is applicable to each affected area. Some solutions include defense structures that would be put close to the coast, raising the elevation of vulnerable roads throughout the state and including pumping stations to drain the water on the surface of the road. With an understanding of all these solutions around the world, the ultimate conclusion came in the form of a six-step plan that Delaware should take in order to best design against sea level rise in these coastal areas.

Arctic Sea Level Variability from Oceanic Reanalysis and Observations

Quantifying the contributions to Arctic sea level(ASL)variability is critical to understand how the Arctic is responsing to ongoing climate change.Here,we use Ocean Reanalysis System 5(ORAS5)reanalysis data and tide gauge and satellite altimetry observations to quantify contributions from different physical processes on the ASL variability.The ORAS5 reanalysis shows that the ASL is rising with a trend of 2.5±0.3 mm yr−1(95%confidence level)over 1979-2018,which can be attributed to four components:(i)the dominant component from the global sea level increase of 1.9±0.5 mm yr−1,explaining 69.7%of the total variance of the ASL time series;(ii)the Arctic Oscillation-induced mass redistribution between the deep central basin and shallow shelves,with no significant trend and explaining 6.3%of the total variance;(iii)the steric sea level increase centering on the Beaufort Gyre region with a trend of 0.5±0.1 mm yr−1 and explaining 29.1%of the total variance of the ASL time series;and(iv)the intrusion of Pacific water into the Arctic Ocean,with no significant trend and contributing 14.2%of the total ASL variability.Furthermore,the dramatic sea ice melting and the larger area of open water changes the impact of the large-scale atmospheric forcing on the ASL variability after 1995,and the ocean dynamic circulation plays a more important role in the ASL variability.

Devonian Sea-Level Change Rhythms in South China and Coupling Relationships Among the Earth-spheres

Twenty orthosequences and their corresponding sea-level change (SLC) cycles have been recognized in the Devonian overlying the Caledonian unconformity, of which 9, 5.5 and 5.5 occurred in the Lower, Middle and Upper Devonian respectively. They can be grouped into 4 orthosequence sets, in which the maximum flooding surfaces lie in the sulcutus Zone (D12), perbonus Zone (D13), Middle and Upper varcus Zone (D22) and gigas Zone (D21) respectively. Four instant palaeogeographical reconstructions of South China have been made in the Emsian and Givetian. Devonian sea-level change rhythms of South China can be divided into 3 categories: the autorhythmic, the worldwide and regional allorhythmic, and the coupling-rhythmic. They developed respectively in the Famennian, Pragian, Eifelian, Lochkovian, Emsian, Givetian, Frasnian and the F / F (between the Frasnian and Famennian) event. The cause of the worldwide allorhythmic SLC of the Pragian and Eifelian under comparatively dry, warm and tranquil conditions may be related to the pulsating expanding and contracting of the oceanic basin volume or the earth volume pulsation, rather than the common glaciation and plate tectonism. The coupling-rhythmic SLC related to the F/ F event is a sensitive indicator of the interaction between terrestrial and extraterrestrial factors, and coordinated action among the earth-spheres.

Terrestrial Temperature, Sea Levels and Ice Area Links with Solar Activity and Solar Orbital Motion

This paper explores the links between terrestrial temperature, sea levels and ice areas in both hemispheres with solar activity indices expressed through averaged sunspot numbers together with the summary curve of eigenvectors of the solar background magnetic field (SBMF) and with changes of Sun-Earth distances caused by solar inertial motion resulting from the gravitation of large planets in the solar system. Using the wavelet analysis of the GLB and HadCRUTS datasets two periods: 21.4 and 36 years in GLB, set and the period of about 19.6 years in the HadCRUTS are discovered. The 21.4-year period is associated with variations in solar activity defined by the summary curve of the largest eigenvectors of the SBMF. A dominant 21.4-year period is also reported in the variations of the sea level, which is linked with the period of 21.4 years detected in the GLB temperature and the summary curve of the SBMF variations. The wavelet analysis of ice and snow areas shows that in the Southern hemisphere, it does not show any links to solar activity periods while in the Northern hemisphere, the ice area reveals a period of 10.7 years equal to a usual solar activity cycle. The TSI in March-August of every year is found to grow with every year following closely the temperature curve, because the Sun moves closer to the Earth orbit owing to gravitation of large planets (solar inertial motion, SIM), while the variations of solar radiation during a whole year have more steady distribution without a sharp TSI increase during the last two centuries. The additional TSI contribution caused by SIM is likely to secure the additional energy input and exchange between the ocean and atmosphere.

Organic Carbon Deposition on the Inner Shelf of the East China Sea Constrained by Sea Level and Climatic Changes Since the Last Deglaciation

The East China Sea(ECS),which is located in the transitional zone between land and ocean,is the main site for the burial of sedimentary organic carbon.Despite good constraints of the modern source to the sinking process of organic carbon,its fate in response to changes in climate and sea level since the last deglaciation remains poorly understood.We aim to fill this gap by presenting a high-resolution sedimentary record of core EC2005 to derive a better understanding of the evolution of the depositional environment and its control on the organic deposition since 17.3 kyr.Our results suggest that sedimentary organic carbon was deposited in a terrestrial environment before the seawater reached the study area around 13.1 kyr.This significant transition from a terrestrial environment to a marine environment is reflected by the decrease in TOC/TN and TOC/TS ratios,which is attributed to deglacial sea level rise.The sea level continued to rise until it reached its highstand at approximately 7.3 kyr when the mud depocenter was developed.Our results further indicate that the deposition of the sedimentary organic carbon could respond quickly to abrupt cold events,including the Heinrich stadial 1 and the Younger Dryas during the last deglaciation,as well as‘Bond events'during the Holocene.We propose that the rapid response of the organic deposition to those cold events in the northern hemisphere is linked to the East Asian winter monsoon.These new findings demonstrate that organic carbon deposition and burial on the inner shelf could effectively document sea level and climatic changes.

A new decomposition model of sea level variability for the sea level anomaly time series prediction

Rising sea level is of great significance to coastal societies;predicting sea level extent in coastal regions is critical.When carrying out predictions,the subsequences obtained using decomposition methods may exhibit a certain regularity and therefore can provide multidimensional information that can be used to improve prediction models.Traditional decomposition methods such as seasonal and trend decomposition using Loess(STL)focus mostly on the fluctuating trend of time series and ignore its impact on prediction.Methods in the signal decomposition domain,such as variational mode decomposition(VMD),have no physical significance.In response to the above problems,a new decomposition method for sea level anomaly time series prediction(DMSLAP)is proposed.With this method,the trend term in a time series can be isolated and the effects of abnormal sea level change behaviors can be attenuated.We decompose multiperiod characteristics using this method while maintaining the smoothness of the analyzed series.Satellite altimetry data from 1993 to 2020 are used in experiments conducted in the study area.The results are then compared with predictions obtained using existing decomposition methods such as the STL and VMD methods and time varying filtering based on empirical mode decomposition(TVF-EMD).The performance of DMSLAP combined with a prediction method resulted in optimal sea level anomaly(SLA)predictions,with a minimum root mean square error(RMSE)of 1.40 cm and a maximum determination coefficient(R^(2))of 0.93 during 2020.The DMSLAP method was more accurate when predicting 1-year data and 3-year data.The TVF-EMD and DMSLAP methods had comparable accuracies,and the periodic term decomposed by the DMSLAP method was more in line with the actual law than that derived using the TVF-EMD method.Thus,DMSLAP can decompose SLA time series better than existing methods and is an effective tool for obtaining short-term SLA prediction.

Improving the understanding of the influencing factors on sea level based on wavelet coherence and partial wavelet coherence

Relationship between sea level change and a single climate indicator has been widely discussed.However,few studies focused on the relationship between monthly mean sea level(MMSL)and several key impact factors,including CO_(2) concentration,sea ice area,and sunspots,on various time scales.In addition,research on the independent relationship between climate factors and sea level on various time scales is lacking,especially when the dependence of climate factors on Nino 3.4 is excluded.Based on this,we use wavelet coherence(WC)and partial wavelet coherence(PWC)to establish a relationship between MMSL and its influencing factors.The WC results show that the influence of climate indices on MMSL has strong regional characteristics.The significant correlation between Southern Hemisphere sea ice area and MMSL is opposite to that between Northern Hemisphere sea ice area and MMSL.The PWC results show that after removing the influence of Nino 3.4,the significant coherent regions of the Pacific Decadal Oscillation(PDO),Dipole Mode Index(DMI),Atlantic Multidecadal Oscillation(AMO),and Southern Oscillation Index(SOI)decrease to varying degrees on different time scales in different regions,demonstrating the influence of Nino 3.4.Our work emphasizes the interrelationship and independent relationship between MMSL and its influencing factors on various time scales and the use of PWC and WC to describe this relationship.The study has an important reference significance for selecting the best predictors of sea level change or climate systems.

Prediction of the joint impacts of sea level rise and land development on distribution patterns of mangrove communities

Mangrove distribution along shorelines shows distinct zonation patterns;thus,different communities may face various influences from sea level rise(SLR)and land use.However,long-term change predictions are usually based only on the total extent of mangroves.Few studies have revealed how SLR and land development such as agriculture,aquaculture,and urbanization jointly affect different intertidal mangrove communities.This study proposed a novel framework combining SLAMM(Sea Level Affecting Marshes Model)and the CLUE-S(Conversion of Land Use and its Effect at Small regional extent)model to assess the potential impacts on upper and lower intertidal mangrove communities.Maoweihai in Guangxi,China,was selected as the study area and the potential impacts from the squeeze effect and mangrove expansion potential were evaluated.We established three scenarios combining SLR and land use patterns to predict mangrove coverage projections by 2070.The results showed that,under a single SLR driver,the upper intertidal mangroves would be more adaptive to rapid SLR than the lower intertidal mangroves.However,under the combined influence of the two drivers,the upper intertidal mangroves would experience larger squeeze effects than the lower intertidal mangroves,with up to 80.5%of suitable habitat lost.Moreover,the expansion potential of upper intertidal mangroves would be considerably more limited than that of lower intertidal mangroves.The length of the expandable habitat patch boundary of upper intertidal mangroves only reached 1.4–1.8 km,while that of the lower intertidal mangroves reached up to99.2–111.2 km.Further,we found that aquaculture ponds and cropland are the top two land development types that could occupy suitable habitat and restrict the mangrove expansion potential.Our results highlight that timely improvement of land use policies to create available landward accommodation space for mangrove migration is essential to maintain the coverage and diversity of mangrove communities under SLR.The proposed method can be a helpful tool for adaptive mangrove conservation and management under climate change.

Potential effects of sea level rise on the soil-atmosphere green-house gas emissions in Kandelia obovata mangrove forests

Mangrove forests are under the stress of sea level rise(SLR)which would affect mangrove soil biogeochemistry.Mangrove soils are important sources of soil-atmosphere greenhouse gas(GHG)emissions,including carbon dioxide(CO_(2)),methane(CH_(4))and nitrous oxide(N_(2)O).Understanding how SLR influences GHG emissions is critical for evaluating mangrove blue carbon capability.In this study,potential effects of SLR on the GHG emissions were quantified through static closed chamber technique among three sites under different intertidal elevations,representing tidal flooding situation of SLR values of 0 cm,40 cm and 80 cm,respectively.Compared with Site SLR 0 cm,annual CO_(2) and N_(2)O fluxes decreased by approximately 75.0%and 27.3%due to higher soil water content,lower salinity and soil nutrient environments at Site SLR 80 cm.However,CH_(4) fluxes increased by approximately 13.7%at Site SLR 40 cm and 8.8%at Site SLR 80 cm because of lower salinity,higher soil water content and soil pH.CO_(2)-equivalent fluxes were 396.61 g/(m^(2)·a),1423.29 g/(m^(2)·a)and 1420.21 g/(m^(2)·a)at Sites SLR 80 cm,SLR 40 cm and SLR 0 cm,respectively.From Site SLR 0 cm to Site SLR 80 cm,contribution rate of N_(2)O and CH_(4) increased by approximately 7.42%and 3.02%,while contribution rate of CO_(2) decreased by approximately 10.44%.The results indicated that warming potential of trace CH_(4) and N_(2)O was non-negligible with SLR.Potential effects of SLR on the mangrove blue carbon capability should warrant attention due to changes of all three greenhouse gas fluxes with SLR.

Evidences of the expanding Earth from space-geodetic data over solid land and sea level rise in recent two decades

According to the space-geodetic data recorded at globally distributed stations over solid land spanning a period of more than 20-years under the International Terrestrial Reference Frame 2008,our previous estimate of the average-weighted vertical variation of the Earth＇s solid surface suggests that the Earth＇s solid part is expanding at a rate of 0.24 ± 0.05 mm/a in recent two decades.In another aspect,the satellite altimetry observations spanning recent two decades demonstrate the sea level rise（SLR） rate 3.2 ± 0.4 mm/a,of which1.8 ± 0.5 mm/a is contributed by the ice melting over land.This study shows that the oceanic thermal expansion is 1.0 ± 0.1 mm/a due to the temperature increase in recent half century,which coincides with the estimate provided by previous authors.The SLR observation by altimetry is not balanced by the ice melting and thermal expansion,which is an open problem before this study.However,in this study we infer that the oceanic part of the Earth is expanding at a rate about 0.4 mm/a.Combining the expansion rates of land part and oceanic part,we conclude that the Earth is expanding at a rate of 0.35 ± 0.47 mm/a in recent two decades.If the Earth expands at this rate,then the altimetry-observed SLR can be well explained.

Multi-mode Multi-frequency GNSS-IR Combination System for Sea Level Retrieval

With the development of Global Navigation Satellite Systems(GNSS),geodetic GNSS receivers have been utilized to monitor sea levels using GNSS-Interferometry Reflectometry(GNSS-IR)technology.The multi-mode,multi-frequency signals of GPS,GLONASS,Galileo,and Beidou can be used for GNSS-IR sea level retrieval,but combining these retrievals remains problematic.To address this issue,a GNSS-IR sea level retrieval combination system has been developed,which begins by analyzing error sources in GNSS-IR sea level retrieval and establishing and solving the GNSS-IR retrieval equation.This paper focuses on two key points:time window selection and equation stability.The stability of the retrieval combination equations is determined by the condition number of the coefficient matrix within the time window.The impact of ill-conditioned coefficient matrices on the retrieval results is demonstrated using an extreme case of SNR data with only ascending or descending trajectories.After determining the time window and removing ill-conditioned equations,the multi-mode,multi-frequency GNSS-IR retrieval is performed.Results from three International GNSS Service(IGS)stations show that the combination method produces high-precision,high-resolution,and high-reliability sea level retrieval combination sequences.

Evaluating the Influence of Sea Level Rise on Beel Kapalia’s Livelihood and Local Adaptation Strategies: Perspectives from the Local Community

Bangladesh is vulnerable to climate change-induced sea level rise due to its location and socioeconomic position. The study examines the Beel Kapalia region in polder no. 24 of the Monirampur upazila of Jessore district, Khulna division. To assess local attitudes on sea level rise-related permanent flooding, Kapalia, Monoharpur, Nehalpur, Balidaha, and Panchakori were polled. This flooding has disrupted residents’ lifestyles, making them vulnerable to increasing sea levels. Viability and adaptability were assessed using livelihood capitals. Participants’ thoughts and knowledge about their resilience in several livelihood factors were gathered using participatory rural appraisal (PRA) instruments and a questionnaire survey in the area. Major discoveries include the impact of permanent floods on Beel Kapalia’s livelihoods, vulnerability and resilience assessments in numerous villages, and community viewpoints on regional adaptation methods to mitigate these consequences. The study found that a sustained 30.5 cm inundation would reduce local human, natural, physical, financial, and social capital resilience to 69.6%, 30.7%, 69.1%, 68.9%, and 69.1%. A constant 61 cm inundation would lower resistance to 40.9%, 8.7%, 42.4%, 45.6%, and 43.8%. Residents believe they can weather a 30.5 cm inundation with local adaptation measures, but if the water level rises to 61 cm, they may be displaced.

Can a global mean sea-level rise reduce the Last Interglacial model-data mismatch in East Asia?

末次间冰期有着丰富的重建和模拟资料,为研究未来温暖气候提供了一个理想的参考.然而,关于末次间冰期的东亚气候,模拟与重建的结果间长期存在着不匹配的情况,模拟结果普遍较重建结果更为冷干。本研究利用挪威地球系统模式(NorESM1-F),探讨了在末次间冰期模拟试验中纳入全球平均海平面上升能否减少模式-数据的不匹配.该试验结果表明,海平面上升情况下东亚地区会产生一定的增温增湿效应,但不足以消除模式-数据不匹配.基于这些结果,作者探讨了其它可能造成不匹配的因素以供进一步研究.

Performance of CMIP6 models in simulating the dynamic sea level:Mean and interannual variance

本研究采用卫星测高数据与第六次国际耦合模式比较计划(CMIP6)海平面动力进行对比,重点针对40S-40N地区的动力海平面(DSL),评估了模式对其平均态与年际变率的综合模拟能力,结果表明,对于DSL平均态的模拟,模式与观测结果非常吻合,模式之间的差异较小.其中,副热带北大西洋是模拟偏差和模式间差异较为显著的区域,对于DSL年际变率的模拟,模式之间保持较高的一致性,但是,模式与观测结果存在明显差异,模式普遍低估了DSL的年际方差;其中,误差大值区域出现在副热带西边界流附近,模式分辨率会影响CMIP6对中小尺度海洋过程的重现能力,这可能是导致CMIP6历史模拟出现误差的原因之一.

Contributors to tidal duration asymmetry with varied coastline configurations on western shelf of Yellow Sea

Coastal management in China is confronted with an urgent choice between natural restoration and maintenance of existing seawalls and reclaimed land for economic development.A key criterion for making this decision is the resilience to coastal flooding,which depends on the ability to predict tidal level.Tidal duration asymmetry(TDA)is a key parameter in determination of the arrival and duration of flood tides.This study selected the western inner shelf of the Yellow Sea(WYS)as the study area and investigated the responses of TDA to different shoreline configurations and relative sea level rise.The responses of TDA to shoreline reconstruction yielded spatial variability locally and remotely.In the nearshore area,the responses of TDA to the complex ocean environment mainly originated from the combined functions of reflection,bottom friction,and advection,which controlled the energy transfer from M2 or S2 constituents to their overtides or compound tides.The sensitivity of TDA to coastline typologies was not limited to coastal waters but could stretch over the entire inner shelf.The vulnerability of tidal responses was due to the displacement of the M2 amphidrome of the Kelvin wave on the WYS,which in turn changed tidal energy fluxes over the regime.The relative sea level rise could intensify the feedback of TDA to seawalls and land reclamation.

Impacts of artificial dams on terrestrial water storage changes and the Earth's elastic load response during 1950-2016: A case study of medium and large reservoirs in Chinese mainland

The construction of dams for intercepting and storing water has altered surface water distributions, landsea water exchanges, and the load response of the solid Earth. The lack of accurate estimation of reservoir properties through the land surface and hydrological models can lead to water storage simulation and extraction errors. This impact is particularly evident in many artificial reservoirs in China. The study aims to comprehensively assess the spatiotemporal distribution and trends of water storage in medium and large reservoirs(MLRs) in Chinese mainland during 1950-2016, and to investigate the gravity,displacement, and strain effects induced by the reservoir mass concentration using the load elasticity theory. In addition, the impoundment contributions of MLRs to the relative sea level changes were assessed using a sea-level equation. The results show impoundment increases in the MLRs during1950-2016, particularly in the Yangtze River(Changjiang) and southern basins, causing significant elastic load effects in the surrounding areas of the reservoirs and increasing the relative sea level in China's offshore. However, long-term groundwater estimation trends are overestimated and underestimated in the Yangtze River and southwestern basins, respectively, due to the neglect of the MLRs impacts or the uncertainty of the hydrological model's output(e.g., soil moisture, etc.). The construction of MLRs may reduce the water mass input from land to the ocean, thus slowing global sea level rise. The results of the impact of human activities on the regional water cycle provide important references and data support for improving the integration of hydrological models, evaluating Earth's viscoelastic responses under longterm reservoir storage, enhancing in-situ and satellite geodetic measurements, and identifying the main factors driving sea level changes.

Earth-Moon Barycentre Excursions and Anomalous Quaternary Sea Level Highstands

Plate tectonics is driven by Earth-Moon barycentre shifts in the lower mantle. The eastern Canary Islands have geographic and geological conditions derived from the movements of the Central American plates. Some features of these islands are influenced by the rotation of the Earth from west to east in the evolution of the marine currents that surround them and the opening of the North Atlantic to the North Pole with little dependence of the glacial isostatic adjustment (GIA). In addition, their position with respect to the Tropic of Cancer and the African continent affect the north-south and east-west climatic change dynamics and their tectonic stability respectively. Dated lavas contain marine and aeolian deposits and some of the Pleistocene marine deposits indicate higher sea level in cooler circumstances, which is anomalous. Relating those marine deposits produced during the warmest interglacial, the last interglacial and the Holocene with their equivalents in the Southern Hemisphere, they reflect shifts in the barycentre. Thanks to Holocene radiocarbon, topographic and day length data and alkenone temperature, we describe a mechanism by which the oscillation of the Moon’s inclination (and declination) reaches extreme values (14ºand 34ºabout 4.9ºmore than current values) approximately every 1450 years. These values occur when there is a harmonic distortion in surface areas of the Earth’s crust as response associated with oscillations in the displacements of the barycentre of the Earth-Moon system. As the declination influences the movement of oceanic waters, there is also a relationship with the Bond Events of the North Atlantic, of unknown cause until now.

Climate-environmental Deteriorations in a Greenhouse Earth System:Causes and Consequences of Short-Term Cretaceous Sea-Level Changes(a Report on IGCP 609)

Anthropogenic global warming and resulting sea-level rise in response to enhanced atmospheric greenhouse gases and melting of the Earth’s continental ice shields have become issues of continuously growing interest for the scientific community as well as the public,pointing to threads of societies in a future greenhouse Earth System.