Community Climate System Model Version Research Articles

Transient simulation of orbital‐scale precipitation variation in monsoonal East Asia and arid central Asia during the last 150 ka

AbstractA long‐term transient simulation is conducted using the Community Climate System Model version 3 and the orbital acceleration technique to analyze the impact of insolation change caused by the Earth's orbital forcing on precipitation in the monsoonal East Asia (EA) and arid central Asia (CA) over the past 150 ka. Our results show that annual precipitation in both EA and CA has strong signals of the 20 ka precessional cycles and varies in phase with the Northern Hemisphere (NH) summer insolation. Similar characteristics can also be observed from previously published oxygen isotope records of stalagmites near EA and CA. Composite analyses based on seven precessional cycles suggest that the increase (decrease) in the NH summer (winter) insolation enhances EA (CA) summer (winter) precipitation by modulation of the Asian monsoon (westerly) circulation in summer (winter). When the precession‐induced NH summer insolation increases, the Asian summer monsoon circulation is enhanced and EA precipitation increases significantly. Meanwhile, the increase in the summer insolation at the precessional scale is accompanied by a decrease in the winter insolation, which causes dramatic cooling of the troposphere in the lower latitudes. Consequently, the CA winter precipitation increases due to the changes in the temperature gradient and the westerly circulation. Therefore, the responses of the Asian monsoon and westerly circulation to summer and winter insolation variations induced by the precessional cycles determine precipitation in the respective rainy seasons and are the primary cause leading to the synchronous variation patterns of annual precipitation in EA and CA at the orbital scale.

The Mean Climate of the Community Atmosphere Model (CAM4) in Forced SST and Fully Coupled Experiments

Abstract The Community Atmosphere Model, version 4 (CAM4), was released as part of the Community Climate System Model, version 4 (CCSM4). The finite volume (FV) dynamical core is now the default because of its superior transport and conservation properties. Deep convection parameterization changes include a dilute plume calculation of convective available potential energy (CAPE) and the introduction of convective momentum transport (CMT). An additional cloud fraction calculation is now performed following macrophysical state updates to provide improved thermodynamic consistency. A freeze-drying modification is further made to the cloud fraction calculation in very dry environments (e.g., the Arctic), where cloud fraction and cloud water values were often inconsistent in CAM3. In CAM4 the FV dynamical core further degrades the excessive trade-wind simulation, but reduces zonal stress errors at higher latitudes. Plume dilution alleviates much of the midtropospheric tropical dry biases and reduces the persistent monsoon precipitation biases over the Arabian Peninsula and the southern Indian Ocean. CMT reduces much of the excessive trade-wind biases in eastern ocean basins. CAM4 shows a global reduction in cloud fraction compared to CAM3, primarily as a result of the freeze-drying and improved cloud fraction equilibrium modifications. Regional climate feature improvements include the propagation of stationary waves from the Pacific into midlatitudes and the seasonal frequency of Northern Hemisphere blocking events. A 1° versus 2° horizontal resolution of the FV dynamical core exhibits superior improvements in regional climate features of precipitation and surface stress. Improvements in the fully coupled mean climate between CAM3 and CAM4 are also more substantial than in forced sea surface temperature (SST) simulations.

The MJO and global warming: a study in CCSM4

The change in Madden–Julian oscillation (MJO) amplitude and variance in response to anthropogenic climate change is assessed in the 1° nominal resolution community climate system model, version 4 (CCSM4), which has a reasonable representation of the MJO characteristics both dynamically and statistically. The twentieth century CCSM4 run is compared with the warmest twenty-first century projection (representative concentration pathway 8.5, or RCP8.5). The last 20 years of each simulation are compared in their MJO characteristics, including spatial variance distributions of winds, precipitation and outgoing longwave radiation, histograms of event amplitude, phase and duration, and composite maps of phases. The RCP8.5 run exhibits increased variance in intraseasonal precipitation, larger-amplitude MJO events, stronger MJO rainfall in the central and eastern tropical Pacific, and a greater frequency of MJO occurrence for phases corresponding to enhanced rainfall in the Indian Ocean sector. These features are consistent with the concept of an increased magnitude for the hydrological cycle under greenhouse warming conditions. Conversely, the number of active MJO days decreases and fewer weak MJO events occur in the future climate state. These results motivate further study of these changes since tropical rainfall variability plays such an important role in the region’s socio-economic well being.

