Request for Help: y-axis label overlapped by x-axis in subplots in big plot

Dear All,
I would like to seek for help on this issue:
1. I set par(mfrow=c(2,2)), hoping to plot 4 subgraphs in a whole graph
2.  Each subgraph has its own x,y axes and each has x-axis label and y-axis label
3. moreover, subgraphs in the left column of the whole graph are all 3D, and have z axes and labels for z axes
4.  subgraphs in the right column of the whole graph are all 2D
5. In each subgraph, x-axis label is at the bottom, y-axis label the left side, z-axis label the right side

5. Issue
Now all subgraphs are plotted successfully, except that
* the y-axis labels for subgraphs in the right column of the whole graph are overlapped by the z-axis labels of the subgraphs in the left column of the whole graph. (meaning that y axis labels for the 2D subplots in the right column of the whole graph are not shown)

6. what I tried
6.1 When plot each graph in its own plot, everything displayed correctly.
6.2 I switched the order of this subgraphs in the whole graph, so that 3D were in the right column of the whole graph, 2D the left column. But in this case, the y axis labels of the 2D graphs are not shown (because I guess they went out of range of graphical area).

Any suggestions?