This paper is in the following e-collection/theme issue:

Public Health (207)

Published on in Vol 2, No 2 (2013): Jul-Dec

Aftershocks Associated With Impaired Health Caused by the Great East Japan Disaster Among Youth Across Japan: A National Cross-Sectional Survey

Aftershocks Associated With Impaired Health Caused by the Great East Japan Disaster Among Youth Across Japan: A National Cross-Sectional Survey

Aftershocks Associated With Impaired Health Caused by the Great East Japan Disaster Among Youth Across Japan: A National Cross-Sectional Survey

Authors of this article:

Takashi Sugimoto1, 2 ;   Tomohiro Shinozaki3 ;   Yuki Miyamoto1
Article Authors Cited by (6) Tweetations (4) Metrics

    Original Paper

    1Graduate School of Medicine, Department of Psychiatric Nursing, The University of Tokyo, Tokyo, Japan

    2Graduate School of Public Policy, Department of Health Technology Assessment and Public Policy, The University of Tokyo, Tokyo, Japan

    3Graduate School of Medicine, Department of Biostatistics, The University of Tokyo, Tokyo, Japan

    Corresponding Author:

    Takashi Sugimoto, BSc

    Graduate School of Medicine

    Department of Psychiatric Nursing

    The University of Tokyo

    7-3-1 Hongo, Bunkyo

    Tokyo, 113-0033

    Japan

    Phone: 81 3 5841 3522

    Fax:81 3 5841 3392

    Email: sugimotot-tky@umin.ac.jp


    Background: The Great East Japan earthquake, subsequent tsunamis and the Fukushima nuclear incident had a tremendous impact on Japanese society. Although small-scale surveys have been conducted in highly affected areas, few have elucidated the disaster’s effect on health from national perspective, which is necessary to prepare national policy and response.

    Objective: The aim of the present study was to describe prefecture-level health status and investigate associations with number of aftershocks, seismic intensity, a closer geographical location to the Fukushima Nuclear Power Plant, or higher reported radiation dose in each prefecture even after adjusting for individual socioeconomic factors, by utilizing individual-level data acquired from a national cross-sectional Internet survey as well as officially reported prefecture-level data.

    Methods: A Japanese government research institute obtained 12,000 participants by quota sampling and 7335 participants were eligible for the analysis in an age range between 17 and 27 years old. We calculated the percentage of people with decreased subjective health in each prefecture after the earthquake. Variability introduced by a small sample size for some prefectures was smoothed using empirical Bayes estimation with a random-intercept logistic model, with and without demographic factors. Multilevel logistic regression was used to calculate adjusted odds ratios (ORs) for change of subjective health associated with prefecture-level and individual-level factors.

    Results: Adjusted empirical Bayes estimates were higher for respondents commuting in the northeast region (Iwate 14%, Miyagi 19%, and Fukushima 28%), which faces the Pacific Ocean, while the values for Akita (10%) and Yamagata (8%) prefectures, which do not face the Pacific Ocean, were lower than those of Tokyo (12%). The values from the central to the western region were clearly lower. The number of aftershocks was coherently associated with decreased health (OR 1.05 per 100 times, 95% CI 1.04-1.06; P<.001) even after adjusting for covariates (OR 1.02 per 100 times, 95% CI 1.00-1.05; 1.32 per 1000 times, 95% CI 1.03-1.71; P=.049). In contrast, seismic intensity of the initial earthquake (OR 0.87, 95% CI 0.65-1.17; P=.36), radiation dose (OR 1.16, 95% CI 0.82-1.64; P=.41), and distance from the Fukushima Nuclear Power Plant (OR 1.00, 95% CI 0.99-1.00; P=.66) were not. Change in job condition (OR 2.05, 95% CI 1.72-2.45; P<.001), female (OR 1.43, 95% CI 1.19-1.69; P<.001), higher age (OR 1.06 per year, 95% CI 1.02-1.11; P=.005), and duration of evacuation longer than 4 weeks (OR 1.44, 95% CI 1.06-1.97; P=.02) seemed to decrease perceived health status.

    Conclusions: We found nationwide differences that show decreased health status because of the Great East Japan disaster according to prefecture. The number of aftershocks, change in work conditions, being female, a higher age, and duration of the evacuation were risk factors for the population after the major earthquake, tsunamis, and nuclear incident.

    Interact J Med Res 2013;2(2):e31

    doi:10.2196/ijmr.2585

    Keywords

    earthquakes; health communication; Internet; nuclear power plants; risk management 


    The Great East Japan earthquake (Tohoku earthquake) on March 11, 2011 and the subsequent tsunami had a tremendous impact on Japanese society [1]. Furthermore, the Fukushima nuclear accident, ranked 7, which is the most severe on the International Nuclear Event Scale, forced large numbers of people and industries to evacuate the area [2,3]. The explosion at the Fukushima Nuclear Power Plant (NPP) released substantial amounts of radioactive materials into the sea [4] and into the atmosphere, resulting in higher levels of radiation reaching as far as the Tokyo metropolitan area, which is 200-300 km away from the Fukushima NPP [5]. A record number of 3945 aftershocks during March after the initial earthquake and nearly 10,000 aftershocks over the following 2 years have been a stark feature of the disaster [6,7]. The persistent aftershocks and lack of credible information provided by authorities of the national government caused distress among people across Japan [8-10]. Despite such concern about health nationwide, there were no rapid assessments of nationwide health status and implementation of strategic announcements, partly because of the logistical challenges involved, but mainly because of the absence of preparedness and coordination between the central and local governments, medical communities, nongovernmental organizations, and volunteer groups [1]. Assessing public health status promptly across a nation is of major relevance to health policy decision makers as well as researchers looking at disasters.

    Unlike impromptu and unsystematic surveys on health, many younger persons tried to empower others in disaster-struck areas promptly through the use of Internet technologies. Physicians and hospital officials in affected areas reported their medical resource status utilizing email lists or social media such as Twitter, which made a definite difference in the disaster response compared with the Great Hanshin earthquake of 1995 [11]. Information on bed availability and so on was collected and became available to the public within the first 48 hours of the earthquake using Google Map technology [11]. Furthermore, use of geographical information system facilitated deployment of medical teams because it provided data on radiation risks in Fukushima promptly [12]. Many younger persons, who were in their early 40s or younger who were not health care professionals, also took the initiative to gather and diffuse relevant information, garner support, and raise money for quake relief [1,13]. Substantially, many people throughout the nation used the Internet to seek accurate information and to also access the mass media. This was in some ways similar to the public health emergency response in the United States after the 9/11 attacks [1,11-14]. The Internet has now been adopted widely enough in Japanese society, especially in the younger generation cohort, to be a valuable assessment tool in a public health emergency.

    Public health assessments as well as constructive advice via the Internet have unique advantages in terms of a more rapid and broader reach [1,11,13,15,16]. Specifically, surveys via the Internet are able to reach more citizens distant from the devastated area, as well as those in highly affected areas, or those with mild health complaints who have not visited a medical facility [15,16]. Most previous studies regarding the effects of the disaster on health have actually been limited to mostly severely affected areas in northeastern Japan, except for a few studies focusing on distance from the epicenter [8,17] and the effect of aftershocks on psychological stress [8]. Most previous studies limited attention to presumed high-risk populations within a highly affected area. For example, previous studies included those with cardiovascular diseases [18,19] and diabetes mellitus [20] or other diseases [21-23], the elderly [24-26], evacuees [24-26], children [27,28], workers [29-31], nonprofessional volunteers [32], caregivers [33], and pregnant women [34] in a disaster area. The limited attention seems to have been common among disaster studies. One study of the previous Niigata-Chuetsu earthquake in 2004 showed that being female had a higher odds ratio of psychological distress [35]. Furthermore, studies showed the importance of preventing mothers from having symptoms of psychological distress caused by anxiety about the health of their children and separation from family members [36-38]. We believed that the inherent restricted nature of surveys after a disaster in Japan is attributed to the scarce knowledge of potential effects on nationwide health status after a massive disaster. This would complicate a coordinated national response to a disaster.

    Although the elderly or children may be the most vulnerable, a previous study indicated disruption of work after natural disaster as being independently associated with decreases in general mental and physical health among university students [39]. Another recent report revealed a significant regional difference in the perception of risk among Japanese university students in a severely hit region, an indirectly affected Tokyo region, and mostly unaffected western regions [40]. Students in the Tokyo region were anticipated to be at increased risk of a future earthquake compared with those in a victimized region as well as western Japan, although the perceived risk of further nuclear risks was approximately the same for these regions [40]. In the case of the Chernobyl disaster, people living relatively far from the disaster site tended to be more concerned about the political and economic situation [41]. To our knowledge, there are few studies on health nationwide that examine geographical factors as well as social factors after large-scale disasters such as that in Japan. Because the Tohoku disaster may have affected nationwide health in every age range due to its nuclear radiation release, persistent aftershocks, widespread concern, and expected differences in these factors among regions [1-10,17,40,41], we wanted to evaluate its impact nationwide, especially among those in their late teens or 20s who typically are seen as less vulnerable. These data are crucial for national planning, programs and recovery.

    We hypothesized that there are positive associations between the decreased subjective health of the young population nationwide and a larger number of aftershocks, a closer geographical location to the Fukushima NPP, or a higher reported radiation dose in the atmosphere in each region even after adjusting for individual socioeconomic factors. Therefore, the aim of the present study is to describe prefecture-level health status and investigate the associations mentioned above after adjusting for individual socioeconomic factors, utilizing individual-level data acquired from a national cross-sectional survey as well as officially reported prefecture-level data.


    Data Acquisition

    The data for this secondary analysis, an Internet survey on the effects of the Great East Japan Disaster on career and wage among a young generation (2012), were provided by the Social Science Japan Data Archive, Centre for Social Research and Data Archives, Institute of Social Science, The University of Tokyo. The Internet survey was conducted in January 2012 to investigate the short-term effect of the Great East Japan earthquake on the wages of college or high-school graduates focusing on the role of the quality of education, by the Economic and Social Research Institute, Cabinet Office, Government of Japan [42]. The survey was conducted according to ethical guidelines for social science research and study participation was voluntary. We did not apply for a research approval from an ethics committee because this secondary analysis used data edited by the Social Science Japan Data Archive that offers archived data to academic researchers and students without any ethical or financial requests. The data had been carefully edited so it was impossible to identify individuals from any analyses.

    Participants

    The survey recruited 12,000 young voluntary participants of a major Internet service in Japan based on quota sampling method [42]. The Economic and Social Research Institute reported that they followed this method rather than random sampling methods due to advantages in rapid assessment, and because random sampling does not always secure valid study participants [42]. The survey subjects graduated from college or high school between March 2009 and March 2011 across Japan, with a resultant age range of 17-27 years old. In this age range, there were 1.2 to 1.5 million peers in each year and 5% more females than males [43]. Consequently, we obtained 7335 participants after we removed participants with data missing for the prefectures in which they commuted. The graduation rates were 23.0% in 2009, 22.1% in 2010, 18.8% in 2011, and 19.1% in 2012. The remaining 17.1% of respondents had graduated from high school and would graduate from college or other schools in 2013 or later. The duration of employment after graduation varied from 0 (25%) to more than 51 months (2.6%), with a triple peak at 10 (10.2%), 22 (7.2%), and 34 months (5.5%) because the survey was conducted in January, taking into account the typical start of the fiscal year for businesses in April.

    Approximately 93% of respondents did not report clearly defined adverse effects caused by the disaster. A small percentage of these participants experienced the loss of second-degree relatives (0.3%) and 0.3% experienced injuries caused by the disaster. Other adverse events included collapse of their house or official evacuation because of the crisis at the Fukushima NPP.

    Perceived Health Status

    The survey asked questions directly related to the disaster. The data included a change of self-perceived health status after the disaster. The question asked was, “Did your health status change because of the Great East Japan Disaster?” There were seven categories for answers: highly improved, improved, relatively improved, unchanged, relatively decreased, decreased, and highly decreased. We categorized these seven categories into two, not-decreased and decreased, because we intended to focus on the binary difference between health and poor health. Furthermore, few respondents answered “highly improved” and “improved,” so the two categories seemed to provide a valid comparison. This binary health status change was set as the outcome variable.

    Prefecture-Level Predictors for Health Status

    We categorized geographic location based on prefectures where respondents commuted according to the values assigned by the Japan Meteorological Agency Seismic Intensity scale (JMA-SI), 0-7, and the radiation dose published by the Ministry of Education, Culture, Sports, Science and Technology (MEXT) of Japan [6,44]. When levels of JMA-SI were recorded in a prefecture, we assigned the largest value as the indicator. The value of the radiation dose was considered a continuous variable, which was determined on March 20, 2011 when it was first made available by MEXT and on January 1, 2012 when the survey was conducted [44]. We calculated the distance in a straight line between each capital city of a prefecture and the Fukushima NPP and included the number of aftershocks from March 11, 2011 to January 31, 2012 recorded in the JMA database [7].

    Six Area Indicators in Japan

    We introduced six area indicators: Northwest region; Iwate, Miyagi, and Fukushima; Kanto region; Central region; Kansai region; and Western region as shown in Figure 1. The Northwest region was an indirectly affected region because it was within the Tohoku area; however, it did not suffer at all from the tsunami and was relatively distant from the Fukushima NPP. Iwate, Miyagi, and Fukushima were directly affected regions because these three areas were most affected by the tsunami, the Fukushima incident, and persistent aftershocks. The Kanto region was an indirectly affected region because it did not bear much of the impact from the tsunami, but suffered from persistent aftershocks as well as anticipation about radiation. The Central and Kansai regions were indirectly affected regions that are relatively distant from Fukushima, but nearer compared to western Japan. The Western region was set as the control group because they were distant from Fukushima and had fewer aftershocks.

    ‎
    Figure 1. Map of Japan divided into six area indicators. The black square indicates the epicentre of the earthquake. The black star indicates the location of the Fukushima NPP.
    View this figure

    Individual-Level Demographic Factors

    Marital status was based on current status, so that widows or divorcées were classified as not married. These participants were categorized in the following manner: never married, 94.1%; married, 5.8%; widows, 0.1%; and divorced, 0.2%. Respondents were asked to answer “No” to the questions about change in employment if the corresponding conditions were caused by intentional career changes or other personal reasons. If a person had lived apart from other family members before the disaster, respondents were asked to answer “No” to questions about separation from their families. A change of working conditions included a reduction of sales (5.1%), reduction and upgrade of graduate recruitment (2.9%), merger and acquisition (0.8%), attrition (1.8%), damage of plants or capital (7.2%), temporary suspension of business (8.1%), and reduction of compensation (2.1%). Participants answered only when the condition was true. Individual economic status was assessed according to income. A difference in income between in 2010 and after the disaster in 2011 was categorized as “minus 2 levels” or “plus 1 level,” where a unit of “level” represented approximately one million yen. In contrast, an expected decrease of income in the 2012 fiscal year was answered yes or no.

    Statistical Analysis

    We calculated the percentage of people with decreased self-perceived health status in each prefecture. Because sampling error from the Internet survey seemed to render the percentages highly variable due to the small sample size in individual prefectures, we smoothed the actual percentages with empirical Bayes predictions from the multilevel (mixed effects) logistic regression by introducing random intercepts for 47 prefectures and an overall constant (the mean across all clusters/population mean) [45]. The empirical Bayes predictions can effectively “shrink” the actual percentages to a population mean, according to each prefecture’s size [46]. To adjust for differences in individual demographic factors between prefectures, we also fitted the model with gender, age, education, marital and employment status, changed job condition, income, death of family member[s], being a parent, family separation, and evacuation. Consequently, the prediction for each prefecture is presented by a regression standardization method [47]. In theory, regression standardization estimates the expectation of an outcome had all individuals been assigned a variable of interest at the specified level (an index level). Because the variable of interest in the analysis was a prefecture membership (incorporated as random intercepts) and index levels were the 47 indicators of prefectures in Japan, we presented the adjusted empirical Bayes predictions for 47 prefectures as the sum of the random intercepts and the mean of 7335 predictions from fixed parts of the model (ie, individual demographic factors) averaged over the whole population, followed by expit transformation: exp(η+ γi)/[1 + exp(η+ γi)] for mean linear predictior η of fixed effects and a random effect γi in prefecture i (i=1,…, 47). Because expit of mean of individual predictiors differs from mean of expit of individual predictiors, our adjusted empirical Bayes predictions approximated the standardized percentages in first order.

    To identify a regional association between the change of self-perceived health status and environmental factors (prefecture-level variables), the odds ratios (ORs) and their 95% CIs were calculated from fixed-effects logistic models. Predictor variables included (a) radiation dose, distance from the Fukushima NPP, JMA-SI, and number of aftershocks from March 11, 2011 to January 31, 2011 with or without (b) the area indicators in Japan. We calculated crude and adjusted ORs from univariable and from multivariable-adjusted logistic models, respectively. Model 1 included prefecture-level variables (a), Model 2 included the area indicator (b), and Model 3 simultaneously introduced all prefecture-level variables (a) and (b). We fitted different models for sensitivity analysis rather than for model building with a rejection of unnecessary covariates, although Akaike's Information Criterion (AIC, which measures prediction error by estimating the mean of Kullback-Leibler divergence over asymptotic sampling distributions) was presented not so as to search for an accurate prediction, but just as a reference for readers.

    To simultaneously estimate the association of prefecture-level variables and demographic variables (individual-level variables), including loss of family, changes in work conditions, individual economic status, experience of evacuation and separation from family, to self-perceived health status, we fitted multilevel logistic regression models that included random-intercepts for 47 prefectures and the above prefecture- and individual-level variables. The sensitivity of individual-level effects to adjustment for prefecture-level variables was analysed by fitting different models conducted in the same manner as above: Model 1 included individual-level variables only; Model 2 adjusted Model 1 with the prefecture-level variables (a); and Model 3 simultaneously adjusted for prefecture-level variables (a) and area indicators (b).

    The variables included in the models were selected from a questionnaire regarding the existing literature, which investigated the effects of the disaster on subsequent distress. Because of unclear previous knowledge on interactions, we did not conduct stratified analyses and did not include interaction terms in the multivariable models. Conformity with a linear gradient in the model was checked graphically before fitting regression models. All statistical and graphical analyses were conducted using R version 3.0.1 for Windows. The lme4 and glmmML packages were primarily used.


    Overview

    Respondents’ demographic characteristics are shown in Table 1. Those with decreased perceived health status tended to include more females, slightly more often changed their job condition, left their family member(s) for a longer duration, and had evacuated longer than those without decreased perceived health status.

    Prefecture-Level Changes in Health Status

    The decrease in self-perceived health status differed significantly among respondents in the prefectures as shown in Table 2.

    Compared with Tokyo, respondents commuting in Miyagi and Fukushima prefectures, which are located adjacent to Fukushima, showed a statistically significant decrease in health status. Respondents commuting in Iwate and Tochigi prefectures also reported an increased reduction in health status, although the differences were not statistically significant. The values were higher for respondents commuting in the Tohoku area (Iwate, Miyagi, and Fukushima), which faces the Pacific Ocean, while the values for Akita and Yamagata prefectures, which do not face the Pacific Ocean, were lower than those of Tokyo. In contrast, there were many prefectures where respondents’ health status was less likely to be reduced, particularly in Hokkaido and the central and western regions of Japan. In addition to those in Tokyo, those commuting in Okinawa, Kochi, and Toyama prefectures also reported a high frequency of decreased health status. In general, young people commuting in Tokyo reported a relatively higher reduction of their perceived health status compared with those commuting in many western regions of Japan.

    Table 1. Demographic characteristics by the change of subjective health.

    		Perceived health status
    VariablesDecreased (n=649)Not decreased (n=6686)

    		n (%) or mean (SD)n (%) or mean (SD)
    Gender, n (%)


    		Male262 (40.37)3097 (46.32)

    		Female387 (59.63)3589 (53.68)
    Age, n (%)


    		17 years old4 (0.62)49 (0.73)

    		18 years old5 (0.77)107 (1.60)

    		19 years old14 (2.16)196 (2.93)

    		20 years old24 (3.70)372 (5.56)

    		21 years old38 (5.86)459 (6.87)

    		22 years old77 (11.86)794 (11.88)

    		23 years old107 (16.49)984 (14.72)

    		24 years old127 (19.57)1363 (20.39)

    		25 years old132 (20.34)1269 (18.97)

    		26 years old67 (10.32)623 (9.32)

    		27 years old54 (8.32)470 (7.03)

    		Mean age (SD)23.74 (2.04)23.48 (2.17)
    Education, n (%)


    		College students462 (71.19)4701 (70.31)

    		Not college students187 (28.81)1985 (29.69)
    Marital status, n (%)


    		Married42 (6.47)390 (5.83)

    		Not married607 (93.53)6296 (94.17)
    Employment status, n (%)


    		Regular employee248 (38.21)2375 (35.52)

    		Not regular employee401 (61.87)4311 (64.48)
    Changed job condition, n (%)


    		Yes245 (37.50)1371 (20.51)

    		No404 (62.25)5315 (79.49)
    Difference of income 2011–2010a, mean (SD)0.27 (1.50)0.20 (1.38)
    Expected income in 2012, n (%)


    		Will be decreased45 (6.93)361 (5.40)

    		Will be increased/stable604 (93.07)6325 (94.60)
    Number of deaths of family members, n (%)


    		> 13 (0.46)16 (0.24)

    		0646 (99.54)6670 (99.76)
    Having a child/children, n (%)


    		Yes24 (3.70)263 (3.93)

    		No625 (96.30)6423 (96.07)
    Duration left for the family member(s), n (%)


    		> 4 weeks66 (10.17)396 (5.92)

    		≤ 4 weeks583 (89.83)6290 (94.08)
    Duration of evacuation, n (%)


    		> 4 weeks91 (14.02)528 (7.90)

    		≤ 4 weeks558 (85.98)6158 (92.10)

    aChange of categorical level, where a unit of “level” represented approximately one million yen. Example: difference is −2 when the level was 5 in 2010 and 3 in 2011.

    Table 2. Commuting location and the decreased self-perceived health.

    		JMA-SI (level) aRadi (µSv/h)Quake (times)Decreased (n)Not decreased (n)Decreased (%)P value
    Miyagi70.111b2841289323.1.002d
    Fukushima62.5c4211264437<.001e
    Tochigi60.1541552114320.09
    Iwate60.027230494018.19
    Ibaraki60.1763422167817.20
    Aomori50.0281883618.25
    Tokyo50.048691232164312.37Reference
    Chiba50.03315262217610.9.73
    Akita50.034644326101.00
    Saitama50.0579462118710.1.37
    Nagano50.06986887310.60
    Kanagawa50.049435403879.4.09
    Yamagata50.048623338.61
    Gunma50.089396757.22
    Niigata50.0478176827.13
    Yamanashi50.0442622307.42
    Gifu40.06221291067.8.18
    Hokkaido40.027395172296.9.01f
    Shizuoka40.037371111726.0.008d
    Aichi40.04179195333.4<.001e
    Toyama30.0477785014.69
    Shiga30.03457756111.00
    Nara30.048425618.34
    Fukui30.045543476.27
    Hyogo30.03747182905.8<.001e
    Kyoto30.03943122085.5.002d
    Mie30.046385875.048f
    Osaka30.04343316934.3<.001e
    Ishikawa30.046882524.056
    Shimane20.036383319.79
    Wakayama20.032924458.51
    Okayama20.0493151094.4.007d
    Totori20.063210280.04f
    Tokushima20.039290410.007d
    Kochi10.0263342613.78
    Kagawa10.053213407.36
    Nagasaki10.029213407.36
    Fukuoka10.03725142245.8.002d
    Oita10.05392395.22
    Hiroshima10.055271474.5.002d
    Kumamoto10.027752494.08
    Kagoshima10.0351301412.053
    Ehime10.047281602.007d
    Saga10.04120190.16
    Okinawa00.0215652915.60
    Miyazaki00.02744324111.00
    Yamaguchi00.094214616.17

    aJMA-SI, Japan Meteorological Agency seismic intensity; Radi, radiation dose on March 20, 2011 (µSv/h); Quakes, total number of aftershocks from March 11, 2011 to January 31, 2012.

    bObtained initially at 19:00 on March 29, 2011.

    cObtained initially at 13:00 on April 6, 2011.

    dMean < .01.

    eMean <.001.

    fMean < .05.

    Empirical Bayes Predictions of the Percentages of Decreased Health Status

    Empirical Bayes estimates from random-effects logistic models of each prefecture’s proportion of respondents with decreased self-perceived health status, as well as the actual percentages from Table 2 are plotted in Figure 2 according to JMA-SI level.

    The percentages of respondents commuting in Tochigi and Ibaraki prefectures as well as Iwate, Miyagi, and Fukushima are the highest. Surprisingly, the percentages for respondents commuting in Tokyo and Chiba are higher compared with the rest of Japan. In contrast, the percentage of reports of decreased health status from central to the western region (Aichi, Osaka, and prefectures located at more western than them) is clearly lower (Figure 3).

    ‎
    Figure 2. Percentage of respondents reporting decreased self-perceived health status and empirical Bayes estimates in each prefecture. JMA-SI, Japan Meteorological Agency Seismic Intensity. For the same JMA-SI levels, we determined the rank order of prefectures based on the values of adjusted empirical Bayes estimates. In adjusted empirical Bayes estimates, percentages were also adjusted according to demographic factors (gender, age, education, marital and employment status, changed job condition, income, death of family member[s], being a parent, family separation, and evacuation).
    View this figure
    ‎
    Figure 3. Map of Japan depicting adjusted empirical Bayes estimates for percentage of people with decreased self-perceived health status. Red (>21%), orange (18%–21%), yellow (15%–18%), chartreuse green (12%–15%), aquamarine (9%–12%), blue (6%–9%), and gray (<6%). The black square indicates the epicenter of the earthquake. The black star indicates the location of the Fukushima NPP.
    View this figure

    Prefecture-Level Factors Associated With Decline in Health Status

    The regional association between the change of self-perceived health status and prefecture-level variables estimated from fixed-effects logistic models are depicted in Table 3. From Model 1, the total number of quakes (OR 1.03, 95% CI 1.01-1.05; P<.001) and the seismic intensity of the primary quake (OR 1.18, 95% CI 1.08-1.29; P<.001) were associated with the prevalence of the decreased self-perceived health, but after adjustment for regional areas (Model 3), the association of seismic intensity vanished (OR 0.90, 95% CI 0.68-1.20; P=.48). On the other hand, the area of Iwate, Miyagi and Fukushima (OR 6.29, 95% CI 4.25-9.32; P<.001 in Model 2) was strongly associated with decreased self-perceived health even after adjustment for the prefecture-level factors in Model 3 (OR 4.63, 95% CI 1.10-19.57; P=.04).

    The results from random-intercept multilevel logistic models for the probability of decreased self-perceived health status including individual-level demographic variables and/or prefecture-level variables are presented in Table 4. Each model provided essentially the same results. Change in job condition (OR 2.05, 95% CI 1.72-2.45; P<.001), female (OR 1.43 [=1/0.70, inverse of OR for male 0.70, 95% CI 0.59-0.84]; P<.001), higher age (OR 1.06 per year, 95% CI 1.02-1.11; P=.005), and duration of evacuation longer than 4 weeks (OR 1.44, 95% CI 1.06-1.97; P=.02) seemed to decrease perceived health status in Model 3. As in Table 3, the total number of quakes (OR 1.02 per 100 times, 95% CI 1.00-1.05; P=.049) was strongly associated with decreased self-perceived health after adjustment for individual-level demographic factors. The OR became 1.32 (95% CI 1.03-1.71) per 1000 times. The area of Iwate, Miyagi, and Fukushima (OR 4.45, 95% CI 0.99-20.01; P=.052) and Kanto region (OR 2.94, 95% CI 0.94-9.18; P=.06) also showed a strong association. As a reviewer pointed out, although Model 3 probably suffers from the unstability of estimates due to multiple collinearity (eg, evacuated seemed obviously a function of distance to epicenter), our models provided coherent results for these variables.

    The radiation levels reported relatively soon, or considerably after the Fukushima NPP crisis and the distance from the nuclear power plant are not significantly associated with decreased health status after adjusting for covariates and potential covariates. The death of family members is not significantly associated with decreased health status because of the small number of events.

    Table 3. Prefecture-level factors associated with decreased health status: regional-level logistic regression analysis.


    		Model 1Model 2Model 3
    Independent variablesCrude ORa (95% CI)Adjusted OR (95% CI)P valueAdjusted OR (95% CI)P valueAdjusted OR (95% CI)P value
    From Fukushima NPPb (km)0.99 (0.99-1.00)1.00 (0.99-1.00).29––1.00 (1.00-1.00).65
    Radiationc (μSv/hr)2.09 (1.72-2.55)1.16 (0.87-1.53).31––1.11 (0.81-1.51).51
    Total quakesd (x10−2)1.05 (1.04-1.06)1.03 (1.01-1.05)<.001––1.03 (1.00-1.05).02
    JMA-SIe1.38 (1.29-1.47)1.18 (1.08-1.29)<.001––0.90 (0.68-1.20).46
    Northwest region0.98 (0.68-1.44)––1.69 (1.07-2.67).022.11 (0.80-5.58).13
    Iwate, Miyagi, and Fukushima3.95 (2.93-5.34)––6.29 (4.25-9.32)<.0014.63 (1.10-19.57).04
    Kanto region1.82 (1.55-2.14)––2.37 (1.77-3.16)<.0012.93 (0.99-8.66).05
    Central region0.57 (0.44-0.74)––1.05 (0.73-1.52).771.34 (0.56-3.23).51
    Kansai region0.53 (0.41-0.67)–　–　1.01 (0.72-1.43).971.23 (0.65-2.32).53


    		AICf value = 4252AIC value = 4244AIC value = 4241

    aOR, odds ratio

    bNPP, nuclear power plant

    cRadiation, radiation dose on March 20, 2011 (μSv/h)

    dTotal quakes, total number of aftershocks from March 11, 2011 to January 31, 2012

    eJMA-SI, Japan Meteorological Agency seismic intensity

    fAIC, Akaike\'s Information Criterion

    Table 4. Individual- and prefecture-level factors associated with decreased health status: multilevel logistic regression analysis with random-intercept for 47 prefectures.


    		Model 1Model 2Model 3
    Independent variablesCrude ORa (95% CI)Adjusted OR (95% CI)P valueAdjusted OR (95% CI)P valueAdjusted OR (95% CI)P value
    Regular employee (ref: not)0.89 (0.75-1.05)0.85 (0.71-1.01).070.85 (0.71-1.01).060.85 (0.71-1.02).07
    Change in job condition (ref: not)2.35 (1.98-2.79)2.12 (1.77-2.52)<.0012.05 (1.72-2.45)<.0012.05 (1.72-2.45)<.001
    Difference in income 2010–2011 (million yen incr.)1.03 (0.98-1.09)1.03 (0.97-1.09).351.03 (0.97-1.09).341.02 (0.97-1.09).35
    Income will decrease
    (ref: not)    		1.31 (0.95-1.80)1.31 (0.94-1.83).111.31 (0.94-1.82).111.31 (0.95-1.84).10
    Sex (male=1, female=0)0.78 (0.67-0.92)0.71 (0.60-0.84)<.0010.71 (0.60-0.84)<.0010.70 (0.59-0.84)<.001
    Age (year)1.06 (1.02-1.10)1.06 (1.02-1.11).0061.06 (1.02-1.11).0051.06 (1.02-1.11).005
    Marital status
    (Married = 1)    		1.12 (0.80-1.55)0.94 (0.66-1.33).720.92 (0.65-1.30).640.92 (0.65-1.30).63
    College student (ref: not)1.04 (0.87-1.25)1.00 (0.82-1.22).990.99 (0.82-1.22).9961.00 (0.83-1.27).95
    Family separation (> 4 weeks; ref: not)1.80 (1.37-2.37)1.22 (0.85-1.74).281.21 (0.85-1.74).281.21 (0.85-1.73).30
    Evacuation (> 4 weeks; ref: not)1.90 (1.50-2.41)1.48 (1.08-2.03).011.45 (1.06-1.98).021.44 (1.06-1.97).02
    Death of family members (ref: not)1.94 (0.56-6.66)2.33 (0.66-8.26).192.28 (0.65-8.03).202.34 (0.66-8.26).19
    Having a child/children (ref: not)0.94 (0.61-1.44)0.93 (0.60-1.44).750.94 (0.61-1.45).770.93 (0.59-1.44).74
    From Fukushima NPPb
    (km)    		0.99 (0.99-1.00)––0.99 (0.99-1.00).571.00 (0.99-1.00).66
    Radiationc (μSv/hr)2.09 (1.72-2.55)––1.17 (0.80-1.69).421.16 (0.82-1.64).41
    Total quakesd ( x10−2)1.05 (1.04-1.06)––1.04 (1.01-1.06).0021.02 (1.00-1.05).049
    JMA-SIe1.38 (1.29-1.47)––1.07 (0.96-1.20).240.87 (0.65-1.17).36
    Northwest region0.98 (0.68-1.44)––––2.21 (0.78-6.26).14
    Iwate, Miyagi, and Fukushima3.95 (2.93-5.34)––––4.45 (0.99-20.01).05
    Kanto region1.82 (1.55-2.14)––––2.94 (0.94-9.18).06
    Central region0.57 (0.44-0.74)––––1.49 (0.59-3.78).40
    Kansai region0.53 (0.41-0.67)–　–　–　–　1.32 (0.66-2.62).43


    		AICf: 4185
    		AIC: 4155
    		AIC: 4158

    aOR, odds ratio

    bNPP, nuclear power plant

    cRadiation, radiation dose on March 20, 2011 (μSv/h)

    dTotal quakes, total number of aftershocks from March 11, 2011 to January 31, 2012

    eJMA-SI, Japan Meteorological Agency seismic intensity

    fAIC, Akaike\'s Information Criterion


    Three Novel Findings

    We first illustrate a prefecture-level difference that shows decreased health status because of the Great East Japan Disaster. Our study found three novel findings. First, the present study shows a coherent association between decreased health status and number of aftershocks. Second, the prefecture-level radiation dose reported after the Fukushima Nuclear Crisis and the distance from the NPP of each prefecture are not significantly associated with decreased health status after adjusting for covariates and potential covariates. Third, we showed coherently that changes in job condition, being female, higher ages in the late teens and 20s, and duration of evacuation longer than 4 weeks were associated with a nationwide decreased health status even after adjusting for regional-level and prefecture-level variables.

    Related Factors for the Three Novel Findings

    People living relatively far from the disaster site tended to be more concerned about the political and economic situation [41]. The same study indicated that concerns about future health may have stimulated greater awareness of physical sensations, according to results indicating that subjective health was partially mediated by perceptions of hazard and risk [41]. Students in the Kanto region felt that they were at higher risk of a future earthquake compared with those in Western Japan [40]. Historically, the significant number of aftershocks in eastern Japan, including the Kanto region, would make subjective health decrease via a perception of hazard and risk for future earthquakes. This assumption might be valid from a nationwide perspective because radiation dose reported at the prefecture-level and at a distance from the Fukushima NPP were not significantly associated with decreased subjective health. Prefecture-level differences in subjective health could not be simply explained by radiation dose reported or distance from the Fukushima NPP. Rather, the health of youth nationwide is affected by aftershocks and sociodemographic factors.

    Consistent with a previous study [39], change in work conditions after the disaster was independently associated with decreases in the subjective health of a young cohort nationwide. Changes in work conditions, whether this was derived from natural disasters, might cause depressive symptoms, a larger burden from work, or economic decline. Being female was robustly indicated as a significantly associated factor consistent with a previous study [35], although we could not determine the reason from the present data. Surprisingly, a higher age within the age range of 17-27 years showed a significant association. We believe that our common sense tends to perceive a higher age as a risk factor for many health problems, but generally we think of this for those aged over 65 or higher. This common thinking about age range, such as those “aged 65 or older” or “children” may result in ignorance on the impact of these events at various ages. Our data regarding an association of decreased subjective health with evacuation from the disaster area must be interpreted with care. To the best of our knowledge, there has been no evidence that evacuation from a disaster area, or from Kanto region to Western region, is a risk factor for those in their late teens or 20s. Previous studies indicated evacuation as a risk factor for high mortality or hospitalization among the elderly because relocation affects living and care quality [25,48], and also the loss of their social network, which results in psychological distress [49].

    Strengths of the Study

    We believe that the present study has several significant strengths. First, we considered self-perceived health, which could include both physical and psychological aspects of health. This allowed us to evaluate general health as a whole rather than focusing on specific diseases. For instance, we could include relatively moderate illnesses that would escape inclusion in studies based on hospital records. A previous study reported that people who were professionally exposed to a disaster reported more physical and mental health complaints even in the absence of abnormal clinical laboratory values [50]. Therefore, it is definitely of value to include subjective health rather than to focus only on objective clinical diagnosis even while assuming that the percentages of self-reported illness are higher than those of clinically verified illnesses as indicated by a previous study [51]. Second, we included all prefectures in Japan rather than focusing only on the devastated area. This study thus illuminates otherwise likely-to-be-missed effects on the health of a population across the nation.

    Limitations

    Our study has several limitations. First, respondents may not represent the entire population of Japanese of the same age range because the data were not collected randomly. Therefore, the data may be biased toward participants possessing higher Internet literacy or health status and who are more likely to answer the questionnaire voluntarily. Despite the relatively limited target population analyzed here, we are confident that the results can be generalized to this entire age group across Japan, because of the common ability to access the Internet among those in this age range. Second, our data was derived from a cross-sectional survey, which does not allow determination of the direction of the relationships between demographic variables and self-perceived health status. Third, there is no information on specific reasons why participants answered that they were experiencing diminished health. However, the questionnaire was sufficiently well-controlled because all questions always included the criterion, “because of the disaster,” and not just “after the disaster.” Therefore, we could assume that we minimized the probability that the reported decrease in subjective health was derived from other causes. Last, one question asked study participants, “In what prefecture do you commute?” Therefore, we could not predetermine the location of residence. Because some employees or students might commute between prefectures, the distribution shown in the map of Japan may change if we specified the location of residence. Despite these limitations, the data presented here on post-earthquake subjective health status of this age group across the nation are worth reporting.

    Implications for Future Disaster Preparedness and Health Policy

    Future major earthquakes may affect health among the broader population, including youth, across the nation via persistent aftershocks and other socioeconomic disruption. Assessing public health status promptly across a nation is of major relevance to health policy decision makers as well as researchers looking at disasters. Assessments via the Internet may be a better measure in public health emergencies and subsequent phases compared to traditional paper-and-pencil-based surveys, especially for subgroups accustomed to Web technologies [52]. Differences in response characteristics and optimizing designs of Web surveys have been examined in some countries [15,16,52,53]. The effects of long-term low-dose radiation exposure on health and useful measurement technologies will be further elucidated in the future [54]. Traditional surveys may suffer from recall bias and low response rates, and impede arriving at reliable results. Rapid and cost-efficient assembly of health information requires the utilization of e-health technologies as well as epidemiological insights to provide better information to all decision makers. Researchers should be cautious in continuing with unconscious old-fashioned strategies for health assessments after a disaster and should better utilize epidemiology and information technologies to further knowledge in this field.

    Conclusions

    We first investigated the extent to which subjective health of participants in each prefecture across Japan decreased as a result of the Great East Japan Disaster. We found that the number of aftershocks was coherently associated with decreased subjective health. In contrast, radiation dose and distance from the Fukushima NPP were not associated. A Web-based survey can provide valuable information on public health issues after a disaster, especially if information technologies are developed that integrate with epidemiology research.

    Acknowledgments

    We appreciate the two reviewers who provided constructive comments on the previous version of our manuscript that we believe considerably improved the quality of the paper. We thank Atsuo Kishimoto, Akio Onishi, Satoko Nishimura, and Shunsuke Yamamoto, Graduate School of Public Policy, University of Tokyo, for fruitful comments on interpretation of the Great East Japan Disaster. We also thank Shigenobu Aoki for sharing R codes to draw a national map of Japan. Asami Matsunaga, Yuki Yonekura, and members of R-lovers let the first author appreciate statistical analysis with R. We thank the participants of the survey; the cabinet office of the japanese government; members of the Centre for Social Research and Data Archives, Institute of Social Science; and The University of Tokyo for their support and contribution to the public benefit. The authors declare that they have no actual or potential competing financial interests. Data for the analyses were obtained without any financial support.

    Conflicts of Interest

    None declared.

    1. Shibuya K, Hashimoto H, Ikegami N, Nishi A, Tanimoto T, Miyata H, et al. Future of Japan's system of good health at low cost with equity: beyond universal coverage. Lancet 2011 Oct 1;378(9798):1265-1273. [CrossRef] [Medline]
    2. Ohnishi T. The disaster at Japan's Fukushima-Daiichi nuclear power plant after the March 11, 2011 earthquake and tsunami, and the resulting spread of radioisotope contamination. Radiat Res 2012 Jan;177(1):1-14. [Medline]
    3. Tanigawa K, Hosoi Y, Hirohashi N, Iwasaki Y, Kamiya K. Loss of life after evacuation: lessons learned from the Fukushima accident. Lancet 2012 Mar 10;379(9819):889-891. [CrossRef] [Medline]
    4. Buesseler KO, Jayne SR, Fisher NS, Rypina II, Baumann H, Baumann Z, et al. Fukushima-derived radionuclides in the ocean and biota off Japan. Proc Natl Acad Sci U S A 2012 Apr 17;109(16):5984-5988 [FREE Full text] [CrossRef] [Medline]
    5. Nagaoka K, Sato S, Araki S, Ohta Y, Ikeuchi Y. Changes of radionuclides in the environment in Chiba, Japan, after the Fukushima nuclear power plant accident. Health Phys 2012 Apr;102(4):437-442. [CrossRef] [Medline]
    6. The official reported table of number of aftershocks over Japan Meteorological Agency seismic intensity level 1 after the main quake on March 11 in 2011.   URL: http://www.seisvol.kishou.go.jp/eq/2011_03_11_tohoku/yukan1.pdf [accessed 2013-05-17] [WebCite Cache]
    7. Database of Japan Meteorological Agency seismic intensity levels.   URL: http://www.seisvol.kishou.go.jp/eq/shindo_db/shindo_index.html [accessed 2013-05-17] [WebCite Cache]
    8. Honma M, Endo N, Osada Y, Kim Y, Kuriyama K. Disturbances in equilibrium function after major earthquake. Sci Rep 2012;2:749 [FREE Full text] [CrossRef] [Medline]
    9. Matsuoka Y, Nishi D, Nakaya N, Sone T, Noguchi H, Hamazaki K, et al. Concern over radiation exposure and psychological distress among rescue workers following the Great East Japan Earthquake. BMC Public Health 2012;12:249 [FREE Full text] [CrossRef] [Medline]
    10. Rubin GJ, Amlôt R, Wessely S, Greenberg N. Anxiety, distress and anger among British nationals in Japan following the Fukushima nuclear accident. Br J Psychiatry 2012 Nov;201(5):400-407 [FREE Full text] [CrossRef] [Medline]
    11. Nagamatsu S, Maekawa T, Ujike Y, Hashimoto S, Fuke N, Japanese Society of Intensive Care Medicine. The earthquake and tsunami-observations by Japanese physicians since the 11 March catastrophe. Crit Care 2011;15(3):167 [FREE Full text] [CrossRef] [Medline]
    12. Nagata T, Kimura Y, Ishii M. Use of a geographic information system (GIS) in the medical response to the Fukushima nuclear disaster in Japan. Prehosp Disaster Med 2012 Apr;27(2):213-215. [CrossRef] [Medline]
    13. Sato N, Yamazaki M. How generation next is rebuilding Japan. Boston: Harvard Business Publishing; 2011.   URL: http://blogs.hbr.org/cs/2011/05/how_generation_next_is_rebuild.html [accessed 2013-05-04] [WebCite Cache]
    14. Kittler AF, Hobbs J, Volk LA, Kreps GL, Bates DW. The Internet as a vehicle to communicate health information during a public health emergency: a survey analysis involving the anthrax scare of 2001. J Med Internet Res 2004 Mar 3;6(1):e8 [FREE Full text] [CrossRef] [Medline]
    15. van den Berg MH, Overbeek A, van der Pal HJ, Versluys AB, Bresters D, van Leeuwen FE, et al. Using Web-based and paper-based questionnaires for collecting data on fertility issues among female childhood cancer survivors: differences in response characteristics. J Med Internet Res 2011;13(3):e76 [FREE Full text] [CrossRef] [Medline]
    16. Ekman A, Dickman PW, Klint A, Weiderpass E, Litton JE. Feasibility of using Web-based questionnaires in large population-based epidemiological studies. Eur J Epidemiol 2006;21(2):103-111. [CrossRef] [Medline]
    17. Sugiura H, Akahane M, Ohkusa Y, Okabe N, Sano T, Jojima N, et al. Prevalence of insomnia among residents of Tokyo and osaka after the great East Japan earthquake: a prospective study. Interact J Med Res 2013;2(1):e2 [FREE Full text] [CrossRef] [Medline]
    18. Aoki T, Fukumoto Y, Yasuda S, Sakata Y, Ito K, Takahashi J, et al. The Great East Japan Earthquake Disaster and cardiovascular diseases. Eur Heart J 2012 Nov;33(22):2796-2803 [FREE Full text] [CrossRef] [Medline]
    19. Omama S, Yoshida Y, Ogasawara K, Ogawa A, Ishibashi Y, Nakamura M, et al. Influence of the great East Japan earthquake and tsunami 2011 on occurrence of cerebrovascular diseases in Iwate, Japan. Stroke 2013 Jun;44(6):1518-1524. [CrossRef] [Medline]
    20. Ogawa S, Ishiki M, Nako K, Okamura M, Senda M, Sakamoto T, et al. Effects of the Great East Japan Earthquake and huge tsunami on glycaemic control and blood pressure in patients with diabetes mellitus. BMJ Open 2012;2(2):e000830 [FREE Full text] [CrossRef] [Medline]
    21. Shiga H, Miyazawa T, Kinouchi Y, Takahashi S, Tominaga G, Takahashi H, et al. Life-event stress induced by the Great East Japan Earthquake was associated with relapse in ulcerative colitis but not Crohn's disease: a retrospective cohort study. BMJ Open 2013;3(2) [FREE Full text] [CrossRef] [Medline]
    22. Kanno T, Iijima K, Abe Y, Koike T, Shimada N, Hoshi T, et al. Hemorrhagic ulcers after Great East Japan Earthquake and Tsunami: features of post-disaster hemorrhagic ulcers. Digestion 2013;87(1):40-46. [CrossRef] [Medline]
    23. Yamanda S, Hanagama M, Kobayashi S, Satou H, Tokuda S, Niu K, et al. The impact of the 2011 Great East Japan Earthquake on hospitalisation for respiratory disease in a rapidly aging society: a retrospective descriptive and cross-sectional study at the disaster base hospital in Ishinomaki. BMJ Open 2013;3(1) [FREE Full text] [CrossRef] [Medline]
    24. Yasumura S, Goto A, Yamazaki S, Reich MR. Excess mortality among relocated institutionalized elderly after the Fukushima nuclear disaster. Public Health 2013 Feb;127(2):186-188. [CrossRef] [Medline]
    25. Nomura S, Gilmour S, Tsubokura M, Yoneoka D, Sugimoto A, Oikawa T, et al. Mortality risk amongst nursing home residents evacuated after the Fukushima nuclear accident: a retrospective cohort study. PLoS One 2013;8(3):e60192 [FREE Full text] [CrossRef] [Medline]
    26. Tsubokura M, Takita M, Matsumura T, Hara K, Tanimoto T, Kobayashi K, et al. Changes in metabolic profiles after the Great East Japan Earthquake: a retrospective observational study. BMC Public Health 2013;13:267 [FREE Full text] [CrossRef] [Medline]
    27. Minoura T, Yanagida N, Watanabe Y, Yamaoka A, Miura K. [The effects of great East Japan earthquake on patients with food allergy in Miyagi Prefecture]. Arerugi 2012 May;61(5):642-651. [Medline]
    28. Usami M, Iwadare Y, Kodaira M, Watanabe K, Aoki M, Katsumi C, et al. Relationships between traumatic symptoms and environmental damage conditions among children 8 months after the 2011 Japan earthquake and tsunami. PLoS One 2012;7(11):e50721 [FREE Full text] [CrossRef] [Medline]
    29. Ohkouchi S, Ebina M, Kamei K, Moriyama H, Tamai T, Shibuya R, et al. Fatal acute interstitial pneumonia in a worker making chips from wooden debris generated by the Great East Japan earthquake and tsunami. Respir Investig 2012 Dec;50(4):129-134. [CrossRef] [Medline]
    30. Shigemura J, Tanigawa T, Saito I, Nomura S. Psychological distress in workers at the Fukushima nuclear power plants. JAMA 2012 Aug 15;308(7):667-669. [CrossRef] [Medline]
    31. Nishi D, Koido Y, Nakaya N, Sone T, Noguchi H, Hamazaki K, et al. Peritraumatic distress, watching television, and posttraumatic stress symptoms among rescue workers after the Great East Japan earthquake. PLoS One 2012;7(4):e35248 [FREE Full text] [CrossRef] [Medline]
    32. Haraoka T, Hayasaka S, Murata C, Ojima T. Prevention of injuries and diseases in non-professional disaster volunteer activities in the Great East Japan Earthquake areas: a preliminary study. Public Health 2013 Jan;127(1):72-75. [CrossRef] [Medline]
    33. Sawa M, Osaki Y, Koishikawa H. Delayed recovery of caregivers from social dysfunction and psychological distress after the Great East Japan Earthquake. J Affect Disord 2013 Jun;148(2-3):413-417. [CrossRef] [Medline]
    34. Sugihara G, Suda S. Need for close watch on children's health after Fukushima disaster. Lancet 2011 Aug 6;378(9790):485-486. [CrossRef] [Medline]
    35. Kuwabara H, Shioiri T, Toyabe S, Kawamura T, Koizumi M, Ito-Sawamura M, et al. Factors impacting on psychological distress and recovery after the 2004 Niigata-Chuetsu earthquake, Japan: community-based study. Psychiatry Clin Neurosci 2008 Oct;62(5):503-507. [CrossRef] [Medline]
    36. Anwar J, Mpofu E, Matthews LR, Shadoul AF, Brock KE. Reproductive health and access to healthcare facilities: risk factors for depression and anxiety in women with an earthquake experience. BMC Public Health 2011;11:523 [FREE Full text] [CrossRef] [Medline]
    37. Heiervang KS, Mednick S, Sundet K, Rund BR. The psychological well-being of Norwegian adolescents exposed in utero to radiation from the Chernobyl accident. Child Adolesc Psychiatry Ment Health 2011;5:12 [FREE Full text] [CrossRef] [Medline]
    38. Hibino Y, Takaki J, Kambayashi Y, Hitomi Y, Sakai A, Sekizuka N, et al. Relationship between the Noto-Peninsula earthquake and maternal postnatal depression and child-rearing. Environ Health Prev Med 2009 Sep;14(5):255-260 [FREE Full text] [CrossRef] [Medline]
    39. Peek-Asa C, Ramirez M, Young T, Cao Y. Flood-related work disruption and poor health outcomes among university students. Prehosp Disaster Med 2012 Dec;27(6):503-508. [CrossRef] [Medline]
    40. Goodwin R, Takahashi M, Sun S, Gaines SO. Modelling psychological responses to the Great East Japan earthquake and nuclear incident. PLoS One 2012;7(5):e37690 [FREE Full text] [CrossRef] [Medline]
    41. Havenaar JM, de Wilde EJ, van den Bout J, Drottz-Sjöberg BM, van den Brink W. Perception of risk and subjective health among victims of the Chernobyl disaster. Soc Sci Med 2003;56(3):569-572. [Medline]
    42. The Short-Term Effect of the Great East Japan Earthquake on the Wage of College Graduates : the Role of the Quality of Education. 287. ESRI Discussion Paper Series No   URL: http://www.esri.go.jp/en/archive/e_dis/abstract/e_dis287-e.html [accessed 2013-05-17] [WebCite Cache]
    43. Population projection data. The Statistics Bureau and the Director-General for Policy Planning of Japan.   URL: http://www.stat.go.jp/data/jinsui/2011np/#a05k23-a [accessed 2013-05-17] [WebCite Cache]
    44. Ministry of Education, Culture, Sports, Science and Technology, Japan. Reading of environmental radioactivity level by prefecture   URL: http://radioactivity.nsr.go.jp/en/list/192/list-1.html [accessed 2013-05-17] [WebCite Cache]
    45. Gelman A, Hill J. Data Analysis Using Regression and Multilevel/Hierarchical Models. New York: Cambridge University Press; 2007.
    46. Greenland S, Robins JM. Empirical-Bayes adjustments for multiple comparisons are sometimes useful. Epidemiology 1991;2(4):244-251. [Medline]
    47. Greenland S. Introduction to regression modelling. In: Rothman KJ, Greenland S, Lash TL, et al , editors. Modern Epidemiology, Third Edition. Philadelphia: Lippincott, Williams and Wilkins; 2008:418-455.
    48. Bernard AM, Hayward RA, Rosevear J, Chun H, McMahon LF. Comparing the hospitalizations of transfer and non-transfer patients in an academic medical center. Acad Med 1996;71(3):262-266. [Medline]
    49. Oyama M, Nakamura K, Suda Y, Someya T. Social network disruption as a major factor associated with psychological distress 3 years after the 2004 Niigata-Chuetsu earthquake in Japan. Environ Health Prev Med 2012;17(2):118-123 [FREE Full text] [CrossRef] [Medline]
    50. Huizink AC, Slottje P, Witteveen AB, Bijlsma JA, Twisk JW, Smidt N, et al. Long term health complaints following the Amsterdam Air Disaster in police officers and fire-fighters. Occup Environ Med 2006;63(10):657-662 [FREE Full text] [CrossRef] [Medline]
    51. Havenaar J, Rumyantzeva G, Kasyanenko A, Kaasjager K, Westermann A, van den Brink W, et al. Health effects of the Chernobyl disaster: illness or illness behavior? A comparative general health survey in two former Soviet regions. Environ Health Perspect 1997;105 Suppl 6:1533-1537 [FREE Full text] [Medline]
    52. Kongsved SM, Basnov M, Holm-Christensen K, Hjollund NH. Response rate and completeness of questionnaires: a randomized study of Internet versus paper-and-pencil versions. J Med Internet Res 2007;9(3):e25 [FREE Full text] [CrossRef] [Medline]
    53. Ekman A, Klint A, Dickman PW, Adami HO, Litton JE. Optimizing the design of Web-based questionnaires--experience from a population-based study among 50,000 women. Eur J Epidemiol 2007;22(5):293-300. [CrossRef] [Medline]
    54. Yasumura S, Hosoya M, Yamashita S, Kamiya K, Abe M, Akashi M, Fukushima Health Management Survey Group. Study protocol for the Fukushima Health Management Survey. J Epidemiol 2012;22(5):375-383 [FREE Full text] [Medline]

    ‎
    AIC: Akaike's Information Criterion
    JMA-SI: Japan Meteorological Agency seismic intensity
    MEXT: Ministry of Education, Culture, Sports, Science and Technology
    NPP: Nuclear Power Plant
    OR: odds ratio

    Edited by G Eysenbach; submitted 26.02.13; peer-reviewed by F Pourmalek, J Reinhardt, J Davis; comments to author 28.04.13; revised version received 28.05.13; accepted 14.11.13; published 20.12.13

    Copyright

    ©Takashi Sugimoto, Tomohiro Shinozaki, Yuki Miyamoto. Originally published in the Interactive Journal of Medical Research (http://www.i-jmr.org/), 20.12.2013.

    This is an open-access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work, first published in the Interactive Journal of Medical Research, is properly cited. The complete bibliographic information, a link to the original publication on http://www.i-jmr.org/, as well as this copyright and license information must be included.