Western U.S. Extreme Precipitation Events and Their Relation to ENSO and PDO in CCSM4

Abstract Water resources and management over the western United States are heavily impacted by both local climate variability and the teleconnected responses of precipitation to the El Niño–Southern Oscillation (ENSO) and Pacific decadal oscillation (PDO). In this work, regional precipitation patterns over the western United States and linkages to ENSO and the PDO are analyzed using output from a Community Climate System Model version 4 (CCSM4) preindustrial control run and observations, with emphasis on extreme precipitation events. CCSM4 produces realistic zonal gradients in precipitation intensity and duration over the western United States, with higher values on the windward side of the Cascade Mountains and Sierra Nevada and lower values on the leeward. Compared to its predecessor CCSM3, CCSM4 shows an improved teleconnected signal of both ENSO and the PDO to large-scale circulation patterns over the Pacific–North America region and also to the spatial pattern and other aspects of western U.S. precipitation. The so-called drizzle problem persists in CCSM4 but is significantly improved compared to CCSM3. In particular, it is found that CCSM4 has substantially less precipitation duration bias than is present in CCSM3. Both the overall and extreme intensity of wintertime precipitation over the western United States show statistically significant linkages with ENSO and PDO in CCSM4. This analysis provides a basis for future studies using greenhouse gas (GHG)-forced CCSM4 runs.

The Atlantic Meridional Heat Transport at 26.5°N and Its Relationship with the MOC in the RAPID Array and the GFDL and NCAR Coupled Models

Abstract The link at 26.5°N between the Atlantic meridional heat transport (MHT) and the Atlantic meridional overturning circulation (MOC) is investigated in two climate models, the GFDL Climate Model version 2.1 (CM2.1) and the NCAR Community Climate System Model version 4 (CCSM4), and compared with the recent observational estimates from the Rapid Climate Change–Meridional Overturning Circulation and Heatflux Array (RAPID–MOCHA) array. Despite a stronger-than-observed MOC magnitude, both models underestimate the mean MHT at 26.5°N because of an overly diffuse thermocline. Biases result from errors in both overturning and gyre components of the MHT. The observed linear relationship between MHT and MOC at 26.5°N is realistically simulated by the two models and is mainly due to the overturning component of the MHT. Fluctuations in overturning MHT are dominated by Ekman transport variability in CM2.1 and CCSM4, whereas baroclinic geostrophic transport variability plays a larger role in RAPID. CCSM4, which has a parameterization of Nordic Sea overflows and thus a more realistic North Atlantic Deep Water (NADW) penetration, shows smaller biases in the overturning heat transport than CM2.1 owing to deeper NADW at colder temperatures. The horizontal gyre heat transport and its sensitivity to the MOC are poorly represented in both models. The wind-driven gyre heat transport is northward in observations at 26.5°N, whereas it is weakly southward in both models, reducing the total MHT. This study emphasizes model biases that are responsible for the too-weak MHT, particularly at the western boundary. The use of direct MHT observations through RAPID allows for identification of the source of the too-weak MHT in the two models, a bias shared by a number of Coupled Model Intercomparison Project phase 5 (CMIP5) coupled models.

Uncertainty in future regional sea level rise due to internal climate variability

