Risk-based cost-benefit analysis of climate adaptation measures for Australian contemporary houses under extreme winds

Climate adaptation measures improve housing resilience to extreme winds, and reduce economic losses associated with wind and rainfall damage under a changing climate. Several adaptation measures are adopted in this study for Australian contemporary houses subjected to non-cyclonic windstorms to either reinforce the building envelope or increase the water resistance of building interior. A risk-based cost-benefit analysis is conducted to evaluate the cost-effectiveness of these adaptation measures that considers the effect of construction defects. It was found that the annual expected losses for houses in Brisbane with construction defects are considerably higher than those without considering construction defects, whereas the influence of construction defects is lower for the Melbourne houses. The cost-benefit analysis reveals that strengthening windows is cost-effective for Brisbane and Melbourne houses. Installing window shutters significantly reduces economic risks associated with extreme winds and is costeffective for houses in Brisbane. Adaptation measures are generally not cost-effective for Melbourne houses due to lower extreme wind speed and associated rainfall.


Introduction
Severe windstorm is one of the major natural disasters in Southeast Australia, which causes millions of annual losses for housing [7]. The climate change impact on the intensity of windstorm and associated rainfall impose more uncertainty in housing damage risks. Climate adaptation measures are thus of significant importance for the improvement of building resilience to windstorms, and the reduction of economic losses associated with wind and rainfall damage under a changing climate. However, the adoption of a specific climate adaptation measure is still in question if its cost-effectiveness is unclear, and therefore there is a need to quantify the costs and benefits of adaptation measures to assess their economic viability.
Many cost-benefit studies have been reported in the literature for houses subjected to tropical cyclones or hurricanes in the United States (e.g. [19,21,41,42]) and Australia (e.g. [35,39]). New South Wales, Victoria and southeastern Queensland of Southeast Australia have nearly 70% of Australia's population, and are classified as non-cyclonic regions in AS/NZS 1170.2 [4]. Only about 5% of Australia's population live in cyclonic regions. Residential construction in Southeast Australia comprises of a large portion of metal-clad contemporary houses, which generally have less wind resistance than houses in cyclonic regions of Australia, and differ from North American houses in materials, construction techniques and building design. Quantitative studies regarding wind risk mitigation and climate adaptation for Australian contemporary houses subjected to noncyclonic windstorms are scarce in the literature. The development of feasible climate adaptation measures and the quantification of their cost-effectiveness are still needed to improve the resilience of residential communities in Australia against wind hazards.
Wind risk mitigation and climate adaptation include measures to either enhance the design during initial construction or retrofit (i.e. upgrade or strengthen) an existing house. Recently, Qin & Stewart [26] developed a probabilistic risk assessment framework for Australian contemporary houses subjected to wind and rainfall damage, and found that most losses result from damage to building interior and contents caused by rainwater intrusion through breaches (e.g. openings due to metal roof sheeting loss and window breakage) and gaps in the building envelope. A set of climate adaptation measures are thus proposed in this study to reinforce the building envelope against wind and improve the water resistance of building interior: (i) strengthening connections for metal roof cladding, (ii) installing shutters for windows, (iii) improving window resistance, and (iv) using waterresistant materials for the building interior.
Most studies regarding the cost-benefit analysis of wind mitigation and climate adaptation measures for housing in Australia are based on empirical vulnerability models to quantify the wind and rainfall damage (e.g. [20,33,36,39]). For example, Stewart [36] employed expert-elicited wind vulnerability curves to assess the housing losses with enhanced building design. Smith & Henderson [33] used insurance data to estimate the cost and benefit of various retrofit and mitigation solutions for houses in cyclonic regions of Australia. However, these empirical vulnerability models have very limited ability to assess the performance of detailed mitigation or adaptation measures applied to specific housing components for wind risk reduction [43]. To overcome this limitation, the cost-benefit analysis in this study is conducted based on the probabilistic risk assessment method developed by Qin & Stewart [26], which systematically integrates engineering judgement, datadriven and engineering-based models for the quantification of extreme wind and associated rainfall, wind damage to housing, rainwater intrusion and economic losses. Such cost-benefit analysis provides a risk-based performance assessment and decision support for the proposed climate adaptation measures.
It has been widely acknowledged that construction defects may significantly increase housing vulnerability under wind hazards [34,45]. Construction defects are common for residential construction, which may also affect the cost-effectiveness of relevant climate adaptation measures. For example, an adaptation measure may provide little reduction in economic losses for houses with good construction quality but it might be worthwhile for houses that have defective components. To the best knowledge of the authors, the influence of construction defects on the assessment of wind damage risks and cost-effectiveness of mitigation or adaptation measures has not yet been taken into account in risk and costbenefit studies. A recent probabilistic model for construction defects in housing by Qin & Stewart [27] is thus incorporated into the present study to examine the effect of construction defects.
The cost-effectiveness of a climate adaptation measure is evaluated in terms of net present value (NPV) that is equal to the benefit minus the cost. An adaptation measure is considered to be economically viable if NPV > 0, which is analogous to the benefit-to-cost ratio or lifecycle cost analysis used in many other studies (e.g. [6,20,21,41]). A break-even economic assessment (e.g. [36]) is also adopted in this study to assess the conditions under which an adaptation measure is costeffective. The break-even analysis is conducted when considerable uncertainties are involved in the assessment of risk reduction and adaptation cost due to either the model limitations and/or a lack of data. As the climate projections for many regions of Australia are under great uncertainty, a scenario-based approach is used in this study to examine the climate change impacts on the cost-benefit analysis.

