Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “GCM”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 559 records · Page 31

A new spatially distributed added value index for regional climate models: the EURO-CORDEX and the CORDEX-CORE highest resolution ensembles

The added value of using regional climate models (RCMs) to downscale data from general circulation models (GCMs) has often been questioned and researched. Although several studies have used different methods to identify (and in some cases quantify) the added value, there is still a need to find a general metric that quantifies the added value of any variable. This paper builds on past studies to propose a new metric of added value in the simulation of present-day climate which measures the difference in the probability density functions (PDFs) at each grid-cell between a model and an observation source, and then compares the results of the RCM and GCM in order to spatially compute the added value index. The same method is also adapted to quantify the climate change downscaling signal in a way that is consistent with the present-day metric. These new metrics are tested on the daily precipitation output from the EURO-CORDEX and CORDEX-CORE projection ensembles and reveal an overall positive added value of RCMs, especially at the tail-end of the distribution. Higher added value is obtained in areas of complex topography and coast-lines, as well as in tropical regions. Areas with large added value in present-day climate are consistent with areas of significant climate change downscaling signal in the RCP 8.5 far future simulations, and when the analysis is repeated at a low-resolution. The use of different resolution observations shows that the added value tends to decrease when models are compared to low-resolution observation datasets.

54 ENVIRONMENTAL SCIENCES↗

Assessment of CMIP6 Performance and Projected Temperature and Precipitation Changes Over South America

We evaluate the performance of a large ensemble of Global Climate Models (GCMs) from the Coupled Model Intercomparison Project Phase 6 (CMIP6) over South America for a recent past reference period and examine their projections of twenty-first century precipitation and temperature changes. The future changes are computed for two time slices (2040–2059 and 2080–2099) relative to the reference period (1995–2014) under four Shared Socioeconomic Pathways (SSPs, SSP1–2.6, SSP2–4.5, SSP3–7.0 and SSP5–8.5). The CMIP6 GCMs successfully capture the main climate characteristics across South America. However, they exhibit varying skill in the spatiotemporal distribution of precipitation and temperature at the sub-regional scale, particularly over high latitudes and altitudes. Future precipitation exhibits a decrease over the east of the northern Andes in tropical South America and the southern Andes in Chile and Amazonia, and an increase over southeastern South America and the northern Andes—a result generally consistent with earlier CMIP (3 and 5) projections. However, most of these changes remain within the range of variability of the reference period. In contrast, temperature increases are robust in terms of magnitude even under the SSP1–2.6. Future changes mostly progress monotonically from the weakest to the strongest forcing scenario, and from the mid-century to late-century projection period. There is an increase in the seasonality of the intra-annual precipitation distribution, as the wetter part of the year contributes relatively more to the annual total. Furthermore, an increasingly heavy-tailed precipitation distribution and a rightward shifted temperature distribution provide strong indications of a more intense hydrological cycle as greenhouse gas emissions increase. The relative distance of an individual GCM from the ensemble mean does not substantially vary across different scenarios. We found no clear systematic linkage between model spread about the mean in the reference period and the magnitude of simulated sub-regional climate change in the future period. Overall, these results could be useful for regional climate change impact assessments across South America.

54 ENVIRONMENTAL SCIENCES↗

Hex-C 558 : A new porous metallic carbon allotrope for lithium-ion battery anode

Due to the limited specific capacity of graphite anodes in lithium-ion batteries (LIBs), it is imperative to find alternatives with better performance. Here in this work, we propose, for the first time, a three dimensional (3D) porous metallic 5-5-8 carbon monolith with hexagonal lattice as a potential candidate. With 5-5-8 carbon nanoribbon as the building block, this material, named Hex-C558, goes beyond graphene-nanoribbon-based 3D porous carbon structures. Using first-principles calculations, we show that Hex-C 558 is not only dynamically and thermally stable, but also is energetically more favorable than many other theoretically predicted carbon allotropes. More importantly, Hex-C 558 is metallic with ordered ionic conducting channels and possesses a low mass density of 1.05 gcm -3 , exhibiting great potential for ion-battery applications. As an anode for LIBs, Hex-C558 possesses a large specific capacity of 591 mAhg -1 , low diffusion energy barrier of 0.27 eV (at low Li concentration) and 0.52 eV (at high Li concentration), a low open-circuit voltage of 0.51 V, and a small volume change of 2.4%. This work provides a new route for the design and synthesis of novel carbon materials for battery applications by using pentagon-based building units.

36 MATERIALS SCIENCE↗