AbstractSea level rise (SLR) is an inescapable consequence of increasing greenhouse gas concentrations, with potentially harmful effects on human populations in coastal and island regions. Observational evidence indicates that global sea level has risen in the 20th century, and climate models project an acceleration of this trend in the coming decades. Here we analyze rates of future SLR on regional scales in a 40‐member ensemble of climate change projections with the Community Climate System Model Version 3. This unique ensemble allows us to assess uncertainty in the magnitude of 21st century SLR due to internal climate variability alone. We find that simulated regional SLR at mid‐century can vary by a factor of 2 depending on location, with the North Atlantic and Pacific showing the greatest range. This uncertainty in regional SLR results primarily from internal variations in the wind‐driven and buoyancy‐driven ocean circulations.

The Norwegian Earth System Model, NorESM1-M – Part 1: Description and basic evaluation of the physical climate

Abstract. The core version of the Norwegian Climate Center's Earth System Model, named NorESM1-M, is presented. The NorESM family of models are based on the Community Climate System Model version 4 (CCSM4) of the University Corporation for Atmospheric Research, but differs from the latter by, in particular, an isopycnic coordinate ocean model and advanced chemistry–aerosol–cloud–radiation interaction schemes. NorESM1-M has a horizontal resolution of approximately 2° for the atmosphere and land components and 1° for the ocean and ice components. NorESM is also available in a lower resolution version (NorESM1-L) and a version that includes prognostic biogeochemical cycling (NorESM1-ME). The latter two model configurations are not part of this paper. Here, a first-order assessment of the model stability, the mean model state and the internal variability based on the model experiments made available to CMIP5 are presented. Further analysis of the model performance is provided in an accompanying paper (Iversen et al., 2013), presenting the corresponding climate response and scenario projections made with NorESM1-M.

On the twenty-first-century wet season projections over the Southeastern United States

This paper reconciles the difference in the projections of the wet season over the Southeastern United States (SEUS) from a global climate model (the Community Climate System Model Version 3 [CCSM3]) and from a regional climate model (the Regional Spectral Model [RSM]) nested in the CCSM3. The CCSM3 projects a dipole in the summer precipitation anomaly: peninsular Florida dries in the future climate, and the remainder of the SEUS region becomes wetter. The RSM forced with CCSM3 projects a universal drying of the SEUS in the late twenty-first century relative to the corresponding twentieth-century summer. The CCSM3 pattern is attributed to the “upped-ante” mechanism, whereby the atmospheric boundary layer moisture required for convection increases in a warm, statically stable global tropical environment. This criterion becomes harder to meet along convective margins, which include peninsular Florida, resulting in its drying. CCSM3 also projects a southwestward expansion of the North Atlantic subtropical high that leads to further stabilizing of the atmosphere above Florida, inhibiting convection. The RSM, because of its high (10-km grid) resolution, simulates diurnal variations in summer rainfall over SEUS reasonably well. The RSM improves upon CCSM3 through the RSM’s depiction of the diurnal variance of precipitation, which according to observations accounts for up to 40 % of total seasonal precipitation variance. In the future climate, the RSM projects a significant reduction in the diurnal variability of convection. The reduction is attributed to large-scale stabilization of the atmosphere in the CCSM3 projections.

Simulating changes in the leaf unfolding time of 20 plant species in China over the twenty-first century

