To enable accurate system-level ESD simulation, the quasi-static I-V model of an IC is enhanced through kernel regression to reflect its circuit board dependency; alternatively, a recurrent neural network may be used to generate a non-quasi-static transient model. Hybrid electromagnetic and circuit simulation is demonstrated for ESD-induced noise coupling analysis.