Natural Hazards and Job Choice: Factors driving job choice and willingness to relocate

We examined all factors related to the job characteristics to understand the drivers of the final job choices by participants. The factors included the salary, the crime risk and hazard risk of the job’s location, the change in hazard risk between where the participant was currently living and the job’s location, the geographic distance (in terms of US Census divisions) the participant would move, and the change in base cost of living between the participant’s current location and job location.

Salary was the strongest driver of job choice, showing the highest statistical significance (p < 0.001) and a positive association, indicating that higher-paying jobs were more likely to be chosen. Crime risk was also a significant factor (p < 0.001), with a negative association, suggesting that jobs located in areas with higher crime risk were less likely to be selected. Hazard risk showed a marginally significant negative effect (p = 0.051), implying that higher hazard risk may slightly reduce the likelihood of job selection, though this effect was weaker than salary or crime risk. Relocation factors, such as the distance between the job and the respondent’s current location (measured by US Census division change), and the change in hazard risk, did not significantly influence final job choice. Similarly, the cost of living, both as an absolute measure and as a change relative to current location, had no meaningful impact on the final decision.

When comparing geoscientists and non-geoscientists separately, salary remained the strongest and most consistent predictor across both cohorts. Crime risk also had a consistently negative and significant effect in both cohorts. Hazard risk was marginally significant only for geoscientists, suggesting that they may be more sensitive to choosing a job located in higher hazard risk areas. No other factors significantly influenced job selection in either group.

Whole Population — Logit Regression Summary

Dep. Variable:finaljob_chosen No. Observations:6,770
Model:Logit Df Residuals:6,762
Method:MLE Df Model:7
Date:Tue, 15 Apr 2025 Pseudo R-squared:0.07682
Time:15:39:52 Log-Likelihood:-1301.2
Converged:True LL-Null:-1409.4
Covariance Type:nonrobust LLR p-value:3.573e-43
variable coef. std. err. z P > |z| [0.025 0.975]
const-3.23540.468-6.9110.000-4.153-2.318
salary1.425e-051.21e-0611.7780.0001.19e-051.66e-05
hazard_risk-0.04910.025-1.9530.051-0.0980.000
crime_risk-0.31400.068-4.5860.000-0.448-0.180
cost_of_living1.053e-059.06e-061.1620.245-7.23e-062.83e-05
division_change0.02320.0380.6080.543-0.0520.098
hazard_risk_change0.00830.0180.4710.637-0.0260.043
cost_of_living_change1.098e-066.32e-060.1740.862-1.13e-051.35e-05
Job-related variables: salary, hazard_risk (cumulative), crime_risk, cost_of_living
Change variables: division_change (number of US Census divisions between current location and job location), hazard_risk_change (difference between cumulative hazard risk at job location and current location), cost_of_living_change (difference between cost of living at job location and current location)

Geoscientists Only — Logit Regression Summary

Dep. Variable:finaljob_chosen No. Observations:3,789
Model:Logit Df Residuals:3,781
Method:MLE Df Model:7
Date:Tue, 15 Apr 2025 Pseudo R-squared:0.09177
Time:16:08:04 Log-Likelihood:-729.33
Converged:True LL-Null:-803.02
Covariance Type:nonrobust LLR p-value:1.440e-28
variable coef. std. err. z P > |z| [0.025 0.975]
const-3.33440.626-5.3250.000-4.562-2.107
salary1.438e-051.49e-069.6270.0001.15e-051.73e-05
hazard_risk-0.06430.033-1.9530.051-0.1290.000
crime_risk-0.31270.091-3.4390.001-0.491-0.134
cost_of_living1.576e-051.20e-051.3170.188-7.7e-063.92e-05
division_change0.06130.0511.2120.226-0.0380.160
hazard_risk_change0.01280.0230.5500.583-0.0330.058
cost_of_living_change8.22e-068.43e-060.4530.650-1.27e-052.03e-05
Job-related variables: salary, hazard_risk (cumulative), crime_risk, cost_of_living
Change variables: division_change (number of US Census divisions between current location and job location), hazard_risk_change (difference between cumulative hazard risk at job location and current location), cost_of_living_change (difference between cost of living at job location and current location)

Non-Geoscientists Only — Logit Regression Summary

