The ever increasing demand for container transportation has led to the congestion of maritime container terminals in the world. In this work, the two interrelated problems of berth and quay crane scheduling are considered in an integrated multiobjective mathematical model. A special character of this model is that the arrival times of vessels and the failure (working) times of quay cranes are not deterministic and can vary based on some scenarios. Hence, a robust model is devised for the problem having three objectives of minimising the deviations from target berthing locations and times as well as departure delays of all vessels. This robust optimisation seeks to minimise the value of the objectives regarding all the scenarios. An exact solution approach based on the 𝜖-constraint method by the Gurobi software is applied. Moreover, regarding the complexity of the problem, two Simulated Annealing (SA) based metaheuristics, namely a Multi-Objective Simulated Annealing (MOSA) and a Pareto Simulated Annealing (PSA) approach are adapted with a novel solution encoding scheme. The three methods are compared based on some multiobjective metrics and a statistical test. The advantage of the integration of berth and quay crane scheduling is examined as well.