Non-stationary precipitation design standards for stormwater infrastructure modernization at USAF installations

The resilience of defense infrastructure systems to a changing climate is critical for national security. Climate induced recurrent flooding is already impacting over 20 U.S. Air Force installations, underscoring the urgency of revisiting precipitation standards and stormwater infrastructure design. Despite growing scientific knowledge and an expanding set of tools for updating outdated precipitation standards based on the assumption of climate stationarity, the adoption of climate informed analyses remain limited in practice. This study utilizes an existing framework to update Intensity (or Depth)-Duration-Frequency (DDF) curves using an ensemble of future climate projections. Change factors in precipitation estimates are derived and applied to six USAF installations across the U.S. The analysis is further extended to evaluate the implications of climate-informed DDFs on stormwater infrastructure performance and flood analysis at Tyndall AFB. Results indicate that the current design precipitation estimates are likely to become obsolete in all six USAF bases by the end of the century. The wide range of change factors across 32 GCM ensembles highlights the need to integrate uncertainty and evolving scientific data into infrastructure planning. The study also finds that the impacts of a changing climate vary spatially and temporally, emphasizing the value of localized analysis for infrastructure decision-making. The work advances ongoing DoD and societal efforts to implement adaptation strategies aimed at enhancing infrastructure resilience.

Intensity-duration-frequency curves↗

Ensemble modeling of global lake evaporation under climate change

Global projections of lake evaporation are typically based on simulations using single mechanistic models. However, because of its complex interactions with various lake physical properties, environmental and anthropogenic drivers, lake evaporation is highly variable and sensitive to the choice of model used. In this study, we present a multi-model analysis to investigate differences across global simulations of lake evaporation during the warm-season using three different lake models driven by outputs from four general circulation models (GCM) (i.e. 12 model combinations in total) for historic and future scenarios. Our results suggest substantial differences among lake-climate model simulations of lake evaporation. These differences varied throughout the 20th and 21st century, with model driver data explaining 74% of the variance in future projections of warm- season lake evaporation. Our projections indicate that, by the end of the 21st century (2070–2099), global annual lake evaporation rates will increase by 10–27% under Representative Concentration Pathways (RCPs) 2.6–8.5. Finally, we highlight the importance of using a multi-model approach for the prediction of future global lake evaporation responses to climate change.

54 ENVIRONMENTAL SCIENCES↗

Water level changes in Lake Erie drive 21st century CO 2 and CH 4 fluxes from a coastal temperate wetland

Wetland water depth influences microbial and plant communities, which can alter the above- and below-ground carbon cycling of a wetland. Wetland water depths are likely to change due to shifting precipitation patterns, which will affect projections of greenhouse gas emissions; however, these effects are rarely incorporated into wetland greenhouse gas models. Seeking to address this gap, we used a mechanistic model, ecosys, to simulate a range of water depth scenarios in a temperate wetland, and analyzed simulated predictions of carbon dioxide (CO 2 ) andmethane (CH 4 ) fluxes over the 21st century. We tested our model using eddy covariance measurements of CO 2 and CH 4 fluxes collected at the Old Woman Creek National Estuarine Research Reserve (OWC) during 2015 and 2016. OWC is a lacustrine, estuarine, freshwater, mineral-soil marsh. An empirical model found that the wetland water depth is highly dependent on the water depth of the nearby Lake Erie. Future wetland surface water depths were modeled based on projection of Lake Erie's water depth using four separate NOAA projections, resulting in four wetland water-depth scenarios. Two of the four 21st century projections for Lake Erie water depths used in this study indicated that the water depth of the wetland would remain nearly steady; however, the other two indicated decreases in the wetland water depth. In our scenario where the wetland dries out, we project the wetland's climatological warming effect will decrease due to smaller CH 4 fluxes to the atmosphere and larger CO 2 uptake by the wetland. We also found that increased water level can lower emissions by shifting the site towards more open water areas, which have lower CH 4 emissions. We found that decreased water depths would cause more widespread colonization of the wetland by macrophyte vegetation. Using an empirical relationship, we also found that further drying could result in other, non-wetland vegetation to emerge, dramatically altering soil carbon cycling. In three of our four projections, we found that in general the magnitude of CO 2 and CH 4 fluxes steadily increase over the next 100 years in response to higher temperatures. However, in our driest simulations, we projected a different response due to increased oxidation of soil carbon, with CH 4 emissions decreasing substantially from an annual cumulative peak of 224.6 to a minimum of 104.7 gC m –2 year –1 . In that same simulation, net cumulative flux of CO 2 changed from being a sink of 56.5 gC m –2 year –1 to a source of 369.6 gCm –2 year –1 over the same period, despite a temperature increase from 13.7 °C to 14.2 °C. This temperature shift in our other three cases with greater water depths increased the source strength of CH 4 and the sink strength of CO 2 . Here, we conclude that the magnitude of wetland greenhouse-gas fluxes depended on the water depth primarily as it affected the areal percentage of the wetland available for plant colonization, but dramatic decreases in water depths could cause significant reductions in the wetland CH 4 fluxes, while simultaneously altering the wetland vegetation.

