Insect physiology is highly dependent on the environmental temperature, and the relationship can be mathematically defined. Thus, many models that aim to predict insect-pest population dynamics, use meteorological data as input to descriptive functions that predict the development rate, survival and reproduction of pest populations. In most cases, however, these functions/models are laboratory-driven and are based on data from constant-temperature experiments. Therefore, they lack an important optimization and validation steps that test their accuracy under field conditions. Here, we developed a realistic and robust regional framework for modeling the field population dynamics of the global insect pest Bemisia tabaci. First, two non-linear functions, development rate (DR) and female reproduction (EN) were fitted to data collected in constant temperature experiments. Next, nine one-generation field experiments were conducted in order to establish a field-derived database of insect performance, representing a variety of growing conditions (different seasons, regions and host plants). Then, sensitivity analyses were performed for identifying the optimal time-scale for which the running-averaged temperatures should be fed to the model. Setting the time to 6 h (i.e., each of the 24-time steps per day represents the last 6 h average) produced the best fit (RMSD score of 1.59 days, 5.7% of the mean) between the field observations and the model simulations. We hypothesize that the 6 h ‘relevant biological time-scale' captures the insect's physiological memory of daily cycling temperature events. Lastly, we evaluated the potential of the developed modeling framework to serve as a decision support tool in pest-management programs by correlating the model predictions with field-observations of three pest control inspectors during 2019. The model successfully predicted the first notable appearance of the insect in the field (completion of the third generation in May). Also, the model correctly identified the sharp rise in abundance (outbreak point) in mid-July (completion of the fifth generation), and the persistent rise in abundance through August and September. Comparing the simulations of the 2018 and 2019 seasons indicated that the model can also serve as a tool for retrospective systematic assessment of major decisions. Taken together, these data demonstrate the model robustness and its potential to provide an excellent decision-making support platform in regional control of pest species.

Insect physiology is highly dependent on the environmental temperature, and the relationship can be mathematically defined. Thus, many models that aim to predict insect-pest population dynamics, use meteorological data as input to descriptive functions that predict the development rate, survival and reproduction of pest populations. In most cases, however, these functions/models are laboratory-driven and are based on data from constant-temperature experiments. Therefore, they lack an important optimization and validation steps that test their accuracy under field conditions. Here, we developed a realistic and robust regional framework for modeling the field population dynamics of the global insect pest Bemisia tabaci. First, two non-linear functions, development rate (DR) and female reproduction (EN) were fitted to data collected in constant temperature experiments. Next, nine one-generation field experiments were conducted in order to establish a field-derived database of insect performance, representing a variety of growing conditions (different seasons, regions and host plants). Then, sensitivity analyses were performed for identifying the optimal time-scale for which the running-averaged temperatures should be fed to the model. Setting the time to 6 h (i.e., each of the 24-time steps per day represents the last 6 h average) produced the best fit (RMSD score of 1.59 days, 5.7% of the mean) between the field observations and the model simulations. We hypothesize that the 6 h ‘relevant biological time-scale' captures the insect's physiological memory of daily cycling temperature events. Lastly, we evaluated the potential of the developed modeling framework to serve as a decision support tool in pest-management programs by correlating the model predictions with field-observations of three pest control inspectors during 2019. The model successfully predicted the first notable appearance of the insect in the field (completion of the third generation in May). Also, the model correctly identified the sharp rise in abundance (outbreak point) in mid-July (completion of the fifth generation), and the persistent rise in abundance through August and September. Comparing the simulations of the 2018 and 2019 seasons indicated that the model can also serve as a tool for retrospective systematic assessment of major decisions. Taken together, these data demonstrate the model robustness and its potential to provide an excellent decision-making support platform in regional control of pest species.