Wind risk
The risk from extreme wind events is expressed as [36] where Pr(C) is the annual probability that a specific climate scenario will occur, Pr(H|C) is the annual probability of a wind hazard conditional on the climate, Pr(DS|H) is the probability of a damage state conditional on the hazard (also known as fragility), Pr(L|DS) is the conditional probability of a loss given occurrence of the damage, and L is the loss or consequence if full damage occurs. The summation sign in Eq. (1) refers to the summation of a number of possible climate scenarios, hazards, damage states and losses. If the probability that a specific climate scenario will occur, Pr(C) is too unreliable, then a scenario-based analysis where climate scenario probability is decoupled from Eq. (1) is typically adopted for risk assessment and decision analysis (e.g. [36,38]).

Cost-benefit analysis
In this study, the cost-effectiveness of a climate adaptation measure is evaluated based on the NPV, which is calculated as where ΔR is the reduction in risk due to the adaptation measure, E(L) is the risk for the house without any adaptation measures (i.e. 'business as usual') given by Eq.
(1), ΔB is the co-benefit of adaptation such as reduced losses to other hazards, increased energy efficiency, etc., and C adapt is the cost for climate adaptation. For an adaptation measure, ΔR can vary from 0 to 100%. If the adaptation measure is applied at the initial design/construction stage, C adapt is the extra money spent on a stronger design/construction. For retrofitting, C adapt is the cost to remove the old housing component, if applicable, plus the cost to upgrade the component. All the climate adaptation measures (see Section "Climate adaptation") considered in this study only require a one-off expense. The present study focuses on the wind risk mitigation, and therefore the quantification of ΔB is not included. The NPV calculated in this study is a dependent variable mainly due to the variability of ΔR.
The confidence bounds of NPV or the probability of NPV > 0 can also be calculated. However, as the decision analysis in this study is based on a maximum expected net benefit/return criterion, only the mean NPV is of interest. Note that the above equations can be generalised for any time period, discounting of future costs and detailed time-dependent cost and damage consequences. In addition to the NPV, a break-even analysis is also employed in this study as a decision tool when significant uncertainties are involved in the estimation of ΔR and C adapt due to the limitation of relevant models and/or a lack of information. The break-even analysis can be viewed as a type of retrospective analysis, the output of which is the condition that enables the mean of NPV to be zero. A climate adaptation measure is not cost-effective if the risk reduction provided by the measure is lower than, or the adaptation cost is higher than the predicted break-even value. Decision-makers can then judge whether an adaptation measure meets these break-even values.

Representative contemporary house
The risk-based cost-benefit analysis is conducted for the representative contemporary housing built in the southeastern Australian suburbs of Melbourne and Brisbane. These two cities account for about 30% of Australia's population. The representative contemporary house is a timber-framed brick-veneer construction with 21.5°timber roof trusses at 600 mm spacings on a complex hipend roof. Windows are generally horizontal sliding aluminum or timber awning with a brick on edge or terracotta tiled window sill. Roof cladding is 762 mm wide corrugated metal sheeting. Metal top-hat battens are attached to timber roof trusses at 900 mm spacings. The dimension, shape and construction type of the representative contemporary house in Brisbane and Melbourne were determined by field surveys completed by the Cyclone Testing Station (CTS), James Cook University (JCU) (see [22] for more details). Figure 1 shows the 3D and plan view of the representative one-storey house [22]. The design of structural members and components for housing subjected to wind loading is based on the design wind classifications specified in AS 4055 [3] for different site conditions. According to AS 4055 [3], most suburban houses in Brisbane have a design wind classification of N2 or N3 excluding those built on the topthird zone of a hill, ridge or escarpment. A design wind classification of N1 or N2 is appropriate for most suburban houses in Melbourne. Thus, the representative contemporary house herein is considered to have design wind classifications ranging from N1 to N3.

