A dataset of historical river discharge into oceans was created using the CaMa-Flood global river routing model and adjusted runoff from the land component of JRA-55. The major rivers were well resolved with a 0.25° horizontal resolution. The total runoff on each drainage basin exhibits a distinctive bias on decadal time scales. The input runoff data were modified using 5-year low-pass-filtered multiplicative factors to fit the annual mean climatology and decadal variations in the reference dataset. The model incorporated data from 1958 to 2016. The yearly and seasonal variations of the major rivers are well represented by the model.