54 ENVIRONMENTAL SCIENCES↗

Theoretical Insights on the Fragmentation of Cellulosic Oligomers to Form Hydroxyacetone and Hydroxyacetaldehyde

The presence of heavy unknown oligomeric sugar products in bio-oil is evidenced in experimental results reported in the literature. In this paper, we study the fragmentation reactions yielding acetol and glycolaldehyde from oligomeric sugars following previous work on dehydration reactions to propose structures of these oligomers. Acetol and glycolaldehyde are primary products of cellulose fast pyrolysis but the fragmentation reaction mechanism of these compounds from oligomers merits further study. The density functional theory (DFT) approach was employed to study this reaction. Results revealed that acetol and glycolaldehyde fragments are favorably removed from the non-reducing end based on their thermodynamic stabilities. Theoretical FTIR and NMR spectra were calculated to aid in understanding the structures of the oligomeric sugars. Also, the thermodynamics and physical properties of these compounds were estimated using the Group Contribution Method (GCM). In conclusion, these properties are essential in the design of processing technologies and the upgrading of products.

10 SYNTHETIC FUELS↗

A Satellite-Based Estimate of Convective Vertical Velocity and Convective Mass Flux: Global Survey and Comparison with Radar Wind Profiler Observations

Convective vertical velocity (w c ) and convective mass flux (M c ) lie at the heart of GCM cumulus parameterizations, but few observations of these critical parameters are available. In this paper, we develop and evaluate a novel, satellite-based method for estimating profiles of w c and M c . Here, comparisons with collocated ground-based radar wind profiler (RWP) observations show that satellite estimated median w c is slightly greater than the RWP estimates, but they show solid agreement when compared at the 95th percentiles (intense updrafts). RWP-derived and satellite estimated M c are broadly comparable in the lower and middle troposphere, with some differences in the upper troposphere due to differences in convective core sampling. A k-means cluster analysis of multiple years of w c data shows that convective characteristics are distinctly different among extratropical convection, tropical land convection, and tropical oceanic convection. Tropical land convection is significantly more intense and more variable than the oceanic counterpart.

54 ENVIRONMENTAL SCIENCES↗

Understanding Differences in California Climate Projections Produced by Dynamical and Statistical Downscaling

Abstract We compare historical and end‐of‐century temperature and precipitation patterns over California from one dynamically downscaled simulation using the Weather Research and Forecast (WRF) model and two simulations statistically downscaled using Localized Constructed Analogs (LOCA). We uniquely separate causes of differences between dynamically and statistically based future climate projections into differences in historical climate (gridded observations versus regional climate model output) and differences in how these downscaling techniques explicitly handle future climate changes (numerical modeling versus analogs). In these methods, solutions between different downscaling techniques differ more in the future compared to the historical period. Changes projected by LOCA are insensitive to the choice of driving data. Only through dynamical downscaling can we simulate physically consistent regional springtime warming patterns across the Sierra Nevada, while the statistical simulations inherit an unphysical signal from their parent Global Climate Model (GCM) or gridded data. The results of our study clarify why these different techniques produce different outcomes and may also provide guidance on which downscaled products to use for certain impact analyses in California and perhaps other Mediterranean regimes.

54 ENVIRONMENTAL SCIENCES↗

Improving Convection Trigger Functions in Deep Convective Parameterization Schemes Using Machine Learning