Probabilistic risk assessment
A probabilistic risk assessment (PRA) method developed by Qin & Stewart [26] provides the basis for risk-based decision-making to mitigate housing damage risks from non-cyclonic windstorms under a changing climate. The PRA framework includes probabilistic models for wind hazard and associated rainfall, reliability-based wind damage assessment, rainwater intrusion evaluation and loss estimation. The economic losses considered for contemporary houses arise from wind damage to metal roof cladding and timber roof framing, windward windows, and rainwater damage to building interior and contents as well as the loss of use. This is based on post-damage surveys (e.g. [18,23]) that the majority of losses to contemporary houses result from wind damage to roof and fenestrations (especially windows), and the subsequent rainwater intrusion. The damage to other housing components (e.g. walls) is rare for contemporary houses in non-cyclonic regions of Australia.
Hazard modelling and climate change impacts Extreme wind speed Extreme value analysis is conventionally used to predict extreme wind speed based on recorded data [15]. The Gumbel distribution is employed to model the annual maximum gust wind speed (the maximum 0.2 s gust velocity at 10 m height in open terrain) in non-cyclonic regions of Australia [44]. Figure 2 shows the relationship between gust wind speed v and return period obtained from the analyses by Wang et al. [44].

Rainfall associated with extreme winds
When assessing the rainwater damage to building interior and contents, the joint probability of wind speed, rainfall intensity and storm duration is needed. The exponential distribution is connected to the Poisson arrival process, and commonly used to model storm duration (e.g. [9,17]). The gamma distribution is a widely used parent distribution for rainfall intensity during a storm (e.g. [16,32]).
According to Qin & Stewart [26], the storm duration (D ur ) is modelled by a two-parameter exponential distribution. The average rainfall intensity R h that is non-zero (i.e. rainfall concurs with strong wind) during an extreme windstorm is modelled by a gamma distribution as a function of the storm duration D ur . The probability of no rain during a windstorm is also taken into account. A twenty-year length (1996 to 2015) of half hourly wind and rainfall data for severe windstorms from two weather stations, i.e. Archerfield airport in Brisbane and Moorabbin airport in Melbourne, were used to estimate the parameters in the hazard model. Figure 3 shows the exponential probability plots for D ur given by Qin & Stewart [26]. It is estimated that the probability of no rain during a windstorm is 25.6% and 45.9% for Brisbane and Melbourne, respectively. Figure 4 shows the mean and quantile values of R h produced by the gamma model as a function of D ur as well as the average rainfall intensity data [26]. Figures 3 and 4 reveal that Brisbane tends to have shorter windstorms with more intense rainfall (e.g. thunderstorms), whereas windstorms in Melbourne are generally longer with lower average rainfall intensity. See Qin & Stewart [26] for more details about the hazard modelling and statistical analysis. Note that the probabilistic and statistical models for extreme wind and associated rainfall described in this paper are data-driven rather than relying on physical processes of the natural phenomena.

Climate change projections
Climate change influences the extreme wind speed and associated rainfall intensity. The latest projections [8,10] for changes in extreme wind speed in Brisbane and Melbourne are summarised in Table 1 for medium and high CO 2 emission scenarios RCP 4.5 and RCP 8.5, respectively, to 2090. Extreme wind projections for Melbourne are only available for RCP 8.5 (changes are not available for RCP 4.5). A drying trend with declining annual average rainfall is predicted for Brisbane and Melbourne, while extreme rainfall is projected to become more intense due to increasing water-holding capacity of the atmosphere under a warmer climate [8,10]. However, the quantitative projections for rainfall concurred with strong winds in a future climate are not available because the climate change impacts considering compound events (e.g. simultaneous occurrence of extreme wind and rainfall) in Australia remains unclear [46]. A high level of uncertainty  is also involved in the projected magnitudes of change in wind speed due to the limitation of general circulation models [37]. To this end, a scenario-based analysis is conducted in this study to examine the climate change impacts on the cost-benefit analysis. While there may be scenarios with no change in climate, and even a reduced climate hazard, climate adaptation measures may still be recommended as there is no surety that current codes of practice are optimised for the current climate.