Insect physiology is highly dependent on the environmental temperature, and the relationship can be mathematically defined. Thus, many models that aim to predict insect-pest population dynamics, use meteorological data as input to descriptive functions that predict the development rate, survival and reproduction of pest populations. In most cases, however, these functions/models are laboratory-driven and are based on data from constant-temperature experiments. Therefore, they lack an important optimization and validation steps that test their accuracy under field conditions. Here, we developed a realistic and robust regional framework for modeling the field population dynamics of the global insect pest Bemisia tabaci. First, two non-linear functions, development rate (DR) and female reproduction (EN) were fitted to data collected in constant temperature experiments. Next, nine one-generation field experiments were conducted in order to establish a field-derived database of insect performance, representing a variety of growing conditions (different seasons, regions and host plants). Then, sensitivity analyses were performed for identifying the optimal time-scale for which the running-averaged temperatures should be fed to the model. Setting the time to 6 h (i.e., each of the 24-time steps per day represents the last 6 h average) produced the best fit (RMSD score of 1.59 days, 5.7% of the mean) between the field observations and the model simulations. We hypothesize that the 6 h ‘relevant biological time-scale' captures the insect's physiological memory of daily cycling temperature events. Lastly, we evaluated the potential of the developed modeling framework to serve as a decision support tool in pest-management programs by correlating the model predictions with field-observations of three pest control inspectors during 2019. The model successfully predicted the first notable appearance of the insect in the field (completion of the third generation in May). Also, the model correctly identified the sharp rise in abundance (outbreak point) in mid-July (completion of the fifth generation), and the persistent rise in abundance through August and September. Comparing the simulations of the 2018 and 2019 seasons indicated that the model can also serve as a tool for retrospective systematic assessment of major decisions. Taken together, these data demonstrate the model robustness and its potential to provide an excellent decision-making support platform in regional control of pest species.

Insect physiology is highly dependent on the environmental temperature, and the relationship can be mathematically defined. Thus, many models that aim to predict insect-pest population dynamics, use meteorological data as input to descriptive functions that predict the development rate, survival and reproduction of pest populations. In most cases, however, these functions/models are laboratory-driven and are based on data from constant-temperature experiments. Therefore, they lack an important optimization and validation steps that test their accuracy under field conditions. Here, we developed a realistic and robust regional framework for modeling the field population dynamics of the global insect pest Bemisia tabaci. First, two non-linear functions, development rate (DR) and female reproduction (EN) were fitted to data collected in constant temperature experiments. Next, nine one-generation field experiments were conducted in order to establish a field-derived database of insect performance, representing a variety of growing conditions (different seasons, regions and host plants). Then, sensitivity analyses were performed for identifying the optimal time-scale for which the running-averaged temperatures should be fed to the model. Setting the time to 6 h (i.e., each of the 24-time steps per day represents the last 6 h average) produced the best fit (RMSD score of 1.59 days, 5.7% of the mean) between the field observations and the model simulations. We hypothesize that the 6 h ‘relevant biological time-scale' captures the insect's physiological memory of daily cycling temperature events. Lastly, we evaluated the potential of the developed modeling framework to serve as a decision support tool in pest-management programs by correlating the model predictions with field-observations of three pest control inspectors during 2019. The model successfully predicted the first notable appearance of the insect in the field (completion of the third generation in May). Also, the model correctly identified the sharp rise in abundance (outbreak point) in mid-July (completion of the fifth generation), and the persistent rise in abundance through August and September. Comparing the simulations of the 2018 and 2019 seasons indicated that the model can also serve as a tool for retrospective systematic assessment of major decisions. Taken together, these data demonstrate the model robustness and its potential to provide an excellent decision-making support platform in regional control of pest species.