Abstract Deficiencies in convection trigger functions, used in deep convection parameterizations in General Circulation Models (GCMs), have critical impacts on climate simulations. A novel convection trigger function is developed using the machine learning (ML) classification model XGBoost. The large‐scale environmental information associated with convective events is obtained from the long‐term constrained variational analysis forcing data from the Atmospheric Radiation Measurement (ARM) program at its Southern Great Plains (SGP) and Manaus (MAO) sites representing, respectively, continental mid‐latitude and tropical convection. The ML trigger is separately trained and evaluated per site, and jointly trained and evaluated at both sites as a unified trigger. The performance of the ML trigger is compared with four convective trigger functions commonly used in GCMs: dilute convective available potential energy (CAPE), undilute CAPE, dilute dynamic CAPE (dCAPE), and undilute dCAPE. The ML trigger substantially outperforms the four CAPE‐based triggers in terms of the F 1 score metric, widely used to estimate the performance of ML methods. The site‐specific ML trigger functions can achieve, respectively, 91% and 93% F 1 scores at SGP and MAO. The unified trigger also has a 91% F 1 score, with virtually no degradation from the site‐specific training, suggesting the potential of a global ML trigger function. The ML trigger alleviates a GCM deficiency regarding the overprediction of convection occurrence, offering a promising improvement to the simulation of the diurnal cycle of precipitation. Furthermore, to overcome the black box issue of the ML methods, insights derived from the ML model are discussed, which may be leveraged to improve traditional CAPE‐based triggers.

54 ENVIRONMENTAL SCIENCES↗

The Influence of Ocean Coupling on Simulated and Projected Tropical Cyclone Precipitation in the HighResMIP–PRIMAVERA Simulations

This study aims to quantify the impacts of atmosphere–ocean coupling on simulated and projected tropical cyclone (TC) precipitation globally. We used global climate model (GCM) simulations over 1950–2050 from the High Resolution Model Intercomparison Project (HighResMIP) and compared its fully coupled atmosphere–ocean GCMs (AOGCMs) with atmosphere-only GCMs (AGCMs). We find that ocean coupling generally leads to decreased TC precipitation over ocean and land. Large-scale sea surface temperature (SST) biases are critical drivers of the precipitation difference, with secondary contributions from local TC–ocean feedbacks via SST cold wakes. The two driving factors, attributed to ocean coupling in the AOGCMs, influence TC precipitation in association with decreased TC intensity and specific humidity. The AOGCMs and AGCMs consistently project TC precipitation increases in 2015–2050 relative to 1950–2014 over ocean for all basins, and for landfalling TCs in the North Atlantic and western North Pacific.

54 ENVIRONMENTAL SCIENCES↗

Factors Governing Cloud Growth and Entrainment Rates in Shallow Cumulus and Cumulus Congestus During GoAmazon2014/5

Shallow cumulus and cumulus congestus clouds play an important role in the large-scale tropical circulation by mixing heat and moisture vertically and preconditioning the environment for deeper convection. Different representations of these shallow clouds account for much of the spread in General Circulation Model (GCMs) climate sensitivity, potentially because of how entrainment is represented in GCM parameterizations. This study uses observations from the Department of Energy's Atmospheric Radiation Measurement (ARM) mobile facility deployed at Manacapuru, Brazil, during the Green Ocean Amazon (GoAmazon2014/5) Campaign. Environmental thermodynamic profiles and observations of cloud top height (CTH) are used to constrain an entraining plume model to estimate bulk entrainment rates (ERs). Estimates of CTH are obtained from a combination of vertically pointing W-band ARM cloud radar and 1,290 MHz Radar Wind Profiler observations. A combination of radiosonde, microwave radiometer profiler, and microwave radiometer observations provides new best estimates of the environmental thermodynamic state. We quantify uncertainty in ERs considering uncertainties in estimated CTH, environmental thermodynamic properties, and assumed initial parcel characteristics. We find ERs ranging from 0.16 to 2.8 km -1 with an average of 0.58 ± 0.10 km -1 over a selected population of 469 shallow cumulus and cumulus congestus clouds. Using the retrieved estimates of ER, we evaluate several entrainment closures that are currently used in atmospheric models or have been proposed based on theory or large eddy simulation. Finally, entrainment rates in cumulus clouds are weakly correlated with low-level buoyancy, cloud depth, and cloud size.

54 ENVIRONMENTAL SCIENCES↗

How May the Choice of Downscaling Techniques and Meteorological Reference Observations Affect Future Hydroclimate Projections?