Wind damage assessment Roof fragility
A reliability-based fragility method developed by Qin & Stewart [28] is used to assess the wind damage to metal roof cladding and timber roof trusses for the representative house described in Section "Representative contemporary house". The fragility of the roof system is defined as the extent of the roof sheeting loss and the roof truss failures as a function of the peak gust wind speed. Since the roof connections are generally the 'weakest links' of the roof system [13], the overloading of cladding-tobatten (CTB), batten-to-rafter/truss (BTR) and rafter/ truss-to-wall (RTW) connections (see Fig. 5) is considered as the limit states, the exceedance of which leads to the failure of roof cladding and trusses. The uplift capacities for roof connections were modelled to follow a lognormal distribution [28]. For CTB and BTR connections, the uplift capacities were taken as the lower of the pull-out and pull-over strengths, and the uplift capacities for triple grip RTW connections were the peak loads in the force-displacement curves obtained from experiment tests [31]. The spatially varying wind uplift pressures acting on the roof surface were calculated from the peak gust wind speed and the wind loading parameters mainly including the terrain and height factor, shielding factor, wind directionality factor, external and internal pressure coefficients, etc. All the wind loading parameters were modelled as random variables with mean-to-nominal ratios and coefficient of variation (COV) values given in Qin & Stewart [28]. The corresponding nominal values of these factors can be obtained from AS/NZS 1170.2 [4] for different site conditions. See Qin & Stewart [28] for more details about roof connection resistances and the probabilistic wind loading modelling. Two typical scenarios were assumed for the internal pressurization in the fragility analysis, i.e. (i) dominant openings existing on windward wall and (ii) without any wall openings. A Monte Carlo Simulation (MCS) analysis in conjunction with a finite element (FE) approach were employed to evaluate the wind fragility for roof cladding and trusses under the two wall opening scenarios, which enables the stochastic characterization of spatially varying wind uplift pressure, structural demands and resistances for roof connections, failure progression and load redistribution, and evolution of internal pressure with increasing sheeting loss. A total of 1646 CTB, 532 BTR and 38 RTW connections were involved in the MCS analysis and FE approach. The damage states for metal roof sheeting and timber roof trusses can then be obtained from the fragility assessment. The MCS and FE approach enables the development of two fragility curves: (i) the extent of roof sheeting loss, and (ii) proportion of roof truss failures. Refer to Qin & Stewart [28] for more details about the MCS analysis, FE approach and fragility analysis.
The occurrence of construction defects reduces the wind resistance of housing components. Damage surveys have indicated that the commonly observed defects in residential construction increase the housing vulnerability   [27] to further investigate the wind damage risks and the costeffectiveness of climate adaptation measures for the representative contemporary house with construction defects. The mean proportions of roof sheeting loss and roof truss failures produced by the fragility assessment for the Brisbane and Melbourne house with/without construction defects are shown in Fig. 6. It is indicated that, the effects of construction defects are considerable for the predicted roof cladding fragility, whereas for roof truss fragility, such effects are lower. Note that the RTW connections are different for houses with design wind classifications of N1, N2 and N3. Refer to Qin & Stewart [27] for more details about the fragility analysis and construction defect model.

Window damage
The windward dominant openings and associated rainwater intrusion are considered to result from the breakage of windward windows by high wind pressure. Window failure caused by windborne debris is not considered in this study because windborne debris is less of a concern in non-cyclonic regions of Australia (AS/NZS 1170.22011). The wind pressure acting on the windward window (W win ) is calculated based on the gust wind speed and the wind loading parameters (e.g. terrain and height factor, shielding factor, wind directionality factor, pressure coefficients, etc.). The ultimate strength (R ult ) and water penetration resistance (R water ) of windows are assumed to follow a normal distribution [12] with the mean and COV estimated from the minimum test pressures specified in AS 2047 [2]. For more details, see Qin & Stewart [26] and Table 2. The limit states used for the windward window are shown in Table 3.
Although Australian standard AS 2047 [2] specifies the design window ratings for housing with different site conditions, the authors found that many window manufacturers do not provide detailed window ratings for their products. Window manufacturers are only required to publish energy data, but the structural performance data is non-mandatory in Australia. This may lead to construction error in practice that windows with unsatisfied or mislabeled window ratings are installed. The effect of this type of construction error in window installation on the cost-benefit analysis is examined in Section "Sensitivity analysis".

Subsequent rainwater intrusion
The volumetric rate of rainwater intrusion is mainly dependent on the wind speed, rainfall intensity, and damage states of the building envelope. A semi-empirical rainwater intrusion model [24,40] is modified by Qin & Stewart [26] and adapted to quantify the amount of rainwater intrusion into the damaged metal roof and windward windows. Empirical parameters involved in the rainwater intrusion model are used to account for the building disturbance to the oblique falling rain (i.e. effect of building geometry and aerodynamics), which are spatially variable, and estimated based on experimental tests and/or field monitoring data as well as engineering judgement.
The quantification of rainwater intrusion is also conducted under the two wall opening scenarios: i) presence of windward wall dominant openings, and ii) absence of any wall openings. For the windward dominant opening scenario, the main source of rainwater intrusion considered is water entering from roof and window breaches. For the scenario without any wall openings, rainwater is considered to enter through roof breaches and gaps around undamaged windward windows. The number and locations of failed roof sheets at a given wind speed obtained from the reliability-based fragility assessment described in Section "Roof fragility" are subsequent used in the MCS analysis by applying the semi-empirical rainwater intrusion model. The total volumetric rate of rainwater intrusion can then be obtained. The volume of rainwater intrusion is calculated by the multiplication of the volumetric rate and the rainfall duration after wind damage. See Qin & Stewart [26] for more details about the quantification of rainwater intrusion. Construction defects in metal roof cladding can cause more roof sheeting loss, which may incur more rainwater intrusion through roof breaches. Figure 7 shows the mean volumetric rate of rainwater intrusion as a function of both the gust wind speed and rainfall intensity for the two wall opening scenarios with/without construction defects. As indicated in this figure, the mean rainwater intrusion rate increases with wind speed and rainfall intensity. The nonlinearity of rainwater intrusion with increasing wind speed is because there is more roof sheeting loss at a higher wind speed allowing for more water ingress. The rainwater intrusion rate is higher when construction defects are considered in the roof fragility assessment.