Insect physiology is highly dependent on the environmental temperature, and the relationship can be mathematically defined. Thus, many models that aim to predict insect-pest population dynamics, use meteorological data as input to descriptive functions that predict the development rate, survival and reproduction of pest populations. In most cases, however, these functions/models are laboratory-driven and are based on data from constant-temperature experiments. Therefore, they lack an important optimization and validation steps that test their accuracy under field conditions. Here, we developed a realistic and robust regional framework for modeling the field population dynamics of the global insect pest Bemisia tabaci. First, two non-linear functions, development rate (DR) and female reproduction (EN) were fitted to data collected in constant temperature experiments. Next, nine one-generation field experiments were conducted in order to establish a field-derived database of insect performance, representing a variety of growing conditions (different seasons, regions and host plants). Then, sensitivity analyses were performed for identifying the optimal time-scale for which the running-averaged temperatures should be fed to the model. Setting the time to 6 h (i.e., each of the 24-time steps per day represents the last 6 h average) produced the best fit (RMSD score of 1.59 days, 5.7% of the mean) between the field observations and the model simulations. We hypothesize that the 6 h ‘relevant biological time-scale' captures the insect's physiological memory of daily cycling temperature events. Lastly, we evaluated the potential of the developed modeling framework to serve as a decision support tool in pest-management programs by correlating the model predictions with field-observations of three pest control inspectors during 2019. The model successfully predicted the first notable appearance of the insect in the field (completion of the third generation in May). Also, the model correctly identified the sharp rise in abundance (outbreak point) in mid-July (completion of the fifth generation), and the persistent rise in abundance through August and September. Comparing the simulations of the 2018 and 2019 seasons indicated that the model can also serve as a tool for retrospective systematic assessment of major decisions. Taken together, these data demonstrate the model robustness and its potential to provide an excellent decision-making support platform in regional control of pest species.

Insect physiology is highly dependent on the environmental temperature, and the relationship can be mathematically defined. Thus, many models that aim to predict insect-pest population dynamics, use meteorological data as input to descriptive functions that predict the development rate, survival and reproduction of pest populations. In most cases, however, these functions/models are laboratory-driven and are based on data from constant-temperature experiments. Therefore, they lack an important optimization and validation steps that test their accuracy under field conditions. Here, we developed a realistic and robust regional framework for modeling the field population dynamics of the global insect pest Bemisia tabaci. First, two non-linear functions, development rate (DR) and female reproduction (EN) were fitted to data collected in constant temperature experiments. Next, nine one-generation field experiments were conducted in order to establish a field-derived database of insect performance, representing a variety of growing conditions (different seasons, regions and host plants). Then, sensitivity analyses were performed for identifying the optimal time-scale for which the running-averaged temperatures should be fed to the model. Setting the time to 6 h (i.e., each of the 24-time steps per day represents the last 6 h average) produced the best fit (RMSD score of 1.59 days, 5.7% of the mean) between the field observations and the model simulations. We hypothesize that the 6 h ‘relevant biological time-scale' captures the insect's physiological memory of daily cycling temperature events. Lastly, we evaluated the potential of the developed modeling framework to serve as a decision support tool in pest-management programs by correlating the model predictions with field-observations of three pest control inspectors during 2019. The model successfully predicted the first notable appearance of the insect in the field (completion of the third generation in May). Also, the model correctly identified the sharp rise in abundance (outbreak point) in mid-July (completion of the fifth generation), and the persistent rise in abundance through August and September. Comparing the simulations of the 2018 and 2019 seasons indicated that the model can also serve as a tool for retrospective systematic assessment of major decisions. Taken together, these data demonstrate the model robustness and its potential to provide an excellent decision-making support platform in regional control of pest species.