Abstract We present an intercomparison of a suite of high‐resolution downscaled climate projections based on a six‐member General Circulation Model (GCM) ensemble from Coupled Models Intercomparison Project (CMIP6). The CMIP6 GCMs have been downscaled using dynamical and statistical downscaling techniques based on two meteorological reference observations over the conterminous United States. We use the regional climate model, RegCM4, for dynamical downscaling, double bias correction constructed analogs method for statistical downscaling, and Daymet and Livneh datasets as the reference observations for statistical training and bias‐correction. We evaluate the performances of downscaled data in both historical and future periods under the SSP585 scenario. While dynamical downscaling improves the simulation of some performance evaluation indices, it adds an extra bias in others, highlighting the need for statistical correction before its use in impact assessments. Downscaled datasets after bias‐correction compare exceptionally well with observations. However, the choice of downscaling techniques and the underlying reference observations influence the hydroclimate characteristics of downscaled data. For instance, the statistical downscaling generally preserves the GCMs climate change signal but overestimates the frequency of hot extremes. Similarly, simulated future changes are sensitive to the choice of reference observations, particularly for precipitation extremes that exhibit a higher projected increase in the ensembles trained and/or corrected by Daymet than Livneh. Overall, these results demonstrate that multiple factors, including downscaling techniques and reference observations, can substantially influence the outcome of downscaled climate projections and stress the need for a comprehensive understanding of such method‐based uncertainties.

54 ENVIRONMENTAL SCIENCES↗

Evaluation of CMIP6 GCMs Over the CONUS for Downscaling Studies

Despite the necessity of Global Climate Models (GCMs) sub-selection in downscaling studies, an objective approach for their selection is currently lacking. Building on the previously established concepts in GCMs evaluation frameworks, we develop a weighted averaging technique to remove the redundancy in the evaluation criteria and rank 37 GCMs from the sixth phase of the Coupled Models Intercomparison Project over the contiguous United States. GCMs are rated based on their average performance across 66 evaluation measures in the historical period (1981–2014) after each metric is weighted between zero and one, depending on its uniqueness. The robustness of the outcome is tested by repeating the process with the empirical orthogonal function analysis in which each GCM is ranked based on its sum of distances from the reference in the principal component space. The two methodologies work in contrasting ways to remove the metrics redundancy but eventually develop similar GCMs rankings. A disparity in GCMs' behavior related to their sensitivity to the size of the evaluation suite is observed, highlighting the need for comprehensive multi-variable GCMs evaluation at varying timescales for determining their skillfulness over a region. The sub-selection goal is to use a representative set of skillful models over the region of interest without substantial overlap in their future climate responses and modeling errors in representing historical climate. Additional analyses of GCMs' independence and spread in their future projections provide the necessary information to objectively select GCMs while keeping all aspects of necessity in view.

54 ENVIRONMENTAL SCIENCES↗

On the Effect of Historical SST Patterns on Radiative Feedback

Here, we investigate the dependence of radiative feedback on the pattern of sea-surface temperature (SST) change in fourteen Atmospheric General Circulation Models (AGCMs) forced with observed variations in SST and sea-ice over the historical record from 1871 to near-present. We find that over 1871-1980, the Earth warmed with feedbacks largely consistent and strongly correlated with long term climate sensitivity feedbacks (diagnosed from corresponding atmosphere-ocean GCM abrupt 4xCO2 simulations). Post 1980 however, the Earth warmed with unusual trends in tropical Pacific SSTs (enhanced warming in the west, cooling in the east) and cooling in the Southern Ocean that drove climate feedback to be uncorrelated with – and indicating much lower climate sensitivity than – that expected for long-term CO 2 increase. We show that these conclusions are not strongly dependent on the AMIP II SST dataset used to force the AGCMs, though the magnitude of feedback post 1980 is generally smaller in nine AGCMs forced with alternative HadISST1 SST boundary conditions. We quantify a ‘pattern effect’ (defined as the difference between historical and longterm CO 2 feedback) equal to 0.48 ± 0.47 [5-95%] W m -2 K -1 for the time-period 1871-2010 when the AGCMs are forced with HadISST1 SSTs, or 0.70 ± 0.47 [5-95%] W m -2 K -1 when forced with AMIP II SSTs. Assessed changes in the Earth’s historical energy budget agree with the AGCM feedback estimates. Furthermore satellite observations of changes in top-of-atmosphere radiative fluxes since 1985 suggest that the pattern effect was particularly strong over recent decades but may be waning post 2014.

54 ENVIRONMENTAL SCIENCES↗

Cloud Phase Simulation at High Latitudes in EAMv2: Evaluation using CALIPSO Observations and Comparison with EAMv1