Recent shifts in phenology reflect the biological response to current climate change. Aiming to enhance our understanding of phenological responses to climate change, we developed, calibrated and validated spatio-temporal models of first leaf date (FLD) for 20 broadleaved deciduous plants in China. Using daily meteorological data from the Chinese Meteorological Administration and the Community Climate System Model, version 3 (CCSM3) created using three IPCC scenarios (A2, A1B and B1), we described the FLD time series of each species over the past 50 years, extrapolating from these results to simulate estimated FLD changes for each species during the twenty-first century. Model validation suggests that our spatio-temporal models can simulate FLD accurately with R² (explained variance) >0.60. Model simulations show that, from 1952 to 2007, the FLD in China advanced at a rate of -1.14 days decade⁻¹) on average. Furthermore, changes in FLD showed noticeable variation between regions, with clearer advances observed in the north than in the south of the country. The model indicates that the advances in FLD observed from 1952-2007 in China will continue over the twenty-first century, although significant differences among species and different climate scenarios are expected. The average trend of FLD advance in China during the twenty-first century is modeled as being -1.92 days decade⁻¹ under the A2 scenario, -1.10 days decade⁻¹ under the A1B scenario and -0.74 days decade⁻¹ under the B2 scenario. The spatial pattern of FLD change for the period 2011-2099 is modeled as being similar but showing some difference from patterns in the 1952-2007 period. At the interspecific level, early-leafing species were found to show a greater advance in FLD, while species with larger distributions tended to show a weaker advance in FLD. These simulated changes in phenology may have significant implications for plant distribution as well as ecosystem structure and function.

Influence of Continental Ice Retreat on Future Global Climate

Abstract Evidence from observations indicates a net loss of global land-based ice and a rise of global sea level. Other than sea level rise, it is not clear how this loss of land-based ice could affect other aspects of global climate in the future. Here, the authors use the Community Climate System Model version 3 (CCSM3) to evaluate the potential influence of shrinking land-based ice on the Atlantic meridional overturning circulation (AMOC) and surface climate in the next two centuries under the Intergovernmental Panel on Climate Change (IPCC) A1B scenario with prescribed rates of melting for the Greenland Ice Sheet, West Antarctic Ice Sheet, and mountain glaciers and ice caps. Results show that the AMOC, in general, is only sensitive to the freshwater discharge directly into the North Atlantic over the next two centuries. If the loss of the West Antarctic Ice Sheet would not significantly increase from its current rate, it would not have much effect on the AMOC. The AMOC slows down further only when the surface freshwater input due to runoff from land-based ice melt becomes large enough to generate a net freshwater gain in the upper North Atlantic. This further-weakened AMOC does not cool the global mean climate, but it does cause less warming, especially in the northern high latitudes and, in particular, in Europe. The projected precipitation increase in North America in the standard run becomes a net reduction in the simulation that includes land ice runoff, but there are precipitation increases in west Australia in the simulations where the AMOC slows down because of the inclusion of land-based ice runoff.

The Impact of Oceanic Near-Inertial Waves on Climate

Abstract The Community Climate System Model, version 4 (CCSM4) is used to assess the climate impact of wind-generated near-inertial waves (NIWs). Even with high-frequency coupling, CCSM4 underestimates the strength of NIWs, so that a parameterization for NIWs is developed and included into CCSM4. Numerous assumptions enter this parameterization, the core of which is that the NIW velocity signal is detected during the model integration, and amplified in the shear computation of the ocean surface boundary layer module. It is found that NIWs deepen the ocean mixed layer by up to 30%, but they contribute little to the ventilation and mixing of the ocean below the thermocline. However, the deepening of the tropical mixed layer by NIWs leads to a change in tropical sea surface temperature and precipitation. Atmospheric teleconnections then change the global sea level pressure fields so that the midlatitude westerlies become weaker. Unfortunately, the magnitude of the real air-sea flux of NIW energy is poorly constrained by observations; this makes the quantitative assessment of their climate impact rather uncertain. Thus, a major result of the present study is that because of its importance for global climate the uncertainty in the observed tropical NIW energy has to be reduced.

Transient Climate Response in Coupled Atmospheric–Ocean General Circulation Models

