Investigation of Climate Change Impacts on Water Resources in the California region
Investigators:
D. Cayan1, M. Dettinger2, R. Hanson 2, T. Brown 3, A. Westerling1
1 Scripps Institution of Oceanography
2 US Geological Survey
3 Desert Research Institute
Last update: 5 November 2001

Click here for ACPI Progress report June 2001

Click here for Comparison of ACPI B06 Simulations February 2001

Click here for ACPI Progress report January 2001

Some useful ACPI links:

The ACPI site: http://www.csm.ornl.gov/ACPI/

The Parallel Climate Model (PCM) web site: http://www.cgd.ucar.edu/pcm/

**********************************************************************
From the 6-hourly average files, the variables available are:
  • LHFLX - surface latent heat flux - W/m^2
  • PRECC - convective precipitation - m/s
  • PRECL - large-scale stable precipitation - m/s
  • QDRAI - LSM sub-surface runoff - mm/s
  • QOVER - LSM surface runoff - mm/s
  • RSW - LSM root zone volumetric soil water - fraction
  • SHFLX - surface sensible heat flux - W/m^2
  • SNOWH - water equivalent snow depth - m
  • SRFRAD - radiative flux absorbed at the surface - W/m^2
  • TAH - LSM "surface" temperature, height - K
  • TGRD - LSM ground "skin" temperature - K
  • TSOI3 - LSM third soil layer temperature - K
  • TSOI4 - LSM fourth soil layer temperature - K
  • TVEG - LSM vegetation "skin" temperature - K
    LSM = Land Surface Model
  • **********************************************************************

    21 June 2001

    Some Recent Precipitation Analyses

    Time series plots of total precipitation for three ACPI simulations are shown below. The three simulations are the control simulation (B06.45, 1995-2048) and two climate change simulations (B06.44 and B06.46, both 1995-2099). Four areas were plotted for land and ocean, for land only and for ocean only: global, 20N-90N, 20S-20N and 90S-20S.
    LAND+OCEAN
    LAND ONLY
    OCEAN ONLY
    Global maps of total precipitation trends for winter (DJF) and summer (JJA) are shown below. These maps show the total seasonal precipitation difference (cm over three months) between the first 10 years and last 10 years of one of the ACPI climate change simualtions (B06.46).
    DJF
    JJA

    **********************************************************************

    Last update: 16 August 2001

    Information about 6-hourly instantaneous data

    Extraction of 6-hourly instataneous data by Scripps has begun. The data will be saved over a section of North America from 72N-20N and from 180 (dateline) to 295E (65W). Five 3-D variables will be saved on all 19 sigma levels: specific humidity, temperature, zonal (u) wind, meridional (v) wind and geopotential height. Four 2-D variables will be saved also: surface pressure, surface temperature, sea level pressure and reference height (2m) air temperature. (Note with the surface pressure, the 3-D variables can be extrapolated to standard pressure levels such as 850mb, 500mb, etc.)

    The 19 latitudes (all North) are:

    71.15775
    68.36775
    65.57761
    62.78735
    59.99702
    57.20663
    54.41620
    51.62573
    48.83524
    46.04472
    43.25419
    40.46365
    37.67309
    34.88252
    32.09194
    29.30136
    26.51077
    23.72017
    20.92957

    The 42 longitudes (all East) are:

    180.000
    182.8125
    185.625
    188.4375
    191.250
    194.0625
    196.875
    199.6875
    202.500
    205.3125
    208.125
    210.9375
    213.750
    216.5625
    219.375
    222.1875
    225.000
    227.8125
    230.625
    233.4375
    236.250
    239.0625
    241.875
    244.6875
    247.500
    250.3125
    253.125
    255.9375
    258.750
    261.5625
    264.375
    267.1875
    270.000
    272.8125
    275.625
    278.4375
    281.250
    284.0625
    286.875
    289.6875
    292.500
    295.3125

    A sample plot of the 2m air temperature (reference height temperature) for 12Z on the first of January, 2030, is shown below. This plot shows the full region for which the 6-hourly instantaneous data is being retrieved. Note the contour interval is 5 Kelvin.


    The following scatter plot shows a sample from the 6-hourly instantaneous data. The lowest model level specific humidity (y-axis, in g/kg) is plotted with the reference height (2m) air temperature for a single PCM grid point: 34.88N 244.7W. The scatter plots are for each of the four time snapshots: 00Z, 06Z, 12Z and 18Z. September and October are shown with a black '+' for 2020 and a red '*' for 2060.


    Lowest Model Level Specific Humidity versus Reference Height Temperature 35N 115W

    **********************************************************************

    June 11, 2001: The Los Angeles Times

    **********************************************************************

    27 September 2001

    Data files for Ruby

    B06.22 2m Air Temperature in Kelvin Mar 1870 thru Nov 1999

    B06.22 Total Precipitation in cm (per month) Mar 1870 thru Nov 1999

    B06.22 Surface Zonal Wind in m/s Mar 1870 thru Nov 1999

    B06.22 Surface Meridional Wind in m/s Mar 1870 thru Nov 1999

    **********************************************************************

    15 May 2001

    Data files for Andrew Wood at UW

    Ascii data on Andy's grid:

    The pcm data is on the T42 grid so the data cells look like:

    233.44 236.25 239.06 241.88 244.69 247.50 250.31 253.13
    54.41620 Land Land Land Land Land Land Land Land 52
    51.62573 Land Land Land Land Land Land Land Land 51
    48.83524 Land Land Land Land Land Land Land Land 50
    46.04472 Sea Land Land Land Land Land Land Land 49
    43.25419 Sea Land Land Land Land Land Land Land 48
    40.46365 Sea Land Land Land Land Land Land Land 47
    37.67309 Sea Land Land Land Land Land Land Land 46
    34.88252 Sea Sea Land Land Land Land Land Land 45
    32.09194 Sea Sea Sea Land Land Land Land Land 45
    84 85 86 87 88 89 90 91

    Four variables are available here for 2 simulations (B06.44 and B06.45):

    B06.44 2m Air Temperature in Kelvin Feb 1995 thru Nov 2099

    B06.45 2m Air Temperature in Kelvin Feb 1995 thru Dec 2048

    B06.44 Total Precipitation in cm (per month) Feb 1995 thru Nov 2099

    B06.45 Total Precipitation in cm (per month) Feb 1995 thru Dec 2048

    B06.44 Surface Zonal Wind in m/s Feb 1995 thru Nov 2099

    B06.45 Surface Zonal Wind in m/s Feb 1995 thru Dec 2048

    B06.44 Surface Meridional Wind in m/s Feb 1995 thru Nov 2099

    B06.45 Surface Meridional Wind in m/s Feb 1995 thru Dec 2048

    Added data for 2 additional simulations (B06.46 and B06.47) on 18 July 2001:

    B06.46 2m Air Temperature in Kelvin Feb 1995 thru Nov 2099

    B06.47 2m Air Temperature in Kelvin Feb 1995 thru Nov 2099

    B06.46 Total Precipitation in cm (per month) Feb 1995 thru Nov 2099

    B06.47 Total Precipitation in cm (per month) Feb 1995 thru Nov 2099

    B06.46 Surface Zonal Wind in m/s Feb 1995 thru Nov 2099

    B06.47 Surface Zonal Wind in m/s Feb 1995 thru Nov 2099

    B06.46 Surface Meridional Wind in m/s Feb 1995 thru Nov 2099

    B06.47 Surface Meridional Wind in m/s Feb 1995 thru Nov 2099

    **********************************************************************

    Plots for Tim Brown; 20 March 2000

    **********************************************************************

    Recent Plots Comparing New "Land Surface Fix" Simulation B06.44 with Old B06 Simulations

    From the first 5 years the precipitation (daily distribution and monthly 5 year average) looks pretty much like another ensemble member to us. We didn't see differences between B06.44 and B06.40 that looked very systematic; rather some months had more and some had less precipitation and this varied from from one latitude to another, when we looked at a transect along the West Coast. We'll be looking at the longer B06.44 run ASAP.

    E-mail discussion on B06.44 versus B06.40 versus Historical

    Global 2m Air Temperature Difference B06.44-B06.40

    JJA 2021 to 2026

    JJA 1995 to 2000

    DJF 2020/21 to 2025/26

    DJF 1995/96 to 2000/01

    Global Total Precipitation Difference B06.44-B06.40

    JJA 1995 to 2000

    DJF 1995/96 to 2000/01

    Land Coastal Strip 2m Air Temperature B06.40 and B06.44

    Feb 1995 to Nov 1995 Land Coastal Strip 2m Air Temperature B06.40 and B06.44

    Feb 2000 to Nov 2000 Land Coastal Strip 2m Air Temperature B06.40 and B06.44

    Land Coastal Strip Total Precipitation B06.40 and B06.44

    Feb 1995 to Nov 1995 Land Coastal Strip Total Precipitation B06.40 and B06.44

    Feb 2000 to Nov 2000 Land Coastal Strip Total Precipitation B06.40 and B06.44

    Binned Total Precipitation

    Sierra (120W 40N) Binned Total Precipitation B06.40 and B06.44


    JJA 2021 to 2026 2m Air Temperature B06.44-B06.40


    JJA 1995 to 2000 2m Air Temperature B06.44-B06.40


    DJF 2020/21 to 2025/26 2m Air Temperature B06.44-B06.40


    DJF 1995/96 to 2000/01 2m Air Temperature B06.44-B06.40


    JJA 1995 to 2000 Total Precipitation B06.44-B06.40


    DJF 1995/96 to 2000/01 Total Precipitation B06.44-B06.40


    Feb 1995 to Nov 1995 Land Coastal Strip 2m Air Temperature B06.40 and B06.44


    Feb 2000 to Nov 2000 Land Coastal Strip 2m Air Temperature B06.40 and B06.44


    Feb 1995 to Nov 1995 Land Coastal Strip Total Precip B06.40 and B06.44


    Feb 2000 to Nov 2000 Land Coastal Strip Total Precip B06.40 and B06.44


    Seasonal Binned Total Precipitation B06.40 (red) and B06.44 (blue)



    Tom - 2 March: I am not surprised that after 5 years B06.44 is not getting "closer" to B06.40, but I anticipate that you will see that it does after 10-15 years. Of course, "closer" is relative to measures of mean+variability.

    Mike - 5 March: Okay, I'll bite. As Dan notes, we haven't yet seen anything that it clearly not "just another ensemble member". Mary just handed me a stack of figures from Tom Bettge's website, and they also look a lot like B0644 is "just another member" by and large. But this same stack confuses me a bit. Why are the "historical (B0628)" traces sooo different from any of the B0640's? Is this a difference between the "free" historical runs and an assimilated initial state? Is there any chance of seeing what the assimilated values also look like? I'm a bit concerned, looking at these, that comparisons between the historical (downscaled) conditions and the b0640 series is likely to be a joke. The good news is that this seems to explain why I have such a large "tear" (towards warmer floodier responses) between the historical runs and the climate-change runs, when I ran b0622 (?) and b0640 through my sierran models and then tried to compare the two. There are a clear and immediate (from first year on) difference between the flood (most notably) statistics from the two runs.

    Ruby - 5 March: I think we should use the first 10-20 years of the B0640 series as the control rather than the historical runs (B0630 series). Even after the initialization problem is fixed for the B0640 series, I suspect they could still be different from the historical runs. The present conditions of each series should be used separately as the control for each set. That's my understanding.

    Mike - 5 March: Yes, I rather gather that. It'll be a major major limitation on my grand objectives. It means that we can't know what the long-term (say, 50-yr recurrence interval extremes) might be in the "historical" simulation, for comparison to the ensemble of climate-change futures (from which such statistics could be extracted). "We're" turning out these lovely multiple long-term futures, with only a (hmmm, multiple) snip at the beginning to compare to. Oh well! More thought is indicated (on my part). My question remains, however: WHY are they so different? I obviously have missed something here. Are the b0640 series runs indeed initiated from data-assimilated versions, different from the "free" historical coupled runs? If so, where is that "assimilated version of history"?

    Tim - 5 March: I also noticed the difference between historical and scenario runs on the PCM page. I bet Dave can shed some light on this. If the SST in the assimilated runs were appreciably different from those of the historical runs in 1995 then that might explain it. I seem to remember this is the case, but Dave will know for sure.

    Dave - 5 March: Hm, I think I am missing the gestalt of this conversation. When I look at the plots I see the following: 1. In volume average temperature, B06.28 is maybe 0.06C cooler than the B06.4x runs. But we changed the T,S field in the B06.4x runs to reflect Detlef's assimilated values. I would expect the volume averaged temperature to be systematically different in the B06.4x cases because of this. (If *something* wasn't different in the B06.4x cases with respect to the historical cases, there would have been no point in doing them!) For a yardstick of what a "big" difference is, note that the difference between the various B06.4x cases spans about 0.05C just due to natural variability. Against that backdrop, the systematic difference due to the assimilated field seems to me to be modest at best. 2. The salinities are different in the B06.4x cases. This is because I changed the salinity field to compensate for the temperature difference -- i.e., it was a density-compensated temperature change. This was done to keep dynamical shocks out of the process, which I believe is important. 3. In surface quantaties, I see no difference between the B06.28 run and the B06.44 (fixed land initialization) run. We know that the B06.40 run (it is specifically mentioned in some of the attached e-mail) had the land initialization problem, so I would advise against using it for anything that is sensitive to land values (especially surface temperature). 4. There are systematic differences in the mean kinetic energy even in the B06.44 run. I'm not sure why, but it seems plausible that holding the T,S values close to the assimilation targets suppresses eddy activity. But if so, there is really no way around this, as it is part of the assimilation process. However, I don't see that this would be a problem for any of the applications or downscaling I have heard about. Maybe if you were looking at eddy fields in the Southern California Bight, but why should the mean K.E. of the ocean affect Sierra precip? Overall, then, it seems to me that the differences I see between the B06.28 run and the B06.44 run are within the range of what I expected to see, are confined to the oceanic interior, and do not yet have a discernable influence on any surface quantity. If the assimilated ocean temperatures never end up having any effect on a surface quantity, then we will conclude that it doesn't matter whether you assimilate 3-D oceanic initial conditions or not, but that remains to be seen.

    Mike - 5 March: Perhaps I should feel chastised, but not quite yet. Yeah, these are small differences in terms of degrees C, but that is not what they represent...surely as a volume average heat content, these must be rather large differences in the scheme of things. (I do like the change that you made to the land initial conditions in b0644 as suggested by the land surface temperatures plot). Yes, the instantaneous influence of these heat (and ice) differences is probably small (since they are deeply buried), but I return to my original question (recast conveniently): Is the historical run really comparable (plus/minus the radiative forcing) to the b0640 series...in the long haul? These are the statistics that caused me to speak up in the first place. I'll take yes for an answer, but not if its only an "in the short term, yes". Sorry for my denseness on this aspect, but I want to know how far to bother with my end of things. If there are (what seem to me to be) large heat content differences in the b0640s initial conditions vs historical run conditions, how do I know what is radiatively forced difference and what is initial conditions?


    For more information please contact Mary Tyree (mtyree@ucsd.edu)