This study performs a comprehensive evaluation of the simulated cloud phase in the U.S. Department of Energy (DOE) Energy Exascale Earth System Model (E3SM) atmosphere model version 2 (EAMv2) and version 1 (EAMv1). Enabled by the CALIPSO (Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation) simulator, EAMv2 and EAMv1 predicted cloud phase is compared against the GCM-Oriented CALIPSO Cloud Product (CALIPSO-GOCCP) at high latitudes where mixed-phase clouds are prevalent. Our results indicate that the underestimation of cloud ice in simulated high-latitude mixed-phase clouds in EAMv1 has been significantly reduced in EAMv2. The increased ice clouds in the Arctic mainly result from the modification on the WBF (Wegner-Bergeron-Findeisen) process in EAMv2. The impact of the modified WBF process is moderately compensated by the low limit of cloud droplet number concentration (CDNC) in cloud microphysics and the new dCAPE_ULL trigger used in deep convection in EAMv2. Moreover, it is found that the new trigger largely contributes to the better cloud phase simulation over the Norwegian Sea and Barents Sea in the Arctic and the Southern Ocean where large errors are found in EAMv1. However, errors in simulated cloud phase in EAMv1, such as the overestimation of supercooled liquid clouds near the surface in both hemispheres and the underestimation of ice clouds over Antarctica, persist in EAMv2. In conclusion, this study highlights the impact of deep convection parameterizations, which has received little attention, on high-latitude mixed-phase clouds, and the importance of continuous improvement of cloud microphysics in climate models for accurately representing mixed-phase clouds.

54 ENVIRONMENTAL SCIENCES↗

Important ice processes are missed by the Community Earth System Model (CESM) in Southern Ocean mixed–phase clouds: Bridging SOCRATES observations to model developments

Global climate models (GCMs) are challenged by difficulties in simulating cloud phase and cloud radiative effect over the Southern Ocean (SO). Some of the new generation GCMs predict too much liquid and too little ice in mixed-phase clouds. This misrepresentation of cloud phase in GCMs results in weaker negative cloud feedback over the SO and a higher climate sensitivity. Based on a model comparison with observational data obtained during the Southern Ocean Cloud Radiation and Aerosol Transport Experimental Study (SOCRATES), this study addresses a key uncertainty in the Community Earth System Model version 2 (CESM2) related to cloud phase, namely ice formation in pristine remote SO clouds. It is found that sea spray organic aerosols (SSOA) are the most important type of ice nucleating particles (INPs) over the SO with concentrations one order of magnitude higher than those of dust INPs based on measurements and CESM2 simulations. Secondary ice production (SIP) which includes riming splintering, rain droplet shattering, and ice-ice collisional fragmentation as implemented in CESM2 is the dominant ice production process in moderately cold clouds with cloud temperatures greater than –20°C. SIP enhances the in-cloud ice number concentrations (Ni) by 1-3 orders of magnitude and predicts more mixed-phase (with percentage occurrence increased from 15% to 21%), in better agreement with the observations. Finally, this study highlights the importance of accurately representing cloud phase over the pristine remote SO by considering the ice nucleation of SSOA and SIP processes, which are currently missing in most GCM cloud microphysics parameterizations.

54 ENVIRONMENTAL SCIENCES↗

Large Ensemble Particle Filter for Spatial Climate Reconstructions Using a Linear Inverse Model

Proxy records that document the last 2000 years of climate provide evidence for the wide range of the natural climate variability from inter-annual to secular timescales not captured by the short window of recent direct observations. Assessing climate models ability to reproduce such natural variations is crucial to understand climate sensitivity and impacts of future climate change. Paleoclimate data assimilation (PDA) offers a powerful way to extend the short instrumental period by optimally combining the physics described by General Circulation Climate Models (GCMs) with information from available proxy records while taking into account their uncertainties. Here we present a new PDA approach based on a sequential importance resampling (SIR) Particle filter (PF) that uses Linear Inverse Modeling (LIM) as an emulator of several CMIP-class GCMs. We examine in a perfect-model framework the skill of the various LIMs to forecast the dynamics of the surface temperatures and provide spatial field reconstructions over the last millennium in a SIR PF. Our results show that the LIMs allow for skillful ensemble forecasts at 1-year lead-time based on GCMs dynamical knowledge with best prediction in the tropics and the North Atlantic. The PDA further provides a set of physically consistent spatial fields allowing robust uncertainty quantification related to climate models biases and proxy spatial sampling. Our results indicate that the LIM yields dynamical memory improving climate variability reconstructions and support the use of the LIM as a GCM-emulator in real reconstruction to propagate large ensembles of particles at low cost in SIR PF.

54 ENVIRONMENTAL SCIENCES↗