To enable power supply in rural areas and to exploit clean energy, fully renewable power systems consisting of cascaded run-of-the-river hydropower and volatile energies such as pv and wind are built around the world. In islanded operation mode, the primary and secondary frequency control, i.e., hydro governors and automatic generation control (AGC), are responsible for the frequency stability. However, due to limited water storage capacity of run-of-the-river hydropower and river dynamics constraints, without coordination between the cascaded plants, the traditional AGC with fixed participation factors cannot fully exploit the adjustability of cascaded hydropower. When imbalances between the volatile energy and load occur, load shedding can be inevitable. To address this issue, this paper proposes a coordinated tertiary control approach by jointly considering power system dynamics and the river dynamics that couples the cascaded hydropower plants. The timescales of the power system and river dynamics are very different. To unify the multi-timescale dynamics to establish a model predictive controller that coordinates the cascaded plants, the relation between AGC parameters and turbine discharge over a time interval is approximated by a data-based second-order polynomial surrogate model. The cascaded plants are coordinated by optimising AGC participation factors in a receding-horizon manner, and load shedding is minimised. Simulation of a real-life system with real-time pv data collected on site shows the proposed method significantly reduces load loss under pv volatility.