Risk analysis Loss estimation
The loss estimation uses an assembly-based approach (e.g. [11,12,25,38]). The representative contemporary house is divided into components/subassemblies based on specific building details as shown in Table 4. The loss estimation takes into account the direct losses from wind damage to windward windows, metal roof cladding and timber roof framing, and the losses to building interior and contents caused by subsequent rainwater intrusion as well as the loss of use.
The losses are estimated in terms of cost ratios, which is defined as the ratio of the cost to complete the subassembly to the building value. The estimated total cost to build a new contemporary house with an approximate floor area of 150 m 2 is L building = $300,000 Australian Dollars [14,30]. Based on cost data provided by Australian housing cost guides [29] and subjective judgement, the subassembly cost ratios for new construction are estimated for a representative contemporary house built to an average standard. The cost ratios are adjusted to account for the additional costs associated with removal, repair and remodelling of an existing house [12]. The adjusted cost ratios are also given in Table 4. The  building interior herein includes internal finishes and fittings, mechanical and electrical systems. The cost ratio of contents is estimated to be 25% of the building value [1]. The empirical loss functions are given by Qin & Stewart [26] to relate economic losses with the damage states of housing components, which are mainly based on engineering judgement and existing loss models given in HAZUS [12]. The losses to roof cladding and framing are directly expressed as a function of the roof damage states, whereas the building interior and contents losses are expressed as a function of the rainwater intrusion, which is indirectly related to roof cladding and window damage. The window losses are estimated based on the window damage model described in Section "Window damage". The loss function for loss of use is taken from HAZUS [12] which is a function of the expected total building loss. See Qin & Stewart [26] for more details about the loss functions.

Annual risk
The annual risk (expressed as the expected loss) for a climate change scenario at time t is where f (v, t) is the probability distribution of the annual maximum gust wind speed at time t, f(R h , t|v, D ur ) is the probability distribution of the average rainfall intensity of a severe windstorm corresponding to a given duration D ur at time t, f (D ur ) is the probabilistic distribution of the windstorm duration which is assumed to be independent of climate change, n d = 8 is the number of cardinal wind directions considered in this study, n c = 6 is the number of components/subassemblies considered in the loss estimation, Pr(DS|v, R h , D ur ) is the probability of damage state (e.g. extent of roof damage, amount of rainwater intrusion) given the gust wind speed, rainfall intensity and storm duration, Pr(L i |DS) is the loss likelihood for the i th component/subassembly given the damage state, L i is the maximum probable loss for the i th component/subassembly, and r is the discount rate. As described in Section "Loss estimation", L i is the cost ratio for the i th component/subassembly given by Table 4, and Pr(L i |DS) is evaluated based on the empirical loss functions given the damage states of the roof and windows, and the amount of rainwater intrusion. More details can be found in Qin & Stewart [26] and (Qin: Risk Assessment and Mitigation for Australian Contemporary Houses Subjected to Non-cyclonic Windstorms, under review). It is noted that Eq. (3) assumes that damage is caused by the largest wind event in any calendar year, which will slightly underestimate damage risks in the event of a lesser damaging windstorm in the same year. The probabilistic risk assessment conducted using a MCS analysis consists of four major components, i.e. i) hazard modelling for extreme wind and associated rainfall, ii) reliability-based wind damage assessment for roof and windows, iii) evaluation of rainwater intrusion, and iv) loss estimation. Figure 8 shows an outline to illustrate the risk analysis method by Qin & Stewart [26] to assess the annual expected economic losses for the representative contemporary house subjected to non-cyclonic windstorms. Table 5 shows the annual expected losses in the present year (normalized by the building value) for Brisbane and Melbourne houses with different design wind classifications with/without construction defects. The average days for loss of use and the expected losses for each housing component/subassembly are also given in this table. Table 5 suggests that losses to building interior and contents caused by rainwater intrusion is the major contributor to the annual expected loss. Houses in Brisbane are generally subjected to higher losses than houses in Melbourne because the extreme wind speed and rainfall intensity are higher in Brisbane as shown in Section "Hazard modelling and climate change impacts". The annual expected losses for houses in Brisbane with construction defects are considerably higher than that without considering construction defects, whereas the effect of construction defects on the calculated risk is slight for Melbourne houses. This is expected because the metal roof sheeting loss is negligible for the Melbourne house even with the consideration of construction defects, and there is not much increase in subsequent rainwater intrusion from roof breaches.