Dep. Variable:finaljob_chosen No. Observations:2,981
Model:Logit Df Residuals:2,973
Method:MLE Df Model:7
Date:Tue, 15 Apr 2025 Pseudo R-squared:0.05999
Time:16:08:04 Log-Likelihood:-569.89
Converged:True LL-Null:-606.26
Covariance Type:nonrobust LLR p-value:4.115e-13
variable coef. std. err. z P > |z| [0.025 0.975]
const-3.11210.710-4.3850.000-4.503-1.721
salary1.39e-052.10e-066.6090.0009.78e-061.80e-05
hazard_risk-0.03360.040-0.8420.400-0.1120.045
crime_risk-0.32730.105-3.1320.002-0.532-0.122
cost_of_living5.706e-061.44e-050.3970.692-2.25e-053.39e-05
division_change-0.01840.058-0.3150.752-0.1330.096
hazard_risk_change0.00420.0280.1500.881-0.0500.059
cost_of_living_change-2.52e-069.75e-06-0.2580.796-2.16e-051.66e-05
Job-related variables: salary, hazard_risk (cumulative), crime_risk, cost_of_living
Change variables: division_change (number of US Census divisions between current location and job location), hazard_risk_change (difference between cumulative hazard risk at job location and current location), cost_of_living_change (difference between cost of living at job location and current location)

Is Job Choice Behavior Different Between Cohorts?

To determine whether final job choice behavior differed between geoscientists and non-geoscientists, we conducted a series of statistical comparisons using three approaches: coefficient-wise z-tests, interaction term testing, and a likelihood ratio test (LRT).

First, we ran separate logistic regression models for each subgroup and compared the resulting coefficients using z-tests. None of the differences in coefficients across key predictors such as salary, hazard risk, crime risk, and cost of living were statistically significant (all p > 0.30).

Coefficient Comparison — Geoscientists vs. Non-Geoscientists

variable coef. (Geo) SE (Geo) coef. (Non-Geo) SE (Non-Geo) Δ coef. z p-value
const-3.33440.6262-3.11210.7097-0.2223-0.2350.8143
salary0.000010.000000.000010.000000.000000.1880.8505
hazard_risk-0.06440.0330-0.03360.0399-0.0308-0.5940.5523
crime_risk-0.31270.0909-0.32740.10450.01470.1060.9156
cost_of_living0.000020.000010.000010.000010.000010.5370.5911
division_change0.06130.0506-0.01840.05840.07981.0320.3022
hazard_risk_change0.01280.02330.00420.02790.00860.2370.8125
cost_of_living_change0.000000.00001-0.000000.000010.000010.4920.6227
Note: Δ coef. is the difference between geoscientist and non-geoscientist coefficients. p-values reflect the significance of the coefficient difference.

Next, we estimated a full interaction model including geoscientist status (geo) and interaction terms between group membership and each predictor. A Wald test on each interaction term revealed no significant group moderation effects (all p > 0.30), indicating that the influence of predictors on job choice did not significantly differ by group.

Interaction Term Tests — Geoscientist × Variable

Interaction Term χ² p-value df constraint
salary × geo0.03550.85051
hazard_risk × geo0.35330.55231
crime_risk × geo0.01120.91561
cost_of_living × geo0.28870.59111
division_change × geo1.06450.30221
hazard_risk_change × geo0.05630.81251
cost_of_living_change × geo0.24200.62271
Note: Results reflect tests of interaction terms between each predictor and the Geoscientist group. χ² is from a nested model comparison; df = 1.

Finally, we conducted a likelihood ratio test comparing the full interaction model to a simpler pooled model with no interaction terms. The result (LR statistic = 3.76, df = 7, p = 0.81) confirmed that the interaction model did not provide a significantly better fit to the data.

Taken together, these results demonstrate that geoscientists and non-geoscientists respond similarly to the job characteristics examined. Therefore, for all subsequent analyses, we use a single pooled model including the full sample. This approach preserves statistical power and simplifies interpretation, without compromising accuracy.

Does Disposable Income Factor into Job Choice?

We calculated disposable income for participant’s final job choice as final job salary minus the cost of living for the location of the job to understand how strongly purchasing power factored into job choice. Disposable income was a strong positive predictor of job selection (p < 0.001), indicating that jobs offering higher earnings relative to cost of living were more likely to be chosen.

Whole Population — Logit Regression

Dep. Variable:finaljob_chosen No. Observations:6,770
Model:Logit Df Residuals:6,763
Method:MLE Df Model:6
Date:Tue, 15 Apr 2025 Pseudo R-squared:0.07944
Time:21:49:22 Log-Likelihood:-1297.5
Converged:True LL-Null:-1409.4
Covariance Type:nonrobust LLR p-value:1.510e-45
variable coef. std. err. z P > |z| [0.025 0.975]
const-4.18530.465-9.0040.000-5.096-3.274
hazard_risk-0.04340.020-2.1460.032-0.083-0.004
crime_risk-0.30060.069-4.3750.000-0.435-0.166
division_change0.01560.0380.4110.681-0.0590.090
disposable_income3.184e-056.42e-064.9630.0001.93e-054.44e-05
cost_of_living4.101e-058.62e-064.7580.0002.41e-055.79e-05
disposable_income × cost_of_living-3.228e-101.16e-10-2.7750.006-5.51e-10-9.48e-11