Insect physiology is highly dependent on the environmental temperature, and the relationship can be mathematically defined. Thus, many models that aim to predict insect-pest population dynamics, use meteorological data as input to descriptive functions that predict the development rate, survival and reproduction of pest populations. In most cases, however, these functions/models are laboratory-driven and are based on data from constant-temperature experiments. Therefore, they lack an important optimization and validation steps that test their accuracy under field conditions. Here, we developed a realistic and robust regional framework for modeling the field population dynamics of the global insect pest Bemisia tabaci. First, two non-linear functions, development rate (DR) and female reproduction (EN) were fitted to data collected in constant temperature experiments. Next, nine one-generation field experiments were conducted in order to establish a field-derived database of insect performance, representing a variety of growing conditions (different seasons, regions and host plants). Then, sensitivity analyses were performed for identifying the optimal time-scale for which the running-averaged temperatures should be fed to the model. Setting the time to 6 h (i.e., each of the 24-time steps per day represents the last 6 h average) produced the best fit (RMSD score of 1.59 days, 5.7% of the mean) between the field observations and the model simulations. We hypothesize that the 6 h ‘relevant biological time-scale' captures the insect's physiological memory of daily cycling temperature events. Lastly, we evaluated the potential of the developed modeling framework to serve as a decision support tool in pest-management programs by correlating the model predictions with field-observations of three pest control inspectors during 2019. The model successfully predicted the first notable appearance of the insect in the field (completion of the third generation in May). Also, the model correctly identified the sharp rise in abundance (outbreak point) in mid-July (completion of the fifth generation), and the persistent rise in abundance through August and September. Comparing the simulations of the 2018 and 2019 seasons indicated that the model can also serve as a tool for retrospective systematic assessment of major decisions. Taken together, these data demonstrate the model robustness and its potential to provide an excellent decision-making support platform in regional control of pest species.

Soil erosion affects agricultural landscapes worldwide, threatening food security and ecosystem viability. In arable environments, soil loss is primarily caused by short, intense rainstorms, typically characterized by high spatiotemporal variability. The complexity of erosive events challenges modeling efforts and explicit inclusion of extreme events in long-term risk assessment is missing. This study is intended to bridge this gap by quantifying the discrete and cumulative impacts of rainstorms on plot-scale soil erosion and providing storm-scale erosion risk analyses for a cropland region in northern Israel. Central to our analyses is the coupling of (1) a stochastic rainfall generator able to reproduce extremes down to 5-minute temporal resolutions; (2) a processes-based event-scale cropland erosion model (Dynamic WEPP, DWEPP); and, (3) a state-of-the-art frequency analysis method that explicitly accounts for rainstorms occurrence and properties. To our knowledge, this is the first study in which DWEPP runoff and soil loss are calibrated at the plot-scale on cropland (NSE is 0.82 and 0.79 for event runoff and sediment, respectively). We generated 300-year stochastic simulations of event runoff and sediment yield based on synthetic precipitation time series. Based on this data, the mean annual soil erosion in the study site is 0.1 kg m−2 [1.1 t ha−1]. Results of the risk analysis indicate that individual extreme rainstorms (>50 return period), characterized by high rainfall intensities (30-minute maximal intensity > $\sim$60 mm h−1) and high rainfall depth (>$\sim$200 mm), can trigger soil losses even one order of magnitude higher than the annual mean. The erosion efficiency of these rainstorms is mainly controlled by the short-duration (≤30 min) maximal intensities. The results demonstrate the importance of incorporating the impact of extreme events into soil conservation and management tools. We expect our methodology to be valuable for investigating future changes in soil erosion with changing climate.