Adaptation measures
The climate adaptation measures proposed in this study mainly aim to reduce the rainwater damage to building interior and contents. The fragility and risk analysis show that wind damage to roof cladding is highly likely to occur for the representative contemporary house in Brisbane. The failure of metal roof sheets can further incur rainwater intrusion through the roof breaches,  though the direct loss from roof cladding damage is limited due to a small cost ratio as shown in Table 4. Window breakage by high wind pressure creates windward dominant openings, which allows for more rainwater intrusion and roof damage due to significantly increased internal pressure. Rainwater intrusion via small gaps around undamaged windows may also occur when the water penetration resistance of window is exceeded. Strengthening roof cladding and windows not only decreases the direct losses from wind damage to the building envelope, but also reduce the amount of rainwater intrusion. In addition to reinforcing the building envelope, the use of water-resistant materials for building interior can also reduce the rainwater damage.
To this end, four adaptation measures are proposed herein: (i) increase of uplift capacities for cladding-to-batten (CTB) connections (RF), (ii) protection of windows by shutters (WS), (iii) improvement of window resistances (WR), and (iv) use of water-resistant materials for ceiling and internal wall finishes (IW). The fragility analysis reveals that the damage proportions of BTR and RTW connections are considerably less than the failures of CTB connections [28], and hence RF is to strengthen CTB connections by (i) increasing the base metal thickness (BMT) of metal roof cladding from 0.42 mm to 0.48 mm, and (ii) upgrading the screw fastener for CTB connections from M6-11 (diameter of 6.0 mm and 11 threads per inch) to 14-12 (diameter of 6.3 mm and 12 threads per inch). This improves the mean pull-over and pull-out capacities of CTB connections by 25% and 7%, respectively. In WS, cyclone-rated steel roller shutters tested for wind pressure up to 2500 Pa are used. The shutter strength against wind Fig. 8 Outline of the risk analysis method [26]  pressure is assumed to follow a normal distribution with a mean of 3000 Pa and a COV of 0.20, which implies that 20% of the shutters cannot satisfy the test pressure to account for the variance of quality in manufacture and installation. A triangular distribution with a lower bound of 50% and an upper bound of 100% as shown in Fig. 9 is assumed for the reduction of wind pressure and winddriven rain on windows due to the installation of shutters, if not damaged. HAZUS [12] assumes the reduction is 50% which is considered as a lower bound for a cyclonerated roller shutter. Window resistances may be increased by using a higher window rating (i.e. N1 to N2, N2 to N3 and N3 to N4). Note that improving window ratings does not necessarily increase the water penetration resistances of windows as shown in Table 2. By adding some additives in manufacture, the water absorption rate of plasterboard can be reduced. Gypsum plasterboard with specified water-resistant grade given by AS/NZS 2588 [5] may be used as the internal linings for ceiling and wall. Waterresistant gypsum plasterboard is required to have an average water absorption of less than 5% of self-weight after two-hour immersion under a minimum of 25 mm of water [5], and a triangular probability distribution given by Fig. 10 is assumed to model the reduction of rainwater damage to internal linings. The internal linings (i.e. ceiling and wall finishes) consist of 30% of the building interior value as shown in Table 4. The proposed climate adaptation measures are summarized in Table 6. The adaptation cost (normalized by building value) and the improved performance due to the adaptation measure (i.e. benefit in Table 6) are also shown in Table 6. The adaptation measures applied at the initial design and construction are RF, WR and IW, whereas window shutters (WS) can be installed at any time during the service life to retrofit an existing house. The adaptation cost is the additional money spent on the adaptation measure, which is considered as a one-off expense. The adaptation costs presented in Table 6 are estimated from an Australian housing construction cost guide [29]. The effect of variations in adaptation costs on the cost-benefit analysis is examined through a break-even analysis in Section "Evaluation of cost-effectiveness".

Cumulative risks
The cumulative expected loss (or risks) over the 50-year service life of the house from 2020 to 2070 is given by where E annual (t) is the annual risk at time t, which is influenced by the climate change and discount rate. In this study, the cumulative economic risks are presented in 2020 Australian dollars and a discount rate of r = 4% is used [38]. A sensitivity analysis is conducted in Section "Sensitivity analysis" using different discount rates. Figure 11 shows the cumulative expected losses to 2070 considering construction defects. A constant climate is assumed in this section (i.e. no climate change impacts on wind and rainfall), and the climate change impacts are examined in Section "Climate change impacts". The cumulative expected losses to 2070 and risk reductions (shown in bracket) provided by the adaptation measures are shown in Table 7. In this section, it is assumed that the adaptation measure WS is applied at the new construction (i.e. Year 2020) to maximise its benefits (i.e. risk reduction). Figure 11 and Table 7 suggest that installing window shutters (WS) is the most effective adaptation measure that provides over 95% risk reduction. Strengthening roof connections (RF) provides up to 8% risk reduction for Brisbane houses but it is not effective for Melbourne houses, which is expected as roof failure is rare and most losses for houses in Melbourne are attributed to rainwater intrusion through windward windows. Improving window ratings (WR) offers more risk reduction when both the ultimate strength and water penetration resistance of windows are increased (i.e. N2 to N3), whereas less risk Fig. 9 Reduction of wind pressure and wind-driven rain on windward windows due to shutters reduction is achieved if only the ultimate strength is increased (i.e. N1 to N2, N3 to N4). Using water-resistant materials for internal linings (IW) is generally effective, and provides about 20% risk reduction. As expected, considerably higher cumulative expected losses are predicted for houses in Brisbane with construction defects compared to those without construction defects, whereas the effect of construction defects on cumulative risks for Melbourne houses is slight.

