Low (web) Quality - BALTEX
Low (web) Quality - BALTEX
Low (web) Quality - BALTEX
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