Soil erosion affects agricultural landscapes worldwide, threatening food security and ecosystem viability. In arable environments, soil loss is primarily caused by short, intense rainstorms, typically characterized by high spatiotemporal variability. The complexity of erosive events challenges modeling efforts and explicit inclusion of extreme events in long-term risk assessment is missing. This study is intended to bridge this gap by quantifying the discrete and cumulative impacts of rainstorms on plot-scale soil erosion and providing storm-scale erosion risk analyses for a cropland region in northern Israel. Central to our analyses is the coupling of (1) a stochastic rainfall generator able to reproduce extremes down to 5-minute temporal resolutions; (2) a processes-based event-scale cropland erosion model (Dynamic WEPP, DWEPP); and, (3) a state-of-the-art frequency analysis method that explicitly accounts for rainstorms occurrence and properties. To our knowledge, this is the first study in which DWEPP runoff and soil loss are calibrated at the plot-scale on cropland (NSE is 0.82 and 0.79 for event runoff and sediment, respectively). We generated 300-year stochastic simulations of event runoff and sediment yield based on synthetic precipitation time series. Based on this data, the mean annual soil erosion in the study site is 0.1 kg m−2 [1.1 t ha−1]. Results of the risk analysis indicate that individual extreme rainstorms (>50 return period), characterized by high rainfall intensities (30-minute maximal intensity > $\sim$60 mm h−1) and high rainfall depth (>$\sim$200 mm), can trigger soil losses even one order of magnitude higher than the annual mean. The erosion efficiency of these rainstorms is mainly controlled by the short-duration (≤30 min) maximal intensities. The results demonstrate the importance of incorporating the impact of extreme events into soil conservation and management tools. We expect our methodology to be valuable for investigating future changes in soil erosion with changing climate.

Soil erosion affects agricultural landscapes worldwide, threatening food security and ecosystem viability. In arable environments, soil loss is primarily caused by short, intense rainstorms, typically characterized by high spatiotemporal variability. The complexity of erosive events challenges modeling efforts and explicit inclusion of extreme events in long-term risk assessment is missing. This study is intended to bridge this gap by quantifying the discrete and cumulative impacts of rainstorms on plot-scale soil erosion and providing storm-scale erosion risk analyses for a cropland region in northern Israel. Central to our analyses is the coupling of (1) a stochastic rainfall generator able to reproduce extremes down to 5-minute temporal resolutions; (2) a processes-based event-scale cropland erosion model (Dynamic WEPP, DWEPP); and, (3) a state-of-the-art frequency analysis method that explicitly accounts for rainstorms occurrence and properties. To our knowledge, this is the first study in which DWEPP runoff and soil loss are calibrated at the plot-scale on cropland (NSE is 0.82 and 0.79 for event runoff and sediment, respectively). We generated 300-year stochastic simulations of event runoff and sediment yield based on synthetic precipitation time series. Based on this data, the mean annual soil erosion in the study site is 0.1 kg m−2 [1.1 t ha−1]. Results of the risk analysis indicate that individual extreme rainstorms (>50 return period), characterized by high rainfall intensities (30-minute maximal intensity > $\sim$60 mm h−1) and high rainfall depth (>$\sim$200 mm), can trigger soil losses even one order of magnitude higher than the annual mean. The erosion efficiency of these rainstorms is mainly controlled by the short-duration (≤30 min) maximal intensities. The results demonstrate the importance of incorporating the impact of extreme events into soil conservation and management tools. We expect our methodology to be valuable for investigating future changes in soil erosion with changing climate.

Soil erosion affects agricultural landscapes worldwide, threatening food security and ecosystem viability. In arable environments, soil loss is primarily caused by short, intense rainstorms, typically characterized by high spatiotemporal variability. The complexity of erosive events challenges modeling efforts and explicit inclusion of extreme events in long-term risk assessment is missing. This study is intended to bridge this gap by quantifying the discrete and cumulative impacts of rainstorms on plot-scale soil erosion and providing storm-scale erosion risk analyses for a cropland region in northern Israel. Central to our analyses is the coupling of (1) a stochastic rainfall generator able to reproduce extremes down to 5-minute temporal resolutions; (2) a processes-based event-scale cropland erosion model (Dynamic WEPP, DWEPP); and, (3) a state-of-the-art frequency analysis method that explicitly accounts for rainstorms occurrence and properties. To our knowledge, this is the first study in which DWEPP runoff and soil loss are calibrated at the plot-scale on cropland (NSE is 0.82 and 0.79 for event runoff and sediment, respectively). We generated 300-year stochastic simulations of event runoff and sediment yield based on synthetic precipitation time series. Based on this data, the mean annual soil erosion in the study site is 0.1 kg m−2 [1.1 t ha−1]. Results of the risk analysis indicate that individual extreme rainstorms (>50 return period), characterized by high rainfall intensities (30-minute maximal intensity > $\sim$60 mm h−1) and high rainfall depth (>$\sim$200 mm), can trigger soil losses even one order of magnitude higher than the annual mean. The erosion efficiency of these rainstorms is mainly controlled by the short-duration (≤30 min) maximal intensities. The results demonstrate the importance of incorporating the impact of extreme events into soil conservation and management tools. We expect our methodology to be valuable for investigating future changes in soil erosion with changing climate.