Evaluation of cost-effectiveness
The cost-benefit analysis is conducted using the annual and cumulative expected losses evaluated based on the PRA method by Qin & Stewart [28] as described in Section "Risk analysis", which provides a risk-based decision support tool for the selection of cost-effective climate adaptation measures for houses subjected to non-cyclonic windstorms. The procedures for the cost-benefit analysis are summarized as: (i) Calculate the annual expected loss in year t for the considered house without adaptation using the PRA method described in Section "Risk analysis", and then calculate the cumulative expected loss for the 50-year building service life using Eq. (4); (ii) Calculate the annual expected loss in year t for the considered house with a given adaptation measure using the PRA method described in Section "Risk analysis", and then calculate the cumulative expected loss for the 50-year building service life using Eq. (4); (iii)Calculate the reduced cumulative expected loss due to the adaptation measure, i.e. E(L)ΔR in Eq.
(2); (iv) Implement (ii) and (iii) for all considered adaptation measures, and then calculate the mean NPVs corresponding to these adaptation measures using Eq. (2); (v) Decision-making based on the mean NPVs corresponding to different adaptation measures. Adaptation measures with the mean NPV greater than zero are deemed as cost-effective, and those with a higher mean NPV is preferred by assuming a risk neutral decision-maker whose primary goal is economic efficiency. Table 8 shows the mean NPVs to 2070 for the representative contemporary house with/without construction defects and no climate change is assumed. Table 8 suggests that, in most cases, the mean NPVs of adaptation measures with the consideration of construction defects are higher than those without considering construction defects. This is expected because construction defects result in higher cumulative risks and a specific adaptation measure can generally offer more reduction in ecnomic losses (i.e. E(L)ΔR) at a given adaptation cost. For the Brisbane house with a design wind classification of N2, the adaptation measures yielding a positive mean NPV are WS and WR. No adaptation measure is cost-effective for the Brisbane house with a design wind classification of N3 without considering construction defects, however, WS turns to a positive mean NPV when construction defects are taken into account. This illustrates the effect of construction defects on the cost-benefit analysis. None of the adaptation measures may be recommended for the Melbourne house with a design wind classification of N1 as the corresponding mean NPVs are negative. Table 8 also shows that WR is the only costeffective adaptation measure for the Melbourne house with a design wind classification of N2.
The adaptation costs given in Table 6 are estimated according to the average cost data given by Australian housing construction cost guide [29]. However, the costs may vary among different locations, contractors and labourers. A break-even analysis is thus conducted to calculate the maximum adaptation cost that enables the mean NPV to be zero with the consideration of construction defects, see 0.7% Improve the uplift capacities of CTB connections with the mean pull-over and pull-out capacity increased by 25% and 7% respectively.
Window WS. Install shutters Cyclone-rated steel roller shutters tested for wind pressure up to 2500 Pa.

2.0%
Reduce wind pressure and the amount of wind-driven rain acting on windows as shown in Fig. 9, and hence decrease the probability of windward dominant openings.

WR. Higher window rating
Increase the window rating to a higher level (e.g. N2 to N3).
0.4% Increase the ultimate strength and water penetration resistance of windows according to Table 2.

Building interior IW. Water resistant materials
Use water-resistant plasterboard for ceiling and internal wall finishes.
1.0% Improve water resistance of building interior and reduce rainwater damage as shown in Fig. 10.  Table 9. It is indicated from Table 9 that an adaptation measure may turn to be cost-effective with reasonable discounts on the cost given by Table 6. For example, only a 13% discount on the adaptation cost makes WR costeffective for the Brisbane house with a design wind classification of N3, whereas a 50% reduction in cost is needed to enable IW to be cost-effective for the Brisbane house with a design wind classification of N2. A few adaptation measures can remain cost-effective even with significantly increased costs. For example, WR is still cost-effective for the Brisbane house with a design wind classification of N2, if the anticipated cost of WR given by Table 6 increases by 300%.

