Representation of bi-parameter singular integrals by dyadic operators

We prove a dyadic representation theorem for bi-parameter singular integrals. That is, we represent certain bi-parameter operators as averages of rapidly decaying sums of what we call bi-parameter shifts. A new version of the product space T1 theorem is established as a consequence.