研究デザイン
この研究では、カナダのケベック州の5つの人口の多い都市での交通事故の結果に対する熱波の影響を評価しました：モントリオール（1,780,000）、ケベック市（542,298）、ラヴァル（437,413）、ロングーイユ（246,855）、およびシェルブーク（167,162人の補給（246,855）） 材料）。分析サンプルには、2015年から2022年まで、6月1日から9月30日までに報告されたすべての交通事故、負傷、および死亡者が含まれていました。
データ
交通事故、怪我、および死亡に関する毎日のデータは、 ケベック自動車保険会社 （saaq） [36]。この情報は、警察の報告に基づいてSAAQによって自動車のcrash落事故が報告されており、処理および割り当てられた交通事故の原因を割り当てられた時点で警察によって収集されます。人口情報はから調達されました ケベック統計研究所 [20]。カナダ環境省が運営する複数の気象局からの気象データ [15] （都市あたりのステーションのリストは、補足資料にあります）。さらに、2020年から2022年にかけて実施された非医薬品COVID-19介入の制御情報は、QC-NP-Covid-19インデックスから考慮されました。 [29]。この研究で使用されたデータは、SAAQデータを除き、法律25に続く具体的な承認を必要とすることを除き、公開されていました。
結果変数
3つの結果変数が分析されました。（i）怪我の種類に関係なく、自動車事故の数。 （ii）軽い交通量の総数。 （iii）殺害され、重傷を負った（KSI）個人の総数。また、クラッシュ場所の管理地域識別コードを使用して、指定された都市コングロマリットの範囲内で発生するすべてのクラッシュにフラグを立てるバイナリインジケーターを作成しました。人口100,000人あたりの割合は、3つの結果すべてについて計算されました。すべてのクラッシュ負傷、歩行者、乗客、自動車の運転手、自転車の運転手、およびすべての衝突負傷と死亡者を含む。
アーバンヒートアイランド
UHIのプロキシとして、4月1日から9月30日までの間に、Landsat 8（30 m解像度）に由来する3年間の平均暖かい季節の陸上表面温度（LST）の3年間の平均最大最大値を使用しました。それは、日中は暑くなり、夜間に過剰な熱を放射するエリアを識別し、都市部の熱島の形成に貢献します [7, 39]。 Canueは、Google Earth Engineを介して適用される公開されているアルゴリズムを使用して、カナダの郵便番号レベルでこの尺度を開発しました [13]。私たちは各都市に関連付けられたUHIを抽出し、時間の経過とともにそれらのエリア全体で地域加重された平均値を表しました [8]。
熱波
この研究では、熱波の2つの定義を採用しました。（i）2日以上連続した90パーセンタイルのしきい値を超える1日の平均温度、（ii）Cheng et al。 [21]。定義は、選択した各都市のデータを使用して適用され、クラッシュ日と一致しました。毎年、これらの値を取得する分析期間は6月1日から9月30日まででした。
制御変数
熱波の影響を適切に評価するために、Québec非薬物介入指数（QCNPI-Index）、降雨、および時間を検討しました。さまざまなソースから派生したQCNPI-Indexは、5つの都市を繰り返して、旅行、物理的距離、閉鎖/開口部の3つの重要な領域に分類された58の介入を網羅しています。 QCNPI-Indexのより高いスコアは、介入のより大きな実装を反映しており、Covid-19のスプレッドを制限するためのより包括的なアプローチを示していますが、スコアの低下は包括的なアプローチを示唆しています。文献は、Covid-19パンデミック中の複数の非医薬品介入の導入後の衝突、交通負傷、および死亡の大幅な減少を一貫して示しているため、このインデックスを既知の交絡効果を制御するために含めました（Nazif-Munoz etal。2024）。さらに、降雨量は、熱イベント中の交通結果に影響を与える時変要因と見なされました。この変数は、1日あたりの沈殿のミリメートルとして報告されています。
統計分析
この研究では、タイム層化されたケースクロスオーバーデザイン [40] 2015年6月から9月にかけて、UHI、熱波、交通事故、負傷、および死亡者との関連を調査するために採用されました。この設計は、急性の健康結果に対する短期暴露効果を分析するために環境健康研究で広く利用されています（Nhung etal。2017）。ケースクロスオーバー設計は、ケースコントロール研究のバリアントであり、各ケースは独自の制御として機能し、時間の経過とともに一定のままである個人レベルの交絡因子を効果的に制御します。このアプローチは、熱波や交通安全の結果に対する即時の影響などの一時的な暴露の調査に特に適しています。具体的には、時間延長されたケースクロスオーバーデザインは、研究期間を通常、数日または数週間の固定長の層に分割します。各ケースの日（つまり、交通事故のある日）について、同じ層から1つ以上の制御日が選択され、通常は季節、月、および年の日に一致して、季節的および長期的な傾向による潜在的な交絡を説明します。 [43]。次に、負の二項回帰を適用して、発生率率（IRR）とそれに対応する信頼区間を推定します。それは、熱波暴露と交通事故、負傷、および死亡のリスクとの関連を定量化します。負の二項回帰モデルは、次の式で表すことができます。 （1）：
$$ begin {aligned} text {log} left（e[Y] right）＆= { beta} _ {0} +{ beta} _ {1} Times Urban Heat Island +{ Beta} _ {2} Times Heatwave \＆ Quad +{ beta} _ {3} end {aligned} $$
（1）
どこ：
y =は、特定の日に発生する交通事故、軽傷、またはKSIの数を表します。
（{ beta} _ {0} ）はインターセプト用語です。
（{ beta} _ {1} ）は、UHI間の関連性の関心係数を表します。
（{ beta} _ {2} ）は、熱波暴露間の関連性の関心係数を表します。
（{ beta} _ {3} ）は、モデルに含まれる追加の共変量（降雨、時間、QCNPI-Index）の係数を表します。
ログ（母集団）は、人口規模を考慮するためのオフセットとして機能します。
感度分析には、ポアソン関数を使用しました。好みのモデル仕様に堅牢になった結果は、補足資料（表S1およびS2）に示されています。逸脱とピアソンのテストは、過剰分散に対処するのに適していることを示しているため、負の二項モデルを選択しました（表S3およびS4）。また、関心のある変数ごとに最大7日間の遅れを含めることにより、分析における潜在的な遅延効果を検討しました。ただし、これらの遅れた変数を含めると、ラグのないモデルよりもベイジアン情報基準が高くなり、モデルの適合性が低いことが示されています。結果はこの仕様にも堅牢でした（リクエストに応じて結果が利用可能です）。
2段階の設計を使用して、5つの都市全体の関心変数の不均一性を調べました。まず、式で説明されているように、各回帰の結果。 （1）各都市について、SeraとGasparriniによって提案された方法論に従って、拡張されたランダム効果メタ分析を使用して、UHIと熱波の全体的な影響のためにプールされました。 [38]。 CochranのQとI2統計は、都市全体の関連性の不均一性を評価するために使用されました。Qは、逆分散を通じて寄与を重み付けすることにより、プールされた推定値から各研究の偏差を測定しますが、I2は都市間の違いに起因する総異種の割合を表します（I2 22）。 [16] 概説して、発生率率（IRR）が都市間で大幅に異なるかどうかを調査しました。
#さまざまな熱景観を持つケベック州の5つの都市の交通安全に対する熱波の影響傷害疫学