Qualitative Assessment of Social Vulnerability to Flood Hazards in Romania

This paper investigates local-scale social vulnerability to flood hazards in Romania, aiming to identify the most vulnerable social and demographic groups across a wide range of geographical locations by considering three dimensions: demographic, socioeconomic, and the built environment. The purpose of the paper is threefold: first, it strives to improve the Social Vulnerability model (SoVI®) by applying a different weighting method adapted to the Romanian context, taking into consideration the municipalities exposed to flood movements. Second, it aims to develop an assessment model for the most vulnerable communities by measuring the heterogeneity according to local indicators related to disaster risks. Third, it aims to facilitate emergency managers to identify community sub-groups that are more susceptible to loss and to increase the resilience of local communities. To perform local-level vulnerability mapping, 28 variables were selected and three aggregated indexes were constructed with the help of the ArcGIS software. Moreover, a model of Geographically Weighted Regression (GWR) between communities directly affected by floods and localities with highand very high values of the Local Social Vulnerability Index (LoSoVI) was used to explore the spatial relationship among them and to compare the appropriateness of Ordinary Least Square (OLS) and GWR for such modelling. The established GWR model has revealed that the negative effects of flood hazards are often associated with communities with a high degree of social vulnerability. Thus, the analysis is able to provide a more comprehensive picture on communities in desperate need of financial resources in order to have the ability to diminish the negative impacts of flood hazards and to provide a more sustainable society.