Soil erosion affects agricultural landscapes worldwide, threatening food security and ecosystem viability. In arable environments, soil loss is primarily caused by short, intense rainstorms, typically characterized by high spatiotemporal variability. The complexity of erosive events challenges modeling efforts and explicit inclusion of extreme events in long-term risk assessment is missing. This study is intended to bridge this gap by quantifying the discrete and cumulative impacts of rainstorms on plot-scale soil erosion and providing storm-scale erosion risk analyses for a cropland region in northern Israel. Central to our analyses is the coupling of (1) a stochastic rainfall generator able to reproduce extremes down to 5-minute temporal resolutions; (2) a processes-based event-scale cropland erosion model (Dynamic WEPP, DWEPP); and, (3) a state-of-the-art frequency analysis method that explicitly accounts for rainstorms occurrence and properties. To our knowledge, this is the first study in which DWEPP runoff and soil loss are calibrated at the plot-scale on cropland (NSE is 0.82 and 0.79 for event runoff and sediment, respectively). We generated 300-year stochastic simulations of event runoff and sediment yield based on synthetic precipitation time series. Based on this data, the mean annual soil erosion in the study site is 0.1 kg m−2 [1.1 t ha−1]. Results of the risk analysis indicate that individual extreme rainstorms (>50 return period), characterized by high rainfall intensities (30-minute maximal intensity > $\sim$60 mm h−1) and high rainfall depth (>$\sim$200 mm), can trigger soil losses even one order of magnitude higher than the annual mean. The erosion efficiency of these rainstorms is mainly controlled by the short-duration (≤30 min) maximal intensities. The results demonstrate the importance of incorporating the impact of extreme events into soil conservation and management tools. We expect our methodology to be valuable for investigating future changes in soil erosion with changing climate.

Soil erosion affects agricultural landscapes worldwide, threatening food security and ecosystem viability. In arable environments, soil loss is primarily caused by short, intense rainstorms, typically characterized by high spatiotemporal variability. The complexity of erosive events challenges modeling efforts and explicit inclusion of extreme events in long-term risk assessment is missing. This study is intended to bridge this gap by quantifying the discrete and cumulative impacts of rainstorms on plot-scale soil erosion and providing storm-scale erosion risk analyses for a cropland region in northern Israel. Central to our analyses is the coupling of (1) a stochastic rainfall generator able to reproduce extremes down to 5-minute temporal resolutions; (2) a processes-based event-scale cropland erosion model (Dynamic WEPP, DWEPP); and, (3) a state-of-the-art frequency analysis method that explicitly accounts for rainstorms occurrence and properties. To our knowledge, this is the first study in which DWEPP runoff and soil loss are calibrated at the plot-scale on cropland (NSE is 0.82 and 0.79 for event runoff and sediment, respectively). We generated 300-year stochastic simulations of event runoff and sediment yield based on synthetic precipitation time series. Based on this data, the mean annual soil erosion in the study site is 0.1 kg m−2 [1.1 t ha−1]. Results of the risk analysis indicate that individual extreme rainstorms (>50 return period), characterized by high rainfall intensities (30-minute maximal intensity > $\sim$60 mm h−1) and high rainfall depth (>$\sim$200 mm), can trigger soil losses even one order of magnitude higher than the annual mean. The erosion efficiency of these rainstorms is mainly controlled by the short-duration (≤30 min) maximal intensities. The results demonstrate the importance of incorporating the impact of extreme events into soil conservation and management tools. We expect our methodology to be valuable for investigating future changes in soil erosion with changing climate.