Abstract The equilibrium climate sensitivity (ECS) has a large uncertainty range among models participating in the Intergovernmental Panel on Climate Change (IPCC) Fourth Assessment Report (AR4) and has recently been presented as “inherently unpredictable.” One way to circumvent this problem is to consider the transient climate response (TCR). However, the TCR among AR4 models also differs by more than a factor of 2. The authors argue that the situation may not necessarily be so pessimistic, because much of the intermodel difference may be due to the fact that the models were run with their oceans at various stages of flux adjustment with their atmosphere. This is shown by comparing multimillennium-long runs of the Goddard Institute for Space Studies model, version E, coupled with the Hybrid Coordinate Ocean Model (GISS-EH) and the Community Climate System Model, version 4 (CCSM4) with what were reported to AR4. The long model runs here reveal the range of variability (~30%) in their TCR within the same model with the same ECS. The commonly adopted remedy of subtracting the “climate drift” is ineffective and adds to the variability. The culprit is the natural variability of the control runs, which exists even at quasi equilibration. Fortunately, for simulations with multidecadal time horizon, robust solutions can be obtained by branching off thousand-year-long control runs that reach “quasi equilibration” using a new protocol, which takes advantage of the fact that forced solutions to radiative forcing forget their initial condition after 30–40 yr and instead depend mostly on the trajectory of the radiative forcing.

Sensitivity to Glacial Forcing in the CCSM4

AbstractResults are presented from the Community Climate System Model, version 4 (CCSM4), simulation of the Last Glacial Maximum (LGM) from phase 5 of the Coupled Model Intercomparison Project (CMIP5) at the standard 1° resolution, the same resolution as the majority of the CCSM4 CMIP5 long-term simulations for the historical and future projection scenarios. The forcings and boundary conditions for this simulation follow the protocols of the Paleoclimate Modeling Intercomparison Project, version 3 (PMIP3). Two additional CCSM4 CO2 sensitivity simulations, in which the concentrations are abruptly changed at the start of the simulation to the low 185 ppm LGM concentrations (LGMCO2) and to a quadrupling of the preindustrial concentration (4×CO2), are also analyzed. For the full LGM simulation, the estimated equilibrium cooling of the global mean annual surface temperature is 5.5°C with an estimated radiative forcing of −6.2 W m−2. The radiative forcing includes the effects of the reduced LGM greenhouse gases, ice sheets, continental distribution with sea level lowered by approximately 120 m from the present, and orbital parameters, but not changes to atmospheric aerosols or vegetation biogeography. The LGM simulation has an equilibrium climate sensitivity (ECS) of 3.1(±0.3)°C, comparable to the CCSM4 4×CO2 result. The LGMCO2 simulation shows a greater ECS of 4.2°C. Other responses found at the LGM in CCSM4 include a global precipitation rate decrease at a rate of ~2% °C−1, similar to climate change simulations in the Intergovernmental Panel on Climate Change Fourth Assessment Report (IPCC AR4); a strengthening of the Atlantic meridional overturning circulation (AMOC) with a shoaling of North Atlantic Deep Water and a filling of the deep basin up to sill depth with Antarctic Bottom Water; and an enhanced seasonal cycle accompanied by reduced ENSO variability in the eastern Pacific Ocean’s SSTs.

A Diagnostic Indicator of the Stability of the Atlantic Meridional Overturning Circulation in CCSM3

AbstractA diagnostic indicator ΔMov is proposed in this paper to monitor the stability of the Atlantic meridional overturning circulation (AMOC). The ΔMov is a diagnostic for a basinwide salt-advection feedback and defined as the difference between the freshwater transport induced by the AMOC across the southern border of the Atlantic Ocean and the overturning liquid freshwater transport from the Arctic Ocean to the North Atlantic. As validated in the Community Climate System Model, version 3 (CCSM3), for an AMOC in the conveyor state, a positive ΔMov (freshwater convergence) in the Atlantic basin indicates a monostable AMOC and a negative ΔMov (freshwater divergence) indicates a bistable AMOC. Based on ΔMov, the authors investigate the AMOC stability in the Last Glacial Maximum (LGM) and analyze the modulation of the AMOC stability by an open/closed Bering Strait. Moreover, the authors estimate that the real AMOC is likely to be bistable in the present day, since some observations suggest a negative ΔMov (freshwater divergence) is currently in the Atlantic basin. However, this estimation is very sensitive to the choice of the observational data.