Climate change impacts
A scenario-based analysis is conducted to examine the costeffectiveness of the climate adaptation measures under a changing climate. The median, 10th and 90th percentile projections for wind speed change under RCP 4.5 and RCP 8.5 as shown in Table 1 are adopted, and the changes in rainfall intensity to 2090 conditional on an extreme wind climate scenario are assumed to be − 20%, − 10%, 0%, + 10% and + 20%. Information is scarce to non-existent on timedependent changes in extreme wind speed and associated rainfall for Australia. A time-dependent linear change [38] is then assumed for both wind speed and rainfall intensity. The mean NPVs for different adaptation measures to 2070 with the consideration of construction defects and climate change impacts are shown in Fig. 12. The mean NPVs in the figure are for Brisbane and Melbourne houses with a design wind classification of N2 as they are more vulnerable than N1 and N3 houses. The error bars in Fig. 12 represent the 10th and 90th percentile projections for wind speed change under RCP 4.5 and RCP 8.5. It is observed that, in general, the larger the wind speed and the associated rainfall intensity in a future climate, the higher the mean NPVs produced by the climate adaptation measures. Climate change tends to have higher influences on the cost-effectiveness of adaptation measures applied to houses in Brisbane.
Installing window shutters (WS) and increasing window ratings (WR) remain cost-effective for the Brisbane house (N2), and WR is cost-effective for the Melbourne house (N2) under all the considered climate scenarios.   For the most adverse climate scenario considered in this study (i.e. 90th percentile projection for extreme wind speed under RCP 8.5 and a 20% increase in rainfall intensity), the mean NPVs increase by approximately 15% and 10% respectively for WS and WR compared to the scenario with no climate change. Therefore, the effect of climate change on the cost-benefit analysis is not significant. Strengthening roof connections (RF) and using water-resistant materials (IW) are not cost-effective under all the considered climate change scenarios. In addition, the co-benefits (ΔB given in Eq. (2)) associated with RF and IW are negligible, whereas considerable cobenefits exist for WS and WR. Using a higher window rating (WR) may provide energy efficiency due to increased glass thickness. Window shutters (WS) can offer more co-benefits such as noise reduction, prevention of heat loss, light control and home security. Although the mean NPV of WR is higher than that corresponding to WS for the Brisbane house (N2) under a changing climate, WS is still a competitive option for climate adaptation due to considerable co-benefits. In terms of implications for insurance premium, an insurer may offer incentives by a discount in annual insurance premium for homeowners who decide to retrofit their houses by installing window shutters.

Sensitivity analysis
The adaptation measure WS is considered as a retrofit that can be applied at any year during the service life. Figure 13 shows the mean NPV to 2070 for houses with no climate change when installing window shutters at different points in time. It is indicated that WS is costeffective for houses in Brisbane if applied no later than 2048 for N2 houses and 2035 for N3 houses. The mean NPVs to 2070 increase for houses in Melbourne if window shutters are installed in a later year.
Two discount rates (i.e. 2% and 7%) are further adopted in the cost-benefit analysis, and the corresponding mean NPVs with the consideration of construction defects and no climate change are given in Table 10. As expected, the mean NPV increases for a lower discount rate, and decreases for a higher discount rate. Installing window shutters (WS) for Brisbane houses with design wind classifications of N2 and N3 is no longer cost-effective under a discount rate of 7%. Increasing window ratings (WR) for the Brisbane house with a design wind classification of N3 turns to be costeffective under a discount rate of 2%.
The construction error of using mislabelled windows as described in Section "Window damage" is not considered in the cost-benefit analysis described above. A scenario-based approach is used herein to examine the effects of mislabelled windows. It is assumed that N1rated windows are incorrectly installed on Brisbane and Melbourne houses with a design wind classification of N2, and N2-rated windows are incorrectly installed on Brisbane houses with a design wind classification of N3.  This increases the cumulative expected losses to 2070 by 32%, 4% and 96% for the Brisbane house (N2), Melbourne house (N2) and Brisbane house (N3), respectively. Table 11 shows the mean NPVs for adaptation measures RF, WS and IW if the mislabelled windows are installed. It is suggested that RF and IW is still not costeffective, and WS is cost-effective for Brisbane houses. When mislabelled windows are installed, there is an approximately twofold and tenfold increase of the mean NPV yielded by WS for Brisbane houses with a design wind classification of N2 and N3, respectively, compared to that without such construction error. This implies that the adverse effects of construction error of using mislabelled windows can be largely counter-balanced by installing window shutters.

Conclusions
In this study, several climate adaptation measures are adopted for Australian contemporary houses subjected to non-cyclonic windstorms to either reinforce the building envelope or increase water resistance of the building interior. A probabilistic risk assessment was conducted considering the effect of construction defects, and a risk-based cost-benefit analysis was then used to evaluate the costeffectiveness of these adaptation measures. The risk analysis results suggest that the annual expected losses for houses in Brisbane with construction defects are considerably higher than those without considering construction defects, whereas the influence of construction defects is slight for the Melbourne houses. The cost-benefit analysis reveals that strengthening windows by increasing window ratings (e.g. N1 to N2, N2 to N3) is cost-effective for Brisbane houses and Melbourne houses. Installing window shutters significantly reduces economic risks incurred by wind and rainfall damage and is cost-effective for houses in Brisbane. The adverse effects of construction error of using mislablled windows can be largely counter-balanced by installing window shutters. Installation of window shutters as a retrofit also provides many co-benefits and has the potential to attract incentives in annual insurance premium, which makes it a strong candidate for climate adaptation.