Introduction
Over the last two decades, the importance of analysing social vulnerability has gained major importance, especially in terms of elaborating different methodologies for quantifying the rate of vulnerability, delimiting territories most exposed to risks as well as elaborating effective development strategies.Quantifying the social dimension of vulnerability allows the identification of the most vulnerable regions along with key factors which, once addressed, could increase the resilience of local communities.By analysing a broad set of studies on topics including methodological development [1][2][3], specific hazard research (climate, floods) [2,4], scale difference analysis (local, international comparison) [1,2], temporal and spatial changes [5], or uncertainty and sensitivity analysis [6], we can see that social vulnerability analysis has come a long way in the last decade.However, due to its multidimensionality, neither a common definition nor a universally accepted measurement has been unanimously accepted so far.Social vulnerability is used, defined, and conceptualised in many different ways [7] and is often linked to associated concepts such as resilience, risk, exposure, sensitivity, and coping capacity [8].In this article we first present the concepts and approaches related to social vulnerability in a critical way, debating afterwards on the issue of developing methods for measuring social vulnerability as well as the advantages and weaknesses of the most frequently used variables.
The evolution of research on vulnerability has been considerably influenced by two distinct approaches: namely, the 'human ecologist school' or 'behavioural paradigm' and the 'structural paradigm'.
The focus of the behavioural paradigm is the human adaptation to natural hazards.Adaptation can be understood as purposeful actions for disaster mitigation, the behavioural paradigm being thus also viewed as a hazard-based approach [9][10][11].The structural paradigm-also called a disaster-based approach-focuses on people's individual socio-economic and demographic characteristics within the specific social, economic, cultural, and political context they live in, emphasizing a range of variables determining vulnerability.The structuralism perspective has also led to important theoretical advancements in the field, resulting in the development of various comprehensive models.
The 'Pressure and Release Model' (PAR) reveals the complex interactions between social processes that create vulnerability and the hazards themselves.Within the model, 'pressure' can be understood as an increasing level of vulnerability and exposure to hazards, while the 'release' aspect encompasses all actions taken to reduce the potential negative impact of a disaster.Although explicitly highlighting vulnerability, the PAR model has been strongly criticized as it does not take into account the role of proximity to the source of the threat and it does not address the physical aspects of the interactions between the social-and natural systems [12,13].The model provides little detail on the hazards' causal sequence and downplays feedback beyond the system of analysis included in the integrative Risk-Hazard models [14].Thus, it represents a specific tool for explaining vulnerability within descriptive analyses, rather than for empirical testing.
An important conceptual development has emerged along with the elaboration of the hazard-of-place approach to vulnerability, which combines the risk/hazard and political ecology perspectives.The model is based on the assumption that people living in a certain area will have an unequal vulnerability to a potential natural hazard due to their social, economic, and political statuses [15].In the model, the geographic context and social structure are linked together, having a strong influence on each other.Thus, the interactions between social vulnerability and biophysical vulnerability create place vulnerability, which is then linked to risk and mitigation elements.The most important deficiencies of this model are its weak capability to explain the root causes of social vulnerability, as well as the inability to include the larger context in which social vulnerability exists [12].However, this approach can be considered a novel contribution with regard to the introduction of GIS technologies not only in risk, hazard, and vulnerability mapping, but also in emergency planning.
Following the internal-external model, which relies on the double structure of vulnerability (exposure/coping), Turner et al. [13] proposed a vulnerability/sustainability framework in order to provide a comprehensive analysis on the complexity and interactions involved in vulnerability analyses, highlighting the factors and linkages that affect the intertwined human and environmental systems exposed to hazards.Since vulnerability is strongly influenced by global processes, the conceptual model stresses the importance of analysing the elements of local vulnerability within the wider context, explicitly considering space, time, and scale.This framework is thus able to show how the interactions of social and environmental forces can result in an increased vulnerability in the event of sudden changes.However, the approach could not escape criticism as it lacks clear differentiation between exposure and sensitivity, making it difficult to identify where vulnerability begins and where it ends [12].
The concept of vulnerability has continuously evolved by including susceptibility, exposure, the capacity to cope and adapt, and by incorporating different thematic areas such as physical, social, economic, environmental and institutional vulnerability [16].Along with the extensive debate on the conceptualisation of social vulnerability, there are also debates on the issue of developing methods for measuring it.The difficulties have mainly been related to two specific aspects: the availability of databases at different territorial levels and the development of the most appropriate methodology which considers all factors of interest, providing necessary information upon which political decisions can be based.Even though there is still no universally agreed-upon set of indicators, methodology or uniquely developed indices for measuring social vulnerability, a series of successful attempts have been published on the topic in recent years.One of the most important pioneering works belongs to Cutter et al. [1] who-by examining social vulnerability in US states-have developed the SoVI ® index.This method has largely been accepted due to its continuous evolution [4,[17][18][19][20], pertinence, and adaptability.The only limitation of the method is that it attributes equal weights to all factors influencing social vulnerability, although there are always certain processes-depending on the analysed region-which can have a stronger impact on shaping social vulnerability.In order to avoid this limitation, most of the research results which have seen the light of day in the last years have focused either on applying different weighting methods [21,22] or on developing completely new methodologies [23,24].
When it comes to assessing social vulnerability, the most frequently used associations of the term are with social processes, economic systems, and relations of power, using variables such as class, ethnicity, gender, age, disability, and health status [25].However, vulnerability is also linked to the lack of access to resources, building stock, and age, type, and density of infrastructure and lifelines [1].The approach of measuring social vulnerability through the use of such variables has certain advantages, but it also has weak points [26].First, by taking into account the demographic characteristics of a society, some authors have argued that women and the elderly are among the most vulnerable, yet information about flood-caused deaths reveals that young and middle-aged men are also vulnerable due to a more risk-taking behaviour [27][28][29].Besides this, not all elderly people are equally vulnerable throughout the entire disaster period [26] as old age itself is relative and depends on several factors (illness, fitness).It has also been proven that gender alone is not an explaining factor of vulnerability, as women's living conditions are greatly influenced by socio-economic status, household structures, and geographic location [30].Education, on the other hand, is a double-faceted variable of social vulnerability.Although well-educated people usually have a higher income and better-positioned properties, the occurrence of natural disasters can potentially result in higher damages.Albeit, the coping capacity of wealthier households is much higher than that of poorer households (having a home insurance can significantly reduce the burden of damages).
The political, environmental and socio-economic context that it is being applied to should be considered.
The most important strength of assessing social vulnerability by using an indicator-and place-based approach is represented by its practical applicability.Thus, once embedded into the wider/national context, it can assess the communities' preparedness for diminishing the negative effects of natural hazards, representing at the same time a valuable input for policy-makers and disaster planners.
In spite of the fact that over the last decades a series of studies have dealt with the issue of social vulnerability in Romania, there is still a strong need to improve the methodology in order to increase its reliability in the context of disaster risk reduction.The most relevant research papers belong to Stângă and Grozavu [31], Armaş and Gavriş [32], Bănică and Muntele [33], Constantin et al. [34], although it is worth mentioning that these studies focus mainly on smaller spatial scales, applying comparative case studies for various parts of regions or for different territorial administrative units (i.e., the Municipality of Bucharest, rural mining settlements in the Apuseni Mountains, Iaşi and Bacău cities and metropolitan areas).Török [35] has applied a modified version of the original SoVI in the Romanian context, proving the workability of the algorithm in a different context of development, using place-specific indicators.
The first and most important result of the study, following the hazard-of-place approach, is the elaboration of place-based indicators along with a specific weighting scheme for measuring local social vulnerability in regions where floods and flood-related hazards are important causes of disasters.The indicators can be well incorporated into other vulnerability models, being relevant and most useful in the vulnerability assessment of other regions characterized by demographic ageing and low to average levels of income as well.Nevertheless, different case studies across countries and regions along with methodological advancements can further enhance the international applicability of the model and the comparison methods of social vulnerability across different areas.Second, analysing social vulnerability across different territories gives emergency managers a meaningful and practical overview of the state of the art, offering the possibility of a much more precise planning and management for emergency situations.The present study found that social status is the most important indicator of vulnerability to disasters in Romanian settlements.Therefore, emergency managers should concentrate their efforts on offering support, advice, and counselling to those specific segments of the population which lack access to basic information on emergency situations and information on emergency management.Third, the present study also reveals the importance of integrating social vulnerability and disaster risk mitigation with the concept of sustainable development.How people live and how society develops has a direct effect on the overall value of potential losses due to disasters.Therefore, the specific results of the study can serve as a starting point for central and local public authorities for drafting and including measures for mitigating existing disaster risks in future development plans and policies, while also striving to develop sustainable communities.
Therefore, the main purpose of the paper is to quantify social vulnerability to natural hazards in Romanian settlements by using a large set of indicators, by applying a new weighting method and by elaborating the Local Social Vulnerability Index (LoSoVI).The results will be used to explore and analyse the spatial relationship between flood hazards and communities with a high degree of social vulnerability, using geographically weighted regression, comparing the results of OLS and GWR modelling at the same time.

The Study Area
With 19.9 million inhabitants according to the latest available census, Romania is one of the largest countries in the EU.
Regarding the territorial-administrative units, the country is organized as follows: counties, cities, towns, communes, and villages.To these, we can add the eight development regions which gained importance after the EU accession, mainly in relation to aspects related to regional development projects and EU funding, although they have neither an administrative status nor a legislative or executive council or government.The names and abbreviations of the counties are listed in the Abbreviations Section.
Hydrographically speaking, Romania is divided into 11 hydrographic units (Figure 1).If we consider the presence of natural hazards, it can be stated that Romania is a country with a high seismic activity and frequent occurrences of extreme phenomena.Romania is also known as one of the most flood-prone countries in Europe with major economic and social consequences, floods usually taking place along the course of internal rivers as well as the Danube [36], with a higher incidence in the north-eastern and north-western parts of the country (Figure 1).

Selection of Variables
While vulnerability is a scale-dependent variable where global, regional, and national scale assessments allow for comparisons across nations, the local scale is more important in disaster risk reduction planning [37].Having this in mind, the analysis of social vulnerability was based on the lowest administrative level which in our case is represented by 3181 study units (including 2861 communes and 320 cities).Statistical data obtained from the National Institute of Statistics, including the 2011 Census, as well as the Tempo Online web database, enabled us to use a large set of variables [38].
The most challenging task in constructing a social vulnerability index is the choice of the major indicators, as the selection depends on both the quality of available variables as well as on the subjectivity of decisions [39].Therefore, the selected indicators used in this study greatly differ from the original SoVI ® model, as the approach was adapted to the Romanian context, taking into account their relevance to the country's situation and reality.The use of different indicators could also be observed in Norway [17], China [18], Portugal [40], Brazil [41], India [42], Vietnam [43], etc.When choosing the indicators, the author has also considered the context and framework drawn up by the main purpose of the study, namely to identify the most vulnerable areas and to reveal the spatial relationship which exists between the flooding and vulnerable communities.After a deep review of the related literature, a total of 27 variables were selected, representing the most relevant factors which have already shown to influence social vulnerability to flood hazards.The variables were then grouped into three major categories, in a completely different manner compared to previous researches [35], comprising also data related to the built infrastructure.Table 1 gives a list of the 27 variables used to develop the social vulnerability index for the settlements, together with their descriptive statistics based on 3181 cases.

Methods
Following the hazard-of-place model approach [15], the methodological framework for assessing the social vulnerability index to flood hazards has included the following steps: first, the socio-economic and built infrastructure related variables were standardized, producing variables with a zero mean and a standard deviation of one.Then, a Principal Component Analysis (PCA) was used to create vulnerability indices from the standardized variables.The PCA technique compresses an original set of variables into a smaller set by identifying patterns in high-dimensional data and revealing the principal factors which best describe variations in the data by finding and clustering variables that measure the same theme [44].While constructing the PCA analysis we have also paid attention to the collinearity analysis and the elimination of redundant data.To check the robustness of the model, two statistical tests-the Kaiser-Meyer-Olkin (KMO) for sampling adequacy and Bartlett's Test of Sphericity-were used.KMO values above 0.6 indicate an acceptable level, while values above 0.8 denote a good compatibility level of variables [45].Bartlett's test of Sphericity verifies whether the correlation matrix is an Identity Matrix, where all the diagonals are 1 and the off-diagonals are 0. If this is the case, none of the variables in the respective matrix are correlated.Finally, Varimax rotation with Kaiser Normalization was used to make the data orthogonal by increasing the interpretability of the factors.This technique minimizes the number of variables which have high loading on a single factor and clarifies the interpretation of the factors.We reduce the number of dimensions by working with only those provisional indexes whose eigenvalues from the PCA are larger than 1.0, provided that we could clearly identify each component.Therefore, components with absolute loading values below 0.5 have been left out.Following the factor rotation stage, a directional adjustment process is applied to the entire factor in order to ensure that the individual component variables all act in the same direction, increasing or decreasing the vulnerability.Thus components enhancing vulnerability were considered positive, while those which reduced it were considered negative.As a consequence, we assign a positive score when the resulting factor increases the total vulnerability, and a negative score when it decreases it.In order to get a more comprehensive picture, we have applied the PCA analysis separately to the built-in, demographic, as well as to the socio-economic dimensions.Figure 2 shows the flowchart guiding the analytical process.The most important/methodological development is represented by the elaboration of an adequate weighting scheme, as not all factors contribute to the overall vulnerability assessment in the same proportion.Since there is no commonly accepted framework defining such a weighting method, several authors used different approaches to assign weights to the specific indices.A brief overview of these can be found in Table 2.As we can observe from the explanation of the cumulative variance provided by the PCA analysis (Table 3), each factor score was assigned according to the number of variables related to each component (e.g., for F 1 , v = 5; for F 2 , v = 2) and multiplied by the weight assigned for each specific factor.The weight was calculated using the total variance explained and the percentage variance explained for each component.
In the first step, we have applied the PCA analysis separately for the built environment resulting in the Built Environment Vulnerability Index (BEVI), separately for the demographic structure forming the Demographic Vulnerability Index (DeVI) as well as for the socio-economic dimension (SEVI), using the following equation.
where, n = the number of territorial units, F = the resulting number of factors, v = the number of variables included in each factor, w = the assigned weight for each factor.Therefore, in the last step, the LoSoVI for all Romanian settlements was calculated with the following equation: The cardinality of each factor was determined in such a way that positive values would increase social vulnerability, while negative values would decrease the overall value.
In order to detect spatial differences in the levels of social vulnerability across the studied area, the LoSoVI scores for each territorial unit were mapped in ArcGIS based on five categories, according to the standard deviation from the mean.Negative values represent low vulnerability while positive values denote a high degree of vulnerability (Table 4).

Standard Deviation
Level of Vulnerability >1.5σ Very high 0.5 to 1.5 High −0.5 to 0.5 Moderate −1.5 to −0.5 Low <−1. 5 Very low After performing the PCA analysis, a model of Geographically Weighted Regression (GWR) was used to explore the spatial relationship between communities directly affected by floods and localities with high and very high LoSoVI values.GWR is a local spatial statistical technique that assumes non-stationarity within relationships-i.e., that is the relationship between the dependent and the explanatory variables changes from location to location [49].Thus the regression coefficient β k takes different values for each location.The equation for the regression model is described below.
where i denotes the coordinates of the points in space and W(i) is a matrix of weights specific to location i constructed in such a way, that observations nearer to i are given a greater weight than others.β represents the vector of global parameters to be estimated, y is a vector of observations on the dependent variable and X is a matrix of independent variables.GWR can be used for uncovering spatial relationships or associations overlooked by OLS.If the relationship does not vary across space, the global model is an appropriate tool for processing the data.Usually, the Akaike's Information Criteria (AICc) is used to compare the overall results from a global OLS linear regression model with those from the local GWR model [50].AICc can be defined as follows: where k is the number of estimated parameters in the model (i.e., the number of variables + 1, to include the intercept) and n is the number of observations in the dataset.If the GWR AICc's value is smaller than the OLS, the value is considered to be the best fit for the model.Besides the AICc information, a set of parameter estimates (local coefficients for each explanatory variable) and associated diagnostics (standard errors, Cook's D statistics, local R 2 statistics) are the major outputs which can be visualized with the help of a GIS platform.By using the GWR4 software to perform the analysis, we have adopted a split-GWR framework using data on just the high and very high social vulnerability communities.This approach has also been adopted by Shoff and Yang [51] while analysing the spatially varying predictors of teenage birth rates in metropolitan and non-metropolitan areas in the USA.GWR was originally developed for spatial point data analysis but it can also be applied to areal data by using the geographical centroids of the analysed areas [51].One of the method's most important strengths is that it allows for the interpolation of values belonging to spatial points not included in the dataset [52].During the computerization, the adaptive bi-square kernel method was applied for geographical weighting in order to estimate local coefficients and a bandwidth size as the observation points of the studied region consist of irregular distances [53].The distinctive feature of this method is that its bandwidth can adapt to the number and density of data points, not requiring that the spatial data points are from a contiguous area, so the spatially continuous character of the observations is not a prerequisite for estimating GWR anymore.The golden-section search was also applied to automatically search for the optimal bandwidth size and the correction of the Akaike Information Criterion (AICc) allowed for an adaptation for selecting the optimum bandwidth [50].Finally, the GWR results were compared to those of the OLS in order to see which model fits the current analysis better.

Measuring Social Vulnerability
By applying the Principal Component Analysis, over a set of 27 variables and a total of 3181 settlements, we have managed to uncover eight latent factors grouped into three categories: Built Environment Vulnerability Index (BEVI), including housing facilities, quality of living, and quality of housing; Demographic Vulnerability Index (DeVI), which incorporates the family structure, demographic vitality and ethnicity, gender, and mobility; and the Socio-Economic Vulnerability Index (SEVI) encompassing education and accessibility, as well as income and employment.The three dimensions of vulnerability have explained 71.7%; 73.2% respectively 71% of the variance within a middling value for the Kaiser-Meyer-Olkin (KMO) Measure of Sampling Adequacy.On the other hand, Bartlett's test also indicates a high probability of significant relationships existing between the analysed variables (p < 0.001) (Table 3).
By analysing the values of the resulted vulnerability scores it must be noted that 34.2% of settlements are characterized by low and very low vulnerability (−1.5 to 0.5 Std.Dev.) while 28.4% fall into categories labelled with high and very high degrees of vulnerability (0.5-1.5 Std.Dev.).We have obtained similar results by adopting the original methodology proposed by Cutter et al. [1], regarding the identification of the least vulnerable and most vulnerable settlements (Figure 3a,b).However, since in the current analysis we have also introduced factors regarding the built environment, specific areas falling into these two extreme categories might differ slightly.According to the factor scores, the most vulnerable areas are to be found-besides the eastern parts of the country-in the southwestern and north-western areas.Within these territories, social vulnerability is mainly impacted by the high proportion of elderly people (with a higher share in the south-western and western parts), by the limited access to resources, and lack of infrastructure.To this, we can add the presence of a high illiteracy rate and a low level of education but also the fact that a relatively high share of the population is working in (subsistence) agriculture-an economic sector highly sensitive to environmental conditions.Low wages and high levels of poverty are just contributing to the social vulnerability of people living in these areas.The only exceptions to the rule are medium-sized cities and suburban localities where job opportunities, the quality of life, and educational opportunities are more favourable, thus lowering the chance of people becoming socially vulnerable (Figure 3b).Looking into the factor scores, we can see that the localities most in need of assistance are Bărbules , ti from Ialomit , a county, Pardos , i and Chiliile from Buzău county, Bulzes , tii de Sus from Hunedoara county, along with Ohaba and Bunila from Alba county.The least vulnerable localities, on the other hand, are concentrated in the northern and central parts of the country, but we could also find some homogenous areas in the south-eastern and western parts, sharing some key characteristics.First of all, the population in these areas has a relatively good access to major public roads and the railway system, disposes of favourable housing facilities and has a relatively high quality of living, in many cases well above the national average.Here we can also add the high demographic vitality (high share of children under 15 years) and the high degree of mobility characterizing mainly the eastern and central parts of the country.According to the factor scores, the least vulnerable areas are Vicovu de Sus, the largest commune from Suceava county, declared a town by Law 83 in 2004, characterized by a dynamic development in recent years; Sângeorz-Băi from Bistrit , a-Năsăud county, a well-known balneo-climateric resort; Cristian from Sibiu county, located in the suburban area of Sibiu Municipality.It must be mentioned that all of the communes characterized by a low SoVI ® index [36] can also be found in this category, like Mos , nit , a Nouă and Remetea Mare from Timis , county; Otopeni, Corbeanca and Mogos , oaia from Ilfov county as well as Gârbău from Cluj county-all of the parts of the greater suburban areas of Timis , oara, Bucharest and Cluj-Napoca urban agglomerations (Figure 3b).These settlements represent 3.9% of the study units (124 communes and municipalities).It is worth mentioning that the level of economic development is not in all cases inversely proportional with the vulnerability to natural hazards since a high level of development can greatly increase the level of internal inequalities, meaning that highly developed urban agglomerations are not necessarily areas with a low level of social vulnerability [54].This can be well observed in large urban areas (of more than 200,000 inhabitants) where the high population density and the concentration of economic activities increase the vulnerability of the built environment.Furthermore, the existence of wealth-poverty dichotomies, social inequalities, and social segregation within cities can amplify their exposure to social risks [36].According to the LoSoVI index, especially Timis , oara, Craiova, Brăila and the capital-Bucharest-are labelled with high values of social vulnerability, while others like Constant , a, Cluj-Napoca and Iaşi belong to intermediate categories, all of them being surrounded by areas with a low level of social vulnerability.As we could observe, some of the most important factors influencing social vulnerability are the levels of economic development, education, and occupation status, but also housing facilities, although these aspects are not only valid for Romania.By integrating these indicators into the vulnerability assessment model, similar results have been achieved by Holand et al. [18] for the municipalities of Norway, by Guillard et al. [41] for Portugal, by Hummell et al. [42] for Brazil, to name just a few in a wider context.

Geographically Weighted Regression Results
The GWR analysis was performed using the GWR4 software.For the dependent variable we have used the number of people affected by flood hazards, the generated eight (especially the high and very high) LoSoVI factors acting as independent variables.The Pearson correlation revealed that a total of five out of eight variables with a p-value less than 0.05 show a significant correlation with frequently flooded communities: housing facilities, quality of housing, family structure, education and accessibility, income, and employment.In order to examine the multicollinearity between the selected variables, the variance inflation factor (VIF) was checked.If the VIF is less than 10, it can be stated that no multicollinearity exists, which means that it will not influence the stability of the parameter estimates significantly [55].As the VIF scores range between 1.348 and 3.005, we used all five LoSoVI scores for the regression model (Table 5).The OLS global model explained about 43% (adjusted R 2 = 0.43) of the variation in the numbers of the population directly affected by flood hazards, with an AICc value of 9139.28.At the same time, the calibrated results obtained by calculating the GWR have shown a significant improvement of the global model.By comparing the AICc values for the two models, we can observe a significant reduction from 9139.28 to 8816.42, the difference of 82.48 implying a strong improvement in adequacy and proving at the same time that the GWR model is much more suitable for analysing the respective data.Since the AICc is a measure of spatial collinearity within the dataset, the lower the value of AICc, the better the model will fit the observed data.Moreover, the local R 2 value of 0.78 also suggests an improvement in the performance of the model (Table 6).Thus, the explanatory power of the local model (GWR) is higher than that of the OLS model.Figure 4 shows the distribution of Local R 2 values (Figure 4a) as well as the territorial distribution of the standard deviation of residuals (Figure 4b).In the case of the Local R 2 , the higher the value, the better the performance of the local model.In our case, 68% of the analysed communities have relatively high R 2 values (between 0.47 and 0.89), which explains why communities characterized by high and very high LoSoVI scores are also the most prone to suffer from flooding.On the other hand, the lower the standard deviation of residuals, the better the GWR model fits the observed data, which, in the case of the analysed regions, explains the relationship between flood hazard and socio-economic aspects.We have found that the relationship between communities directly affected by floods and localities with high and very high LoSoVI values can be investigated over multiple spatial scales by incorporating geographic information into GWR analyses.The improved performance of the GWR over the global regression model was demonstrated by significantly lower values both for the AICc as well as deviance.An additional strength of the GWR method was the ability to explore spatial variability in the relationship between the population directly affected by flooding and explanatory variables like housing facilities and quality of housing, family structure, education and accessibility as well as income and employment.The value of local R 2 varies between 0.01 and 0.89, indicating that the GWR model fits the observed y value well.As we have mentioned, the higher the value of the indicator, the better the model fits the data, while lower values indicate that the fit of the model is not unequivocal.
With a few exceptions, most of the settlements included in the analysis have a high explanatory power (R 2 = 0.47, denoting a value above 47%), while the low values of the residuals representing the difference between the observed and estimated values (Figure 4b) denote a proper fit of the GWR model to the observed data.

Discussion
Addressing social vulnerability to disaster risks requires the understanding that the social context itself can increase both vulnerability (socio-economic variables), as well as the building of capacities to deal effectively with flood hazards [56].Eliminating variables which lead to vulnerability (such as a low level of education, high poverty rate and low accessibility to basic infrastructure as well as the imbalanced family structure) and capacity building which include those activities that help people mitigate the impact of a disaster, should represent the most important policy implications.Disaster risk prevention and reduction are crucial to enhance the socio-economic resilience of communities.In the context of climate change, where floods are expected to become more frequent, flood risk mitigation is a major concern of water policies and sustainable development strategies at different spatial scales.
Managing climate change-related risks requires decisions and strategies for adaptation, with implications for future generations, economies, and environments [57].The most important national document regulating flood risk management in Romania is the National Strategy for Flood Risk Management for the medium and long-term, which was adopted in 2010.This strategy defines the technical, institutional and legal framework for diminishing the negative consequences of floods on socio-economic activities, life and population health, as well as on the environment, for the 2010-2035 period [58].Romania has a central water management authority, the Romanian Waters National Administration (RWNA) which is responsible for the implementation of the national water management strategy and policies, aspects related to quantitative and qualitative water management as well as the operation of the different management structures.This Authority has 11 regional branches (see Figure 1) organized according to the hydrographic basins of Romania, each with its own Flood Risk Management Plan.These plans address all aspects of flood risk management focusing on prevention, protection, preparedness, including flood forecasts and early warning systems, considering the characteristics of the particular river basin or sub-basin.Flood risk management plans also include actions to promote sustainable land use practices, the improvement of water retention as well as the controlled flooding of certain areas in the case of a flood event [59].
As it could be observed, government agencies-at least on paper-have elaborated and implemented all necessary measures associated with flood risks (prevention).Related to our analysis, the central question is related to policy or planning mechanisms aiming to cope with existing social vulnerabilities, increasing public participation, public relations/communication, and education on flood risk.According to the results of local surveys conducted in Niraj Walley, local authorities are generally aware of their settlements being exposed to a flood event, but usually have no information on multiple problems involved which can be related to social and economic aspects of the built infrastructure, as they have been revealed by the present analysis.As this is the first attempt ever to measure social vulnerability for all settlements in Romania, inserting the main findings into the revised National Strategy for Flood Risk Management, followed by an update of the Flood Risk Management Plans for all the 11 hydrographic districts would have a great added value in increasing the effectiveness of the mentioned documents in diminishing the impact of future floods.Even so, the most important activities for flood risk prevention depend on the affected population, as the responsibilities of local authorities are mainly related to transmitting relevant information in due time and to as many people as possible.With this regard, it must be mentioned that Romania has introduced a compulsory insurance system enforced by law no.260/2008 and law 191/2015, making it compulsory for all Romanian real-estate owners to purchase multi-hazard insurance.The amount of insurance depends on the quality of the building-mostly reflected by the construction type-and not on the probability of exposure to the specific hazard.The mentioned legislation also states that no compensation will be paid from the State Budget to uninsured households.Even though this provision is compulsory for all Romanian real-estate owners, according to the latest available statistics (31 August 2018), the rate of ensured households has barely reached 20%, with higher rates in urban than in rural areas [60].The low penetration in rural areas is mostly due to a lack of trust regarding insurers, missing information on how to buy insurance, difficulties proving the ownership of properties but also due to the self-assurance of people believing that such a hazard event could not happen to them.According to the survey performed in July 2015, when 741 households were queried in the Niraj Walley (Central Romania) about their perception on natural hazards, only 8.1 percent of the respondents mentioned that they are prepared for all hazardous events, 36.6 percent claimed they can only take the most basic measures, 43.2 percent considered themselves unable to take any action, while 8.5 percent believed that their household could not suffer any hazardous event.These values are all the more surprising when we consider that most settlements have strategic documents elaborated like Defence Plans in case of emergencies, including information on the Local Committee for Emergency Situations, schemes for transmitting information in case of emergency situations at the local scale, defence measures against floods, limits of flooded areas, etc.This makes us realize that it is not enough to publish or disseminate such information on the City Hall website or on posters with textual information or maps regarding flood risk in buildings of the local council, since in this way information does not reach the elderly or people from isolated households which are in the greatest need of help.In this case, the most important task should be overtaken not only by local authorities, but also by other public figures like religious leaders or teachers, by informing people and raising awareness regarding the flood risk.According to the results of the above-mentioned survey, we consider that the best solutions are actions which are held in schools in order to raise awareness and educate the children at an early age.In the case of the elderly, personal contact and giving help in increasing the level of preparedness are of utmost importance.Flood risk adaptation is key for the existence of people in certain areas, enabling them to survive in extreme conditions, such as floods.Therefore, exercises and simulations of flood crises should be organised at regular time intervals in order to train the population at risk.As it could be observed, adaptation to flood risk can be achieved by various hard and soft measures, all contributing to flood risk mitigation.At the same time, LoSoVI can provide the basis for discussing the main causes of risk and vulnerability, mitigation measures, preparedness, response planning as well as recovery.

Conclusions
In recent years, there has been an increased interest towards initiatives aiming to reduce the impact of hazards, not just from a physical perspective, but also from a social perspective.Despite the absence of a generally accepted definition and an agreed-upon methodology for the measurement of social vulnerability, we can conclude that an important number of initiatives providing different methodological approaches and indicators offer a valuable contribution to its assessment.It is a universally accepted fact that there is neither a single common approach nor a universal catalogue of vulnerability indicators: vulnerability is highly context-specific in terms of the socio-economic, demographic, political, and cultural contexts.Furthermore, considering the temporal and spatial variations of social vulnerability, a one-size-fits-all approach for increasing preparedness, improving response, and facilitating recovery may be the least effective for reducing it [61].
This paper was aimed to provide an assessment of social vulnerability in Romanian settlements by applying the SoVI method in a local context.Using GIS-based local models and global statistics to explore the relationship between populations exposed to flood hazards and settlements characterized by high vulnerability scores, we were able to detect and extract certain key information pieces concerning non-stationarity in spatial data.By comparing the results regarding the adequacy of global OLS and local GWR models with certain parameter estimates, this paper has managed to explicitly demonstrate that within spatial data, the relationship is non-static across geographic space.The results from our analysis complement the main findings and existing debates in contemporary literature in many ways.
First, the conducted research represents the first in-depth attempt to quantify community-level vulnerability to flood hazards in Romania.The path dependency of some localities, the specific historical, social and economic factors have put their mark on the evolution of social vulnerability, especially on variables like the ageing of the population, occupation and income levels.These factors also represent the main origin of socio-spatial inequalities which have now been reflected in the territorial distribution of the LoSoVI model.In general, a high level of vulnerability associated with a low level of resilience characterizes peripheral rural areas, which are further affected by natural hazards, especially floods (eastern and north-western parts of rural Romania).What is more, due to the high rate of poverty and limited access to resources and services, these areas are often associated with low levels of education and lack of sanitation facilities, decreasing the ability to cope with the impact of natural hazards.In this sense, their present status may very well give serious reason for concern.As the study has confirmed, even the most developed urban areas can be vulnerable in the face of natural disasters (generally due to the high population density), indicating that large urban agglomerations do not necessarily have low levels of social vulnerability and that economic development is not always inversely proportional to vulnerability, although as mentioned earlier, the respective communities have much better capabilities to cope with-, to resist-, and to recover from losses.
Second, the developed GWR model can be used to carry out predictions for flood hazard damages among highly vulnerable population groups in Romania.
Third, from the point of view of sustainable spatial planning, the present study is particularly important because by capturing the spatial characteristics of social vulnerability, the above in-depth analysis may offer a viable and integrative practical guide to authorities and policy-makers for taking the necessary measures in order to diminish social vulnerability in specific areas.Thus, the analysis of vulnerability at the local level enables the identification of settlements in need of financial resources in order to diminish the negative impact of natural hazards.By assessing the social vulnerability of a certain area we can determine its unique position in relation to flood hazards.This, in turn, can serve as the starting point for local planners and decision-makers in identifying weak points to be included into future development strategies, thus, improving the communities' own ability to cope with the effects of flood hazards.Further on, the collaboration between researchers of flood hazards and local practitioners regarding the effective vulnerability assessment can contribute not only to the sharing of knowledge, but can also help facilitate the insertion of specific local needs into local policy decisions.
Although the analysis provides the possibility to explore fundamental social-, environmentaland economic phenomena using an innovative approach, due to the limited access to reliable and up-to-date economic and environmental data at the local level, it is still not able to provide us with an overall assessment on social vulnerability.At the same time, the final results of such analyses are also very much dependent not only on the used variables, selected methods and weighting schemes, but also on the different interpretations of the concept of vulnerability.The present analysis has managed to capture some of the social aspects which can be considered essential building blocks in the overall assessment of physical vulnerability, therefore-as Holand et al. mentioned [18]-differing results should be seen as complementarities rather than contradictions in applying the different models.

Figure 1 .
Figure 1.The administrative units affected by floods.

Figure 2 .
Figure 2. The framework of the analysis.
Note: The symbol indicates whether the used variable increases or decreases social vulnerability (+ = increases vulnerability, − = decreases social vulnerability).

Figure 3 .
Figure 3.The comparison between the Social Vulnerability Index scores in Romania calculated without weighting (a) and using a specific weighting method (b) for 2011.

Figure 4 .
Figure 4. (a) The spatial distributions of the local R 2 and (b) the standard deviation of the residuals.

Table 1 .
Descriptive statistics of variables considered for social vulnerability assessment.

Table 2 .
The relevant studies on measuring social vulnerability.

Table 3 .
The significant components and their adjusted loadings used to assess local social vulnerability in Romania.

Table 4 .
The vulnerability level classification based on the standard deviation.

Table 5 .
The selection of variables using correlation and the multicollinearity test.

Table 6 .
The results of the Ordinary Least Square (OLS) and Geographically Weighted Regression (GWR) analysis models.