09.11.2014 Views

Low (web) Quality - BALTEX

Low (web) Quality - BALTEX

Low (web) Quality - BALTEX

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

162<br />

JP10: 59-year 10 km dynamical downscaling of re-analysis over Japan<br />

Hideki Kanamaru, Kei Yoshimura, Wataru Ohfuchi, Kozo Ninomiya and Masao Kanamitsu<br />

Climate Change and Bioenergy Unit, Viale delle Terme di Caracalla 00153 Rome, Italy. Hideki.Kanamaru@fao.org<br />

1. Introduction<br />

To date, multi-decadal dynamical downscaling simulations<br />

driven by the Reanalysis data have been completed over<br />

many regions (e.g., Vidale et al., 2003), but only a few are<br />

in very high resolution, which is required for application<br />

studies such as agricultural study or energy and water<br />

resources management. As one such very high resolution<br />

long-term simulation, Kanamitsu and Kanamaru (2007;<br />

KK07 hereafter) conducted a 10 km resolution run over<br />

California for a half-century. KK07 revealed that these<br />

simulations had better skill than the regional Reanalysis<br />

product, NARR, whose spatial resolution (32 km) is coarser<br />

than the dynamical downscaling simulations (10 km),<br />

because the current data assimilation system is incapable of<br />

effectively utilizing high-density near-surface observations,<br />

and places more weight on the initial guess produced by the<br />

regional high-resolution numerical model.<br />

In this study, we conduct a similar long-term high resolution<br />

dynamical downscaling as CaRD10 over Japan. The product<br />

is hereafter called JP10. The major objective is to<br />

demonstrate that the dynamical downscaling is capable of<br />

reproducing small scale detail which agrees better with<br />

station observations than the coarse resolution analysis,<br />

without injecting small scale observation. To do so, the<br />

accuracy of JP10 will be compared with independent high<br />

density in-situ observation data, and how JP10 reproduces<br />

some historical extreme events will be also investigated.<br />

2. Model and Observation<br />

The Regional Spectral Model (RSM; Juang and<br />

Kanamitsu 1994) originates from the one used at the<br />

National Centers for Environmental Prediction (NCEP), but<br />

the code was updated with greater flexibility and much<br />

higher efficiency (Kanamitsu et al. 2005) at the Scripps<br />

Institution of Oceanography. The RSM utilizes a spectral<br />

method (with sine and cosine series) in two dimensions. A<br />

unique aspect of the model is that the spectral decomposition<br />

is applied to the difference between the full field and the<br />

time-evolving background global analysis field. The model<br />

configuration and the downscaling method in this study are<br />

basically the same as that of CaRD10 (10 km California<br />

Reanalysis Downscaling; Kanamitsu and Kanamaru, 2007)<br />

but for a domain covering Japan Islands (22.123°–49.163°N<br />

and 119.960°–151.577°E; shown in Figure 1a with<br />

topography) for 1948 to 2006, and with narrower lateral<br />

boundary nudging zones that extends only 2.5% of the total<br />

width in each of four lateral boundaries instead of 11.5% in<br />

CaRD10 to increase the useable domain.<br />

As same as CaRD10, a spectral nudging scheme, i.e.,<br />

scale selective bias correction (SSBC, Kanamaru and<br />

Kanamitsu 2007), is applied to the Reanalysis large scale<br />

thermodynamic fields for a 10 km horizontal resolution<br />

downscaling simulation, to reduce the growth of large-scale<br />

error spanning the regional domain. The scheme consists of<br />

three components: 1) dampening the large-scale (more than<br />

1000 km scale) part of the wind perturbation toward zero,<br />

with dampening coefficient of 0.9, 2) removing the area<br />

average perturbation of temperature and moisture at every<br />

model level, and 3) adjusting the area mean perturbation<br />

logarithm of surface pressure to the corresponding<br />

difference of logarithm of surface pressure due to the area<br />

mean difference in the global and regional topography.<br />

For verification of the JP10 dataset, this study uses<br />

hourly precipitation, surface temperature, and surface<br />

wind direction and speed from the AMeDAS re-statistic<br />

dataset (JMBSC, 2007), which is a surface meteorological<br />

observation dataset of 29 years (1976-2004) from over<br />

1,300 Japanese automated meteorological data-acquisition<br />

system (AMeDAS) sites. These observation sites are<br />

distributed across Japan at intervals of about 20 km on<br />

average. The locations of all the observatories are shown<br />

in Figure 1b. These about 1300 sites are divided into 9<br />

groups in accordance to the geographical divisions in<br />

Japan Meteorology Agency (JMA).<br />

Figure 1. (a) Domain and topography and (b)<br />

AMeDAS observatories location. Nine regions are<br />

classified as N-Hokkaido (red), S-Hokkaido (green),<br />

Tohoku (blue), Kanto (yellow), Chubu (pink), Kinki<br />

(sky blue), Shikoku (orange), Kyushu (purple), and<br />

Okinawa (light green).<br />

3. Results<br />

Climatology<br />

Figure 2 shows monthly climatology variations averaged<br />

over each of the nine AMeDAS regions. It is notable that<br />

there is significant systematic over estimation of<br />

precipitation all over Japan whereas air temperature is<br />

almost perfectly agreed. By intensive investigations, we<br />

found that this is due to the over correction of humidity<br />

field by this version of SSBC.<br />

Daily variations<br />

As shown in Figure 2, downscaled daily variations of<br />

surface meteorology at one of the AMeDAS observatories<br />

in Okinawa Island are well reproduced, particularly for<br />

wind and air temperature. Precipitation includes errors in<br />

amount and timing. In Figure 3, daily air temperature in a<br />

month at each observatory is compared with those of JP10,<br />

and the distributions of their correlation coefficients are<br />

shown for four months (Feb., May, Aug., and Nov.) in<br />

1976. This figure indicates that high frequency variations<br />

of the air temperature associated with the synoptic scale<br />

weather changes are accurately captured in the<br />

downscaled product (in average, more than 0.8 of<br />

correlation coefficient), whereas the original coarse scale<br />

Reanalysis could not reproduce them in this degree of<br />

accuracy. Note that the accuracy is relatively low in

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!