Model support for forcing of the 8.2 ka event by meltwater from the Hudson Bay ice dome

Previous model experiments of the 8.2 ka event forced by the drainage of Lake Agassiz often do not produce climate anomalies as long as those inferred from proxies. In addition to the Agassiz forcing, there is new evidence for significant amounts of freshwater entering the ocean at 8.2 ka from the disintegration of the Laurentide ice sheet (LIS). We use the Community Climate System Model version 3 (CCSM3) to test the contribution of this additional meltwater flux. Similar to previous model experiments, we find that the estimated freshwater forcing from Lake Agassiz is capable of sustaining ocean and climate anomalies for only two to three decades, much shorter than the event duration of ~150 years in proxies. Using new estimates of the LIS freshwater flux (~0.13 Sv for 100 years) from the collapse of the Hudson Bay ice dome in addition to the Agassiz drainage, the CCSM3 generates climate anomalies with a magnitude and duration that match within error those from proxies. This result is insensitive to the duration of freshwater release, a major uncertainty, if the total volume remains the same. An analysis of the modeled North Atlantic freshwater budget indicates that the Agassiz drainage is rapidly transported out of the North Atlantic while the LIS contribution generates longer-lasting freshwater anomalies that are also subject to recirculation by the subtropical gyre back into the North Atlantic. Thus, the meltwater flux originating from the LIS appears to be more important than the Agassiz drainage in generating 8.2 ka climate anomalies and is one way to reconcile some model-data discrepancies.

Last Millennium Climate and Its Variability in CCSM4

Abstract An overview of a simulation referred to as the “Last Millennium” (LM) simulation of the Community Climate System Model, version 4 (CCSM4), is presented. The CCSM4 LM simulation reproduces many large-scale climate patterns suggested by historical and proxy-data records, with Northern Hemisphere (NH) and Southern Hemisphere (SH) surface temperatures cooling to the early 1800s Common Era by ~0.5°C (NH) and ~0.3°C (SH), followed by warming to the present. High latitudes of both hemispheres show polar amplification of the cooling from the Medieval Climate Anomaly (MCA) to the Little Ice Age (LIA) associated with sea ice increases. The LM simulation does not reproduce La Niña–like cooling in the eastern Pacific Ocean during the MCA relative to the LIA, as has been suggested by proxy reconstructions. Still, dry medieval conditions over the southwestern and central United States are simulated in agreement with proxy indicators for these regions. Strong global cooling is associated with large volcanic eruptions, with indications of multidecadal colder climate in response to larger eruptions. The CCSM4’s response to large volcanic eruptions captures some reconstructed patterns of temperature changes over Europe and North America, but not those of precipitation in the Asian monsoon region. The Atlantic multidecadal oscillation (AMO) has higher variance at centennial periods in the LM simulation compared to the 1850 nontransient run, suggesting a long-term Atlantic Ocean response to natural forcings. The North Atlantic Oscillation (NAO), Pacific decadal oscillation (PDO), and El Niño–Southern Oscillation (ENSO) variability modes show little or no change. CCSM4 does not simulate a persistent positive NAO or a prolonged period of negative PDO during the MCA, as suggested by some proxy reconstructions.

The Sensitivity of the Arctic Ocean Sea Ice Thickness and Its Dependence on the Surface Albedo Parameterization

