Event-based methods are used in flood estimation to obtain the entire flood hydrograph. Previously, such methods adopted in the UK have relied on pre-determined values of the input variables (e.g. rainfall and antecedent conditions) to a rainfall-runoff model, which is expected to result in an output flood of a particular return period. In contrast, this paper presents a method that allows all the input variables to take on values across the full range of their individual distributions. These values are then brought together in all possible combinations as input to an event-based rainfall-runoff model in a Monte Carlo simulation approach. Further, this simulation strategy produces a long string of events (on average 10 per year), where dependencies from one event to the next, as well as between different variables within a single event, are accounted for. Frequency analysis is then applied to the annual maximum peak flows and flow volumes.