The aim of the paper is twofold. Firstly, to analyze the historical data of the earthquakes in the boarder area of Greece and then to produce a reliable model for the risk dynamics of the magnitude of the earthquakes, using advanced techniques from the Extreme Value Theory. Secondly, to discuss briefly the relevant theory of incomplete markets and price earthquake catastrophe bonds, combining the model found for the earthquake risk and an appropriate model for the interest rate dynamics in an incomplete market framework. The paper ends by providing some numerical results using Monte Carlo simulation techniques and stochastic iterative equations.