Abstract In this study, the response of sea ice thickness to changes in the external forcing is investigated and particularly how this response depends on the surface albedo formulation by means of a one-dimensional coupled ocean–ice–atmosphere model. The main focus is on the thickness response to the atmospheric heat advection Fwall, solar radiation FSW, and amount of snow precipitation Sprec. Different albedo parameterization schemes [ECHAM5, CSIRO, and Community Climate System Model, version 3 (CCSM3)] representing albedos commonly used in global climate models are compared together with more simplified schemes. Using different albedo schemes with the same external forcing produces large differences in ice thickness. The ice thickness response is similar for all realistic albedo schemes with a nearly linear decrease with increasing Fwall in the perennial ice regime and with a steplike transition into seasonal ice when Fwall exceeds a certain threshold. This transition occurs at an annual-mean ice thickness of 1.7–2.0 m. Latitudinal differences in solar insolation generally leads to increasing ice thickness toward the North Pole. The snow response varies significantly depending on which albedo scheme is used. The ECHAM5 scheme yields thinner ice with Sprec, the CSIRO scheme gives ice thickness nearly independent of Sprec, and with the CCSM3 scheme the ice thickness decreases with Sprec. A general result is that the modeled ice cover is rather sensitive to positive perturbations of the external heat supply when it is close to the transition such that just a small increase of, for example, Fwall can force the ice cover into the seasonal regime.

Pacific Decadal Variability: Paced by Rossby Basin Modes?

Abstract A systematic study is presented of decadal climate variability in the North Pacific. In particular, the hypothesis is addressed that oceanic Rossby basin modes are responsible for enhanced energy at decadal and bidecadal time scales. To this end, a series of statistical analyses are performed on a 500-yr control integration of the Community Climate System Model, version 3 (CCSM3). In particular, a principal oscillation pattern (POP) analysis is performed to identify modal behavior in the subsurface pressure field. It is found that the dominant energy of sea surface temperature (SST) variability at 25 yr (the model equivalent of the Pacific decadal oscillation) cannot be explained by the resonant excitation of an oceanic basin mode. However, significant energy in the subsurface pressure field at time scales of 17 and 10 yr appears to be related to internal ocean oscillations. However, these oscillations lack the characteristics of the classical basin modes, and must either be deformed beyond recognition by the background circulation and inhomogeneous stratification or have another dynamical origin altogether. The 17-yr oscillation projects onto the Pacific decadal oscillation and, if present in the real ocean, has the potential to enhance the predictability of low-frequency climate variability in the North Pacific.

Short-term climatic impacts of afforestation in the East Asian monsoon region

The influence of changes in vegetation cover on short-term climate over the East Asian monsoon region is simulated using the Community Climate System Model Version 3.5. The results show the annual mean surface air temperature significantly decreases by 0.93°C in response to afforestation over the East Asian monsoon region. Also, surface air temperature decreases by 1.46 and 0.40°C in summer and winter, respectively. The cooling is caused by enhanced evapotranspiration (ET) produced by increased forest cover. Evapotranspiration is greater in summer than in winter, so summer cooling is greater than winter cooling. The annual mean precipitation increases in response to afforestation, with a maximum of 7% in April. Water vapor increases significantly because of greater latent heat flux release. Meanwhile, afforestation leads to higher surface roughness, which decreases surface wind speed and induces an ascending air motion. These factors can produce more clouds and precipitation. Moreover, the surface albedo and the reflective solar radiation are reduced in response to afforestation.

A Pseudoproxy Evaluation of Bayesian Hierarchical Modeling and Canonical Correlation Analysis for Climate Field Reconstructions over Europe*

Abstract A pseudoproxy comparison is presented for two statistical methods used to derive annual climate field reconstructions (CFRs) for Europe. The employed methods use the canonical correlation analysis (CCA) procedure presented by Smerdon et al. and the Bayesian hierarchical model (BHM) method adopted from Tingley and Huybers. Pseudoproxy experiments (PPEs) are constructed from modeled temperature data sampled from the 1250-yr paleo-run of the NCAR Community Climate System Model (CCSM) version 1.4 model by Ammann et al. Pseudoproxies approximate the distribution of the multiproxy network used by Mann et al. over the European region of interest. Gaussian white noise is added to the temperature data to mimic the combined signal and noise properties of real-world proxies. Results indicate that, while both methods perform well in areas with good proxy coverage, the BHM method outperforms the CCA method across the entire field and additionally returns objective error estimates.