You can describe the system with an operator acting over an input $x(t)$ transforming it into $z(t)$. If $L$ is the operator, $z(t)=L[x(t)]$.
Remember that the system and the operator are linear and time-invariant, and consider the use of the impulse function to write $x(t)$ in the following way:
$x(t)=\int_\limits{-\infty}^{\infty} x(\zeta)\delta(t-\zeta)d\zeta$.
Then:
$z(t)=L\left[\int_\limits{-\infty}^{\infty} x(\zeta)\delta(t-\zeta)d\zeta\right]$
as $L$ is a linear operator that acts over a function of $t$:
$z(t)=\int_\limits{-\infty}^{\infty} x(\zeta)L[\delta(t-\zeta)]d\zeta$
We can define $h(t,\zeta)=L[\delta(t-\zeta)]$ as the response of the system to the impulse $\delta(t-\zeta)$. This is useful for a time-invariant system (and operator) where $h(t,\zeta)=h(t-\zeta)$ and $z(t)$ can be written as:
$z(t)=\int_\limits{-\infty}^{\infty} x(\zeta)h(t-\zeta)d\zeta$
This equation is known as the convolution integral, then in a shorter form can be written as:
$z(t)=x(t)*h(t)$.
This result is interesting because it shows that a linear and time-invariant system can be described just by its impulse response $h(t)$.
If you want more detail about this, I recommend you: Probability, random variables, and random signal principles, P. Peebles, chapter 8 in 2nd ed.