Does Willingness to Relocate Influence Job Choice?

We examined whether a participant’s willingness to relocate influenced job choice behavior. The willingness to move variable represented individuals preferred geographic job search scope (i.e., only within their current US Census division, both within and outside of their current US Census division, only outside of their current US Census division). The coefficient for willingness to move was negative but not statistically significant, indicating that individuals’ openness to relocating did not significantly affect the likelihood of choosing a particular job once job characteristics were accounted for.

Whole Population — Logit Regression

Dep. Variable:finaljob_chosen No. Observations:6,770
Model:Logit Df Residuals:6,764
Method:MLE Df Model:5
Date:Tue, 15 Apr 2025 Pseudo R-squared:0.07706
Time:22:15:06 Log-Likelihood:-1300.8
Converged:True LL-Null:-1409.4
Covariance Type:nonrobust LLR p-value:5.845e-45
variable coef. std. err. z P > |z| [0.025 0.975]
const-3.30080.374-8.8170.000-4.035-2.567
willingness_to_move-0.09850.085-1.1630.245-0.2640.068
salary1.422e-051.21e-0611.7410.0001.18e-051.66e-05
cost_of_living1.34e-056.98e-061.9190.055-2.89e-072.71e-05
hazard_risk-0.04530.020-2.2330.026-0.085-0.006
crime_risk-0.31790.068-4.6470.000-0.452-0.184

We also examined whether an individual’s willingness to move affected the likelihood that they would choose a job that’s in a different US Census division than their current location of residence (division change). Neither willingness to move nor the interaction with division change were statistically significant, suggesting that individuals’ openness to relocating did not significantly moderate their likelihood of selecting jobs located outside their current region, after accounting for job attributes.

Whole Population — Logit Regression Summary

Dep. Variable:finaljob_chosen No. Observations:6,770
Model:Logit Df Residuals:6,762
Method:MLE Df Model:7
Date:Tue, 15 Apr 2025 Pseudo R-squared:0.07775
Time:22:22:11 Log-Likelihood:-1299.8
Converged:True LL-Null:-1409.4
Covariance Type:nonrobust LLR p-value:9.860e-44
variable coef. std. err. z P > |z| [0.025 0.975]
const-3.31830.384-8.6380.000-4.071-2.565
willingness_to_move-0.09990.123-0.8150.415-0.3400.140
salary1.421e-051.21e-0611.7350.0001.18e-051.66e-05
cost_of_living1.19e-057.09e-061.6800.093-1.98e-062.58e-05
hazard_risk-0.04350.020-2.1400.032-0.083-0.004
crime_risk-0.31010.069-4.5200.000-0.445-0.176
division_change0.10510.0891.1780.239-0.0700.280
wtm × division_change-0.04860.072-0.6700.503-0.1910.093
Note: This model includes an interaction term: wtm × division_change (willingness to move × number of Census divisions moved).

To evaluate whether the effect of salary on job selection varies by willingness to relocate, we included an interaction term between willingness to move and salary in the logistic regression model. The interaction between salary and willingness to move was not statistically, nor was the main effect of willingness to move. These findings suggest that individuals’ stated openness to relocation does not meaningfully alter how salary influences their final job choice. There was no evidence that more mobile individuals were more responsive to higher-paying jobs. As such, stated willingness to move had no independent or moderating effect on job selection once job-level features were accounted for.  

Whole Population — Logit Regression

Dep. Variable:finaljob_chosen No. Observations:6,770
Model:Logit Df Residuals:6,763
Method:MLE Df Model:6
Date:Tue, 15 Apr 2025 Pseudo R-squared:0.07736
Time:22:30:29 Log-Likelihood:-1300.4
Converged:True LL-Null:-1409.4
Covariance Type:nonrobust LLR p-value:2.671e-44
variable coef. std. err. z P > |z| [0.025 0.975]
const-3.44650.407-8.4590.000-4.245-2.648
willingness_to_move0.05020.1820.2760.783-0.3070.407
salary1.579e-052.08e-067.5720.0001.17e-051.99e-05
cost_of_living1.361e-056.99e-061.9470.052-8.91e-082.73e-05
hazard_risk-0.04550.020-2.2430.025-0.085-0.006
crime_risk-0.31920.069-4.6580.000-0.454-0.185
wtm × salary-1.65e-061.79e-06-0.9210.357-5.16e-061.86e-06
Note: This model includes an interaction term: wtm × salary (willingness to move × salary).

Continue reading