The derivation of the Poisson distribution indicated that P_k =P^k/k!. If we compare the ratio of the next P_k to P_k we get λ/(k+1). The ratio of succeeding numbers in the fit is the same. This ratio is not constant but steadily decreases with k (or M). At the peak the ratio is 1 and prior to that it is greater than 1 which explains the growth in numbers. So the Poisson fit suggests a mechanism in which there are changes in what happens for each state of the Markov process. In the simple mechanism considered the ratio between succeeding states was 1-r. If r is a fraction this the ratio does not exceed 1 and so there is no peak. Any mechanism would also have to explain the magnitudes as well as how the number of earthquakes depend on the occupancy of the states. Perhaps a delay in advancement or some blockage in the advancement would give a more realistic fit.
The computed distribution gave equilibrium values. Blockage at some magnitude would prevent advancement from one state to the next and one would expect the occupancy of the states below the point of blockage to increase relative to the equilibrium values and those above to decrease somewhat. One needs to consider processes which might compete with the occurrance of earthquakes. The goal is a better explaination of the process behind earthquakes based on evidence. It is a search for knowledge.