Riza Andrian IbrahimSukonoHerlina NapitupuluRose Irnawaty Ibrahim2025-04-212025-04-2120252025-4-21Riza Andrian Ibrahim , None Sukono, Herlina Napitupulu, & Rose Irnawaty Ibrahim. (2025). Modeling earthquake bond prices with correlated dual trigger indices and the approximate solution using the Monte Carlo algorithm. AIMS Mathematics, 10(2), 2223–2253. https://doi.org/10.3934/math.20251032473-69881524-5410.3934/math.2025103https://oarep.usim.edu.my/handle/123456789/26589https://www.aimspress.com/article/doi/10.3934/math.2025103?viewType=HTMLhttps://www.scopus.com/record/display.uri?eid=2-s2.0-85219067933&origin=resultslist&sort=plf-f&src=s&sot=b&sdt=b&s=TITLE-ABS-KEY%28Modeling+earthquake+bond+prices+with+correlated+dual+trigger+indices+and+the+approximate+solution+using+the+Monte+Carlo+algorithm%29AIMS Mathematics Volume 10, Issue 2 (page:2223–2253)Countries prone to earthquakes face increasing seismic activity, often resulting in losses that exceed national budgets. To mitigate these losses, earthquake bonds present a promising alternative funding source; however, pricing them is complex, requiring simultaneous accounting for financial and seismic risks. Therefore, this study aimed to model earthquake bond pricing. The model incorporates earthquake intensity to account for rising seismic activity. It also includes depth and maximum magnitude as correlated dual trigger indices, making the bonds more attractive to investors, as claims are generated if both events occur. These three factors were modeled together as a compound stochastic process. The bond price was then formulated using a risk-neutral pricing measure with a stochastic interest rate under the Cox-Ingersoll-Ross model. Since the model lacks a closed-form solution, we employed an algorithm based on the Monte Carlo method for estimation. Through this algorithm, we showed that bond prices for terms of one to six years follow a normal distribution. The use of stochastic interest rates becomes significant as the bond term increases. We also found that earthquake intensity and bond terms negatively correlate with bond prices, while annual coupons positively correlate. Additionally, including dual triggers lowers claim probability and increases the bond demand, but is compensated by higher prices. This study can assist issuers in pricing earthquake bonds based on earthquake severity-maximum magnitude, depth, and intensity-and aid geological institutions in estimating earthquake risk in observed areas.en-USearthquake bondpriceintensitydepthmaximum magnituderisk-neutral pricing measureMonte Carlo methodtext::journal::journal article22232253102