Soil erosion affects agricultural landscapes worldwide, threatening food security and ecosystem viability. In arable environments, soil loss is primarily caused by short, intense rainstorms, typically characterized by high spatiotemporal variability. The complexity of erosive events challenges modeling efforts and explicit inclusion of extreme events in long-term risk assessment is missing. This study is intended to bridge this gap by quantifying the discrete and cumulative impacts of rainstorms on plot-scale soil erosion and providing storm-scale erosion risk analyses for a cropland region in northern Israel. Central to our analyses is the coupling of (1) a stochastic rainfall generator able to reproduce extremes down to 5-minute temporal resolutions; (2) a processes-based event-scale cropland erosion model (Dynamic WEPP, DWEPP); and, (3) a state-of-the-art frequency analysis method that explicitly accounts for rainstorms occurrence and properties. To our knowledge, this is the first study in which DWEPP runoff and soil loss are calibrated at the plot-scale on cropland (NSE is 0.82 and 0.79 for event runoff and sediment, respectively). We generated 300-year stochastic simulations of event runoff and sediment yield based on synthetic precipitation time series. Based on this data, the mean annual soil erosion in the study site is 0.1 kg m−2 [1.1 t ha−1]. Results of the risk analysis indicate that individual extreme rainstorms (>50 return period), characterized by high rainfall intensities (30-minute maximal intensity > $\sim$60 mm h−1) and high rainfall depth (>$\sim$200 mm), can trigger soil losses even one order of magnitude higher than the annual mean. The erosion efficiency of these rainstorms is mainly controlled by the short-duration (≤30 min) maximal intensities. The results demonstrate the importance of incorporating the impact of extreme events into soil conservation and management tools. We expect our methodology to be valuable for investigating future changes in soil erosion with changing climate.

Soil erosion affects agricultural landscapes worldwide, threatening food security and ecosystem viability. In arable environments, soil loss is primarily caused by short, intense rainstorms, typically characterized by high spatiotemporal variability. The complexity of erosive events challenges modeling efforts and explicit inclusion of extreme events in long-term risk assessment is missing. This study is intended to bridge this gap by quantifying the discrete and cumulative impacts of rainstorms on plot-scale soil erosion and providing storm-scale erosion risk analyses for a cropland region in northern Israel. Central to our analyses is the coupling of (1) a stochastic rainfall generator able to reproduce extremes down to 5-minute temporal resolutions; (2) a processes-based event-scale cropland erosion model (Dynamic WEPP, DWEPP); and, (3) a state-of-the-art frequency analysis method that explicitly accounts for rainstorms occurrence and properties. To our knowledge, this is the first study in which DWEPP runoff and soil loss are calibrated at the plot-scale on cropland (NSE is 0.82 and 0.79 for event runoff and sediment, respectively). We generated 300-year stochastic simulations of event runoff and sediment yield based on synthetic precipitation time series. Based on this data, the mean annual soil erosion in the study site is 0.1 kg m−2 [1.1 t ha−1]. Results of the risk analysis indicate that individual extreme rainstorms (>50 return period), characterized by high rainfall intensities (30-minute maximal intensity > $\sim$60 mm h−1) and high rainfall depth (>$\sim$200 mm), can trigger soil losses even one order of magnitude higher than the annual mean. The erosion efficiency of these rainstorms is mainly controlled by the short-duration (≤30 min) maximal intensities. The results demonstrate the importance of incorporating the impact of extreme events into soil conservation and management tools. We expect our methodology to be valuable for investigating future changes in soil erosion with changing climate.